Compressive advection and multi‐component methods for interface‐capturing

This paper develops methods for interface‐capturing in multiphase flows. The main novelties of these methods are as follows: (a) multi‐component modelling that embeds interface structures into the continuity equation; (b) a new family of triangle/tetrahedron finite elements, in particular, the P₁DG‐P₂(linear discontinuous between elements velocity and quadratic continuous pressure); (c) an interface‐capturing scheme based on compressive control volume advection methods and high‐order finite element interpolation methods; (d) a time stepping method that allows use of relatively large time step sizes; and (e) application of anisotropic mesh adaptivity to focus the numerical resolution around the interfaces and other areas of important dynamics. This modelling approach is applied to a series of pure advection problems with interfaces as well as to the simulation of the standard computational fluid dynamics benchmark test cases of a collapsing water column under gravitational forces (in two and three dimensions) and sloshing water in a tank. Two more test cases are undertaken in order to demonstrate the many‐material and compressibility modelling capabilities of the approach. Numerical simulations are performed on coarse unstructured meshes to demonstrate the potential of the methods described here to capture complex dynamics in multiphase flows. Copyright © 2015 John Wiley & Sons, Ltd.     

Exact Solution of a(2+1)-Dimensional Anisotropic Star in Finch and Skea Spacetime

We provide a new class of interior solution of a(2+1)-dimensional anisotropic star in Finch and Skea spacetime corresponding to the BTZ black hole. We develop the model by considering the MIT bag model EOS and a particular ansatz for the metric function grrproposed by Finch and Skea [M.R. Finch and J.E.F. Skea, Class. Quantum.Grav. 6(1989) 467]. Our model is free from central singularity and satisfies all the physical requirements for the acceptability of the model.     

Automated chemical resonance generation and structure filtration for kinetic modeling

This work discusses efficient and automated methods for constructing a set of representative resonance structures for arbitrary chemical species, including radicals and biradicals, consisting of the elements H, C, O, N, and S. Determining the representative reactive structures of chemical species is crucial for identification of reactive sites and consequently applying the correct reaction templates to generate the set of important reactions during automated chemical kinetic model generation. We describe a fundamental set of resonance pathway types, accounting for simple resonating structures, as well as global approaches for polycyclic aromatic species. Automatically discovering potential localized structures along with filtration to identify the representative structures was shown to be robust and relatively fast. The algorithms discussed here were recently implemented in the Reaction Mechanism Generator (RMG) software. The final structures proposed by this method were found to be in reasonable agreement with quantum chemical computation results of localized structure contributions to the resonance hybrid.     

Keller–Osserman estimates and removability result for the anisotropic porous medium equation with gradient absorption term

We study “large” nonnegative solutions for a class of quasilinear equations model of which is We give a sufficient condition on the exponents and for the removability of isolated singularities.     

A higher resolution edge‐based finite volume method for the simulation of the oil–water displacement in heterogeneous and anisotropic porous media using a modified IMPES method

In this article, we present a higher‐order finite volume method with a ‘Modified Implicit Pressure Explicit Saturation’ (MIMPES) formulation to model the 2D incompressible and immiscible two‐phase flow of oil and water in heterogeneous and anisotropic porous media. We used a median‐dual vertex‐centered finite volume method with an edge‐based data structure to discretize both, the elliptic pressure and the hyperbolic saturation equations. In the classical IMPES approach, first, the pressure equation is solved implicitly from an initial saturation distribution; then, the velocity field is computed explicitly from the pressure field, and finally, the saturation equation is solved explicitly. This saturation field is then used to re‐compute the pressure field, and the process follows until the end of the simulation is reached. Because of the explicit solution of the saturation equation, severe time restrictions are imposed on the simulation. In order to circumvent this problem, an edge‐based implementation of the MIMPES method of Hurtado and co‐workers was developed. In the MIMPES approach, the pressure equation is solved, and the velocity field is computed less frequently than the saturation field, using the fact that, usually, the velocity field varies slowly throughout the simulation. The solution of the pressure equation is performed using a modification of Crumpton's two‐step approach, which was designed to handle material discontinuity properly. The saturation equation is solved explicitly using an edge‐based implementation of a modified second‐order monotonic upstream scheme for conservation laws type method. Some examples are presented in order to validate the proposed formulation. Our results match quite well with others found in literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Brake Orbits of First Order Convex Hamiltonian Systems with Particular Anisotropic Growth

Combining the dual least action principle with Mountain-pass lemma,we obtain the existence of brake orbits for first-order convex Hamiltonian systems with particular anisotropic growth.     

Microfunnel‐filter‐based emulsification microextraction followed by gas chromatography for simple determination of organophosphorus pesticides in environmental water samples

A simple and fast method named microfunnel‐filter‐based emulsification microextraction is introduced for an efficient determination of some organophosphorus pesticides including diazinon, malathion, and chlorpyrifos in the environmental samples including the river, sea, and well water. This method is based upon the dispersion of a low‐toxicity organic solvent (dihexyl ether), as the extractant, in a high volume of an aqueous sample solution (45 mL). It is implemented without a centrifugation step, and using a syringe filter and a micro‐funnel, the phase separation and transfer of the enriched analytes to the gas chromatograph are simply achieved. By filtration of the extractant phase, a suitable sample clean‐up is obtained, and the total extraction time is just a few minutes. The factors influencing the extraction efficiency are optimized, and under the optimal conditions, the proposed method provides a good linearity (in the range of 15–1500 ng/mL (R² > 0.996). A high enrichment factor is obtained (in the range of 306–342), and the method provides low limits of detection and quantification (in the ranges of 4–8 and 15–25 ng/mL, respectively).     

Azimuthal orientational correlations due to excluded volume epitaxy in growing anisotropic grains

The microstructure of anisotropically shaped grains can strongly influence a range of material properties, including transport, mechanical and electro-optical. A grain-structure-related phenomenon, known as excluded volume epitaxy (EVE), is reported in this study. EVE is a local, inter-grain orientational correlations effect, which results from a combination of continuous nucleation of anisotropic grains and impingement of growing grains. Due to EVE, anisotropically shaped grains have a tendency to be similarly aligned in a local neighbourhood, despite the absence of any forced global orientation in the sample. The effect has been repeatedly observed by the authors in block copolymers, as illustrated by a representative TEM image. Optical microscopy of anisotropically shaped non-polymeric crystals revealed the generality of this effect. The simulation study revealed a tendency for azimuthal, inter-grain orientational correlation and re-confirmed the experimental observation of EVE.     

A comparison of numerical and analytical radiative-transfer solutions for plane albedo of natural waters

Several numerical and analytical solutions of the radiative transfer equation (RTE) were compared for plane albedo in a problem of solar light reflection by sea water. The study incorporated the simplest case—a semi-infinite one-dimensional plane—parallel absorbing and scattering homogeneous layer illuminated by a monodirectional light beam. Inelastic processes (such as Raman scattering and fluorescence), polarization and air-water surface refraction-reflection effects, were not considered. Algorithms were based on the invariant imbedding method and two different variants of the discrete ordinate method (DOM). Calculations were performed using parameters across all possible ranges (single-scattering albedo ω₀ and refracted solar zenith angle θ₁), but with a special emphasis on natural waters. All computations were made for two scattering phase functions, which included an almost isotropic Rayleigh phase function and strongly anisotropic double-peaked Fournier-Forand-Mobley phase function. Models were validated using quasi-single-scattering (QSSA) and exponential approximations, which represent the extreme cases of ω₀→0 and ω₀→1, respectively. All methods yielded relative differences within 1.8% for modeled natural waters. An analysis of plane albedo behavior resulted in the development of a new extended QSSA approximation, which when applied in conjunction with the extended Hapke approximation developed earlier, resulted in a maximum relative error of 2.7%. The study results demonstrated that for practical applications, the estimation of inherent optical properties from observed reflectance can best be achieved using an extended Hapke approximation.