Toward assessing the effects of aerosols on deep convection: a numerical study using the WRF-Chemistry model

As the formative agents of cloud droplets, aerosols play an undeniably important role in the development of clouds and precipitation. Few meteorological models have been developed or adapted to simulate aerosols and their contribution to cloud and precipitation processes. The Weather Research and Forecasting model (WRF) has recently been coupled with an atmospheric chemistry suite and is jointly referred to as WRF-Chem, allowing atmospheric chemistry and meteorology to influence each other’s evolution within a mesoscale modeling framework. Provided that the model physics are robust, this framework allows the feedbacks between aerosol chemistry, cloud physics, and dynamics to be investigated. This study focuses on the effects of aerosols on meteorology, specifically, the interaction of aerosol chemical species with microphysical processes represented within the framework of the WRF-Chem. Aerosols are represented by eight size bins using the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) sectional parameterization, which is linked to the Purdue Lin bulk microphysics scheme. The aim of this study is to examine the sensitivity of deep convective precipitation modeled by the 2D WRF-Chem to varying aerosol number concentration and aerosol type. A systematic study has been performed regarding the effects of aerosols on parameters such as total precipitation, updraft/downdraft speed, distribution of hydrometeor species, and organizational features, within idealized maritime and continental thermodynamic environments. Initial results were obtained using WRFv3.0.1, and a second series of tests were run using WRFv3.2 after several changes to the activation, autoconversion, and Lin et al. microphysics schemes added by the WRF community, as well as the implementation of prescribed vertical levels by the author. The results of WRFv3.2 runs contrasted starkly with WRFv3.0.1 runs. The WRFv3.0.1 runs produced a propagating system resembling a developing squall line, whereas the WRFv3.2 runs did not. The response of total precipitation, updraft/downdraft speeds, and system organization to increasing aerosol concentrations were opposite between runs with different versions of WRF. Results of the WRFv3.2 runs, however, were in better agreement in timing and magnitude of vertical velocity and hydrometeor content with a WRFv3.0.1 run using single-moment Lin et al. microphysics, than WRFv3.0.1 runs with chemistry. One result consistent throughout all simulations was an inhibition in warm-rain processes due to enhanced aerosol concentrations, which resulted in a delay of precipitation onset that ranged from 2-3 minutes in WRFv3.2 runs, and up to 15 minutes in WRFv.3.0.1 runs. This result was not observed in a previous study by Ntelekos et al. (2009) using the WRF-Chem, perhaps due to their use of coarser horizontal and vertical resolution within their experiment. The changes to microphysical processes such as activation and autoconversion from WRFv3.0.1 to WRFv3.2, along with changes in the packing of vertical levels, had more impact than the varying aerosol concentrations even though the range of aerosol tested was greater than that observed in field studies. In order to take full advantage of the input of aerosols now offered by the chemistry module in WRF, the author recommends that a fully double-moment microphysics scheme be linked, rather than the limited double-moment Lin et al. scheme that currently exists. With this modification, the WRF-Chem will be a powerful tool for studying aerosol-cloud interactions and allow comparison of results with other studies using more modern and complex microphysical parameterizations.

Physical modeling of processes of sediment transport in inlets

During the transformation of the low tide to the high tide, an exactly inverse phenomenon is occurred and the high tidal delta is formed at the mouth upstream. Increasing the tidal range does not affect the nature of this phenomenon and just change its intensity. In this situation, the inlet will be balance over time. A new relationship between equilibrium cross section and tidal prism for different tidal levels as well as sediment grading has been provided which its results are corresponded with results of numerical modeling. In the combination state, the wave height significantly affects the current and sedimentary pattern such that the wave height dimensionless index (Hw/Ht) determines the dominant parameter (the short period wave or tide) in the inlet. It is notable that in this state, the inlet will be balanced over the time. In order to calculate sedimentary phenomena, each of which are individually determined under solely wave and only tide conditions and then they are added. Estimated values are similar to numerical modeling results of the combination state considering nonlinear terms. Also, it is clear that the wave and tide performance is of meaning in the direct relationship with the water level. The water level change causes variations of the position of the breaking line and sedimentary active area. It changes the current and sedimentary pattern coastward while does not change anything seaward. Based on modeling results of sediment transport due to the wave, tide and their combination, it could be said that the erosion at the mouth due to the wave is less than that due to the wave and tide combination. In these situations, tide and wave-tide combination increase the low tidal and high tidal delta volume, respectively. Hence, tide plays an effective role in changing sedimentary phenomena at the channel and mouth downstream. Whereas, short period and combined waves have a crucial role in varying the morphology and sediment transport coast ward.

Geotechnical Data and Numerical Analysis of Edifice Collapse and Related Hazards at Pacaya Volcano, Guatemala

The continual eruptive activity, occurrence of an ancestral catastrophic collapse, and inherent geologic features of Pacaya volcano (Guatemala) demands an evaluation of potential collapse hazards. This thesis merges techniques in the field and laboratory for a better rock mass characterization of volcanic slopes and slope stability evaluation. New field geological, structural, rock mechanical and geotechnical data on Pacaya is reported and is integrated with laboratory tests to better define the physical-mechanical rock mass properties. Additionally, this data is used in numerical models for the quantitative evaluation of lateral instability of large sector collapses and shallow landslides. Regional tectonics and local structures indicate that the local stress regime is transtensional, with an ENE-WSW sigma 3 stress component. Aligned features trending NNW-SSE can be considered as an expression of this weakness zone that favors magma upwelling to the surface. Numerical modeling suggests that a large-scale collapse could be triggered by reasonable ranges of magma pressure (greater than or equal to 7.7 MPa if constant along a central dyke) and seismic acceleration (greater than or equal to 460 cm/s2), and that a layer of pyroclastic deposits beneath the edifice could have been a factor which controlled the ancestral collapse. Finally, the formation of shear cracks within zones of maximum shear strain could provide conduits for lateral flow, which would account for long lava flows erupted at lower elevations.

Groundwater flow model of the Logan river alluvial aquifer system Josephville, South East Queensland

The study focuses on an alluvial plain situated within a large meander of the Logan River at Josephville near Beaudesert which supports a factory that processes gelatine. The plant draws water from on site bores, as well as the Logan River, for its production processes and produces approximately 1.5 ML per day (Douglas Partners, 2004) of waste water containing high levels of dissolved ions. At present a series of treatment ponds are used to aerate the waste water reducing the level of organic matter; the water is then used to irrigate grazing land around the site. Within the study the hydrogeology is investigated, a conceptual groundwater model is produced and a numerical groundwater flow model is developed from this. On the site are several bores that access groundwater, plus a network of monitoring bores. Assessment of drilling logs shows the area is formed from a mixture of poorly sorted Quaternary alluvial sediments with a laterally continuous aquifer comprised of coarse sands and fine gravels that is in contact with the river. This aquifer occurs at a depth of between 11 and 15 metres and is overlain by a heterogeneous mixture of silts, sands and clays. The study investigates the degree of interaction between the river and the groundwater within the fluvially derived sediments for reasons of both environmental monitoring and sustainability of the potential local groundwater resource. A conceptual hydrogeological model of the site proposes two hydrostratigraphic units, a basal aquifer of coarse-grained materials overlain by a thick semi-confining unit of finer materials. From this, a two-layer groundwater flow model and hydraulic conductivity distribution was developed based on bore monitoring and rainfall data using MODFLOW (McDonald and Harbaugh, 1988) and PEST (Doherty, 2004) based on GMS 6.5 software (EMSI, 2008). A second model was also considered with the alluvium represented as a single hydrogeological unit. Both models were calibrated to steady state conditions and sensitivity analyses of the parameters has demonstrated that both models are very stable for changes in the range of ± 10% for all parameters and still reasonably stable for changes up to ± 20% with RMS errors in the model always less that 10%. The preferred two-layer model was found to give the more realistic representation of the site, where water level variations and the numerical modeling showed that the basal layer of coarse sands and fine gravels is hydraulically connected to the river and the upper layer comprising a poorly sorted mixture of silt-rich clays and sands of very low permeability limits infiltration from the surface to the lower layer. The paucity of historical data has limited the numerical modelling to a steady state one based on groundwater levels during a drought period and forecasts for varying hydrological conditions (e.g. short term as well as prolonged dry and wet conditions) cannot reasonably be made from such a model. If future modelling is to be undertaken it is necessary to establish a regular program of groundwater monitoring and maintain a long term database of water levels to enable a transient model to be developed at a later stage. This will require a valid monitoring network to be designed with additional bores required for adequate coverage of the hydrogeological conditions at the Josephville site. Further investigations would also be enhanced by undertaking pump testing to investigate hydrogeological properties in the aquifer.

Determination of elastic modulus of thin coating materials using nanoindentation and finite element analysis

A deconvolution method that combines nanoindentation and finite element analysis was developed to determine elastic modulus of thin coating layer in a coating-substrate bilayer system. In this method, the nanoindentation experiments were conducted to obtain the modulus of both the bilayer system and the substrate. The finite element analysis was then applied to deconvolve the elastic modulus of the coating. The results demonstrated that the elastic modulus obtained using the developed method was in good agreement with that reported in literature.

Meshless TL and UL approaches for large deformation analysis

To accurately and effectively simulate large deformation is one of the major challenges in numerical modeling of metal forming. In this paper, an adaptive local meshless formulation based on the meshless shape functions and the local weak-form is developed for the large deformation analysis. Total Lagrangian (TL) and the Updated Lagrangian (UL) approaches are used and thoroughly compared each other in computational efficiency and accuracy. It has been found that the developed meshless technique provides a superior performance to the conventional FEM in dealing with large deformation problems for metal forming. In addition, the TL has better computational efficiency than the UL. However, the adaptive analysis is much more efficient using the UL approach than using in the TL approach.

Anomalous subdiffusion equations by the meshless collocation method

This paper aims to develop an implicit meshless collocation technique based on the moving least squares approximation for numerical simulation of the anomalous subdiffusion equation(ASDE). The discrete system of equations is obtained by using the MLS meshless shape functions and the meshless collocation formulation. The stability and convergence of this meshless approach related to the time discretization are investigated theoretically and numerically. The numerical examples with regular and irregular nodal distributions are used to the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling of ASDEs.

Nonlinear analysis of plastic ball grid array solder joints

A nonlinear finite element analysis was carried out to investigate the viscoplastic deformation of solder joints in a ball grid array (BGA) package under temperature cycle. The effects of constraint on print circuit board (PCB) and stiffness of substrate on the deformation behaviour of the solder joints were also studied. A relative damage stress was adopted to analyze the potential failure sites in the solder joints. The results indicated that high inelastic strain and strain energy density were developed in the joints close to the package center. On the other hand, high constraint and high relative damage stress were associated with the joint closest to the edge of the silicon chip. The joint closest to the edge of the silicon chip was regarded as the most susceptible failure site if cavitation instability is the dominant failure mechanism. Increase the external constraint on the print circuit board (PCB) causes a slight increase in stress triaxiality (m/eq) and relative damage stress in the joint closest to the edge of silicon die. The relative damage stress is not sensitive to the Young’s modulus of the substrate.

Simulation of plant cell shrinkage during drying : A SPH–DEM approach

Plant based dried food products are popular commodities in global market where much research is focused to improve the products and processing techniques. In this regard, numerical modelling is highly applicable and in this work, a coupled meshfree particle-based two-dimensional (2-D) model was developed to simulate micro-scale deformations of plant cells during drying. Smoothed Particle Hydrodynamics (SPH) was used to model the viscous cell protoplasm (cell fluid) by approximating it to an incompressible Newtonian fluid. The visco-elastic characteristic of the cell wall was approximated to a Neo-Hookean solid material augmented with a viscous term and modelled with a Discrete Element Method (DEM). Compared to a previous work [H. C. P. Karunasena, W. Senadeera, Y. T. Gu and R. J. Brown, Appl. Math. Model., 2014], this study proposes three model improvements: linearly decreasing positive cell turgor pressure during drying, cell wall contraction forces and cell wall drying. The improvements made the model more comparable with experimental findings on dried cell morphology and geometric properties such as cell area, diameter, perimeter, roundness, elongation and compactness. This single cell model could be used as a building block for advanced tissue models which are highly applicable for product and process optimizations in Food Engineering.

A particle based model to simulate microscale morphological changes of plant tissues during drying

Fundamental understanding on microscopic physical changes of plant materials is vital to optimize product quality and processing techniques, particularly in food engineering. Although grid-based numerical modelling can assist in this regard, it becomes quite challenging to overcome the inherited complexities of these biological materials especially when such materials undergo critical processing conditions such as drying, where the cellular structure undergoes extreme deformations. In this context, a meshfree particle based model was developed which is fundamentally capable of handling extreme deformations of plant tissues during drying. The model is built by coupling a particle based meshfree technique: Smoothed Particle Hydrodynamics (SPH) and a Discrete Element Method (DEM). Plant cells were initiated as hexagons and aggregated to form a tissue which also accounts for the characteristics of the middle lamella. In each cell, SPH was used to model cell protoplasm and DEM was used to model the cell wall. Drying was incorporated by varying the moisture content, the turgor pressure, and cell wall contraction effects. Compared to the state of the art grid-based microscale plant tissue drying models, the proposed model can be used to simulate tissues under excessive moisture content reductions incorporating cell wall wrinkling. Also, compared to the state of the art SPH-DEM tissue models, the proposed model better replicates real tissues and the cell-cell interactions used ensure efficient computations. Model predictions showed good agreement both qualitatively and quantitatively with experimental findings on dried plant tissues. The proposed modelling approach is fundamentally flexible to study different cellular structures for their microscale morphological changes at dehydration.

Thin single-walled carbon nanotubes with narrow chirality distribution : constructive interplay of plasma and Gibbs–Thomson effects

Multiscale, multiphase numerical modeling is used to explain the mechanisms of effective control of chirality distributions of single-walled carbon nanotubes in direct plasma growth and suggest effective approaches to further improvement. The model includes an unprecedented combination of the plasma sheath, ion/radical transport, species creation/loss, plasma–surface interaction, heat transfer, surface/bulk diffusion, graphene layer nucleation, and bending/lift-off modules. It is shown that the constructive interplay between the plasma and the Gibbs–Thomson effect can lead to the effective nucleation and lift-off of small graphene layers on small metal catalyst nanoparticles. As a result, much thinner nanotubes with narrower chirality distributions can nucleate at much lower process temperatures and pressures compared to thermal CVD. This approach is validated by a host of experimental results, substantially reduces the amounts of energy and atomic matter required for the nanotube growth, and can be extended to other nanoscale structures and materials systems, thereby nearing the ultimate goal of energy- and matter-efficient nanotechnology.

Rapid, simultaneous activation of thin nanowire growth in low-temperature, low-pressure chemically active plasmas

Multiscale numerical modeling of the species balance and transport in the ionized gas phase and on the nanostructured solid surface complemented by the heat exchange model is used to demonstrate the possibility of minimizing the Gibbs-Thompson effect in low-temperature, low-pressure chemically active plasma-assisted growth of uniform arrays of very thin Si nanowires, impossible otherwise. It is shown that plasma-specific effects drastically shorten and decrease the dispersion of the incubation times for the nucleation of nanowires on non-uniform Au catalyst nanoparticle arrays. The fast nucleation makes it possible to avoid a common problem of small catalyst nanoparticle burying by amorphous silicon. These results explain a multitude of experimental observations on chemically active plasma-assisted Si nanowire growth and can be used for the synthesis of a range of inorganic nanowires for environmental, biomedical, energy conversion, and optoelectronic applications.

Nanopowder management and control of plasma parameters in electronegative SiH4 plasmas

The nanopowder management and control of plasma parameters in electronegative SiH4 plasmas were discussed. The spatial profiles of electron and positive/negative ion number densities, electron temperature and charge of the fine particles were obtained. It was found that management of powder charge distribution is also possible through control of the external parameters.

Diagnostics and two-dimensional simulation of low-frequency inductively coupled plasmas with neutral gas heating and electron heat fluxes

This article presents the results on the diagnostics and numerical modeling of low-frequency (∼460 KHz) inductively coupled plasmas generated in a cylindrical metal chamber by an external flat spiral coil. Experimental data on the electron number densities and temperatures, electron energy distribution functions, and optical emission intensities of the abundant plasma species in low/intermediate pressure argon discharges are included. The spatial profiles of the plasma density, electron temperature, and excited argon species are computed, for different rf powers and working gas pressures, using the two-dimensional fluid approach. The model allows one to achieve a reasonable agreement between the computed and experimental data. The effect of the neutral gas temperature on the plasma parameters is also investigated. It is shown that neutral gas heating (at rf powers≥0.55kW) is one of the key factors that control the electron number density and temperature. The dependence of the average rf power loss, per electron-ion pair created, on the working gas pressure shows that the electron heat flux to the walls appears to be a critical factor in the total power loss in the discharge.

Magnetic fields and uniformity of radio frequency power deposition in low-frequency inductively coupled plasmas with crossed internal oscillating currents

Radial and axial distributions of magnetic fields in a low-frequency (∼460 kHz)inductively coupled plasmasource with two internal crossed planar rf current sheets are reported. The internal antenna configuration comprises two orthogonal sets of eight alternately reconnected parallel and equidistant copper litz wires in quartz enclosures and generates three magnetic (H z, H r, and H φ) and two electric (E φ and E r) field components at the fundamental frequency. The measurements have been performed in rarefied and dense plasmas generated in the electrostatic(E) and electromagnetic (H)discharge modes using two miniature magnetic probes. It is shown that the radial uniformity and depth of the rf power deposition can be improved as compared with conventional sources of inductively coupled plasmas with external flat spiral (“pancake”) antennas. Relatively deeper rf power deposition in the plasma source results in more uniform profiles of the optical emission intensity, which indicates on the improvement of the plasma uniformity over large chamber volumes. The results of the numerical modeling of the radial magnetic field profiles are found in a reasonable agreement with the experimental data.