Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shear correction”

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 361 records · Page 20

Lidar Measurements and High-Resolution Mesoscale Modeling of Coastally Trapped Disturbances off the Coast of California

Coastally Trapped disturbances (CTDs) are shifts in wind direction from the pre-dominant direction to equatorward to poleward for a period of time. These CTDs occur during the warm season off the California coast and impact coastal weather conditions and planned offshore wind plants. This study assesses the characteristics of CTD events as observed by lidar and other offshore buoys, then evaluates the ability of modeling systems to capture the correct characteristics, leveraging model output from the High-Resolution Rapid Refresh (HRRR) operational modeling system and the NOW-23 (National Offshore Wind) model dataset. CTDs were analyzed for October 2020 and May through to October of 2021, identifying 18 unique CTD events, confirmed by a nearby National Data Buoy Center (NDBC) buoy. The HRRR model captured most of these events, but the NOW-23 model output contained only 12 events. Composites of the wind, temperature, and pressure perturbations pre-, during, and post-event demonstrated the diminishment in wind speed, particularly for the alongshore component. Although the NOW-23 model captured the alongshore wind component and pressure perturbations well, the cross-shore wind component and temperature perturbations varied substantially. When the turbulent kinetic energy deviation and wind shear was positive across all levels pre-event, the NOW-23 modeling system was less likely to capture the CTD event. In contrast, the events that were captured by the model tended to have negative wind shear aloft pre-event.

coastal meteorology↗

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗

Modeling the Non-Linear Response of Fiber-Reinforced Laminates Using a Combined Damage/Plasticity Model

The present work is concerned with modeling the non-linear response of fiber reinforced polymer laminates. Recent experimental data suggests that the non-linearity is not only caused by matrix cracking but also by matrix plasticity due to shear stresses. To capture the effects of those two mechanisms, a model combining a plasticity formulation with continuum damage has been developed to simulate the non-linear response of laminates under plane stress states. The model is used to compare the predicted behavior of various laminate lay-ups to experimental data from the literature by looking at the degradation of axial modulus and Poisson s ratio of the laminates. The influence of residual curing stresses and in-situ effect on the predicted response is also investigated. It is shown that predictions of the combined damage/plasticity model, in general, correlate well with the experimental data. The test data shows that there are two different mechanisms that can have opposite effects on the degradation of the laminate Poisson s ratio which is captured correctly by the damage/plasticity model. Residual curing stresses are found to have a minor influence on the predicted response for the cases considered here. Some open questions remain regarding the prediction of damage onset.

Schuecker, Clara↗

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗

Data Assimilation Cycling for Weather Analysis

This software package runs the atmospheric model MM5 in data assimilation cycling mode to produce an optimized weather analysis, including the ability to insert or adjust a hurricane vortex. The program runs MM5 through a cycle of short forecasts every three hours where the vortex is adjusted to match the observed hurricane location and storm intensity. This technique adjusts the surrounding environment so that the proper steering current and environmental shear are achieved. MM5cycle uses a Cressman analysis to blend observation into model fields to get a more accurate weather analysis. Quality control of observations is also done in every cycle to remove bad data that may contaminate the analysis. This technique can assimilate and propagate data in time from intermittent and infrequent observations while maintaining the atmospheric field in a dynamically balanced state. The software consists of a C-shell script (MM5cycle.driver) and three FORTRAN programs (splitMM5files.F, comRegrid.F, and insert_vortex.F), and are contained in the pre-processor component of MM5 called "Regridder." The model is first initialized with data from a global model such as the Global Forecast System (GFS), which also provides lateral boundary conditions. These data are separated into single-time files using splitMM5.F. The hurricane vortex is then bogussed in the correct location and with the correct wind field using insert_vortex.F. The modified initial and boundary conditions are then recombined into the model fields using comRegrid.F. The model then makes a three-hour forecast. The three-hour forecast data from MM5 now become the analysis for the next short forecast run, where the vortex will again be adjusted. The process repeats itself until the desired time of analysis is achieved. This code can also assimilate observations if desired.

Tran, Nam↗

Signal processing for airborne doppler radar detection of hazardous wind shear as applied to NASA 1991 radar flight experiment data

Radar data collected during the 1991 NASA flight tests have been selectively analyzed to support research directed at developing both improved as well as new algorithms for detecting hazardous low-altitude windshear. Analysis of aircraft attitude data from several flights indicated that platform stability bandwidths were small compared to the data rate bandwidths which should support an assumption that radar returns can be treated as short time stationary. Various approaches at detection of weather returns in the presence of ground clutter are being investigated. Non-coventional clutter rejection through spectrum mode tracking and classification algorithms is a subject of continuing research. Based upon autoregressive modeling of the radar return time sequence, this approach may offer an alternative to overcome errors in conventional pulse-pair estimates. Adaptive filtering is being evaluated as a means of rejecting clutter with emphasis on low signal-to-clutter ratio situations, particularly in the presence of discrete clutter interference. An analysis of out-of-range clutter returns is included to illustrate effects of ground clutter interference due to range aliasing for aircraft on final approach. Data are presented to indicate how aircraft groundspeed might be corrected from the radar data as well as point to an observed problem of groundspeed estimate bias variation with radar antenna scan angle. A description of how recorded clutter return data are mixed with simulated weather returns is included. This enables the researcher to run controlled experiments to test signal processing algorithms. In the summary research efforts involving improved modelling of radar ground clutter returns and a Bayesian approach at hazard factor estimation are mentioned.

Baxa, Ernest G., Jr.↗

Gamma-ray bursts from sheared Alfven waves

The physical process by which sheared Alfven waves can accelerate electrons to a Lorentz factor of 10,000 to 100,000 within 5 km of the stellar surface is applied to a study of gamma-ray bursts, taking both resonant and nonresonant scattering into account. Several very encouraging features of the model are discussed. Although the field is oscillatory, virtually all the charges are ejected from the system, resulting in very little backheating of the stellar surface. The particle number density is accounted for naturally in terms of BA0 and m, which in principle are known from the physical manifestation of the agent causing the crustal disturbance. The resulting gamma-ray spectrum compares very favorably with the observation. The model restricts the geometry of the emission region, in the sense that only the Compton upscattering of soft photons from a warm polar cap can produce the correct spectral shape.

Melia, Fulvio↗

Thermocapillary Bubble Migration: Thermal Boundary Layers for Large Marangoni Numbers

The migration of an isolated gas bubble in an immiscible liquid possessing a temperature gradient is analyzed in the absence of gravity. The driving force for the bubble motion is the shear stress at the interface which is a consequence of the temperature dependence of the surface tension. The analysis is performed under conditions for which the Marangoni number is large, i.e. energy is transferred predominantly by convection. Velocity fields in the limit of both small and large Reynolds numbers are used. The thermal problem is treated by standard boundary layer theory. The outer temperature field is obtained in the vicinity of the bubble. A similarity solution is obtained for the inner temperature field. For both small and large Reynolds numbers, the asymptotic values of the scaled migration velocity of the bubble in the limit of large Marangoni numbers are calculated. The results show that the migration velocity has the same scaling for both low and large Reynolds numbers, but with a different coefficient. Higher order thermal boundary layers are analyzed for the large Reynolds number flow field and the higher order corrections to the migration velocity are obtained. Results are also presented for the momentum boundary layer and the thermal wake behind the bubble, for large Reynolds number conditions.

Balasubramaniam, R.↗

Euler and Navier-Stokes solutions for supersonic shear flow past a circular cylinder

Euler and Navier-Stokes solutions of the supersonic shear flow past a circular cylinder are obtained. These solutions are used to study the basic flow structure around the cylinder. Both the inviscid and viscous calculations show the formation of a large recirculating flow region around the front stagnation point. The calculations further show that the overall size of the recirculating region is approximately the same for the Euler and Navier-Stokes solutions but the inside structure is quite different. The inviscid flow shows only one vortex whereas the viscous flow shows two vortices inside the recirculating flow region. The inner vortex in the Navier-Stokes solution is formed primarily due to the viscous effects near the body surface and its size depends upon the Reynolds number. It is found that with increasing Reynolds number, the inner vortex diminishes in size and the Navier-Stokes solution asymptotically approaches the Euler solution. These results indicate that the Euler equations may correctly predict certain high Reynolds number separation phenomenon in flows with natural inviscid vorticity source.

Kumar, A.↗

Correction to: Investigating the rheology of fluidized and non-fluidized gas-particle beds: implications for the dynamics of geophysical flows and substrate entrainment

We report natural geophysical mass flows are among the most complex granular systems and their dynamics are often modified by the presence of an interstitial fluid. Prediction of their runout requires the development of models estimating the solid stresses in these hazardous currents wherein excess pore-fluid pressure can develop. We use discrete element modelling (DEM-CFD) with a Coarse-Graining post-processing technique (CG) to investigate the rheology of unsteady gas-particle fluidized to non-fluidized granular beds placed on horizontal and inclined planes. Similar to fluidized beds immersed in viscous fluids, the effective friction coefficient of air-fluidized beds can be defined as a function of the classic μ(I)-rheology and the non-dimensional fluid or solid pressure to explain the failure and dynamics of granular flows with excess pore pressure on inclines. However, dilation imposed by fluid drag and particle collisions in gas-particle fluidized beds can drastically change its effective frictional properties. In contrast with the common assumption in water-particle flows that granular temperature is negligible, in our gas-particle simulations, the contribution of the velocity fluctuations to the stress tensor is significant. Hence, the shear stress is found to be non-zero even when the flow is fully fluidized in the inertial regime. These results suggest the need to better understand velocity fluctuations to predict the effective viscosity of sheared fluidized mixtures and are relevant for many applications. Notably, a unified approach is useful for many geophysical flows that encompass a range of fluidization conditions in a single flow such as pyroclastic density currents and snow avalanches.

58 GEOSCIENCES↗

Preliminary Assessment of Mars Exploration Rover Landing Site Predictions

Selection of the Mars Exploration Rover (MER) landing sites took place over a three year period in which engineering constraints were identified, 155 possible sites were downselected to the final two, surface environments and safety considerations were developed, and the potential science return at the sites was considered. Landing sites in Gusev crater and Meridiani Planum were selected because they appeared acceptably safe for MER landing and roving and had strong morphologic and mineralogical indicators of liquid water in their past and thus appeared capable of addressing the science objectives of the MER missions, which are to determine the aqueous, climatic, and geologic history of sites on Mars where conditions may have been favorable to the preservation of evidence of possible pre-biotic or biotic processes. Engineering constraints important to the selection included: latitude (10 N-15 S) for maximum solar power; elevation (<-1.3 km) for sufficient atmosphere to slow the lander; low horizontal winds, shear and turbulence in the last few kilometers to minimize horizontal velocity; low 10-m scale slopes to reduce airbag spinup and bounce; moderate rock abundance to reduce abrasion or stroke-out of the airbags; and a radar-reflective, load-bearing and trafficable surface safe for landing and roving that is not dominated by fine-grained dust. In selecting the MER landing sites these engineering constraints were addressed via comprehensive evaluation of surface and atmospheric characteristics from existing remote sensing data and models as well as targeted orbital information acquired from Mars Global Surveyor and Mars Odyssey. This evaluation resulted in a number of predictions of the surface characteristics of the sites, which are tested in this abstract. Relating remote sensing signatures to surface characteristics at landing sites allows these sites to be used as ground truth for the orbital data, is essential for selecting and validating landing sites for future missions, and is required for correctly interpreting the surfaces and materials globally present on Mars.

Golombek, M.↗

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Excitation of neutron star normal modes during binary inspiral

As a compact binary inspirals due to the emission of gravitational waves, its orbital period decreases continuously down to approximately 1 ms, its value at coalescence. During the last part of the inspiral, the two stars are close together, and their tidal interactions become strong. Neutron stars have many normal modes (core g-modes, crustal discontinuity modes, shear modes, etc.) whose periods lie in the range (approximately several ms) swept by the orbital period. Some of these modes are resonantly excited by the tidal force. The amount of energy a mode absorbs is proportional to the square of the overlap integral between its displacement field and the tidal force field. For all modes of interest, this overlap is poor, resulting in relatively weak excitation. For the best case, the absorbed energy is only a small fraction (approximately 10(exp -6)) of the orbital energy, so the orbital phase shift is too weak to be detected by observations of the gravitational wave signal emitted by the inspiraling binary. However, with displacement amplitudes of excited quadrupole modes ranging up to 0.5% of the stellar radius, the possibility of a detectable electromagnetic signature cannot be dismissed. Both the periods of the modes and the energy they absorb depend quite strongly on the internal structure of the star. Their observation could shed light on the correct high-density equation of state.

Reisenegger, Andreas↗

Study of nonlinear MHD equations governing the wave propagation in twisted coronal loops

The solar corona, modelled by a low beta, resistive plasma slab, sustains MHD wave propagations due to shearing footpoint motions in the photosphere. By using a numerical algorithm the excitation and nonlinear development of MHD waves in twisted coronal loops are studied. The plasma responds to the footpoint motion by sausage waves if there is no twist. The twist in the magnetic field of the loop destroys initially developed sausage-like wave modes and they become kinks. The transition from sausage to kink modes is analyzed. The twist brings about mode degradation producing high harmonics and this generates more complex fine structures. This can be attributed to several local extrema in the perturbed velocity profiles. The Alfven wave produces remnants of the ideal 1/x singularity both for zero and non-zero twist and this pseudo-singularity becomes less pronounced for larger twist. The effect of nonlinearity is clearly observed by changing the amplitude of the driver by one order of magnitude. The magnetosonic waves also exhibit smoothed remnants of ideal logarithmic singularities when the frequency of the driver is correctly chosen. This pseudo-singularity for fast waves is absent when the coronal loop does not undergo any twist but becomes pronounced when twist is included. On the contrary, it is observed for slow waves even if there is no twist. Increasing the twist leads to a higher heating rate of the loop. The larger twist shifts somewhat uniformly distributed heating to layers inside the slab corresponding to peaks in the magnetic field strength.

Parhi, S.↗

Characterization of the NASA Langley Arc Heated Scramjet Test Facility Using NO PLIF

The nitric oxide planar laser-induced fluorescence (NO PLIF) imaging was used to characterize the air flow of the NASA Langley Arc Heated Scramjet Test Facility (AHSTF) configured with a Mach 6 nozzle. The arc raises the enthalpy of the test gas in AHSTF, producing nitric oxide. Nitric oxide persists as the temperature drops through the nozzle into the test section. NO PLIF was used to qualitatively visualize the flowfield at different experimental conditions, measure the temperature of the gas flow exiting the facility nozzle, and visualize the wave structure downstream of the nozzle at different operating conditions. Uniformity and repeatability of the nozzle flow were assessed. Expansion and compression waves on the free-jet shear layer as the nozzle flow expands into the test section were visualized. The main purpose of these experiments was to assess the uniformity of the NO in the freestream gas for planned experiments, in which NO PLIF will be used for qualitative fuel-mole-fraction sensitive imaging. The shot-to-shot fluctuations in the PLIF signal, caused by variations in the overall laser intensity as well as NO concentration and temperature variations in the flow was 20-25% of the mean signal, as determined by taking the standard deviation of a set of images obtained at constant conditions and dividing by the mean. The fluctuations within individual images, caused by laser sheet spatial variations as well as NO concentration and temperature variations in the flow, were about 28% of the mean in images, determined by taking standard deviation within individual images, dividing by the mean in the same image and averaged over the set of images. Applying an averaged laser sheet intensity correction reduced the within-image intensity fluctuations to about 10% suggesting that the NO concentration is uniform to within 10%. There was no significant difference in flow uniformity between the low and high enthalpy settings. While not strictly quantitative, the temperature maps show qualitative agreement with the computations of the flow.

Kidd, F. Gray, III↗

Mitigation of the Brighter-fatter Effect in the LSST Camera

Abstract Thick, fully depleted charge-coupled devices are known to exhibit nonlinear behavior at high signal levels due to the dynamic behavior of charges collecting in the potential wells of pixels, called the brighter-fatter effect (BFE). The effect results in distorted images of bright calibration stars, creating a flux-dependent point-spread function that if left unmitigated, could make up a large fraction of the error budget in Stage IV weak-lensing (WL) surveys such as the Legacy Survey of Space and Time (LSST). In this paper, we analyze image measurements of flat fields and artificial stars taken at different illumination levels with the LSST Camera (LSSTCam) at SLAC National Accelerator Laboratory in order to quantify this effect in the LSSTCam before and after a previously introduced correction technique. We observe that the BFE evolves anisotropically as a function of flux due to higher-order BFEs, which violates the fundamental assumption of this correction method. We then introduce a new method based on a physically motivated model to account for these higher-order terms in the correction, and then we test the modified correction on both data sets. We find that the new method corrects the effect in flat fields better than it corrects the effect in artificial stars, which we suggest is the result of sub-pixel physics not included in this correction model. We use these results to define a new metric for the full-well capacity of our sensors and advise image processing strategies to further limit the impact of the effect on LSST WL science pathways.

47 OTHER INSTRUMENTATION↗

Design and Calibration of Autonomous Coherent Doppler Lidar for Space Missions

Developed a new algorithm for the simulation of three dimensional homogeneous turbulent velocity fields. For typical atmospheric conditions it is impossible to produce a simulated velocity field that simultaneously satisfy a given spatial correlation and the corresponding spatial spectrum because of spectral aliasing. The new algorithms produce a turbulent velocity field which has accurate spatial correlations which is required for performance predictions from space-based systems. Developed a new algorithm for extracting the spatial statistics of the atmospheric velocity field using coherent Doppler lidar. The performance of the algorithm was compared with past methods and the new algorithm produces useful results for space-based data, which was not possible before. Developed new methods for verification of the errors in ground-based and space-based Doppler lidar wind measurements. These new methods do not require independent in situ data. This is an important issue for the verification of space-based Doppler lidar measurements of the global wind field. The performance of the new algorithm was compared with past results for both space-based and ground-based operation. The new algorithm has the best performance and is the only algorithm that performed satisfactory for spacebased operation. The performance of coherent Doppler lidar for a space missions with various scanning geometries was determined using computer simulation which contained the effects of random instrumental velocity errors, wind shear, wind variability along the range-gate and from shot-to-shot, and random variations in atmospheric aerosol backscatter over the measurement volume. The bias in the velocity estimates was small and the accuracy in the is typically less than 0.5 m/s for high signal conditions. For a large number of shot per velocity estimate, the threshold signal level for acceptable estimates is proportional to the number of shots to the minus one half power. This agrees with previous results determined for ground-based measurements. The use of multi-element optical detectors for autonomous operation of coherent Doppler lidar was shown to be a very promising technique. Optimal detector geometries were determined by computer simulation of performance: for ground-based testing with a fixed calibration target and for space-based operation using the random surface returns. The effects of refractive turbulence on ground-based calibration of coherent Doppler lidar was determined by computer simulations and compared with theoretical predictions. New techniques were required to correctly predict performance for the focused beam geometry commonly used for verification of space-based operation. An improved velocity estimator was evaluated for space-based applications were signal shot measurements are used to produce vector wind measurements. This permits more accurate measurements when the signal level is not known a priori or not available from multiple shot measurements. The average Doppler lidar signal spectrum including the effects of velocity turbulence was derived and calculated. This permits new estimation algorithms for turbulence based on spectral estimates. In situ atmospheric measurements were conducted and analyzed using an instrumented kite-platform. This work helps provide the required in situ data for verification of Doppler lidar velocity statistics.

Frehlich, Rod G.↗

The Interaction Between Boundary Layer and Convection Schemes in a WRF Simulation of Post‐Cold‐Frontal Clouds over the ARM East North Atlantic Site

The correct representation of lowlevel midlatitude clouds found in the wake of cold fronts strongly relies on the representation of planetary boundary layer (PBL) and convection processes, which are typically parameterized separately in numerical models. Using the Weather Research and Forecasting Model (WRF), this study investigates how distinct pairs of PBL and convection parameterization schemes represent cloud fraction in the postcoldfrontal region. The simulations focus on the region of the DOEARM Eastern North Atlantic observation site in the Azores Islands in the wake of a cold front that passed on December 25, 2015. Different PBL and convection schemes are combined to create 12 distinct configurations. The main differences between the selected physical parameterizations are the strength of vertical mixing and the entrainment. The simulations produce a wide range of cloud fractions, where some configurations significantly underestimate while others clearly overestimate satellite and surfacebased estimates of cloud fraction. A skill score is used to quantitatively assess the performance of each configuration with respect to groundbased radar data. The key processes that are found to significantly impact the cloud fraction distribution are the strength of the PBL decoupling, the vertical wind shear, entrainment and detrainment rates in shallow convection and the occurrence of drizzle. This indicates that to successfully simulate postcoldfrontal clouds, modeled physics must balance strong internal vertical mixing and weak exchange with the free troposphere. For this case study, cloud fraction was more sensitive to the choice of convection scheme than PBL scheme.

moist convection interaction between parameterizat↗