Flux difference splitting for inviscid, real gases with non-equilibrium chemistry

A flux-difference splitting method is presented for the inviscid terms of the compressible flow equations for chemical non-equilibrium gases

An efficient finite difference scheme for three-dimensional, non-equilibrium flows using operator splitting

An efficient finite difference scheme is presented for the inviscid terms of the three-dimensional, compressible flow equations for chemical non-equilibrium gases. This scheme represents an extension and an improvement of one proposed by the author, and includes operator splitting.

Response theory for equilibrium and non-equilibrium statistical mechanics: causality and generalized Kramers-Kronig relations

We consider the general response theory recently proposed by Ruelle for describing the impact of small perturbations to the non-equilibrium steady states resulting from Axiom A dynamical systems. We show that the causality of the response functions entails the possibility of writing a set of Kramers-Kronig (K-K) relations for the corresponding susceptibilities at all orders of nonlinearity. Nonetheless, only a special class of directly observable susceptibilities obey K-K relations. Specific results are provided for the case of arbitrary order harmonic response, which allows for a very comprehensive K-K analysis and the establishment of sum rules connecting the asymptotic behavior of the harmonic generation susceptibility to the short-time response of the perturbed system. These results set in a more general theoretical framework previous findings obtained for optical systems and simple mechanical models, and shed light on the very general impact of considering the principle of causality for testing self-consistency: the described dispersion relations constitute unavoidable benchmarks that any experimental and model generated dataset must obey. The theory exposed in the present paper is dual to the time-dependent theory of perturbations to equilibrium states and to non-equilibrium steady states, and has in principle similar range of applicability and limitations. In order to connect the equilibrium and the non equilibrium steady state case, we show how to rewrite the classical response theory by Kubo so that response functions formally identical to those proposed by Ruelle, apart from the measure involved in the phase space integration, are obtained. These results, taking into account the chaotic hypothesis by Gallavotti and Cohen, might be relevant in several fields, including climate research. In particular, whereas the fluctuation-dissipation theorem does not work for non-equilibrium systems, because of the non-equivalence between internal and external fluctuations, K-K relations might be robust tools for the definition of a self-consistent theory of climate change.

Mesoscale simulations of organized convection: importance of convective equilibrium

The validity of convective parametrization breaks down at the resolution of mesoscale models, and the success of parametrized versus explicit treatments of convection is likely to depend on the large-scale environment. In this paper we examine the hypothesis that a key feature determining the sensitivity to the environment is whether the forcing of convection is sufficiently homogeneous and slowly varying that the convection can be considered to be in equilibrium. Two case studies of mesoscale convective systems over the UK, one where equilibrium conditions are expected and one where equilibrium is unlikely, are simulated using a mesoscale forecasting model. The time evolution of area-average convective available potential energy and the time evolution and magnitude of the timescale of convective adjustment are consistent with the hypothesis of equilibrium for case 1 and non-equilibrium for case 2. For each case, three experiments are performed with different partitionings between parametrized and explicit convection: fully parametrized convection, fully explicit convection and a simulation with significant amounts of both. In the equilibrium case, bulk properties of the convection such as area-integrated rain rates are insensitive to the treatment of convection. However, the detailed structure of the precipitation field changes; the simulation with parametrized convection behaves well and produces a smooth field that follows the forcing region, and the simulation with explicit convection has a small number of localized intense regions of precipitation that track with the mid-levelflow. For the non-equilibrium case, bulk properties of the convection such as area-integrated rain rates are sensitive to the treatment of convection. The simulation with explicit convection behaves similarly to the equilibrium case with a few localized precipitation regions. In contrast, the cumulus parametrization fails dramatically and develops intense propagating bows of precipitation that were not observed. The simulations with both parametrized and explicit convection follow the pattern seen in the other experiments, with a transition over the duration of the run from parametrized to explicit precipitation. The impact of convection on the large-scaleflow, as measured by upper-level wind and potential-vorticity perturbations, is very sensitive to the partitioning of convection for both cases. © Royal Meteorological Society, 2006. Contributions by P. A. Clark and M. E. B. Gray are Crown Copyright.

A finite difference scheme for inviscid flows with non-eequilibrium chemistry and internal energy

A finite difference scheme is presented for the inviscid terms of the equations of compressible fluid dynamics with general non-equilibrium chemistry and internal energy.

Influence of a non-ionic amphiphilic copolymer on the self-assembly of a peptide amphiphile that forms nanotapes

The influence of a non-ionic polymeric surfactant on the self-assembly of a peptide amphiphile (PA) that forms nanotapes is investigated using a combination of microscopic, scattering and spectroscopic techniques. Mixtures of Pluronic copolymer P123 with the PA C16-KTTKS in aqueous solution were studied at a fixed concentration of the PA at which it is known to self-assemble into extended nanotapes, but varying P123 concentration. We find that P123 can disrupt the formation of C16- KTTKS nanotapes, leading instead to cylindrical nanofibril structures. The spherical micelles formed by P123 at room temperature are disrupted in the presence of the PA. There is a loss of cloudiness in the solutions as the large nanotape aggregates formed by C16-KTTKS are broken up, by P123 solubilization. At least locally, b-sheet structure is retained, as confirmed by XRD and FTIR spectroscopy, even for solutions containing 20 wt% P123. This indicates, unexpectedly, that peptide secondary structure can be retained in solutions with high concentration of non-ionic surfactant. Selfassembly in this system exhibits slow kinetics towards equilibrium, the initial self-assembly being dependent on the order of mixing. Heating above the lipid chain melting temperature assists in disrupting trapped non-equilibrium states.

Thermal physics of the atmosphere

Thermal Physics of the Atmosphere offers a concise and thorough introduction on how basic thermodynamics naturally leads on to advanced topics in atmospheric physics. The book starts by covering the basics of thermodynamics and its applications in atmospheric science. The later chapters describe major applications, specific to more specialized areas of atmospheric physics, including vertical structure and stability, cloud formation, and radiative processes. The book concludes with a discussion of non-equilibrium thermodynamics as applied to the atmosphere. This book provides a thorough introduction and invaluable grounding for specialised literature on the subject. Introduces a wide range of areas associated with atmospheric physics Starts from basic level thermal physics Ideally suited for readers with a general physics background Self-assessment questions included for each chapter Supplementary website to accompany the book

A global climatology of wind–wave interaction

Generally, ocean waves are thought to act as a drag on the surface wind so that momentum is transferred downwards, from the atmosphere into the waves. Recent observations have suggested that when long wavelength waves, characteristic of remotely generated swell, propagate faster than the surface wind momentum can also be transferred upwards. This upward momentum transfer acts to accelerate the near-surface wind, resulting in a low-level wave-driven wind jet. Previous studies have suggested that the sign reversal of the momentum flux is well predicted by the inverse wave age, the ratio of the surface wind speed to the speed of the waves at the peak of the spectrum. ECMWF ERA-40 data has been used here to calculate the global distribution of the inverse wave age to determine whether there are regions of the ocean that are usually in the wind-driven wave regime and others that are generally in the wave-driven wind regime. The wind-driven wave regime is found to occur most often in the mid-latitude storm tracks where wind speeds are generally high. The wave-driven wind regime is found to be prevalent in the tropics where wind speeds are generally light and swell can propagate from storms at higher latitudes. The inverse wave age is also a useful indicator of the degree of coupling between the local wind and wave fields. The climatologies presented emphasise the non-equilibrium that exists between the local wind and wave fields and highlight the importance of swell in the global oceans.

A statistical mechanical approach for the computation of the climatic response to general forcings

The climate belongs to the class of non-equilibrium forced and dissipative systems, for which most results of quasi-equilibrium statistical mechanics, including the fluctuation-dissipation theorem, do not apply. In this paper we show for the first time how the Ruelle linear response theory, developed for studying rigorously the impact of perturbations on general observables of non-equilibrium statistical mechanical systems, can be applied with great success to analyze the climatic response to general forcings. The crucial value of the Ruelle theory lies in the fact that it allows to compute the response of the system in terms of expectation values of explicit and computable functions of the phase space averaged over the invariant measure of the unperturbed state. We choose as test bed a classical version of the Lorenz 96 model, which, in spite of its simplicity, has a well-recognized prototypical value as it is a spatially extended one-dimensional model and presents the basic ingredients, such as dissipation, advection and the presence of an external forcing, of the actual atmosphere. We recapitulate the main aspects of the general response theory and propose some new general results. We then analyze the frequency dependence of the response of both local and global observables to perturbations having localized as well as global spatial patterns. We derive analytically several properties of the corresponding susceptibilities, such as asymptotic behavior, validity of Kramers-Kronig relations, and sum rules, whose main ingredient is the causality principle. We show that all the coefficients of the leading asymptotic expansions as well as the integral constraints can be written as linear function of parameters that describe the unperturbed properties of the system, such as its average energy. Some newly obtained empirical closure equations for such parameters allow to define such properties as an explicit function of the unperturbed forcing parameter alone for a general class of chaotic Lorenz 96 models. We then verify the theoretical predictions from the outputs of the simulations up to a high degree of precision. The theory is used to explain differences in the response of local and global observables, to define the intensive properties of the system, which do not depend on the spatial resolution of the Lorenz 96 model, and to generalize the concept of climate sensitivity to all time scales. We also show how to reconstruct the linear Green function, which maps perturbations of general time patterns into changes in the expectation value of the considered observable for finite as well as infinite time. Finally, we propose a simple yet general methodology to study general Climate Change problems on virtually any time scale by resorting to only well selected simulations, and by taking full advantage of ensemble methods. The specific case of globally averaged surface temperature response to a general pattern of change of the CO2 concentration is discussed. We believe that the proposed approach may constitute a mathematically rigorous and practically very effective way to approach the problem of climate sensitivity, climate prediction, and climate change from a radically new perspective.

Evidence of dispersion relations for the nonlinear response of the Lorenz 63 system

Along the lines of the nonlinear response theory developed by Ruelle, in a previous paper we have proved under rather general conditions that Kramers-Kronig dispersion relations and sum rules apply for a class of susceptibilities describing at any order of perturbation the response of Axiom A non equilibrium steady state systems to weak monochromatic forcings. We present here the first evidence of the validity of these integral relations for the linear and the second harmonic response for the perturbed Lorenz 63 system, by showing that numerical simulations agree up to high degree of accuracy with the theoretical predictions. Some new theoretical results, showing how to derive asymptotic behaviors and how to obtain recursively harmonic generation susceptibilities for general observables, are also presented. Our findings confirm the conceptual validity of the nonlinear response theory, suggest that the theory can be extended for more general non equilibrium steady state systems, and shed new light on the applicability of very general tools, based only upon the principle of causality, for diagnosing the behavior of perturbed chaotic systems and reconstructing their output signals, in situations where the fluctuation-dissipation relation is not of great help.

Relevance of sampling schemes in light of Ruelle's linear response theory

We reconsider the theory of the linear response of non-equilibrium steady states to perturbations. We �rst show that by using a general functional decomposition for space-time dependent forcings, we can de�ne elementary susceptibilities that allow to construct the response of the system to general perturbations. Starting from the de�nition of SRB measure, we then study the consequence of taking di�erent sampling schemes for analysing the response of the system. We show that only a speci�c choice of the time horizon for evaluating the response of the system to a general time-dependent perturbation allows to obtain the formula �rst presented by Ruelle. We also discuss the special case of periodic perturbations, showing that when they are taken into consideration the sampling can be �ne-tuned to make the de�nition of the correct time horizon immaterial. Finally, we discuss the implications of our results in terms of strategies for analyzing the outputs of numerical experiments by providing a critical review of a formula proposed by Reick.

Observations of wind speed profiles over Greater London, UK, using a Doppler lidar

To calculate the potential wind loading on a tall building in an urban area, an accurate representation of the wind speed profile is required. However, due to a lack of observations, wind engineers typically estimate the characteristics of the urban boundary layer by translating the measurements from a nearby reference rural site. This study presents wind speed profile data obtained from a Doppler lidar in central London, UK, during an 8 month observation period. Used in conjunction with wind speed data measured at a nearby airport, the data have been used to assess the accuracy of the predictions made by the wind engineering tools currently available. When applied to multiple changes in surface roughness identified from morphological parameters, the non-equilibrium wind speed profile model developed by Deaves (1981) provides a good representation of the urban wind speed profile. For heights below 500 m, the predicted wind speed remains within the 95% confidence interval of the measured data. However, when the surface roughness is estimated using land use as a proxy, the model tends to overestimate the wind speed, particularly for very high wind speed periods. These results highlight the importance of a detailed assessment of the nature of the surface when estimating the wind speed above an urban surface.

Characterization and antimicrobial efficacy against E. coli of a helium/air plasma at atmospheric pressure created in a plastic package

A plasma source, sustained by the application of a floating high voltage (±15 kV) to parallel-plate electrodes at 50 Hz, has been achieved in a helium/air mixture at atmospheric pressure (P = 105 Pa) contained in a zip-locked plastic package placed in the electrode gap. Some of the physical and antimicrobial properties of this apparatus were established with a view to ascertain its performance as a prototype for the disinfection of fresh produce. The current–voltage (I–V) and charge–voltage (Q–V) characteristics of the system were measured as a function of gap distance d, in the range (3 × 103 ≤ Pd ≤ 1.0 × 104 Pa m). The electrical measurements showed this plasma source to exhibit the characteristic behaviour of a dielectric barrier discharge in the filamentary mode and its properties could be accurately interpreted by the two-capacitance in series model. The power consumed by the discharge and the reduced field strength were found to decrease quadratically from 12.0 W to 4.5 W and linearly from 140 Td to 50 Td, respectively, in the range studied. Emission spectra of the discharge were recorded on a relative intensity scale and the dominant spectral features could be assigned to strong vibrational bands in the 2+ and 1− systems of N2 and ${\rm N}_2^+$ , respectively, with other weak signatures from the NO and OH radicals and the N+, He and O atomic species. Absolute spectral intensities were also recorded and interpreted by comparison with the non-equilibrium synthetic spectra generated by the computer code SPECAIR. At an inter-electrode gap of 0.04 m, this comparison yielded typical values for the electron, vibrational and translational (gas) temperatures of (4980 ± 100) K, (2700 ± 200) K and (300 ± 100) K, respectively and an electron density of 1.0 × 1017 m−3. A Boltzmann plot also provided a value of (3200 ± 200 K) for the vibrational temperature. The antimicrobial efficacy was assessed by studying the resistance of both Escherichia coli K12 its isogenic mutants in soxR, soxS, oxyR, rpoS and dnaK selected to identify possible cellular responses and targets related with 5 min exposure to the active gas in proximity of, but not directly in, the path of the discharge filaments. Both the parent strain and mutants populations were significantly reduced by more than 1.5 log cycles in these conditions, showing the potential of the system. Post-treatment storage studies showed that some transcription regulators and specific genes related to oxidative stress play an important role in the E. coli repair mechanism and that plasma exposure affects specific cell regulator systems.

Understanding and representing atmospheric convection across scales: recommendations from the meeting held at Dartington Hall, Devon, UK, 28–30 January 2013

As weather and climate models move toward higher resolution, there is growing excitement about potential future improvements in the understanding and prediction of atmospheric convection and its interaction with larger-scale phenomena. A meeting in January 2013 in Dartington, Devon was convened to address the best way to maximise these improvements, specifically in a UK context but with international relevance. Specific recommendations included increased convective-scale observations, high-resolution virtual laboratories, and a system of parameterization test beds with a range of complexities. The main recommendation was to facilitate the development of physically based convective parameterizations that are scale-aware, non-local, non-equilibrium, and stochastic.

Characterisation of convective regimes over the British Isles

Convection-permitting modelling has led to a step change in forecasting convective events. However, convection occurs within different regimes which exhibit different forecast behaviour. A convective adjustment timescale can be used to distinguish between these regimes and examine their associated predictability. The convective adjustment timescale is calculated from radiosonde ascents and found to be consistent with that derived from convection-permitting model forecasts. The model-derived convective adjustment timescale is then examined for three summers in the British Isles to determine characteristics of the convective regimes for this maritime region. Convection in the British Isles is predominantly in convective quasi-equilibrium with 85%of convection having a timescale less than or equal to three hours. This percentage varies spatially with more non-equilibriumevents occurring in the south and southwest. The convective adjustment timescale exhibits a diurnal cycle over land. The nonequilibrium regime occurs more frequently at mid-range wind speeds and with winds from southerly to westerly sectors. Most non-equilibrium convective events in the British Isles are initiated near large coastal orographic gradients or on the European continent. Thus, the convective adjustment timescale is greatest when the location being examined is immediately downstream of large orographic gradients and decreases with distance from the convective initiation region. The dominance of convective quasiequilibrium conditions over the British Isles argues for the use of large-member ensembles in probabilistic forecasts for this region.