基于C#的二维核磁共振测井观测模式设计软件的开发

二维核磁共振测井观测模式是以获取特定地层信息为目标的数据采集方式，它直接决定着核磁共振谱仪对不同类型储层的适应性以及获取原始数据的可信度．本文从观测模式的组成元素以及工作机理出发，采用C#语言开发了一套具有可视化功能的二维核磁共振测井观测模式设计软件．该软件支持二维核磁共振测井以及核磁共振岩心分析仪的观测模式设计，同时提供了三种采集参数编辑方式、观测模式优化机制以及图形化显示等功能，实现了观测模式的灵活调整，解决了当前主流二维核磁共振测井观测模式中存在的采集参数固定、不宜灵活调整从而导致所采集的原始数据信噪比低下、储层中流体弛豫特性加载不完全等问题．此外，软件还提供了观测模式的采集参数信息以及采集时序信息的输出功能等．     

Computing communities in complex networks using the Dirichlet processing Gaussian mixture model with spectral clustering

Community detection becomes a significant tool for the complex network analysis. The study of the community detection algorithms has received an enormous amount of attention. It is still an open question whether a highly accurate and efficient algorithm is found in most data sets. We propose the Dirichlet Processing Gaussian Mixture Model with Spectral Clustering algorithm for detecting the community structures. The combination of traditional spectral algorithm and new non-parametric Bayesian model provides high accuracy and quality. We compare the proposed algorithm with other existing community detecting algorithms using different real-world data sets and computer-generated synthetic data sets. We show that the proposed algorithm results in high modularity, and better accuracy in a wide range of networks. We find that the proposed algorithm works best for the large size of the data sets.     

Design and fabrication of fluid controlled dynamic optical lens system

A novel small fluid controlled optical lens system that is capable of displaying dynamic variation of its focal length and field-of-view (FOV) is designed and fabricated. In this active lens system, appropriate volume of the optical fluid can be pumped into or out of the lens chamber to provide double-convex (DCX) or double-concave (DCV) lens effect. Simple optical imaging experiments were performed using different sets of glass lenses with fixed focal lengths to determine the optimum lens configuration required for designing a dynamic optical lens system. The experimental results obtained from the glass lenses demonstrate that a combination of a single DCX lens with three DCV lenses provides a wider FOV. The flexible membranes for fluid controlled lenses were fabricated using polydimethylsiloxane polymer material, which has good optical transparency and elasticity. A simple fluid injection system is used to vary the radius of curvature of the lenses, and thereby to change the focal length. A dynamic optical lens system with a combination of one DCX and multiple variable focal length DCV lenses as designed here can image an object with a wide range of focal length and FOV. With this fluid controlled optical system, the FOV and focal length could be continuously varied and a maximum FOV of 118.3° could be achieved. The smallest f-number (f_/#) for this fluid controlled single lens system was found to be 1.3, which corresponds to the numerical aperture value of 0.35.     

Dynamical van der Waals Model of Glassy Behavior

A model fluid system whose intermolecular forces consist of the short range part characterizing usual stable fluid and the long range part of Kac-type is considered. When the ratio of the force ranges is large enough, the system can be described by regarding the system having only short range interactions as a reference fluid treated as a continuum, on which long range forces acting among tiny fluid elements are superimposed. We discuss the glassy behavior of this model relating it to the mode coupling theory and using real replica theory. These theories lead to the two equations for non-ergodicity parameter which are totally different from each other. We argue that our model can be a basis for examining nature of the drastic approximations entering derivations of the mode coupling equations. We further explore the possibility of developing the dynamical real replica approach for our model system with the hope of providing a framework to cope with different time scales characterizing complex glassy behavior.     

方阱链流体在固液界面分布的密度泛函理论研究

应用Yethiraj的加权密度近似泛函理论研究平板狭缝中方阱链流体的密度分布，系统的Helm holtz自由能泛函分为理想气体的贡献利剩余贡献两部分，其中剩余贡献部分分别采用刘洪 来等人建立的基于空穴相关函数的方阱链流体状态方程和Gil-Villegas等人提出的统计缔合 流体理论状态方程(SAFT-VR)结合简单加权密度近似计算.考察了不同链长、温度、系统密度 和壁面吸引强度下平板狭缝中方阱链流体的密度分布，并与Monte Carlo(MC)模拟结果进行 了比较.结果表明采用不同的状态方程对密度分布的计算有明显的影响，对于受限于硬壁狭 缝中的方阱链流体，温度和密度比较高时，两种状态方程计算的结果均与MC模拟符合得比较 好，在低温和低密度下效果变差，SAFT-VR方程的计算结果更接近于MC模拟结果.对于受限于 方阱壁狭缝中的方阱链流体，由于系统密度分布的非均匀性加强，采用两种状态方程计算的 结果均与MC模拟结果有一定偏差，寻找更合适的权重函数是进一步改进的关键. 关键词： 密度泛函理论 非均匀流体 密度分布 固液界面 方阱链     

Scattering from a partially fluid-filled, elastic-shelled sphere

In this article a sphere is taken to be partially filled with fluid so that its interior is part fluid and part air. A set of basis of functions, based upon an origin at the fluid/air interface, is used for the interior and a set of basis functions based upon the center of the sphere is used for the shell and exterior of the sphere. These sets of basis functions are coupled at the shell/interior interface and the resulting coupled system of equations solved to yield the scattered field. Numerical computations using this approach are presented for varying amounts of fluid-fill and for varying incident plane waves.     

Ultrasoft pseudopotentials and projector augmented-wave data sets: application to diatomic molecules

We test several ultrasoft pseudopotentials (US-PPs) and projector augmented-wave (PAW) data sets, calculating the bond lengths, the atomization energies and the frequencies of the vibrational stretch modes of various diatomic molecules. The US-PPs and the PAW data sets are constructed with the same recipe and using the local density approximation or the Perdew, Burke and Ernzerhof generalized gradient approximation for the exchange and correlation energies. We study the dimers H(2), N(2), O(2), F(2), Al(2), Si(2), P(2), S(2) and Cl(2) and several monohydrides, carbides, nitrides and oxides of boron, carbon, nitrogen, oxygen, fluorine, aluminum, silicon, phosphorus, sulfur, chlorine, iron and nickel. We find that US-PPs and PAW data sets constructed with the same parameters provide almost equivalent results for the bond lengths and the vibrational stretch frequencies while, for some molecules, the PAW method is superior to the US-PP method for the calculation of the atomization energies. Our geometries and vibrational frequencies are compared with the results present in the literature and obtained by localized basis sets. It is found that the agreement is very good, with discrepancies comparable to those due to the use of different localized basis sets.     

Hausdorff dimension of attractors for the two-dimensional navier-stokes equations with boundary conditions

We consider a viscous incompressible fluid enclosed in a region of IR², and subject to boundary conditions. We obtain explicit bounds (depending only on the data) for the entropy (Kolmogorov-Sinaï invariant) and the Hausdorff dimension of attracting sets.This work has been mostly performed at Centre de Physique Théorique, Ecole Polytechnique, F-91128 Palaiseau Cédex     

A new approach to design a control system for a FGR furnace using the combination of the CFD and linear system identification techniques

A new approach to design control systems for an industrial furnace with flue gas recirculation (FGR) is presented. To facilitate the control system design, a linear dynamic model is needed for the furnace. Full-scale computational fluid dynamics (CFD) simulations are used to generate the required small signal input and output data sets. Subsequently, a least squares based system identification technique is used to obtained the linear dynamic models. After model validation, feedback controller is designed based on these linear dynamic models. Finally, the performance of the designed closed-loop control system is also evaluated using both linear dynamic model and full-scale nonlinear CFD model. The comparison shows that the control system designed using the proposed approach can minimize the deviation of nitric oxides (NO) emission from the design point by minimize the dynamic NO formation, hence to prevent any excessive NO formation in the combustion process when the system subjects to disturbances.     

Oscillatory rise of bubbles in wormlike micellar fluids with different microstructures

Previous observations of the nontransient oscillations of rising bubbles and falling spheres in wormlike micellar fluids were limited to a single surfactant system. We present an extensive survey of rising bubbles in another system, an aqueous solution of cetylpyridinium chloride and sodium salicylate, with and without NaCl, across a range of concentrations and temperatures. Two different types of oscillations are seen in different concentration ranges, each with its own temperature dependence. Rheological data identify these different hydrodynamic states with different fluid microstructures.     

Automated assessment of the quality of diffusion tensor imaging data using color cast of color-encoded fractional anisotropy images

Diffusion tensor imaging (DTI) data often suffer from artifacts caused by motion. These artifacts are especially severe in DTI data from infants, and implementing tight quality controls is therefore imperative for DTI studies of infants. Currently, routine procedures for quality assurance of DTI data involve the slice-wise visual inspection of color-encoded, fractional anisotropy (CFA) images. Such procedures often yield inconsistent results across different data sets, across different operators who are examining those data sets, and sometimes even across time when the same operator inspects the same data set on two different occasions. We propose a more consistent, reliable, and effective method to evaluate the quality of CFA images automatically using their color cast, which is calculated on the distribution statistics of the 2D histogram in the color space as defined by the International Commission on Illumination (CIE) on lightness and a and b (LAB) for the color-opponent dimensions (also known as the CIELAB color space) of the images. Experimental results using DTI data acquired from neonates verified that this proposed method is rapid and accurate. The method thus provides a new tool for real-time quality assurance for DTI data.     

Benchmark calculations for Monte Carlo simulations of electrontransport

Benchmark calculations have been performed for electron transport coefficients with an aim to produce a body of data required to verify the codes used in plasma modeling. The present code for the time resolved Monte Carlo simulation (MCS) was shown to represent properly DC transport coefficients in a purely electric field, in crossed electric and magnetic fields, and in the presence of nonconservative collisions, ionization, and attachment. In addition, we have suggested tests of the time dependent solutions. Relaxation of the initial transport coefficient may serve as an accurate test of the code as well as the input data for some fluid codes. In this paper, we show only one example, but several different sets of conditions and cross sections should be used as well. Finally, we propose application of the quasi-steady state results in RF fields. As an example we suggest calculation of the components of diffusion tensor showing anomalous longitudinal diffusion and calculations made with nonconservative collisions (ionization in this case). We also check the application of approximate formulas to determine drift velocity on the basis of total collision frequency and to determine a diffusion coefficient by using the Einstein relation. Other tests required to verify the transport data calculations are discussed as well     

The inversion of NMR log data sets with different measurement errors.

We present a composite-data processing method which simultaneously processes two or more data sets with different measurement errors. We examine the role of the noise level of the data in the singular value decomposition inversion process, the criteria for a proper cutoff, and its effect on the uncertainty of the solution. Examples of processed logs using the composite-data processing method are presented and discussed. The possible usefulness of the apparent T(1)/T(2) ratio extracted from the logs is illustrated.     

Calculating the dimension of attractors from small data sets

Calculations of the order-2 information dimension for data sets generated from numerical simulations as well as from an experiment show that it is possible to distinguish between periodic, chaotic and random signals, as well as to characterize different kinds of chaotic behavior using relatively small, noisy data sets.     

Nuclear corrections in neutrino-nucleus deep inelastic scattering and their compatibility with global nuclear parton-distribution-function analyses

We perform a global χ2 analysis of nuclear parton distribution functions using data from charged current neutrino-nucleus (νA) deep-inelastic scattering (DIS), charged-lepton-nucleus (?(±)A) DIS, and the Drell-Yan (DY) process. We show that the nuclear corrections in νA DIS are not compatible with the predictions derived from ?(±)A DIS and DY data. We quantify this result using a hypothesis-testing criterion based on the χ2 distribution which we apply to the total χ2 as well as to the χ2 of the individual data sets. We find that it is not possible to accommodate the data from νA and ?(±)A DIS by an acceptable combined fit. Our result has strong implications for the extraction of both nuclear and proton parton distribution functions using combined neutrino and charged-lepton data sets.     

Spatial and temporal resolution effects on dynamic contrast-enhanced magnetic resonance mammography

We tested the hypothesis that partial volume effects due to poor in-plane resolution and/or low temporal resolution used in clinical dynamic contrast-enhanced magnetic resonance imaging results in erroneous diagnostic information based on inaccurate estimates of tumor contrast agent extravasation and tested whether reduced encoding techniques can correct for dynamic data volume averaging. Image spatial resolution was reduced from 469 x 469 microm2 to those reported below by selecting a subset of k-space data. We then compared the top five K(trans)/V(T) "hot spots" obtained from the original data set, 469 x 469-microm in-plane spatial resolution and an 18-s temporal resolution processed by fast Fourier transform (FFT), with values obtained from data sets having in-plane spatial resolutions of 938 x 938, 1875 x 1875 and 2500 x 2500 microm2 and a temporal resolution of 18 s, or data sets with temporal resolutions of 36, 54 and 72 and a spatial resolution of 469 x 469 microm2, and found them to statistically differ from the parent data sets. We then tested four different post processing methods for improving the spatial resolution without sacrificing temporal resolution: zero-filled FFT, keyhole, reduced-encoding imaging by generalized-series reconstruction (RIGR) and two-reference RIGR (TRIGR). The top five values of K(trans)/V(T) obtained from data sets, the in-plane spatial resolutions of which were improved to 469 x 469 microm2 by zero-filling FFT, Keyhole and RIGR, statistically differed from those obtained from the original 469 x 469 microm2 FFT parent image data set. Only the 938 x 938 and 1875 x 1875 microm2 data sets reconstructed to 469 x 469 microm2 with TRIGR reconstruction method yielded values of the top five K(trans)/V(T) hot spots statistically the same as the original parent data set, 469 x 469 microm2 in-plane spatial and 18-s temporal-resolution FFT. That is, partial volume effects from data sets of different in-plane spatial resolution resulted in statistically different values of the top five K(trans)/V(T) hot spots relative to a high spatial and temporal resolution data set, and TRIGR reconstruction of these low resolution data sets to high resolution images provided statistically similar values with a savings in temporal resolution of 2 to 4 times.     

Adhesion and fingering in the lifting Hele-Shaw cell: Role of the substrate

The lifting Hele-Shaw cell (LHSC) is used to study adhesion as well as viscous fingering. In the present paper we report a series of observations of development of the interface for different viscous fluids, both Newtonian and non-Newtonian, in a LHSC operated at a constant lifting force. Glass and perspex are used as the plates in two different sets of experiments. The objectives are 1) to measure the time required to separate the plates as a function of the lifting force and 2) to note the force above which viscous fingering appears. We find that for the Newtonian fluids, the plate separation time follows a universal power law with the lifting force, irrespective of fluid and substrate. The non-Newtonian fluids too, with proper scaling obey the same power law. The appearance of fingering, however, depends on the properties of the fluid as well as the substrate. We suggest a modified form of the capillary number which controls the onset of fingering; this new quantity, termed the “fingering parameter” involves the dielectric constants of the substrate and fluid in addition to the viscosity and surface tension.     

Further studies in aesthetic field theory IV

We are able to obtain a bounded particle, with no indication of a singularity appearing, in several ways different from our previous papers. For one set of data we find slightly greater structure (more turnabout points) than previously. We discuss some of the properties of ten different sets of data.     

Interpolating Strange Attractors via Fractional Brownian Bridges

We present a novel method for interpolating univariate time series data. The proposed method combines multi-point fractional Brownian bridges, a genetic algorithm, and Takens’ theorem for reconstructing a phase space from univariate time series data. The basic idea is to first generate a population of different stochastically-interpolated time series data, and secondly, to use a genetic algorithm to find the pieces in the population which generate the smoothest reconstructed phase space trajectory. A smooth trajectory curve is hereby found to have a low variance of second derivatives along the curve. For simplicity, we refer to the developed method as PhaSpaSto-interpolation, which is an abbreviation for phase-space-trajectory-smoothing stochastic interpolation. The proposed approach is tested and validated with a univariate time series of the Lorenz system, five non-model data sets and compared to a cubic spline interpolation and a linear interpolation. We find that the criterion for smoothness guarantees low errors on known model and non-model data. Finally, we interpolate the discussed non-model data sets, and show the corresponding improved phase space portraits. The proposed method is useful for interpolating low-sampled time series data sets for, e.g., machine learning, regression analysis, or time series prediction approaches. Further, the results suggest that the variance of second derivatives along a given phase space trajectory is a valuable tool for phase space analysis of non-model time series data, and we expect it to be useful for future research.     

Validation subset selections for extrapolation oriented QSPAR models

One of the most important features of QSPAR models is their predictive ability. The predictive ability of QSPAR models should be checked by external validation. In this work we examined three different types of external validation set selection methods for their usefulness in in-silico screening. The usefulness of the selection methods was studied in such a way that: 1) We generated thousands of QSPR models and stored them in 'model banks'. 2) We selected a final top model from the model banks based on three different validation set selection methods. 3) We predicted large data sets, which we called 'chemical universe sets', and calculated the corresponding SEPs. The models were generated from small fractions of the available water solubility data during a GA Variable Subset Selection procedure. The external validation sets were constructed by random selections, uniformly distributed selections or by perimeter-oriented selections. We found that the best performing models on the perimeter-oriented external validation sets usually gave the best validation results when the remaining part of the available data was overwhelmingly large, i.e., when the model had to make a lot of extrapolations. We also compared the top final models obtained from external validation set selection methods in three independent and different sizes of 'chemical universe sets'.