Modeling the co-transport of viruses and colloids in unsaturated porous media

A mathematical model is developed to simulate the co-transport of viruses and colloids in unsaturated porous media under steady-state flow conditions. The virus attachment to the mobile and immobile colloids is described using a linear reversible kinetic model. Colloid transport is assumed to be decoupled from virus transport; that is, we assume that colloids are not affected by the presence of attached viruses on their surface. The governing equations,are solved numerically using an alternating three-step operator splitting approach. The model is verified by fitting three sets of experimental data published in the literature: (1) Syngouna and Chrysikopoulos (2013) and (2) Walshe et al. (2010), both on the co-transport of viruses and clay colloids under saturated conditions, and (3) Syngouna and Cluysikopoulos (2015) for the co-transport of viruses and clay colloids under unsaturated conditions. We found a good agreement between observed and fitted breakthrough curves (BTCs) under both saturated and unsaturated conditions. Then, the developed model was used to simulate the co-transport of viruses and colloids in porous media under unsaturated conditions, with the aim of understanding the relative importance of various processes on the co-transport of viruses and colloids in unsaturated porous media. The virus retention in porous media in the presence of colloids is greater during unsaturated conditions as compared to the saturated conditions due to: (1) virus attachment to the air-water interface (AWI), and (2) co-deposition of colloids with attached viruses on its surface to the AWL A sensitivity analysis of the model to various parameters showed that the virus attachment to AWI is the most sensitive parameter affecting the BTCs of both free viruses and total mobile viruses and has a significant effect on all parts of the BTC. The free and the total mobile viruses BTCs are mainly influenced by parameters describing virus attachment to the AIM, virus interaction with mobile and immobile colloids, virus attachment to solid-water interface (SWI), and colloid interaction with SWI and AWL The virus BTC is relatively insensitive to parameters describing the maximum adsorption capacity of the AWI for colloids, inlet colloid concentration, virus detachment rate coefficient from the SW!, maximum adsorption capacity of the AWI for viruses and inlet virus concentration. (C) 2015 Elsevier B.V. All rights reserved.

PLUTO plus : Near-Complete Modeling of Affine Transformations for Parallelism and Locality

Affine transformations have proven to be very powerful for loop restructuring due to their ability to model a very wide range of transformations. A single multi-dimensional affine function can represent a long and complex sequence of simpler transformations. Existing affine transformation frameworks like the Pluto algorithm, that include a cost function for modern multicore architectures where coarse-grained parallelism and locality are crucial, consider only a sub-space of transformations to avoid a combinatorial explosion in finding the transformations. The ensuing practical tradeoffs lead to the exclusion of certain useful transformations, in particular, transformation compositions involving loop reversals and loop skewing by negative factors. In this paper, we propose an approach to address this limitation by modeling a much larger space of affine transformations in conjunction with the Pluto algorithm's cost function. We perform an experimental evaluation of both, the effect on compilation time, and performance of generated codes. The evaluation shows that our new framework, Pluto+, provides no degradation in performance in any of the Polybench benchmarks. For Lattice Boltzmann Method (LBM) codes with periodic boundary conditions, it provides a mean speedup of 1.33x over Pluto. We also show that Pluto+ does not increase compile times significantly. Experimental results on Polybench show that Pluto+ increases overall polyhedral source-to-source optimization time only by 15%. In cases where it improves execution time significantly, it increased polyhedral optimization time only by 2.04x.

Microkinetic Modeling of CO Oxidation on Ionic Palladium-Substituted Ceria

This study presents a plausible dual-site mechanism and microkinetic model for CO oxidation over palladium-substituted ceria incorporating the theoretical oxygen storage capacity of different-catalysts into the kinetic model. A rate expression without prior assumption of rate-determining steps has been developed for the proposed microkinetic model using reaction route analysis. Experiments were conducted using various percentages of palladium in ceria that were synthesized by solution combustion. Obtained catalysts were characterized by X-ray diffraction, X-ray photoelectron spectra, and Brunauer-Emmett-Teller surface area measurements. A detailed mechanism was, developed, and the kinetic parameters and rate expression were validated with the conversion data obtained in the presence of the catalysts. Furthermore, a reduced rate expression based on rate-determining step and most abundant reactive intermediate approximation was obtained and tested against the original rate expression for different experimental conditions. From the results obtained it was concluded that the simulated rate predictions matched the experimental trend with reasonable accuracy, validating the kinetic parameters proposed it this study.

Numerical prediction of nozzle flow separation: Issue of turbulence modeling

Numerical simulation of separated flows in rocket nozzles is challenging because existing turbulence models are unable to predict it correctly. This paper addresses this issue with the Spalart-Allmaras and Shear Stress Transport (SST) eddy-viscosity models, which predict flow separation with moderate success. Their performances have been compared against experimental data for a conical and two contoured subscale nozzles. It is found that they fail to predict the separation location correctly, exhibiting sensitivity to the nozzle pressure ratio (NPR) and nozzle type. A careful assessment indicated how the model had to be tuned for better, consistent prediction. It is learnt that SST model's failure is caused by limiting of the shear stress inside boundary layer according to Bradshaw's assumption, and by over prediction of jet spreading rate. Accordingly, SST's coefficients were empirically modified to match the experimental wall pressure data. Results confirm that accurate RANS prediction of separation depends on the correct capture of the jet spreading rate, and that it is feasible over a wide range of NPRs by modified values of the diffusion coefficients in the turbulence model. (C) 2015 Elsevier Masson SAS. All rights reserved.

Modeling and mesoscopic damage constitutive relation of brittle short-fiber-reinforced composites

Aimed at brittle composites reinforced by randomly distributed short-fibers with a relatively large aspect ratio, stiffness modulus and strength, a mesoscopic material model was proposed. Based on the statistical description, damage mechanisms, damage-induced anisotropy, damage rate effect and stress redistribution, the constitutive relation were derived. By taking glass fiber reinforced polypropylene polymers as an example, the effect of initial orientation distribution of fibers, damage-induced anisotropy, and damage-rate effect on macro-behaviors of composites were quantitatively analyzed. The theoretical predictions compared favorably with the experimental results.

A Dynamic Loading Device for Suction Foundations in Centrifuge Modeling

Suction bucket foundations are widely used in the offshore platform for the exploitation of the offshore petroleum and natural gas resources. During winter seasons, ice sheets formed in Bohai Bay will impose strong impact and result in strong vibration on the platform. This paper describes a dynamic loading device developed on the geotechnical centrifuge and its application in modeling suction bucket foundation under the equivalent ice-induced vibration loadings. Some experimental results are presented. It is shown that when the loading amplitude is over a critical value, the sand at the upper part around the bucket softens or even liquefies. The excess pore pressure decreases from the upper part to the lower part of the sand foundation in vertical direction while decreases from near to far away from the bucket's side wall in the horizontal direction. Large settlements of the bucket and the sand around the bucket occur under the horizontal dynamic loading. The dynamic responses of the bucket with smaller size are heavier.

A Discussion on Modeling Shape Memory Alloy Embedded in A Composite Laminate as Axial Force and Elastic Foundation

In this paper, the possible error sources of the composite natural frequencies due to modeling the shape memory alloy (SMA) wire as an axial force or an elastic foundation and anisotropy are discussed. The great benefit of modeling the SMA wire as an axial force and an elastic foundation is that the complex constitutive relation of SMA can be avoided. But as the SMA wire and graphite-epoxy are rigidly bonded together, such constraint causes the re-distribution of the stress in the composite. This, together with anisotropy, which also reduces the structural stiffness can cause the relatively large error between the experimental data and theoretical results.

Three-Dimensional Modeling of Heat Transfer and Fluid Flow in Laminar-Plasma Material Re-Melting Processing

Modeling study is performed concerning the heat transfer and fluid flow for a laminar argon plasma jet impinging normally upon a flat workpiece exposed to the ambient air. The diffusion of the air into the plasma jet is handled by using the combined-diffusion-coefficient approach. The heat flux density and jet shear stress distributions at the workpiece surface obtained from the plasma jet modeling are then used to study the re-melting process of a carbon steel workpiece. Besides the heat conduction within the workpiece, the effects of the plasma-jet inlet parameters (temperature and velocity), workpiece moving speed, Marangoni convection, natural convection etc. on the re-melting process are considered. The modeling results demonstrate that the shapes and sizes of the molten pool in the workpiece are influenced appreciably by the plasma-jet inlet parameters, workpiece moving speed and Marangoni convection. The jet shear stress manifests its effect at higher plasma-jet inlet velocities, while the natural convection effect can be ignored. The modeling results of the molten pool sizes agree reasonably with available experimental data.

Modeling and Simulation of Liquid Pulsed Particulate Fluidized Beds

Pulsed fluidization is of considerable interest in process engineering for improving fluidization quality. Quantitative understanding of the pulsed two-phase flow behaviors is very important for proper design and optimum operation of such contactors. The

Physical modeling of untrenched submarine pipeline instability

Wave-induced instability of untrenched pipeline on sandy seabed is a `wave-soil-pipeline' coupling dynamic problem. To explore the mechanism of the pipeline instability, the hydrodynamic loading with U-shaped oscillatory flow tunnel is adopted, which is quite different from the previous experiment system. Based on dimensional analysis, the critical conditions for pipeline instability are investigated by altering pipeline submerged weight, diameter, soil parameters, etc. Based on the experimental results, different linear relationships between Froude number (Fr) and non-dimensional pipeline weight (G) are obtained for two constraint conditions. Moreover, the effects of loading history on the pipeline stability are also studied. Unlike previous experiments, sand scouring during the process of pipe's losing stability is detected in the present experiments. In addition, the experiment results are compared with the previous experiments, based on Wake II model for the calculation of wave-induced forces upon pipeline. It shows that the results of two kinds of experiments are comparable, but the present experiments provide better physical insight of the wave-soil-pipeline coupling effects.

Modeling on Runoff Concentration Caused by Rainfall on Hillslopes and Application in Maoping Slope

Based on the fact that the concentration flowlines of overland flow depend on the surface landform of hillslope, a kinematic wave model was developed for simulating runoff generation and flow concentration caused by rainfall on hillslopes. The model-simulated results agree well with experimental observations. Applying the model to the practical case of Maoping slope, we obtained the characteristics of runoff generation and infiltration on the slope. Especially, the simulated results adequately reflected the confluent pattern of surface runoff, which offers a scientific foundation for designing the drainage engineering on the Maoping slope.

The Head-On Colliding Process Of Binary Liquid Droplets At Low Velocity: High-Speed Photography Experiments And Modeling

The experimental and theoretical studies are reported in this paper for the head-on collisions of a liquid droplet with another of the same fluid resting on a solid substrate. The droplet on the hydrophobic polydimethylsiloxane (PDMS) substrate remains in a shape of an approximately spherical segment and is isometric to an incoming droplet. The colliding process of the binary droplets was recorded with high-speed photography. Head-on collisions saw four different types of response in our experiments: complete rebound, coalescence, partial rebound With conglutination, and coalescence accompanied by conglutination. For a complete rebound, both droplets exhibited remarkable elasticity and the contact time of the two colliding droplets was found to be in the range of 10-20 ms. With both droplets approximately considered as elastic bodies, Hertz contact theory was introduced to estimate the contact time for the complete rebound case. The estimated result Was found to be on the same order of magnitude as the experimental data, which indicates that the present model is reasonable. (C) 2008 Elsevier Inc. All rights reserved.

Micromechanics Modeling Of Strength For Nanocrystalline Copper

We present a model in this paper for predicting the inverse Hall-Petch phenomenon in nanocrystalline (NC) materials which are assumed to consist of two phases: grain phase of spherical or spheroidal shapes and grain boundary phase. The deformation of the grain phase has an elasto-viscoplastic behavior, which includes dislocation glide mechanism, Coble creep and Nabarro-Herring creep. However the deformation of grain boundary phase is assumed to be the mechanism of grain boundary diffusion. A Hill self-consistent method is used to describe the behavior of nanocrystalline pure copper subjected to uniaxial tension. Finally, the effects of grain size and its distribution, grain shape and strain rate on the yield strength and stress-strain curve of the pure copper are investigated. The obtained results are compared with relevant experimental data in the literature.

Microstructure Observation And Mechanical Behavior Modeling For Limnetic Nacre

In the present research, microstructure of a kind of limnetic shell (Hyriopsis cumingii) is observed and measured by using the scanning electron microscopy, and mechanical behavior experiments of the shell nacre are carried out by using bending and tensile tests. The dependence of mechanical properties of the shell nacre on its microstructure is analyzed by using a modified shear-lag model, and the overall stress-strain relation is obtained. The experimental results reveal that the mechanical properties of shell nacre strongly depend on the water contents of the limnetic shell. Dry nacre shows a brittle behavior, whereas wetting nacre displays a strong ductility. Compared to the tensile test, the bending test overestimates the strength and underestimates the Young's modulus. The modified shear-lag model can characterize the deformation features of nacre effectively.

Semi-gas kinetics model for performance modeling of flowing chemical oxygen-iodine lasers (COIL)

A semi-gas kinetics (SGK) model for performance analyses of flowing chemical oxygen-iodine laser (COIL) is presented. In this model, the oxygen-iodine reaction gas flow is treated as a continuous medium, and the effect of thermal motions of particles of different laser energy levels on the performances of the COIL is included and the velocity distribution function equations are solved by using the double-parameter perturbational method. For a premixed flow, effects of different chemical reaction systems, different gain saturation models and temperature, pressure, yield of excited oxygen, iodine concentration and frequency-shift on the performances of the COIL are computed, and the calculated output power agrees well with the experimental data. The results indicate that the power extraction of the SGK model considering 21 reactions is close to those when only the reversible pumping reaction is considered, while different gain saturation models and adjustable parameters greatly affect the output power, the optimal threshold gain range, and the length of power extraction.