Laboratory simulations show diabatic heating drives cumulus-cloud evolution and entrainment

Clouds are the largest source of uncertainty in climate science, and remain a weak link in modeling tropical circulation. A major challenge is to establish connections between particulate microphysics and macroscale turbulent dynamics in cumulus clouds. Here we address the issue from the latter standpoint. First we show how to create bench-scale flows that reproduce a variety of cumulus-cloud forms (including two genera and three species), and track complete cloud life cycles-e.g., from a ``cauliflower'' congestus to a dissipating fractus. The flow model used is a transient plume with volumetric diabatic heating scaled dynamically to simulate latent-heat release from phase changes in clouds. Laser-based diagnostics of steady plumes reveal Riehl-Malkus type protected cores. They also show that, unlike the constancy implied by early self-similar plume models, the diabatic heating raises the Taylor entrainment coefficient just above cloud base, depressing it at higher levels. This behavior is consistent with cloud-dilution rates found in recent numerical simulations of steady deep convection, and with aircraft-based observations of homogeneous mixing in clouds. In-cloud diabatic heating thus emerges as the key driver in cloud development, and could well provide a major link between microphysics and cloud- scale dynamics.

Filtered density function for large eddy simulation of turbulent reacting flows

A methodology termed the “filtered density function” (FDF) is developed and implemented for large eddy simulation (LES) of chemically reacting turbulent flows. In this methodology, the effects of the unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities. A transport equation is derived for the FDF in which the effect of chemical reactions appears in a closed form. The influences of scalar mixing and convection within the subgrid are modeled. The FDF transport equation is solved numerically via a Lagrangian Monte Carlo scheme in which the solutions of the equivalent stochastic differential equations (SDEs) are obtained. These solutions preserve the Itô-Gikhman nature of the SDEs. The consistency of the FDF approach, the convergence of its Monte Carlo solution and the performance of the closures employed in the FDF transport equation are assessed by comparisons with results obtained by direct numerical simulation (DNS) and by conventional LES procedures in which the first two SGS scalar moments are obtained by a finite difference method (LES-FD). These comparative assessments are conducted by implementations of all three schemes (FDF, DNS and LES-FD) in a temporally developing mixing layer and a spatially developing planar jet under both non-reacting and reacting conditions. In non-reacting flows, the Monte Carlo solution of the FDF yields results similar to those via LES-FD. The advantage of the FDF is demonstrated by its use in reacting flows. In the absence of a closure for the SGS scalar fluctuations, the LES-FD results are significantly different from those based on DNS. The FDF results show a much closer agreement with filtered DNS results. © 1998 American Institute of Physics.

Command Tracking in High Performance Aircrafts: A New Dynamic Inversion Design

This paper proposes a new straight forward technique based on dynamic inversion, which is applied for tracking the pilot commands in high performance aircrafts.Pilot commands assumed in longitudinal mode are normal acceleration and total velocity(while roll angle and lateral acceleration are maintained at zero). In lateral mode, roll rate and total velocity are used as pilot commands (while climb rate and lateral acceleration are maintained at zero). Ensuring zero lateral acceleration leads to a better turn co-ordination. A six degree-of-freedom model of F-16 aircraft is used for both control design as well as simulation studies. Promising results are obtained which are found to be superior as compared to an existing approach (which is also based on dynamic inversion). The new approach has two potential benefits, namely reduced oscillatory response and reduced control magnitude. Another advantage of this approach is that it leads to a significant reduction of tuning parameters in the control design process.

Model-following Neuro-adaptive approach for robust control of high performance aircrafts

Based on dynamic inversion, a relatively straightforward approach is presented in this paper for nonlinear flight control design of high performance aircrafts, which does not require the normal and lateral acceleration commands to be first transferred to body rates before computing the required control inputs. This leads to substantial improvement of the tracking response. Promising results are obtained from six degree-offreedom simulation studies of F-16 aircraft, which are found to be superior as compared to an existing approach (which is also based on dynamic inversion). The new approach has two potential benefits, namely reduced oscillatory response (including elimination of non-minimum phase behavior) and reduced control magnitude. Next, a model-following neuron-adaptive design is augmented the nominal design in order to assure robust performance in the presence of parameter inaccuracies in the model. Note that in the approach the model update takes place adaptively online and hence it is philosophically similar to indirect adaptive control. However, unlike a typical indirect adaptive control approach, there is no need to update the individual parameters explicitly. Instead the inaccuracy in the system output dynamics is captured directly and then used in modifying the control. This leads to faster adaptation, which helps in stabilizing the unstable plant quicker. The robustness study from a large number of simulations shows that the adaptive design has good amount of robustness with respect to the expected parameter inaccuracies in the model.

Single Particle and Packed Bed Combustion in Modern Gasifier Stoves—Density Effects

This article is concerned with a study of an unusual effect due to density of biomass pellets in modern stoves based on close-coupled gasification-combustion process. The two processes, namely, flaming with volatiles and glowing of the char show different effects. The mass flux of the fuel bears a constant ratio with the air flow rate of gasification during the flaming process and is independent of particle density; char glowing process shows a distinct effect of density. The bed temperatures also have similar features: during flaming, they are identical, but distinct in the char burn (gasification) regime. For the cases, wood char and pellet char, the densities are 350, 990 kg/m(3), and the burn rates are 2.5 and 3.5 g/min with the bed temperatures being 1380 and 1502 K, respectively. A number of experiments on practical stoves showed wood char combustion rates of 2.5 +/- 0.5 g/min and pellet char burn rates of 3.5 +/- 0.5 g/min. In pursuit of the resolution of the differences, experimental data on single particle combustion for forced convection and ambient temperatures effects have been obtained. Single particle char combustion rate with air show a near-d(2) law and surface and core temperatures are identical for both wood and pellet char. A model based on diffusion controlled heat release-radiation-convection balance is set up. Explanation of the observed results needs to include the ash build-up over the char. This model is then used to explain observed behavior in the packed bed; the different packing densities of the biomass chars leading to different heat release rates per unit bed volume are deduced as the cause of the differences in burn rate and bed temperatures.

An experimental study of boundary layer transition induced by a cylinder wake

Boundary layer transition induced by the wake of a circular cylinder in the free stream has been investigated using the particle image velocimetry technique. Some differences between simulation and experimental studies have been reported in the literature, and these have motivated the present study. The appearance of spanwise vortices in the early stage is further confirmed here. Lambda spanwise vortex appears to evolve into a Lambda/hairpin vortex; the flow statistics also confirm such vortices. With increasing Reynolds number, based on the cylinder diameter, and with decreasing cylinder height from the plate, the physical size of these hairpin-like structures is found to decrease. Some mean flow characteristics, including the streamwise growth of the disturbance energy, in a wake-induced transition resemble those in bypass transition induced by free stream turbulence. Streamwise velocity streaks that are eventually generated in the late stage often undergo sinuous-type oscillations. Similar to other transitional flows, an inclined shear layer in the wall-normal plane is often seen to oscillate and shed vortices. The normalized shedding frequency of these vortices, estimated from the spatial spacing and the convection velocity of these vortices, is found to be independent of the Reynolds number, similar to that in ribbon-induced transition. Although the nature of free stream disturbance in a wake-induced transition and that in a bypass transition are different, the late-stage features including the flow breakdown characteristics of these two transitions appear to be similar.

Nature of Electronic States in Atomically Thin MoS(2) Field-Effect Transistors

We present low-temperature electrical transport experiments in five field-effect transistor devices consisting of monolayer, bilayer, and trilayer MoS(2) films, mechanically exfoliated onto Si/SiO(2) substrate. Our experiments reveal that the electronic states In all films are localized well up to room temperature over the experimentally accessible range of gate voltage. This manifests in two-dimensional (2D) variable range hopping (VRH) at high temperatures, while below similar to 30 K, the conductivity displays oscillatory structures In gate voltage arising from resonant tunneling at the localized sites. From the correlation energy (T(0)) of VRH and gate voltage dependence of conductivity, we suggest that Coulomb potential from trapped charges In the substrate is the dominant source of disorder in MoS(2) field-effect devices, which leads to carrier localization, as well.

Modeling of growth and motion of equiaxed dendrites in a convecting melt

A numerical micro-scale model is developed to study the behavior of dendrite growth in presence of melt convection. In this method, an explicit, coupled enthalpy model is used to simulate the growth of an equiaxed dendrite, while a Volume of Fluid (VOF) method is used to track the movement of the dendrite in the convecting melt in a two-dimensional Eulerian framework. Numerical results demonstrate the effectiveness of the enthalpy model in simulating the dendritic growth involving complex shape, and the accuracy of VOF method in conserving mass and preserving the complex dendritic shape during motion. Simulations are performed in presence of uniform melt flow for both fixed and moving dendrites, and the difference in dendrite morphology is shown.

Filling Characteristics for an Activated Carbon Based Adsorbed Natural Gas Storage System

The storage capacity of an activated carbon bed is studied using a 2D transport model with constant inlet flow conditions. The predicted filling times and variation in bed pressure and temperature are in good agreement with experimental observations obtained using a 1.82 L prototype ANG storage cylinder. Storage efficiencies based on the maximum achievable V/V (volume of gas/volume of container) and filling times are used to quantify the performance of the charging process. For the high permeability beds used in the experiments, storage efficiencies are controlled by the rate of heat removal. Filling times, defined as the time at which the bed pressure reaches 3.5 MPa, range from 120 to 3.4 min for inlet flow rates of 1.0 L min(-1) and 30.0 L min(-1), respectively. The corresponding storage efficiencies, eta(s), vary from 90% to 76%, respectively. Simulations with L/D ratios ranging from 0.35 to 7.8 indicate that the storage efficiencies can be improved with an increase in the LID ratios and/or with water cooled convection. Thus for an inlet flow rate of 30.0 L min(-1), an eta(s) value of 90% can be obtained with water cooling for an L/D ratio of 7.8 and a filling time of a few minutes. In the absence of water cooling the eta(s) value reduces to 83% at the same L/D ratio. Our study suggests that with an appropriate choice of cylinder dimensions, solutions based on convective cooling during adsorptive storage are possible with some compromise in the storage capacity.

A Mechanism for the Occurrence of Ions in Distilled Water

Solar distillation can be used to produce potable water from contaminated water. However, studies show that ions such as F(-) and NO(3)(-) occur in distillates from solar stills. In order to understand the reasons for this behavior, imaging and distillation experiments were conducted. White dots were seen in the vapor space above the interface of hot water poured into containers. The concentrations of various ions such as F(-) and SO(4)(2-) in the distillates from thermal and solar distillation experiments were roughly comparable when the feed consisted of deionized water and also solutions having fluoride concentrations of 100 and 10 000 mg/L. These observations suggest that aerosols enter the distillation setup through leaks and provide nuclei for the condensation of water vapor. The water-soluble component of aerosols dissolves in the drops formed, and some of the drops are transferred to the distillate by buoyancy-driven convection.

Computer-aided analysis of high-speed protective relays

Most of the modern distance relays are designed to avoid overreaching due to the transient d.c. component of the fault current, whereas a more likely source of transients in e.h.v. systems is the oscillatory discharge of the system charging current into the fault. Until now attempts have not been made to reproduce these transients in the laboratory. This paper describes an analogue and an accurate digital simulation of these harmonic transients. The dynamic behaviour of a typical polarised mho-type relay is analysed, and results are presented. The paper also advocates the use of active filters for filtering the harmonics associated with e.h.v. system, and hence, to improve the speed of response and accuracy of the protective relays.

Validation of a Modified Eddy Dissipation Concept Model for Stationary and Nonstationary Turbulent Diffusion Flames

A transient flame simulation tool based on unsteady Reynolds average Navier Stokes (RANS) is characterized for stationary and nonstationary flame applications with a motivation of performing computationally affordable flame stability studies. Specifically, the KIVA-3V code is utilized with incorporation of a recently proposed modified eddy dissipation concept for simulating turbulence-chemistry interaction along with a model for radiation loss. Detailed comparison of velocities, turbulent kinetic energies, temperature, and species are made with the experimental data of the turbulent, non-premixed DLR_A CH4/H-2/N-2 jet flame. The comparison shows that the model is able to predict flame structure very well. The effect of some of the modeling assumptions is assessed, and strategies to model a stationary diffusion flame are recommended. Unsteady flame simulation capabilities of the numerical model are assessed by simulating an acoustically excited, experimental, oscillatory H-2-air diffusion flame. Comparisons are made with oscillatory velocity field and OH plots, and the numerical code is observed to predict transient flame structure well.

Simple theoretical analysis of the photoemission from quantum confined effective mass superlattices of optoelectronic materials

The photoemission from quantum wires and dots of effective mass superlattices of optoelectronic materials was investigated on the basis of newly formulated electron energy spectra, in the presence of external light waves, which controls the transport properties of ultra-small electronic devices under intense radiation. The effect of magnetic quantization on the photoemission from the aforementioned superlattices, together with quantum well superlattices under magnetic quantization, has also been investigated in this regard. It appears, taking HgTe/Hg1-xCdxTe and InxGa1-xAs/InP effective mass superlattices, that the photoemission from these quantized structures is enhanced with increasing photon energy in quantized steps and shows oscillatory dependences with the increasing carrier concentration. In addition, the photoemission decreases with increasing light intensity and wavelength as well as with increasing thickness exhibiting oscillatory spikes. The strong dependence of the photoemission on the light intensity reflects the direct signature of light waves on the carrier energy spectra. The content of this paper finds six different applications in the fields of low dimensional systems in general.

Dynamical behavior in the nonlinear rheology of surfactant solutions

Several surfactant molecules self-assemble in solution to form long, flexible wormlike micelles which get entangled with each other, leading to viscoelastic gel phases. We discuss our recent work on the rheology of such a gel formed in the dilute aqueous solutions of a surfactant CTAT. In the linear rheology regime, the storage modulus G′(ω) and loss modulus G″(ω) have been measured over a wide frequency range. In the nonlinear regime, the shear stress σ shows a plateau as a function of the shear rate math above a certain cutoff shear rate mathc. Under controlled shear rate conditions in the plateau regime, the shear stress and the first normal stress difference show oscillatory time-dependence. The analysis of the measured time series of shear stress and normal stress has been done using several methods incorporating state space reconstruction by embedding of time delay vectors. The analysis shows the existence of a finite correlation dimension and a positive Lyapunov exponent, unambiguously implying that the dynamics of the observed mechanical instability can be described by that of a dynamical system with a strange attractor of dimension varying from 2.4 to 2.9.

Down-to-earth temperatures: The mechanics of the thermal environment

This lecture describes some recent attempts at unravelling the mechanics of the temperature distribution near ground, especially during calm, clear nights. In particular, a resolution is offered of the so-called Ramdas paradox, connected with observations of a temperature minimum some decimetres above bare soil on calm clear nights, in apparent defiance of the Rayleigh criterion for instability due to thermal convection. The dynamics of the associated temperature distribution is governed by radiative and convective transport and by thermal conduction, and is characterised by two time constants, involving respectively quick radiative adjustments and slow diffusive relaxation. The theory underlying the work described here suggests that surface parameters like ground emissivity and soil thermal conductivity can exert appreciable influence on the development of nocturnal inversions.