GENERIC-compliant simulations of Brownian multi-particle systems: modeling stochastic lubrication

A stochastic Lagrangian model for simulating the dynamics and rheology of a Brownian multi-particle system interacting with a simple liquid medium is presented. The discrete particle model is formulated within the GENERIC framework for Non-Equilibrium Thermodynamics and therefore it satisfies discretely the First/Second Laws of Thermodynamics and the Fluctuation Dissipation Theorem (FDT). Long-range fluctuating hydrodynamics interactions between suspended particles are described by an explicit solvent model. To this purpose, the Smoothed Dissipative Particle Dynamics method is adopted, which is a GENERIC-compliant Lagrangian meshless discretization of the fluctuating Navier–Stokes equations. In dense multi-particle systems, the average inter-particle distance is typically small compared to the particle size and short-range hydrodynamics interactions play a major role. In order to bypass an explicit—computationally costly—solution for these forces, a lubrication correction is introduced based on semi-analytical expressions for spheres under Stokes flow conditions. We generalize here the lubrication formalism to Brownian conditions, where an additional thermal-lubrication contribution needs to be taken into account in a way that discretely satisfies FDT. The coupled lubrication dynamics is integrated in time using a generalized semi-implicit splitting scheme for stochastic differential equations. The model is finally validated for a single particle diffusion as well as for a Brownian multi-particle system under homogeneous shear flow. Results for the diffusional properties as well as the rheological behavior of the whole suspension are presented and discussed.

Consideration of the effects of increasing spray cone angle and turbulence intensity on heavy-duty diesel engine pollution and specific outputs using CFD

In this article, the effects of increasing spray cone angle and turbulence intensity on the performance and emission of heavy-duty diesel engine has been examined in two separate stages using AVL-Fire CFD code. First, the injector and its spray have been simulated with various geometries. In this step, the Eulerian-Eulerian model has been applied for injector simulation and the Eulerian -Lagrangian model has been applied for spray simulation. The numerical results of this step indicate that creating swirly flow inside the nozzle decreasing penetration length while, fuel spray cone angle increasing during the injection process. In the subsequent step, the heavy-duty diesel engine has been simulated with its conventional and different nozzle hole geometries. In this step, the Eulerian-Lagrangian model has been applied to simulate the engine cycle. The numerical results of this step show that the nozzle with spiral rifling like guides has better performance and lower emission compared to other nozzle geometries. In this case, the fuel consumption is decreasing 32% than cylindrical nozzle hole, while the engine power and its torque increasing 63%. In addition, the amount of nitrogen oxide (NOx) and carbon monoxide (CO) for the spiral convergent conical nozzle geometry reducing 15% and 30% respectively than cylindrical nozzle hole while engine has no soot emission problem. Diesel injector and engine CFD results and experimental data have been validated from previous researches.

Impact of parameterized convection on a numerically simulated tropical Cyclone structure

ABSTRACT. Condensational heating is a primary source of energy for disturbances like a tropical storm. The resolvable scale condensation and the parameterized convection, in many fine mesh numerical models, are evaluated at intervals greater than the time step, order of a minute, used for computing dynamical processes. The latent heating may depend on the model resolution and the interval at which the precipitation physics is evaluated. Numerical results from a series of short range forecasts are compared to study the impact of varying the horizontal resolution and the interval for evaluating condensation physics, and of excluding the parameterized convective heating. A horizontal grid spacing of 40 km (coarse mesh) or 20 km (fine mesh) in National Centers for Environmental Prediction's Quasi-Lagrangian Model (QLM), and the initial data for a tropical storm case, are utilized. Resolvable scale condensation is invoked only at supersaturated grid points, and a Kuo-type convective parameterization procedure is employed. Significant structural differences are produced when the interval for computing both parameterized convection and resolvable scale heating is changed, and these differences broaden when the horizontal resolution is increased. The central warm cote structure and storm intensity are simulated better when both condensational processes are evaluated at an interval of twelve time steps than at each time step. Vertical columns in central storm area rapidly become convectively stable, and the maximum in vertical motion and strongest horizontal winds shift in the outer storm area, when both condensational processes are invoked at each time step. The central storm area remains conditionally unstable, and strongest winds develop close to the center, when both condensational processes are evaluated at intervals of twelve time steps. The central storm area remains conditionally unstable also in the fine mesh experiment in which the parameterized convective heating is excluded and the resolvable scale heating is evaluated at each time step. Intense vertical motion and vigorous heating develop in deep vertical columns, indicating that the heating on the convective scale is simulated as the resolvable scale heating. The vertical distribution of heating and the storm structure, during the first six hours in this case, are similar to those in the fine mesh run in which both condensational processes are evaluated at intervals of twelve time steps. However, the storm intensifies more rapidly after 6 h in the former than in the later case. Numerical results from additional experiments are presented to show that predicted storm structure is modified with a change in interval for invoking either or both condensational processes, and these circulation differences are not due to the initial spin up. Transfer of moisture and heat from low levels into the higher troposphere in cumulonimbus clouds takes place in several minutes. Above cited and other predictions from the QLM suggest that storm structure. intensity and motion in a mesoscale model are likely to, be improved when parameterized convective heating is included; however, a parameterization scheme that concurrently produces alterations in the entire model cloud depth should be invoked at intervals of several minutes.

Estimating Brewer-Dobson circulation trends from changes in stratospheric water vapour and methane

The stratospheric meridional overturning circulation, also referred to as the Brewer-Dobson circulation (BDC), controls the composition of the stratosphere, which, in turn, affects radiation and climate. As the BDC cannot be directly measured, one has to infer its strength and trends indirectly. For instance, trace gas measurements allow the calculation of average transit times. Satellite measurements provide information on the distributions of trace gases for the entire stratosphere, with measurements of particularly long and dense coverage available for stratospheric water vapour (H2O). Although chemical processes and boundary conditions confound interpretation, the influence of CH4 oxidation on H2O is relatively straightforward, and thus H2O is an appealing tracer for transport analysis despite these caveats. In this work, we explore how mean age of air trends can be estimated from the combination of stratospheric H2O and CH4 data. We carry out different sensitivity studies with the Chemical Lagrangian Model of the Stratosphere (CLaMS) and focus on the analysis of the periods of 1990–2006 and 1990–2017. In particular, we assess the methodological uncertainties related to the two commonly-used approximations of (i) instantaneous stratospheric entry mixing ratio propagation, and (ii) constant correlation between mean age and the fractional release factor of methane. Our results show that the estimated mean age of air trends from the combination of observed stratospheric H2O and CH4 changes may be significantly affected by the assumed approximations. Depending on the investigated stratospheric region and the considered period, the error in estimated mean age of air decadal trends can be large – the discrepancies are up to 10 % per decade or even more at the lower stratosphere. For particular periods, the errors from the two approximations can lead to opposite effects, which may even cancel out. Finally, we propose an improvement to the approximation method by using an idealised age spectrum to propagate stratospheric entry mixing ratios. The findings of this work can be used for improving and assessing the uncertainties in stratospheric BDC trend estimation from global satellite measurements.

Tropical cyclones prediction by numerical models in India Meteorological Department

lkj & bl 'kks/k&Ik= esa Hkkjr ekSle foKku foHkkx ¼Hkk- ekS- fo- fo-½ esa viukbZ xbZ pØokr izfr:fir djus dh dfYir rduhdksa ij ppkZ dh xbZ gSA vDrwcj 1999 esa mM+hlk esa vk, egkpØokr ds izkjfEHkd {ks=ksa esa dkYifud Hkzfeyrk dk mi;ksx djds] pØokr ds fof’k"V ekWMy] Doklh ySaxjfx;u ekWMy ¼D;w- ,y- ,e-½ ls 72 ?kaVs ds iwokZuqeku vkSj Hkkjr ekSle foKku foHkkx ds lhfer {ks= fun’kZ ¼,y- ,- ,e-½ ls 36 ?kaVs ds iwokZuqeku izfr:fir fd, x,A bl 'kks/k esa] 26 ls 28 vDrwcj rd dh izkjafHkd fLFkfr;ksa ds vk/kkj ij D;w- ,y- ,e- ls pØokr ds ekxZ ds iwokZuqeku dh vkSlr =qfV;k¡ 24 ?kaVs ds fy, 21 fd-eh-] 48 ?kaVs ds fy, 91 fd-eh- vkSj 72 ?kaVs ds fy, 179 fd-eh- jghA 1998&2004 rd ds fiNys lkr o"kksZa ds nkSjku D;w- ,y- ,e- ls pØokr ds ekxZ ds iwokZuqeku dh =qfV;ksa ds vk¡dM+ksa ij Hkh blesa ppkZ dh xbZ gSA blds vykok] ,y- ,- ,e- ls fd, x, iwokZuqeku ij izkjafHkd fLFkfr;ksa ds izHkko dh Hkh tk¡p dh xbZA fofHkUu izkjafHkd fLFkfr;ksa ls rS;kj fd, x, vkSlr ¼lesfdr½ iwokZuqeku ls 24 ?kaVs ds iwokZuqeku esa 123 fd-eh- vkSj 36 ?kaVs ds iwokZuqeku esa 81 fd-eh- dh =qfV;k¡ ikbZ xbZ] tks ,dek= iwokZuqeku dh rqyuk esa de jghA bu iz;ksxksa ls ;g irk pyk fd dkYifud Hkzfeyrk okys D;w- ,y- ,e- ekWMy ls pØokr ds ekxZ dk lVhd iwokZuqeku izkIr fd;k tk ldrk gS tks vHkh rd la[;kRed ekWMyksa ls miyC/k gks ikrk FkkA In the present paper, the cyclone bogusing techniques followed in India Meteorological Department (IMD) were discussed. Using the idealized vortex in the initial fields for Orissa super cyclone October 1999, the specialized cyclone model, Quasi-Lagrangian Model (QLM) 72 hours track forecast and also 36 hours forecast with IMD limited area model (LAM) were simulated. In this case, the QLM average track forecast errors based on 26-28 October initial conditions were 21 km for 24 hours, 91 km for 48 hours and 179 km for 72 hours. Also the QLM track forecast error statistics during the last 7 years 1998-2004 are discussed. In addition, the impact of initial conditions on the LAM forecast was examined. It was observed that the mean (ensemble) forecast generated from different initial conditions was shown track error of 123 km in 24 hours and 81 km in 36 hours forecast which is less than individual forecast. These experiments have established that the QLM model, with idealized vortex, provides track forecast within an accuracy level that was currently available from numerical models.

A Lagrangian Tool for Simulating the Transport of Chemical Pollutants in the Arabian/Persian Gulf

A rapid-response Lagrangian model for the use in simulation of the transport of a chemical pollutant in the Arabian/Persian Gulf is described. The model is well suited to the provision of a fast response after an emergency due to an accident or a deliberate spill. It is easy to set up for any situation since only requires the modification of a few input files specifying the pollutant properties and release characteristics. Running times are short, even on a desktop PC, which makes it appropriate for a rapid assessment of a hypothetical accident occurring in the region. Baroclinic circulation was obtained from an HYCOM ocean model, and tides were calculated using a barotropic model. The interactions of pollutants with sediments (uptake/release processes) were described using a dynamic approach based on kinetic transfer coefficients and a stochastic numerical method. Some examples of model applications are shown, showing the influence of the geochemical behaviour of the pollutant in its distribution patterns.

Characterization of transport from the Asian summer monsoon anticyclone into the UTLS via shedding of low-potential vorticity cutoffs

Air mass transport within the summertime Asian monsoon circulation provides a major source of anthropogenic pollution for the upper troposphere and lower stratosphere (UTLS). Here, we investigate the quasi-horizontal transport of air masses from the Asian summer monsoon anticyclone (ASMA) into the extratropical lower stratosphere and their chemical evolution. For that reason, we developed a method to identify and track the air masses exported from the monsoon. This method is based on the anomalously low potential vorticity (PV) of these air masses (tropospheric low–PV cutoffs) compared to the lower-stratosphere, and uses trajectory calculations and chemical fields from the Chemical Lagrangian Model of the Stratosphere (CLaMS). The results show evidence for frequent summertime transport from the monsoon anticyclone to mid-latitudes over the North Pacific, even reaching high latitude regions of Siberia and Alaska. Most of the low–PV cutoffs related to air masses exported from the ASMA have lifetimes shorter than one week (about 90 %) and sizes smaller than 1 percent of the northern hemisphere (NH) area. The chemical composition of these air masses is characterised by carbon monoxide, ozone and water vapour mixing ratios at an intermediate range between values typical for the monsoon anticyclone and the lower-stratosphere. The chemical evolution during transport within these low–PV cutoffs shows a gradual change from characteristics of the monsoon anticyclone to characteristics of the lower stratospheric background during about one week, indicating continuous mixing with the background atmosphere.

Numerical simulation of dropwise condensation over hydrophobic surfaces using vapor-diffusion model

Dropwise condensation of humid air over hydrophilic and hydrophobic surfaces is numerically investigated using a phenomenological, Lagrangian model. Mass flux through droplets free surface is predicted via a vapor-diffusion model. Validation with literature experimental data is successfully conducted at different air humidities and air velocities. The accuracy of the implemented condensation model is compared with a standard analogy between convective heat and mass transfer, showing that the latter is not able to predict heat transfer performances in the investigated air velocity range.

Numerical Simulation on Double-nozzle Spray Evaporation of Desulfurization Wastewater

To achieve the near zero emission of wastewater in the flue gas desulfurization (FGD) system in coal-fired power plant and better utilize the exhaust heat from flue gas, a feasible technology of spraying FGD wastewater in the flue duct for evaporation is discussed in the present study. A full-scale influencing factor investigation on the wastewater droplet evaporation performance is established under the Eulerian-Lagrangian model numerically. The dominant factors, including the characters of wastewater droplets, flue gas and the spray nozzles were analyzed under different conditions, respectively. Considering the multiple factors and conditions in the process, a Least-Square support vector machine (LSSVM) model is introduced to predict the evaporation rate based on the numerical results. Conclusions are made that the flue gas temperature and droplet diameter are of great importance in the evaporation process. The spray direction of droplet parallel with the flue gas flow direction is profitable for the dispersion of droplet, resulting the maximal evaporation rate. A double-nozzle arrangement optimized with relatively small flow rate is recommended. The LSSVM model can accurately predict the evaporation rate using the numerical results with different conditions.