Geological realism in hydrogeological and geophysical inverse modeling: A review

Scientific curiosity, exploration of georesources and environmental concerns are pushing the geoscientific research community toward subsurface investigations of ever-increasing complexity. This review explores various approaches to formulate and solve inverse problems in ways that effectively integrate geological concepts with geophysical and hydrogeological data. Modern geostatistical simulation algorithms can produce multiple subsurface realizations that are in agreement with conceptual geological models and statistical rock physics can be used to map these realizations into physical properties that are sensed by the geophysical or hydrogeological data. The inverse problem consists of finding one or an ensemble of such subsurface realizations that are in agreement with the data. The most general inversion frameworks are presently often computationally intractable when applied to large-scale problems and it is necessary to better understand the implications of simplifying (1) the conceptual geological model (e.g., using model compression); (2) the physical forward problem (e.g., using proxy models); and (3) the algorithm used to solve the inverse problem (e.g., Markov chain Monte Carlo or local optimization methods) to reach practical and robust solutions given today's computer resources and knowledge. We also highlight the need to not only use geophysical and hydrogeological data for parameter estimation purposes, but also to use them to falsify or corroborate alternative geological scenarios.

Effect of hamstring injuries on kicking in soccer using inverse dynamics

Las lesiones musculares del muslo tienen una gran incidencia en el fútbol. El objetivo del estudio ha sido desarrollar un nuevo procedimiento para evaluar el efecto de las lesiones de isquiotibiales en los golpeos con el pie s en el fútbol utilizando los principios de la dinámica inversa. El trabajo se ha centrado en la evaluación de la diferencia entre sujetos que habían sufrido la lesión en los últimos 5 años y los que no. Se analizaron 17 jugadores de fútbol profesionales realizando cinco tiros con el empeine y cinco con el interior del pie. Los movimientos se registraron mediante una plataforma de fuerza y un sistema de captura de movimiento Vicon funcionando a 500Hz. Los participantes también tomaron parte en una prueba de isocinético en la que se midió el torque isocinético en 60 º/s y 120 º/s. Se observaron diferencias significativas en los parámetros cinemáticos y cinéticos entre los dos grupos (lesionados y no lesionados) en la fase posterior del golpeo y en el instante de máxima flexión de cadera. No se encontraron diferencias significativas entre los dos grupos en la prueba isocinética tradicional. Estos resultados indican que el procedimiento empleado probablemente podría ser muy útil en la evaluación del efecto de las lesiones de isquiotibiales en el fútbol.

Inverse Spectral Methods in Acoustic Normal Mode Velocimetry of High Reynolds Number Spherical Couette Flows

Experimental geophysical fluid dynamics often examines regimes of fluid flow infeasible for computer simulations. Velocimetry of zonal flows present in these regimes brings many challenges when the fluid is opaque and vigorously rotating; spherical Couette flows with molten metals are one such example. The fine structure of the acoustic spectrum can be related to the fluid’s velocity field, and inverse spectral methods can be used to predict and, with sufficient acoustic data, mathematically reconstruct the velocity field. The methods are to some extent inherited from helioseismology. This work develops a Finite Element Method suitable to matching the geometries of experimental setups, as well as modelling the acoustics based on that geometry and zonal flows therein. As an application, this work uses the 60-cm setup Dynamo 3.5 at the University of Maryland Nonlinear Dynamics Laboratory. Additionally, results obtained using a small acoustic data set from recent experiments in air are provided.

Distance and overcurrent relay coordination considering non standardized inverse time curves

Protective relaying comprehends several procedures and techniques focused on maintaining the power system working safely during and after undesired and abnormal network conditions, mostly caused by faulty events. Overcurrent relay is one of the oldest protective relays, its operation principle is straightforward: when the measured current is greater than a specified magnitude the protection trips; less variables are required from the system in comparison with other protections, causing the overcurrent relay to be the simplest and also the most difficult protection to coordinate; its simplicity is reflected in low implementation, operation, and maintenance cost. The counterpart consists in the increased tripping times offered by this kind of relays mostly before faults located far from their location; this problem can be particularly accentuated when standardized inverse-time curves are used or when only maximum faults are considered to carry out relay coordination. These limitations have caused overcurrent relay to be slowly relegated and replaced by more sophisticated protection principles, it is still widely applied in subtransmission, distribution, and industrial systems. In this work, the use of non standardized inverse-time curves, the model and implementation of optimization algorithms capable to carry out the coordination process, the use of different levels of short circuit currents, and the inclusion of distance relays to replace insensitive overcurrent ones are proposed methodologies focused on the overcurrent relay performance improvement. These techniques may transform the typical overcurrent relay into a more sophisticated one without changing its fundamental principles and advantages. Consequently a more secure and still economical alternative can be obtained, increasing its implementation area

Ziziphus mauritiana leaf extract emulsion for skin rejuvenation

Purpose: To formulate stable water in oil (W/O) emulsion containing hydroalcoholic crude extract of Ziziphus mauritiana leaves for skin rejuvenation. Methods: Placebo (base) without any plant extract and formulation with 4 % Ziziphus mauritiana extract were prepared by mixing. Samples of the emulsions were subjected to varying storage conditions, i.e., 8, 25, 40 oC and 40 oC + 75 % relative humidity for a period of 4 weeks to predict their stability. During this period, stability parameters, including liquefaction, phase separation, color, electrical conductivity, centrifugation and pH were monitored at specified time intervals. Skin rejuvenation was evaluated using 13 healthy human volunteers over a period of 8 weeks. During this period, various skin parameters such as erythema, melanin level, moisture content, elasticity and sebum content of the skin were evaluated at specified intervals. Results: Both the active formulation and placebo were stable in terms of liquifaction, phase separation and color at all the storage conditions of temperature and humidity. Active formulation showed statistically significant (p < 0.05) improvement in skin melanin as well as in skin moisture and sebum levels, whereas these properties were reduced or even absent in the placebo formulation (p > 0.05). Both active and placebo formulations changed skin elasticity and erythema significantly (p < 0.05). Conclusion: It is evident from the findings that the leaf extract of Ziziphus mauritiana possesses antiaging properties as well as exert skin lightening, moisturizing and viscoelastic effects on human skin.

Formulation and stability of topical water in oil emulsion containing corn silk extract

Purpose: To formulate the water in oil (W/O) emulsion of corn silk (CS) extract and to evaluate its stability at various storage conditions. Methods: Ethanol CS extract was prepared using maceration (cold) technique. A 4 % CS emulsion was prepared using varying concentrations of liquid paraffin, ABIL EM90 and water. The formulations were kept at 40 oC for 28 days and to screen out the less stable formulations. The remaining formulations were further stressed at 50 oC to choose the most stable formulation. The optimized formulation was evaluated for physical characteristics including phase separation, rheology and mean droplet size. The physical stability of the formulation was evaluated by monitoring these parameters over a period of 12 weeks at 8, 25, 40 and 40 oC, and 75 % RH. Results: The chosen formulation showed good resistance to phase separation on centrifugation under all storage conditions. Rheological behavior followed non-Newtonian pseudoplastic pattern at various storage conditions. Mean droplet size of freshly prepared formulation was 2.98 ± 1.32 µm and did not show significant (p < 0.05) changes at normal storage conditions (8 and 25 oC). Conclusion: The findings indicate that the developed CS extract W/O emulsion is stable and therefore may be suitable for topical use on skin as an antioxidant preparation.

Fouling in emulsion polymerization reactors

160 p.

STUDY OF A DIRECT SAMPLING METHOD FOR THE INVERSE MEDIUM SCATTERING PROBLEM

Direct sampling methods are increasingly being used to solve the inverse medium scattering problem to estimate the shape of the scattering object. A simple direct method using one incident wave and multiple measurements was proposed by Ito, Jin and Zou. In this report, we performed some analytic and numerical studies of the direct sampling method. The method was found to be effective in general. However, there are a few exceptions exposed in the investigation. Analytic solutions in different situations were studied to verify the viability of the method while numerical tests were used to validate the effectiveness of the method.

Minor and major consolidations in inverse DEA:definition and determination

Many production systems have acquisition and merge operations to increase productivity. This paper proposes a novel method to anticipate whether a merger in a market is generating a major or a minor consolidation, using InvDEA model. A merger between two or more decision making units (DMUs) producing a single merged DMU that affects the efficiency frontier, defined by the pre-consolidation market conditions, is called a major consolidation. The corresponding alternative case is called a minor consolidation. A necessary and sufficient condition to distinguish the two types of consolidations is proven and two numerical illustrations in banking and supply chain management are discussed. The crucial importance of anticipating the magnitude of a consolidation in a market is outlined.

An approximate method for the solution of an influence function foil rolling model

Fleck and Johnson (Int. J. Mech. Sci. 29 (1987) 507) and Fleck et al. (Proc. Inst. Mech. Eng. 206 (1992) 119) have developed foil rolling models which allow for large deformations in the roll profile, including the possibility that the rolls flatten completely. However, these models require computationally expensive iterative solution techniques. A new approach to the approximate solution of the Fleck et al. (1992) Influence Function Model has been developed using both analytic and approximation techniques. The numerical difficulties arising from solving an integral equation in the flattened region have been reduced by applying an Inverse Hilbert Transform to get an analytic expression for the pressure. The method described in this paper is applicable to cases where there is or there is not a flat region.

Micropolar flow in a porous channel with high mass transfer

Two dimensional flow of a micropolar fluid in a porous channel is investigated. The flow is driven by suction or injection at the channel walls, and the micropolar model due to Eringen is used to describe the working fluid. An extension of Berman's similarity transform is used to reduce the governing equations to a set of non-linear coupled ordinary differential equations. The latter are solved for large mass transfer via a perturbation analysis where the inverse of the cross-flow Reynolds number is used as the perturbing parameter. Complementary numerical solutions for strong injection are also obtained using a quasilinearisation scheme, and good agreement is observed between the solutions obtained from the perturbation analysis and the computations.

Adaptive unsupervised learning of human actions

Automatic detection of suspicious activities in CCTV camera feeds is crucial to the success of video surveillance systems. Such a capability can help transform the dumb CCTV cameras into smart surveillance tools for fighting crime and terror. Learning and classification of basic human actions is a precursor to detecting suspicious activities. Most of the current approaches rely on a non-realistic assumption that a complete dataset of normal human actions is available. This paper presents a different approach to deal with the problem of understanding human actions in video when no prior information is available. This is achieved by working with an incomplete dataset of basic actions which are continuously updated. Initially, all video segments are represented by Bags-Of-Words (BOW) method using only Term Frequency-Inverse Document Frequency (TF-IDF) features. Then, a data-stream clustering algorithm is applied for updating the system's knowledge from the incoming video feeds. Finally, all the actions are classified into different sets. Experiments and comparisons are conducted on the well known Weizmann and KTH datasets to show the efficacy of the proposed approach.

Fast reconstruction of aerodynamic shapes using evolutionary algorithms and virtual Nash strategies in a CFD design environment

This paper compares the performances of two different optimisation techniques for solving inverse problems; the first one deals with the Hierarchical Asynchronous Parallel Evolutionary Algorithms software (HAPEA) and the second is implemented with a game strategy named Nash-EA. The HAPEA software is based on a hierarchical topology and asynchronous parallel computation. The Nash-EA methodology is introduced as a distributed virtual game and consists of splitting the wing design variables - aerofoil sections - supervised by players optimising their own strategy. The HAPEA and Nash-EA software methodologies are applied to a single objective aerodynamic ONERA M6 wing reconstruction. Numerical results from the two approaches are compared in terms of the quality of model and computational expense and demonstrate the superiority of the distributed Nash-EA methodology in a parallel environment for a similar design quality.

Krylov subspace methods for approximating functions of symmetric positive definite matrices with applications to applied statistics and anomalous diffusion

Matrix function approximation is a current focus of worldwide interest and finds application in a variety of areas of applied mathematics and statistics. In this thesis we focus on the approximation of A^(-α/2)b, where A ∈ ℝ^(n×n) is a large, sparse symmetric positive definite matrix and b ∈ ℝ^n is a vector. In particular, we will focus on matrix function techniques for sampling from Gaussian Markov random fields in applied statistics and the solution of fractional-in-space partial differential equations. Gaussian Markov random fields (GMRFs) are multivariate normal random variables characterised by a sparse precision (inverse covariance) matrix. GMRFs are popular models in computational spatial statistics as the sparse structure can be exploited, typically through the use of the sparse Cholesky decomposition, to construct fast sampling methods. It is well known, however, that for sufficiently large problems, iterative methods for solving linear systems outperform direct methods. Fractional-in-space partial differential equations arise in models of processes undergoing anomalous diffusion. Unfortunately, as the fractional Laplacian is a non-local operator, numerical methods based on the direct discretisation of these equations typically requires the solution of dense linear systems, which is impractical for fine discretisations. In this thesis, novel applications of Krylov subspace approximations to matrix functions for both of these problems are investigated. Matrix functions arise when sampling from a GMRF by noting that the Cholesky decomposition A = LL^T is, essentially, a `square root' of the precision matrix A. Therefore, we can replace the usual sampling method, which forms x = L^(-T)z, with x = A^(-1/2)z, where z is a vector of independent and identically distributed standard normal random variables. Similarly, the matrix transfer technique can be used to build solutions to the fractional Poisson equation of the form ϕn = A^(-α/2)b, where A is the finite difference approximation to the Laplacian. Hence both applications require the approximation of f(A)b, where f(t) = t^(-α/2) and A is sparse. In this thesis we will compare the Lanczos approximation, the shift-and-invert Lanczos approximation, the extended Krylov subspace method, rational approximations and the restarted Lanczos approximation for approximating matrix functions of this form. A number of new and novel results are presented in this thesis. Firstly, we prove the convergence of the matrix transfer technique for the solution of the fractional Poisson equation and we give conditions by which the finite difference discretisation can be replaced by other methods for discretising the Laplacian. We then investigate a number of methods for approximating matrix functions of the form A^(-α/2)b and investigate stopping criteria for these methods. In particular, we derive a new method for restarting the Lanczos approximation to f(A)b. We then apply these techniques to the problem of sampling from a GMRF and construct a full suite of methods for sampling conditioned on linear constraints and approximating the likelihood. Finally, we consider the problem of sampling from a generalised Matern random field, which combines our techniques for solving fractional-in-space partial differential equations with our method for sampling from GMRFs.