Análise matemática de soluções descontínuas de leis de conservação hiperbólicas e resoluções numéricas para a captura de ondas de choque em escoamentos multifásicos em meios porosos

O processo de recuperação secundária de petróleo é comumente realizado com a injeção de água ou gás no reservatório a fim de manter a pressão necessária para sua extração. Para que o investimento seja viável, os gastos com a extração precisam ser menores do que o retorno financeiro obtido com a produção de petróleo. Objetivando-se estudar possíveis cenários para o processo de exploração, costuma-se utilizar simulações dos processos de extração. As equações que modelam esse processo de recuperação são de caráter hiperbólico e não lineares, as quais podem ser interpretadas como Leis de Conservação, cujas soluções são complexas por suas naturezas descontínuas. Essas descontinuidades ou saltos são conhecidas como ondas de choque. Neste trabalho foi abordada uma análise matemática para os fenômenos oriundos de leis de conservação, para em seguida utilizá-la no entendimento do referido problema. Foram estudadas soluções fracas que, fisicamente, podem ser interpretadas como ondas de choque ou rarefação, então, para que fossem distinguidas as fisicamente admissíveis, foi considerado o princípio de entropia, nas suas diversas formas. As simulações foram realizadas nos âmbitos dos escoamentos bifásicos e trifásicos, em que os fluidos são imiscíveis e os efeitos gravitacionais e difusivos, devido à pressão capilar, foram desprezados. Inicialmente, foi feito um estudo comparativo de resoluções numéricas na captura de ondas de choque em escoamento bifásico água-óleo. Neste estudo destacam-se o método Composto LWLF-k, o esquema NonStandard e a introdução da nova função de renormalização para o esquema NonStandard, onde obteve resultados satisfatórios, principalmente em regiões onde a viscosidade do óleo é muito maior do que a viscosidade da água. No escoamento bidimensional, um novo método é proposto, partindo de uma generalização do esquema NonStandard unidimensional. Por fim, é feita uma adaptação dos métodos LWLF-4 e NonStandard para a simulação em escoamentos trifásicos em domínios unidimensional. O esquema NonStandard foi considerado mais eficiente nos problemas abordados, uma vez que sua versão bidimensional mostrou-se satisfatória na captura de ondas de choque em escoamentos bifásicos em meios porosos.

Some characteristics of thin shear layers in homogeneous turbulent flow.

Tomographic particle image velocimetry measurements of homogeneous isotropic turbulence that have been made in a large mixing tank facility at Cambridge are analysed in order to characterize thin highly sheared regions that have been observed. The results indicate that such regions coincide with regions of high enstrophy, dissipation and stretching. Large velocity jumps are observed across the width of these regions. The thickness of the shear layers seems to scale with the Taylor microscale, as has been suggested previously. The results discussed here concentrate on examining individual realizations rather than statistics of these regions.

Eliminating discharge imbalances in predicting shallow water flows with shocks

The shallow water equations are widely used in modelling environmental flows. Being a hyperbolic system of differential equations, they admit shocks that represent hydraulic jumps and bores. Although the water surface can be solved satisfactorily with the modern shock-capturing schemes, the predicted flow rate often suffers from imbalances where shocks occur, eg the mass conservation is violated by failing to maintain a constant discharge rate at every cross-section in a steady open channel flow. A total-variation-diminishing Lax-Wendroff scheme is developed, and used to demonstrate how to achieve an exact flux balance. The performance of the proposed methods is inspected through some test cases, which include 1- and 2-dimensional, flat and irregular bed scenarios. The proposed methods are shown to preserve the mass exactly, and can be easily extended to other shock-capturing models.

Local equilibrium of the Gibbs interface in two-phase systems

We analyze the local equilibrium assumption for interfaces from the perspective of gauge transformations, which are the small displacements of Gibbs' dividing surface. The gauge invariance of thermodynamic properties turns out to be equivalent to conditions for jumps of bulk densities across the interface. This insight strengthens the foundations of the local equilibrium assumption for interfaces and can be used to characterize nonequilibrium interfaces in a compact and consistent way, with a clear focus on gauge-invariant properties. Using the principle of gauge invariance, we show that the validity of Clapeyron equations can be extended to nonequilibrium interfaces, and an additional jump condition for the momentum density is recognized to be of the Clapeyron type. © 2012 Europhysics Letters Association.

The temporal evolution of feedback gains rapidly update to task demands.

Recent theoretical frameworks such as optimal feedback control suggest that feedback gains should modulate throughout a movement and be tuned to task demands. Here we measured the visuomotor feedback gain throughout the course of movements made to "near" or "far" targets in human subjects. The visuomotor gain showed a systematic modulation over the time course of the reach, with the gain peaking at the middle of the movement and dropping rapidly as the target is approached. This modulation depends primarily on the proportion of the movement remaining, rather than hand position, suggesting that the modulation is sensitive to task demands. Model-predictive control suggests that the gains should be continuously recomputed throughout a movement. To test this, we investigated whether feedback gains update when the task goal is altered during a movement, that is when the target of the reach jumped. We measured the visuomotor gain either simultaneously with the jump or 100 ms after the jump. The visuomotor gain nonspecifically reduced for all target jumps when measured synchronously with the jump. However, the visuomotor gain 100 ms later showed an appropriate modulation for the revised task goal by increasing for jumps that increased the distance to the target and reducing for jumps that decreased the distance. We conclude that visuomotor feedback gain shows a temporal evolution related to task demands and that this evolution can be flexibly recomputed within 100 ms to accommodate online modifications to task goals.

Investigation on multiple-port microcylinder lasers based on coupled modes

Microcylinder resonators with multiple ports connected to waveguides are investigated by 2D finite-difference time-domain (FDTD) simulation for realizing microlasers with multiple outputs. For a 10 mu m radius microcylinder with a refractive index of 3.2 and three 2 mu m wide waveguides, confined mode at the wavelength of 1542.3 nm can have a mode Q factor of 6.7 x 10(4) and an output coupling efficiency of 0.76. AlGaInAs/InP microcylinder lasers with a radius of 10 mu m and a 2 mu m wide output waveguide are fabricated by planar processing techniques. Continuous-wave electrically injected operation is realized with a threshold current of 4 mA at room temperature, and the jumps of output power are observed accompanying a lasing mode transformation.

Peculiar photocurrent response due to Gamma-X coupling in a GaAs/AlAs heterostructure

Peculiar current jumps and hysteresis in current-voltage curves are reported in an illuminated heterostructure consisting basically of a thick AlAs layer and a narrow GaAs quantum well. These novel features come from the photon-assisted transfer of electron-hole pairs and the resultant charge polarization in the structure, mainly caused by the resonant Gamma-X coupling at the heterointerfaces. Using the transfer-matrix method, the simulated current density-voltage curve reproduces the main features of the experimental observations in the case where the influence of resonant Gamma-X coupling at the heterointerfaces is included, further confirming the physical mechanism involved. The structure presented here may be used as a new type of photonic memory cell and also as an optically controlled switch.

Highly conductive, methanol resistant fuel cell membranes fabricated by layer-by-layer self-assembly of inorganic heteropolyacid

Layer-by-layer (LBL) self-assembly is a simple and elegant method of constructing organic-inorganic composite thin films from environmentally benign aqueous solutions. In this paper, we utilize this method to develop proton-exchange membranes for fuel cells. The multilayer film is constructed onto the surface of sulfonated poly(arylene ether ketone) (SPAEK-COOH) membrane by LBL self-assembly of polycation chitosan (CTS) and negatively charged inorganic particle phosphotungstic acid (VIA). The highly conductive inorganic nanoparticles ensure SPAEK-COOH-(CTS/PTA)(n) membranes to maintain high proton conductivity values up to 0.086 S cm(-1) at 25 degrees C and 0.24S cm(-1) at 80 degrees C, which are superior than previous LBL assembled electrolyte systems.

基于三重震相波形模拟研究亚洲东北部上地幔过渡带速度结构

Both the global and regional P wave tomographic studies have revealed significant deep structural heterogeneities in subduction zone regions. In particular, low-velocity anomalies have been observed beneath the descending high-velocity slabs in a number of subduction zones. The limited resolution at large depths and possible trade-off between the high and low velocities, however, make it difficult to substantiate this feature and evaluate the vertical extent of the low-velocity structure. From broadband waveform modeling of triplicated phases near the 660-km discontinuity for three deep events, we constrained both the P and SH wave velocity structures around the base of the upper mantle in northeast Asia. For the two events beneath the southern Kurile, the rays traveled through the lowermost transition zone and uppermost lower mantle under the descending Pacific slab. Our preferred models consistently suggest normal-to-lower P and significantly low SH velocities above and below the 660-km discontinuity extending to about 760-km depth compared with the global IASP91 model, corroborating previous observations for a slow structure underneath the slab. In contrast, both high P and SH velocity anomalies are shown in our preferred models for the Japan subduction zone region, likely reflecting the structural feature of a slab stagnant above the 660-km discontinuity. The velocity jumps across the 660-km discontinuity were found to be on average 4.5% and 7% for P and S waves under the south Kurile, and 3% and 6% under the Japan subduction zone. The respective velocity contrasts in the two regions are consistent with mineralogical models for colder slab interior and hotter under-slab areas. Based on mineral physics data, the depth-averaged ~1.5% P and ~2.5% SH velocity differences in the depth range of 560-760 km between the two regions could be primarily explained by a 350~450K temperature variation, although the presence of about 0.5wt%~1wt% water might also contribute to the subtle velocity variations near the base of the transition zone in the southern Kurile. From our modeling results, we speculate that the slow structure in the southern Kurile may be correlated to the low velocity zone observed previously around the 410-km discontinuity under Northern Honshu. Both are probably associated with a thermal anomaly rooted in the lower mantle beneath the subduction zone in northeast Asia.

中国大中型气田天然气成藏物理化学模拟研究

Natural gas pays more important role in the society as clean fuel. Natural gas exploration has been enhanced in recent years in many countries. It also has prospective future in our country through "85" and "95" national research. Many big size gas fields have been discovered in different formations in different basins such as lower and upper Paleozoic in Erdos basin, Tertiary system in Kuche depression in Tarim basin, Triassic system in east of Sichuan basin. Because gas bearing basins had been experienced multiple tectogenesis. The characteristics of natural gases usually in one gas field are that they have multiple source rocks and are multiple maturities and formed in different ages. There has most difficult to research on the gas-rock correlation and mechanism of gas formation. Develop advanced techniques and methods and apply them to solve above problems is necessary. The research is focused on the critical techniques of geochemistry and physical simulation of gas-rock correlation and gas formation. The lists in the following are conclusions through research and lots of experiments. I 8 advanced techniques have been developed or improved about gas-rock correlation and gas migration, accumulation and formation. A series of geochemistry techniques has been developed about analyzing inclusion enclave. They are analyzing gas and liquid composition and biomarker and on-line individual carbon isotope composition in inclusion enclave. These techniques combing the inclusion homogeneous temperature can be applied to study on gas-rock correlation directly and gas migration, filling and formation ages. Technique of on-line determination individual gas carbon isotope composition in kerogen and bitumen thermal pyrolysis is developed. It is applied to determine the source of natural is kerogen thermal degradation or oil pyrolysis. Method of on-line determination individual gas carbon isotope composition in rock thermal simulation has being improved. Based on the "95"former research, on-line determination individual gas carbon isotope composition in different type of maceral and rocks thermal pyrolys is has been determined. The conclusion is that carbon isotope composition of benzene and toluene in homogenous texture kerogen thermal degradation is almost same at different maturity. By comparison, that in mixture type kerogen thermal pyrolysis jumps from step to step with the changes of maturity. This conclusion is a good proof of gas-rock dynamic correlation. 3. Biomarker of rock can be determined directly through research. It solves the problems such as long period preparing sample, light composition losing and sample contamination etc. It can be applied to research the character of source rock and mechanism of source rock expulsion and the path of hydrocarbon migration etc. 4. The process of hydrocarbon dynamic generation in source rock can be seen at every stage applying locating observation and thermal simulation of ESEM. The mechanism of hydrocarbon generation and expulsion in source rock is discussed according to the experiments. This technique is advanced in the world. 5. A sample injection system whose character is higher vacuum, lower leaks and lower blank has been built up to analyze inert gas. He,Ar,Kr and Xe can be determined continuously on one instrument and one injection. This is advanced in domestic. 7. Quality and quantity analysis of benzene ring compounds and phenolic compounds and determination of organic acid and aqueous gas analysis are applied to research the relationship between compounds in formation water and gas formation. This is another new idea to study the gas-rock correlation and gas formation. 8. Inclusion analysis data can be used to calculate the Paleo-fluid density, Paleo-geothermal gradient and Paleo-geopressure gradient and then to calculate the Paleo-fluid potential. It's also a new method to research the direction of hydrocarbon migration and accumulation. 9. Equipment of natural gas formation simulation is produced during the research to probe how the physical properties of rock affect the gas migration and accumulation and what efficiency of gas migrate and factors of gas formation and the models of different type of migration are. II study is focused on that if the source rocks of lower Paleozoic generated hydrocarbon and what the source rocks of weathered formation gas pool and the mechanism of gas formation are though many advanced techniques application. There are four conclusions. 1.The maturity of Majiagou formation source rocks is higher in south than that in north. There also have parts of the higher maturity in middle and east. Anomalous thermal pays important role in big size field formation in middle of basin. 2. The amount of gas generation in high-over maturity source rocks in lower Paleozoic is lager than that of most absorption of source rocks. Lower Paleozoic source rocks are effective source rocks. Universal bitumen exists in Ordovician source rocks to prove that Ordovician source rocks had generated hydrocarbon. Bitumen has some attribution to the middle gas pool formation. 3. Comprehensive gas-rock correlation says that natural gases of north, west, south of middle gas field of basin mainly come from lower Paleozoic source rocks. The attribution ratio of lower Paleozoic source rocks is 60%-70%. Natural gases of other areas mainly come from upper Paleozoic. The attribution ratio of upper Paleozoic source rocks is 70%. 4. Paleozoic gases migration phase of Erdos basin are also interesting. The relative abundance of gasoline aromatic is quite low especially toluene that of which is divided by that of methyl-cyclohexane is less than 0.2 in upper Paleozoic gas pool. The migration phase of upper Paleozoic gas may be aqueous phase. By comparison, the relative abundance of gasoline aromatic is higher in lower Paleozoic gas. The distribution character of gasoline gas is similar with that in source rock thermal simulation. The migration phase of it may be free phase. IH Comprehensive gas-rock correlation is also processed in Kuche depression Tarim basin. The mechanism of gas formation is probed and the gas formation model has been built up. Four conclusions list below. 1. Gases in Kuche depression come from Triassic-Jurassic coal-measure source rocks. They are high-over maturity. Comparatively, the highest maturity area is Kelasu, next is Dabei area, Yinan area. 2. Kerogen thermal degradation is main reason of the dry gas in Kuche depression. Small part of dry gas comes from oil pyrolysis. VI 3.The K12 natural gas lays out some of hydro-gas character. Oil dissolved in the gas. Hydro-gas is also a factor making the gas drier and carbon isotope composition heavier. 4. The mechanism and genesis of KL2 gas pool list as below. Overpressure has being existed in Triassic-Jurassic source rocks since Keche period. Natural gases were expulsed by episode style from overpressure source rocks. Hetero-face was main migration style of gas, oil and water at that time. The fluids transferred the pressure of source rocks when they migrated and then separated when they got in reservoir. After that, natural gas migrated up and accumulated and formed with the techno-genesis. Tectonic extrusion made the natural gas overpressure continuously. When the pressure was up to the critical pressure, the C6-C7 composition in natural gas changed. The results were that relative abundance of alkane and aromatic decreased while cycloalkane and isoparaffin increased. There was lots of natural gas filling during every tectonic. The main factors of overpressure of natural gas were tectonic extrusion and fluid transferring pressure of source rocks. Well preservation was also important in the KL2 gas pool formation. The reserves of gas can satisfy the need of pipeline where is from west to east. IV A good idea of natural gas migration and accumulation modeling whose apparent character is real core and formation condition is suggested to model the physical process of gas formation. Following is the modeling results. 1. Modeling results prove that the gas accumulation rule under cap layer and gas fraction on migration path. 2. Natural gas migration as free phase is difficult in dense rock. 3. Natural gases accumulated easily in good physical properties reservoirs where are under the plugging layer. Under the condition of that permeability of rock is more than 1 * 10~(-3)μm~(-1), the more better the physical properties and the more bigger pore of rock, the more easier the gas accumulation in there. On the contrary, natural gas canonly migrate further to accumulate in good physical properties of rock. 4. Natural gas migrate up is different from that down. Under the same situation, the amount of gas migration up is lager than that of gas migration down and the distance of migration up is 3 times as that of migration down. 5. After gas leaks from dense confining layer, the ability of its dynamic plug-back decreased apparently. Gas lost from these arils easily. These confining layer can confine again only after geology condition changes. 6. Water-wetted and capillary-blocking rocks can't block water but gases generally. The result is that water can migrate continuously through blocking rocks but the gases stay under the blocking rocks then form in there. The experiments have proved the formation model of deep basin gas.

A Network Charge-Orineted MOS Transistor Model

The MOS transistor physical model as described in [3] is presented here as a network model. The goal is to obtain an accurate model, suitable for simulation, free from certain problems reported in the literature [13], and conceptually as simple as possible. To achieve this goal the original model had to be extended and modified. The paper presents the derivation of the network model from physical equations, including the corrections which are required for simulation and which compensate for simplifications introduced in the original physical model. Our intrinsic MOS model consists of three nonlinear voltage-controlled capacitors and a dependent current source. The charges of the capacitors and the current of the current source are functions of the voltages $V_{gs}$, $V_{bs}$, and $V_{ds}$. The complete model consists of the intrinsic model plus the parasitics. The apparent simplicity of the model is a result of hiding information in the characteristics of the nonlinear components. The resulted network model has been checked by simulation and analysis. It is shown that the network model is suitable for simulation: It is defined for any value of the voltages; the functions involved are continuous and satisfy Lipschitz conditions with no jumps at region boundaries; Derivatives have been computed symbolically and are available for use by the Newton-Raphson method. The model"s functions can be measured from the terminals. It is also shown that small channel effects can be included in the model. Higher frequency effects can be modeled by using a network consisting of several sections of the basic lumped model. Future plans include a detailed comparison of the network model with models such as SPICE level 3 and a comparison of the multi- section higher frequency model with experiments.

Estimating stochastic volatility diffusion using conditional moments of integrated volatility

We exploit the distributional information contained in high-frequency intraday data in constructing a simple conditional moment estimator for stochastic volatility diffusions. The estimator is based on the analytical solutions of the first two conditional moments for the latent integrated volatility, the realization of which is effectively approximated by the sum of the squared high-frequency increments of the process. Our simulation evidence indicates that the resulting GMM estimator is highly reliable and accurate. Our empirical implementation based on high-frequency five-minute foreign exchange returns suggests the presence of multiple latent stochastic volatility factors and possible jumps. © 2002 Elsevier Science B.V. All rights reserved.

Effect of different jump distributions on the dynamics of jump processes

The paper investigates stochastic processes forced by independent and identically distributed jumps occurring according to a Poisson process. The impact of different distributions of the jump amplitudes are analyzed for processes with linear drift. Exact expressions of the probability density functions are derived when jump amplitudes are distributed as exponential, gamma, and mixture of exponential distributions for both natural and reflecting boundary conditions. The mean level-crossing properties are studied in relation to the different jump amplitudes. As an example of application of the previous theoretical derivations, the role of different rainfall-depth distributions on an existing stochastic soil water balance model is analyzed. It is shown how the shape of distribution of daily rainfall depths plays a more relevant role on the soil moisture probability distribution as the rainfall frequency decreases, as predicted by future climatic scenarios. © 2010 The American Physical Society.

Exploring the effects of image persistence in low frame rate virtual environments

© 2015 IEEE.In virtual reality applications, there is an aim to provide real time graphics which run at high refresh rates. However, there are many situations in which this is not possible due to simulation or rendering issues. When running at low frame rates, several aspects of the user experience are affected. For example, each frame is displayed for an extended period of time, causing a high persistence image artifact. The effect of this artifact is that movement may lose continuity, and the image jumps from one frame to another. In this paper, we discuss our initial exploration of the effects of high persistence frames caused by low refresh rates and compare it to high frame rates and to a technique we developed to mitigate the effects of low frame rates. In this technique, the low frame rate simulation images are displayed with low persistence by blanking out the display during the extra time such image would be displayed. In order to isolate the visual effects, we constructed a simulator for low and high persistence displays that does not affect input latency. A controlled user study comparing the three conditions for the tasks of 3D selection and navigation was conducted. Results indicate that the low persistence display technique may not negatively impact user experience or performance as compared to the high persistence case. Directions for future work on the use of low persistence displays for low frame rate situations are discussed.

Observation of stable phase jump lines in convection of a twisted nematic liquid crystal

We report observations of stable, localized, linelike structures in the spatially periodic pattern formed by nematic electroconvection, along which the phase of the pattern jumps by pi. With increasing electric voltage, these lines form a gridlike structure that goes over into a structure indistinguishable from the well-known grid pattern. We present theoretical arguments that suggest that the twisted cell geometry we are using is indirectly stabilizing the phase jump lines, and that the phase jump lines lattice is caused by an interaction of phase jump lines and a zig-zag instability of the surrounding pattern.