Modelling of atmospheric pollutant transport and deposition in lakeshore environments

In order to study the effects of lakes on weather modification, pollutant transport and deposition, the hydrodynamic equations were solved numerically on a variable horizontal and vertical grid. Time splitting was applied to couple an explicit pseudospectral and an implicit finite difference scheme. The resulted model difference equations were used to predict the local wind, temperature, and pressure fields as well as the evolution of the atmospheric planetary boundary layer parameters. The modeling domain was centered around a lake to investigate the potential impact of nearby industrial sources. The situation studied was characteristic of summertime conditions, driven by a synoptic weather system with air moving from land, over a cooler lake surface. This flow is shown to become highly modified by a mesoscale meteorological phenomenon known as lake-breeze. It is also shown that the presence of a lake limits pollutant deposition over its length, due to stabilization and plume entrainment in the lake-breeze. On the hand, deposition was enhanced over the upwind land due to the lake-breeze effect, and the downwind shore due to intense mxing.     

Modeling of data center airflow and heat transfer: State of the art and future trends

An assessment of the current thermal modeling methodologies for data centers is presented, with focus on the use of computational fluid dynamics and heat transfer as analysis tools, and model validation. Future trends in reduced or compact modeling of data center airflow and heat transfer are presented to serve as an overview of integrating rack-level compact models into full-scale facility level numerical computations. Compact models can be used to efficiently model data centers through varying model fidelity across length scales. Dynamic effects can be included to develop next-generation control schemes to maximize data center energy efficiency. Recommended by: Monem Beitelmal     

Towards coupling a 3D hydrodynamic lake model with the Canadian Regional Climate Model: Simulation on Great Slave Lake

Recently, it has been recognized that large lakes exert considerable influence on regional climate systems and vice versa and that the Canadian Regional Climate Model (CRCM), which does not currently have a lake component, requires the development of a coupled lake sub-model. Prior to a full effort for this model development, however, studies are needed to select and assess the suitability of a lake hydrodynamic model in terms of its capability to couple with the CRCM. This paper evaluates the performance of the 3-dimensional hydrodynamic model ELCOM on Great Slave Lake, one of Canada's largest lakes in the northern climatic system. Model simulations showed dominant circulation patterns that can create relatively large spatial and temporal gradients in temperature. Simulated temperatures compared well with cross-lake temperature observations both at the surface and vertically. Sensitivity analysis was applied to determine the critical meteorological variables affecting simulations of temperature and surface heat fluxes. For example, a 10% increase in air temperature and solar radiation was found to result in a 3.1% and 8.3% increase in water surface temperature and 8.5% increase in latent heat flux. Knowledge of the model sensitivity is crucial for future research in which the hydrodynamic model coupled with the atmosphere will be forced from the CRCM output.     

Coupled model of surface water flow, sediment transport and morphological evolution

This paper presents a mathematical model coupling water flow and sediment transport dynamics that enables calculating the changing surface morphology through time and space. The model is based on the shallow water equations for flow, conservation of sediment concentration, and empirical functions for bed friction, substrate erosion and deposition. The sediment transport model is a non-capacity formulation whereby erosion and deposition are treated independently and influence the sediment flux by exchanging mass across the bottom boundary of the flow. The resulting hyperbolic system is solved using a finite volume, Godunov-type method with a first-order approximate Riemann solver. The model can be applied both to short time scales, where the flow, sediment transport and morphological evolution are strongly coupled and the rate of bed evolution is comparable to the rate of flow evolution, or to relatively long time scales, where the time scale of bed evolution associated with erosion and/or deposition is slow relative to the response of the flow to the changing surface and, therefore, the classical quasi-steady approximation can be invoked. The model is verified by comparing computed results with documented solutions. The developed model can be used to investigate a variety of problems involving coupled flow and sediment transport including channel initiation and drainage basin evolution associated with overland flow and morphological changes induced by extreme events such as tsunami.     

Quantifying sources of error in multitemporal multisensor lake mapping

Regional- to global-scale lake maps can now be produced using existing technology and freely available data and serve as powerful tools for a variety of lake- and water-related studies. The accuracy of these studies depends in part on the accuracy of the lake map that they use. Mapping lakes using remote sensing requires a careful study of error and uncertainty. Errors in lake maps are caused by sensor-specific, lake-specific, and processing-specific factors. These can be further broken down to spatial, spectral/radiometric, and temporal factors. In this study, we analyse and compare these factors using modern and historical Landsat images along with intensive ground surveys of lakes in northern Alaska. Percentage error in lake area (relative to lake size) decreases for larger and more circular lakes, making a minimum size threshold an effective error mitigation practice. Image resampling involved in image transformation significantly increased error in lake area and is easily avoided by performing co-registration in the vector domain. Spectral properties varied for individual lakes due to depth, suspended sediment, vegetation, and other in situ factors, necessitating a normalized water index and independently derived threshold values for each lake. For lake change detection studies, spatially degrading a finer resolution image to the resolution of the coarser image (a common practice) does not significantly affect the difference in observed lake area. Due to the large numbers of lakes, particularly in the climatologically sensitive Arctic region, small errors in individual lake areas can compound to significantly impact results on regional to global scales. This study is intended to inform future static and multitemporal lake remote-sensing studies by evaluating errors and uncertainties in lake area, as measured by remote sensing.     

Automated calculation of surface energy fluxes with high-frequency lake buoy data

Lake Heat Flux Analyzer is a program used for calculating the surface energy fluxes in lakes according to established literature methodologies. The program was developed in MATLAB for the rapid analysis of high-frequency data from instrumented lake buoys in support of the emerging field of aquatic sensor network science. To calculate the surface energy fluxes, the program requires a number of input variables, such as air and water temperature, relative humidity, wind speed, and short-wave radiation. Available outputs for Lake Heat Flux Analyzer include the surface fluxes of momentum, sensible heat and latent heat and their corresponding transfer coefficients, incoming and outgoing long-wave radiation. Lake Heat Flux Analyzer is open source and can be used to process data from multiple lakes rapidly. It provides a means of calculating the surface fluxes using a consistent method, thereby facilitating global comparisons of high-frequency data from lake buoys.     

基于WebGIS的湖北省湖泊时空演变系统设计与实现

位于长江中游的湖北省,分布有丰富广阔的湖泊资源,有着"千湖之省"的美誉.近年来水系自然演变和人类活动的影响使得湖泊面积变化显著.为充分了解湖北省湖泊变迁,文章利用WebGIS的空间信息处理功能、数据库等关键技术,设计建立湖北省湖泊信息数据库,实现对湖北省湖泊基本情况的时空演变的展示.系统可提供对湖北省湖泊信息的存储管理...     

Modeling research on the sorption kinetics of pentachlorophenol (PCP) to sediments based on neural networks and neuro-fuzzy systems

The sorption kinetics of pentachlorophenol (PCP) to sediments (from 8 different lakes, south of China) was studied in batch experiments. The comparison of a radial basis function neural network (RBFN) and an adaptive neuro-fuzzy inference system (ANFIS) applied for modeling the sorption behaviors was presented. Although the physical and chemical characteristics were different, the modeling results showed that the sorption behaviors of 8 different sediments were similar. Both RBFN and ANFIS could model the sorption behaviors and make predictions in high accuracy, which illustrated that the two models reflected the internal principle of the sorption better than the traditional model Fick's second law (FSL). Especially, RBFN held the promise of being able to work under noisy conditions to obtain high accuracy. In conclusion, RBFN was the valid options for modeling the sorption kinetics of PCP to the lake sediments. The individual changes of the three different inputs (the concentration of PCP in the aqueous phase, the reversible fraction and the irreversible fraction) affected the modeling results to the similar extent, from which we could infer that the sorption kinetics of PCP to these sediments were affected by all the three factors rather than only by the concentration of PCP in the aqueous phase.     

Influence of lake morphology and clarity on water surface temperature as measured by EOS ASTER

Relationships between lake morphometric parameters and nighttime lake surface temperatures were investigated in North American temperate lakes using the ASTER kinetic temperature (AST08) product. Nighttime ASTER kinetic temperature measurements were found to be a good analogue for nighttime surface temperatures. Linear regression between ASTER and buoy-measured temperatures in a test lake were better during the evening (R² = 0.98) than the day (R² = 0.90), presumably due to the greater influence of radiation and latent heat fluxes during daylight hours. Nighttime lake surface temperatures measured in three ASTER scenes were significantly correlated to logarithm of lake area, maximum lake depth, Secchi depth (a measure of lake clarity) and lake order (a measure of lake connection with surface drainage), during October and November. Nighttime lake surface temperatures were significantly correlated only with lake area in July. We hypothesize that morphology was more strongly related to surface temperature in the fall months due to lake turnover during that season. This study suggests that satellite derived thermal data may be useful for calculation of lake heat budgets and evaporation rates, provided surface temperatures are measured in well-mixed lakes.     

Thin film evaporation in microchannels with interfacial slip

We investigate the role of interfacial slip on evaporation of a thin liquid film in a microfluidic channel. The effective slip mechanism is attributed to the formation of a depleted layer adhering to the substrate–fluid interface, either in a continuum or in a rarefied gas regime, as a consequence of intricate hydrophobic interactions in the narrow confinement. We appeal to the fundamental principles of conservation in relating the evaporation mechanisms with fluid flow and heat transfer over interfacial scales. We obtain semi-analytical solutions of the pertinent governing equations, with coupled heat and mass transfer boundary conditions at the liquid–vapor interface. We observe that a general consequence of interfacial slip is to elongate the liquid film, thereby leading to a film thickening effect. Thicker liquid films, in turn, result in lower heat transfer rates from the wall to liquid film, and consequently lower mass transfer rates from the liquid film to the vapor phase. Nevertheless, the total mass of evaporation (or equivalently, the net heat transfer) turns out to be higher in case of interfacial slip due to the longer film length. We also develop significant physical insights on the implications of the relative thickness of the depleted layer with reference to characteristic length scales of the microfluidic channel on the evaporation process, under combined influences of the capillary pressure, disjoining pressure, and the driving temperature differential for the interfacial transport.     

A reduced-order model for whole-chip thermal analysis of microfluidic lab-on-a-chip systems

This paper presents a Krylov subspace projection-based reduced-order model (ROM) for whole microfluidic chip thermal analysis, including conjugate heat transfer. Two key steps in the reduced-order modeling procedure are described in detail: (1) the acquisition of a 3D full-scale computational model in the state-space form to capture the dynamic thermal behavior of the entire microfluidic chip; and (2) the model order reduction using the block Arnoldi algorithm to markedly lower the dimension of the full-scale model. Case studies using practically relevant thermal microfluidic chip are undertaken to establish the capability and to evaluate the computational performance of the reduced-order modeling technique. The ROM is compared against the full-scale model and exhibits good agreement in spatiotemporal thermal profiles (<0.5 % relative error in pertinent time scales) and over three-orders-of-magnitude acceleration in computational speed. The salient model reusability and real-time simulation capability render it amenable for operational optimization and in-line thermal control and management of microfluidic systems and devices.     

Multi-time-scale heat transfer modeling of turbid tissues exposed to short-pulsed irradiations

A combined hyperbolic radiation and conduction heat transfer model is developed to simulate multi-time-scale heat transfer in turbid tissues exposed to short-pulsed irradiations. An initial temperature response of a tissue to an ultrashort pulse irradiation is analyzed by the volume-average method in combination with the transient discrete ordinates method for modeling the ultrafast radiation heat transfer. This response is found to reach pseudo steady state within 1 ns for the considered tissues. The single pulse result is then utilized to obtain the temperature response to pulse train irradiation at the microsecond/millisecond time scales. After that, the temperature field is predicted by the hyperbolic heat conduction model which is solved by the MacCormack's scheme with error terms correction. Finally, the hyperbolic conduction is compared with the traditional parabolic heat diffusion model. It is found that the maximum local temperatures are larger in the hyperbolic prediction than the parabolic prediction. In the modeled dermis tissue, a 7% non-dimensional temperature increase is found. After about 10 thermal relaxation times, thermal waves fade away and the predictions between the hyperbolic and parabolic models are consistent.     

Remote sensing of the colour and temperature of volcanic lakes

Crater lakes on active volcanoes act as heat and chemical traps, and are amenable to surveillance from space. By use of all seven spectral bands, the Landsat Thematic Mapper (TM) can simultaneously measure: (i) lake surface temperature, (ii) lake surface area, and (iii) lake colour which is related to water chemistry. Use of the total instrument in this way enhances its utility for volcano surveillance. This work examines TM data for crater lakes at the following volcanoes: Ruapehu (New Zealand), Taal (Philippines), Kawah Ijen and Kelut (Indonesia), Poas (Costa Rica), and Apoyeque and Jiloa (Nicaragua). Observatory data indicate that lake surface temperatures derived by TM band 6 are typically 1-4°C less than contemporaneously measured bulk temperatures, probably due to the skin effect, the difference between water bulk and surface temperature. For Ruapehu, TM band 6 hotspots coincide approximately with known upwelling sites above volcanic vents on the lakebed. Field observations at Kawah Ijen show that the skin effect (< 3°C) is strongly correlated with windspeed: wind gusts peaking at ≈ 5m s^?1 caused rapid decreases in surface temperature of ≈ 0.5-1.0°C. These fluctuations are small compared with the magnitude of volcanogenic changes in lake temperature and do not reduce the utility of infrared surveillance. TM-derived water surface spectral reflectances indicate high concentrations of suspended chemical sediment in the most active crater lakes: Ruapehu, Poás and Ijen. For Ruapehu, imaged on two dates, the later scene reveals an upwelling slick, bright in bands 5 and 7, possibly composed of hollow sulphur spherules emitted from a subaqueous vent. The Advanced Spaceborne Thermal Emission and Reflectance Radiometer (ASTER) and Enhanced Thematic Mapper (ETM +) both due for launch in 1998, will offer improved capabilities for remote surveillance of crater lakes.     

Correcting for the influence of frozen lakes in satellite microwave radiometer observations through application of a microwave emission model

The spatial resolution of passive microwave observations from space is of the order of tens of kilometers with currently available instruments, such as the Special Sensor Microwave/Imager (SSM/I) and Advanced Microwave Scanning Radiometer (AMSR-E). The large field of view of these instruments dictates that the observed brightness temperature can originate from heterogeneous land cover, with different vegetation and surface properties.In this study, we assess the influence of freshwater lakes on the observed brightness temperature of AMSR-E in winter conditions. The study focuses on the geographic region of Finland, where lakes account for 10% of the total terrestrial area. We present a method to mitigate for the influence of lakes through forward modeling of snow covered lakes, as a part of a microwave emission simulation scheme of space-borne observations. We apply a forward model to predict brightness temperatures of snow covered sceneries over several winter seasons, using available data on snow cover, vegetation and lake ice cover to set the forward model input parameters. Comparison of model estimates with space-borne observations shows that the modeling accuracy improves in the majority of examined cases when lakes are accounted for, with respect to the case where lakes are not included in the simulation. Moreover, we present a method for applying the correction to the retrieval of Snow Water Equivalent (SWE) in lake-rich areas, using a numerical inversion method of the forward model. In a comparison to available independent validation data on SWE, also the retrieval accuracy is seen to improve when applying the influence of snow covered lakes in the emission model.     

Understanding fire and smoke flow through modeling and visualization

Computer modeling and visualization are important tools for understanding the processes of fire behavior. Fire models range in complexity from simple correlations for predicting quantities such as flame heights or flow velocities to moderately complex zone fire models for predicting time-dependent smoke layer temperatures and heights. Zone fire model calculations can run on today's computers within minutes because they solve only four differential equations per room. Zone models approximate the entire upper layer with just one temperature. This approximation works remarkably well but breaks down for complicated flows or geometries. For such cases, computational fluid dynamics (CFD) techniques are required.     

混合记录气膜温度场仿真

混合记录是实现超高密度存储的重要途径之一,是当前国际超高密度存储研究的热点。该技术利用高矫顽力的磁记录材料进行热辅助写入,并在常温下读出和保存,可克服超顺磁极限实现超高密度存储。基于对流换热理论建立了混合记录介质与浮动气膜间的对流换热模型。磁头在不同的飞行姿态下分别可以用层流换热模型和紊流换热模型来求解气膜温度场。注意到磁记录磁头飞高在纳米量级的特点,运用边界层理论求解了浮动气膜的温度场,得到了近场光、磁混合记录中由近场馈热导致气膜温升的影响规律及相关的数据,     

Numerical Simulation of growth and collapse of a bubble induced by a pulsed microheater

A three-dimensional numerical analysis of the growth and collapse of a bubble on a microheater is presented. SIMULENT code, which solves the full Navier-Stokes equations with surface tension effects, is used in these simulations. A volume of fluid (VOF) interface tracking algorithm is used to track the evolution of the free surface flow. A one-dimensional heat conduction model is used to consider the energy transfer between the bubble and the surrounding liquid, as well as the temperature distribution in the liquid layer. Details of the velocity and pressure distribution in the liquid during the growth and collapse of the vapor bubble are obtained. Numerical results for the growth and the collapse of the bubbles are compared with those of experiments under similar conditions. Comparisons show that the volume evolution of the vapor bubble is well predicted by the numerical model.     

Pseudo bond graph model of coupled heat and mass transfers in a plastic tunnel greenhouse

In greenhouses, models built are classified into two kinds; heterogeneous approaches based on computational fluid dynamics (CFD) codes and homogeneous one based on heat and mass balance equations. Bond graph modelling was not seriously introduced in greenhouse modelling in spite of the concordance of the energetic based approach of bond graphs with the nature of the greenhouse plant.In this research work, a pseudo bond graph model of a greenhouse was elaborated to simulate temperature and relative humidity inside. The model proposed is an energetic lumped approach which describes coupled heat and mass transfers in a plastic tunnel greenhouse. This model includes convection, evaporation/condensation phenomena, air change flow and soil heat and mass transfer. Multiport (3 ports) bond graph elements are introduced to describe the state of the two elements (dry air and water vapor) fluid. New bond graph schemes are used to characterize the coupling effect between heat and mass transfer, and for modelling free evaporation/condensation mechanism. Big leaf assumption was used to model the canopy. New boundary layer elements are added in the model, these elements allow a separation of the different phenomena inside greenhouse and thus a simplification of the modelling task. The practical results obtained from an experimental tunnel greenhouse are used here as validation elements for the greenhouse bond graph model. A good correlation is observed between measured and predicted samples.     

Implementation of a discontinuous Galerkin morphological model on two-dimensional unstructured meshes

The shallow water equations are used to model large-scale surface flow in the ocean, coastal rivers, estuaries, salt marshes, bays, and channels. They can describe tidal flows as well as storm surges associated with extreme storm events, such as hurricanes. The resulting currents can transport bed load and suspended sediment and result in morphological changes to the seabed. Modeling these processes requires tightly coupling a bed morphology equation to the shallow water equations. Discontinuous Galerkin finite element methods are a natural choice for modeling this coupled system, given the need to solve these problems on unstructured computational meshes, as well as the desire to implement hp-adaptivity for capturing the dynamic features of the solution. However, because of the presence of non-conservative products in the momentum equations, the standard DG method cannot be applied in a straightforward manner. To rectify this situation, we summarize and follow an extended approach described by Rhebergen et al., which uses theoretical results due to Dal Maso et al. appearing in earlier work. In this paper, we focus on aspects of the implementation of the morphological model for bed evolution within the Advanced Circulation (ADCIRC) modeling framework, as well as the verification of the RKDG method in both h (mesh spacing) and p (polynomial order). This morphological model is applied to a number of coastal engineering problems, and numerical results are presented, with attention paid to the effects of h- and p-refinement in these applications. In particular, it is observed that for sediment transport, piecewise constant (i.e., finite volume) approximations of the bed are very over-diffusive and lead to poor sediment solutions.