Inter-comparison and performance evaluation of chemistry transport models over Indian region

Aerosol loading over the South Asian region has the potential to affect the monsoon rainfall, Himalayan glaciers and regional air-quality, with implications for the billions in this region. While field campaigns and network observations provide primary data, they tend to be location/season specific. Numerical models are useful to regionalize such location-specific data. Studies have shown that numerical models underestimate the aerosol scenario over the Indian region, mainly due to shortcomings related to meteorology and the emission inventories used. In this context, we have evaluated the performance of two such chemistry-transport models: WRF-Chem and SPRINTARS over an India-centric domain. The models differ in many aspects including physical domain, horizontal resolution, meteorological forcing and so on etc. Despite these differences, both the models simulated similar spatial patterns of Black Carbon (BC) mass concentration, (with a spatial correlation of 0.9 with each other), and a reasonable estimates of its concentration, though both of them under-estimated vis-a-vis the observations. While the emissions are lower (higher) in SPRINTARS (WRF-Chem), overestimation of wind parameters in WRF-Chem caused the concentration to be similar in both models. Additionally, we quantified the under-estimations of anthropogenic BC emissions in the inventories used these two models and three other widely used emission inventories. Our analysis indicates that all these emission inventories underestimate the emissions of BC over India by a factor that ranges from 1.5 to 2.9. We have also studied the model simulations of aerosol optical depth over the Indian region. The models differ significantly in simulations of AOD, with WRF-Chem having a better agreement with satellite observations of AOD as far as the spatial pattern is concerned. It is important to note that in addition to BC, dust can also contribute significantly to AOD. The models differ in simulations of the spatial pattern of mineral dust over the Indian region. We find that both meteorological forcing and emission formulation contribute to these differences. Since AOD is column integrated parameter, description of vertical profiles in both models, especially since elevated aerosol layers are often observed over Indian region, could be also a contributing factor. Additionally, differences in the prescription of the optical properties of BC between the models appear to affect the AOD simulations. We also compared simulation of sea-salt concentration in the two models and found that WRF-Chem underestimated its concentration vis-a-vis SPRINTARS. The differences in near-surface oceanic wind speeds appear to be the main source of this difference. In-spite of these differences, we note that there are similarities in their simulation of spatial patterns of various aerosol species (with each other and with observations) and hence models could be valuable tools for aerosol-related studies over the Indian region. Better estimation of emission inventories could improve aerosol-related simulations. (C) 2015 Elsevier Ltd. All rights reserved.

Transmission loss analysis of rectangular expansion chamber with arbitrary location of inlet/outlet by means of Green's functions

Transmission loss of a rectangular expansion chamber, the inlet and outlet of which are situated at arbitrary locations of the chamber, i.e., the side wall or the face of the chamber, are analyzed here based on the Green's function of a rectangular cavity with homogeneous boundary conditions. The rectangular chamber Green's function is expressed in terms of a finite number of rigid rectangular cavity mode shapes. The inlet and outlet ports are modeled as uniform velocity pistons. If the size of the piston is small compared to wavelength, then the plane wave excitation is a valid assumption. The velocity potential inside the chamber is expressed by superimposing the velocity potentials of two different configurations. The first configuration is a piston source at the inlet port and a rigid termination at the outlet, and the second one is a piston at the outlet with a rigid termination at the inlet. Pressure inside the chamber is derived from velocity potentials using linear momentum equation. The average pressure acting on the pistons at the inlet and outlet locations is estimated by integrating the acoustic pressure over the piston area in the two constituent configurations. The transfer matrix is derived from the average pressure values and thence the transmission loss is calculated. The results are verified against those in the literature where use has been made of modal expansions and also numerical models (FEM fluid). The transfer matrix formulation for yielding wall rectangular chambers has been derived incorporating the structural–acoustic coupling. Parametric studies are conducted for different inlet and outlet configurations, and the various phenomena occurring in the TL curves that cannot be explained by the classical plane wave theory, are discussed.

Orographic effects on the southwest monsoon: A review

An overview of the problem of orographic effects on the southwest monsoon using the contributions of all the available analytical and numerical models is attempted. A quasi-geostrophic model is applied to deduce the effect of the topographic complex on the Indian peninsula. This model suggests that the southward bending of the low-level isobars on the peninsula can be ascribed to the topographically-induced southward velocity. This southward velocity triggers a Rossby wave to the east of the peninsula which is manifested as a trough on the southern Bay of Bengal.

Studies on design, fabrication and reliability assessment of embedded passives on a high-density interconnect (HDI) organic substrate using a sequential build-up process

One of the foremost design considerations in microelectronics miniaturization is the use of embedded passives which provide practical solution. In a typical circuit, over 80 percent of the electronic components are passives such as resistors, inductors, and capacitors that could take up to almost 50 percent of the entire printed circuit board area. By integrating passive components within the substrate instead of being on the surface, embedded passives reduce the system real estate, eliminate the need for discrete and assembly, enhance electrical performance and reliability, and potentially reduce the overall cost. Moreover, it is lead free. Even with these advantages, embedded passive technology is at a relatively immature stage and more characterization and optimization are needed for practical applications leading to its commercialization.This paper presents an entire process from design and fabrication to electrical characterization and reliability test of embedded passives on multilayered microvia organic substrate. Two test vehicles focusing on resistors and capacitors have been designed and fabricated. Embedded capacitors in this study are made with polymer/ceramic nanocomposite (BaTiO3) material to take advantage of low processing temperature of polymers and relatively high dielectric constant of ceramics and the values of these capacitors range from 50 pF to 1.5 nF with capacitance per area of approximately 1.5 nF/cm(2). Limited high frequency measurement of these capacitors was performed. Furthermore, reliability assessments of thermal shock and temperature humidity tests based on JEDEC standards were carried out. Resistors used in this work have been of three types: 1) carbon ink based polymer thick film (PTF), 2) resistor foils with known sheet resistivities which are laminated to printed wiring board (PWB) during a sequential build-up (SBU) process and 3) thin-film resistor plating by electroless method. Realization of embedded resistors on conventional board-level high-loss epoxy (similar to 0.015 at 1 GHz) and proposed low-loss BCB dielectric (similar to 0.0008 at > 40 GHz) has been explored in this study. Ni-P and Ni-W-P alloys were plated using conventional electroless plating, and NiCr and NiCrAlSi foils were used for the foil transfer process. For the first time, Benzocyclobutene (BCB) has been proposed as a board level dielectric for advanced System-on-Package (SOP) module primarily due to its attractive low-loss (for RF application) and thin film (for high density wiring) properties.Although embedded passives are more reliable by eliminating solder joint interconnects, they also introduce other concerns such as cracks, delamination and component instability. More layers may be needed to accommodate the embedded passives, and various materials within the substrate may cause significant thermo -mechanical stress due to coefficient of thermal expansion (CTE) mismatch. In this work, numerical models of embedded capacitors have been developed to qualitatively examine the effects of process conditions and electrical performance due to thermo-mechanical deformations.Also, a prototype working product with the board level design including features of embedded resistors and capacitors are underway. Preliminary results of these are presented.

Analytical prediction of break-out noise from a reactive rectangular plenum with four flexible walls

This paper describes an analytical calculation of break-out noise from a rectangular plenum with four flexible walls by incorporating three-dimensional effects along with the acoustical and structural wave coupling phenomena. The breakout noise from rectangular plenums is important and the coupling between acoustic waves within the plenum and structural waves in the flexible plenum walls plays a critical role in prediction of the transverse transmission loss. The first step in breakout noise prediction is to calculate the inside plenum pressure field and the normal flexible plenum wall vibration by using an impedance-mobility approach, which results in a compact matrix formulation. In the impedance-mobility compact matrix (IMCM) approach, it is presumed that the coupled response can be described in terms of finite sets of the uncoupled acoustic subsystem and the structural subsystem. The flexible walls of the plenum are modeled as an unfolded plate to calculate natural frequencies and mode shapes of the uncoupled structural subsystem. The second step is to calculate the radiated sound power from the flexible walls using Kirchhoff-Helmholtz (KH) integral formulation. Analytical results are validated with finite element and boundary element (FEM-BEM) numerical models. (C) 2010 Acoustical Society of America. DOI: 10.1121/1.3463801]

Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

Accelerating frequency-domain diffuse optical tomographic image reconstruction using graphics processing units

Diffuse optical tomographic image reconstruction uses advanced numerical models that are computationally costly to be implemented in the real time. The graphics processing units (GPUs) offer desktop massive parallelization that can accelerate these computations. An open-source GPU-accelerated linear algebra library package is used to compute the most intensive matrix-matrix calculations and matrix decompositions that are used in solving the system of linear equations. These open-source functions were integrated into the existing frequency-domain diffuse optical image reconstruction algorithms to evaluate the acceleration capability of the GPUs (NVIDIA Tesla C 1060) with increasing reconstruction problem sizes. These studies indicate that single precision computations are sufficient for diffuse optical tomographic image reconstruction. The acceleration per iteration can be up to 40, using GPUs compared to traditional CPUs in case of three-dimensional reconstruction, where the reconstruction problem is more underdetermined, making the GPUs more attractive in the clinical settings. The current limitation of these GPUs in the available onboard memory (4 GB) that restricts the reconstruction of a large set of optical parameters, more than 13, 377. (C) 2010 Society of Photo-Optical Instrumentation Engineers. DOI: 10.1117/1.3506216]

Particle filters for structural system identification using multiple test and sensor data: A combined computational and experimental study

The problem of structural system identification when measurements originate from multiple tests and multiple sensors is considered. An offline solution to this problem using bootstrap particle filtering is proposed. The central idea of the proposed method is the introduction of a dummy independent variable that allows for simultaneous assimilation of multiple measurements in a sequential manner. The method can treat linear/nonlinear structural models and allows for measurements on strains and displacements under static/dynamic loads. Illustrative examples consider measurement data from numerical models and also from laboratory experiments. The results from the proposed method are compared with those from a Kalman filter-based approach and the superior performance of the proposed method is demonstrated. Copyright (C) 2009 John Wiley & Sons, Ltd.

Nonaxisymmetric convection in stationary Gas Tungsten Arc Weld Pools

Observations of surface flow patterns of steel and aluminum GTAW pools have been made using a pulsed laser visualization system. The weld pool convection is found to be three-dimensional with the azimuthal circulation depending on the location of the clamp with respect to the torch. Oscillation of steel pools and undulating motion in aluminum weld pools are also observed even with steady process parameters. Current axisymmetric numerical models are unable to explain such phenomena. A three-dimensional computational study is carried out in this study to explain the rotational flow in aluminum weld pools.

Prediction of Breakout Noise from a Rectangular Duct with Compliant Walls

Breakout noise from HVAC ducts is important at low frequencies, and the coupling between the acoustic waves and the structural waves plays a critical role in the prediction of the transverse transmission loss. This paper describes the analytical calculation of breakout noise by incorporating three-dimensional effects along with the acoustical and structural wave coupling phenomena. The first step in the breakout noise prediction is to calculate the inside duct pressure field and the normal duct wall vibration by using the solution of the governing differential equations in terms of Green's function. The resultant equations are rearranged in terms of impedance and mobility, which results in a compact matrix formulation. The Green's function selected for the current problem is the cavity Green's function with modification of wave number in the longitudinal direction in order to incorporate the terminal impedance. The second step is to calculate the radiated sound power from the compliant duct walls by means of an ``equivalent unfolded plate'' model. The transverse transmission loss from the duct walls is calculated using the ratio of the incident power due to surface source inside the duct to the acoustic power radiated from the compliant duct walls. Analytical results are validated with the FE-BE numerical models.

Macromolecular Transport Through Porous Arterial Walls

Arteries are heterogeneous, composite structures that undergo large cyclic deformations during blood transport. Presence, build-up and consequent rupture of blockages in blood vessels, called atherosclerotic plaques, lead to disruption in the blood flow that can eventually be fatal. Abnormal lipid profile and hypertension are the main risk factors for plaque progression. Treatments span from pharmacological methods, to minimally invasive balloon angioplasty and stent procedures, and finally to surgical alternatives. There is a need to understand arterial disease progression and devise methods to detect, control, treat and manage arterial disease through early intervention. Local delivery through drug eluting stents also provide an attractive option for maintaining vessel integrity and restoring blood flow while releasing controlled amount of drug to reduce and alleviate symptoms. Development of drug eluting stents is hence interesting albeit challenging because it requires an integration of knowledge of mechanical properties with material transport of drug through the arterial wall to produce a desired biochemical effect. Although experimental models are useful in studying such complex multivariate phenomena, numerical models of mass transport in the vessel have proved immensely useful to understand and delineate complex interactions between chemical species, physical parameters and biological variables. The goals of this review are to summarize literature based on studies of mass transport involving low density lipoproteins in the arterial wall. We also discuss numerical models of drug elution from stents in layered and porous arterial walls that provide a unique platform that can be exploited for the design of novel drug eluting stents.

Characteristics of Arabian Sea mini warm pool and Indian summer monsoon

Arabian Sea Mini Warm Pool (ASMWP) is a part of the Indian Ocean Warm Pool and formed in the eastern Arabian Sea prior to the onset of the summer monsoon season. This warm pool attained its maximum intensity during the pre-monsoon season and dissipated with the commencement of summer monsoon. The main focus of the present work was on the triggering of the dissipation of this warm pool and its relation to the onset of summer monsoon over Kerala. This phenomenon was studied utilizing NCEP/NCAR (National Center for Environmental Prediction/National Center for Atmospheric and Research) re-analysis data, TRMM Micro wave Imager (TMI) and observational data. To define the ASMWP, sea surface temperature exceeding 30.25A degrees C was taken as the criteria. The warm pool attained its maximum dimension and intensity nearly 2 weeks prior to the onset of summer monsoon over Kerala. Interestingly, the warm pool started its dissipation immediately after attaining its maximum core temperature. This information can be included in the present numerical models to enhance the prediction capability. It was also found that the extent and intensity of the ASMWP varied depending on the type of monsoon i.e., excess, normal, and deficient monsoon. Maximum core temperature and wide coverage of the warm pool observed during the excess monsoon years compared to normal and deficient monsoon years. The study also revealed a strong relationship between the salinity in the eastern Arabian Sea and the nature of the monsoon.

Geothermal reservoirs - A brief review

A brief discussion and review of the geothermal reservoir systems, geothermal energy and modeling and simulation of the geothermal reservoirs has been presented here. Different types of geothermal reservoirs and their governing equations have been discussed first. The conceptual and numerical modeling along with the representation of flow though fractured media, some issues related to non isothermal flow through fractured media, the efficiency of the geothermal reservoir, structure of the numerical models, boundary conditions and calibration procedures have been illustrated. A brief picture of the Indian scenario and some barriers related with geothermal power are discussed and presented thereafter. Finally some gaps of the existing knowledge and recent focuses of research are discussed.