Perturbational finite volume scheme for the one-dimensional Navier-Stokes equations

Starting from the second-order finite volume scheme,though numerical value perturbation of the cell facial fluxes, the perturbational finite volume (PFV) scheme of the Navier-Stokes (NS) equations for compressible flow is developed in this paper. The central PFV scheme is used to compute the one-dimensional NS equations with shock wave.Numerical results show that the PFV scheme can obtain essentially non-oscillatory solution.

On the contact width in generalized plain strain JKR model for isotropic solids

Adhesive contact model between an elastic cylinder and an elastic half space is studied in the present paper, in which an external pulling force is acted on the above cylinder with an arbitrary direction and the contact width is assumed to be asymmetric with respect to the structure. Solutions to the asymmetric model are obtained and the effect of the asymmetric contact width on the whole pulling process is mainly discussed. It is found that the smaller the absolute value of Dundurs' parameter beta or the larger the pulling angle theta, the more reasonable the symmetric model would be to approximate the asymmetric one.

Extraction kinetics of Sc(III), Y(III), La(III) and Gd(III) from chloride medium by Cyanex 302 in heptane using the constant interfacial cell with laminar flow

The extraction kinetics of Sc, Y, La and Gd(III) from the hydrochloric acid medium using Cyanex 302 (hereafter HL) in heptane solution have been measured by the constant interfacial cell with laminar flow. Reaction regions are explored at liquid-liquid interface. Extraction regimes are deduced to be diffusion-controlled for Sc(Ill) and mixed controlled for Y, La and Gd(Ill). Extraction mechanisms are discussed according to the dimeric model of Cyanex 302 in non-polar solution. From the temperature dependence of rate measurement, the values of E-a, Delta H-+/-, Delta S-+/- and Delta G(300)(+/-) are calculated and it is found that the absolute values of these parameters keep crescent trend for Sc, Y, La and Gd(III). At the same time, it is found that it can easily achieve the mutual separation among the Sc, Y and La(III) with kinetics extraction methods.

电阻率测井反演方法应用研究

Formation resistivity is one of the most important parameters to be evaluated in the evaluation of reservoir. In order to acquire the true value of virginal formation, various types of resistivity logging tools have been developed. However, with the increment of the proved reserves, the thickness of interest pay zone is becoming thinner and thinner, especially in the terrestrial deposit oilfield, so that electrical logging tools, limited by the contradictory requirements of resolution and investigation depth of this kinds of tools, can not provide the true value of the formation resistivity. Therefore, resitivity inversion techniques have been popular in the determination of true formation resistivity based on the improving logging data from new tools. In geophysical inverse problems, non-unique solution is inevitable due to the noisy data and deficient measurement information. I address this problem in my dissertation from three aspects, data acquisition, data processing/inversion and applications of the results/ uncertainty evaluation of the non-unique solution. Some other problems in the traditional inversion methods such as slowness speed of the convergence and the initial-correlation results. Firstly, I deal with the uncertainties in the data to be processed. The combination of micro-spherically focused log (MSFL) and dual laterolog(DLL) is the standard program to determine formation resistivity. During the inversion, the readings of MSFL are regarded as the resistivity of invasion zone of the formation after being corrected. However, the errors can be as large as 30 percent due to mud cake influence even if the rugose borehole effects on the readings of MSFL can be ignored. Furthermore, there still are argues about whether the two logs can be quantitatively used to determine formation resisitivities due to the different measurement principles. Thus, anew type of laterolog tool is designed theoretically. The new tool can provide three curves with different investigation depths and the nearly same resolution. The resolution is about 0.4meter. Secondly, because the popular iterative inversion method based on the least-square estimation can not solve problems more than two parameters simultaneously and the new laterolog logging tool is not applied to practice, my work is focused on two parameters inversion (radius of the invasion and the resistivty of virgin information ) of traditional dual laterolog logging data. An unequal weighted damp factors- revised method is developed to instead of the parameter-revised techniques used in the traditional inversion method. In this new method, the parameter is revised not only dependency on the damp its self but also dependency on the difference between the measurement data and the fitting data in different layers. At least 2 iterative numbers are reduced than the older method, the computation cost of inversion is reduced. The damp least-squares inversion method is the realization of Tikhonov's tradeoff theory on the smooth solution and stability of inversion process. This method is realized through linearity of non-linear inversion problem which must lead to the dependency of solution on the initial value of parameters. Thus, severe debates on efficiency of this kinds of methods are getting popular with the developments of non-linear processing methods. The artificial neural net method is proposed in this dissertation. The database of tool's response to formation parameters is built through the modeling of the laterolog tool and then is used to training the neural nets. A unit model is put forward to simplify the dada space and an additional physical limitation is applied to optimize the net after the cross-validation method is done. Results show that the neural net inversion method could replace the traditional inversion method in a single formation and can be used a method to determine the initial value of the traditional method. No matter what method is developed, the non-uniqueness and uncertainties of the solution could be inevitable. Thus, it is wise to evaluate the non-uniqueness and uncertainties of the solution in the application of inversion results. Bayes theorem provides a way to solve such problems. This method is illustrately discussed in a single formation and achieve plausible results. In the end, the traditional least squares inversion method is used to process raw logging data, the calculated oil saturation increased 20 percent than that not be proceed compared to core analysis.

Effects of the orientation of microgravity on the concentration distribution in crystallization from solution

A quasi-steady state growth and dissolution in a 2-D rectangular enclosure is numerically investigated. This paper is an extension to indicate the effects of the orientation of gravity on the concentration field in crystallization from solution under microgravity, especially on the lateral non-uniformity of concentration distribution at the growth surface. The thermal and solute convection are included in this model.

Optical Diagnostic and Modeling Solution Growth Process of Sodium Chlorate Crystals

Both a real time optical interferometric experiment and a numerical simulation of two-dimension non-steady state model were employed to study the growth process of aqueous sodium chlorate crystals. The parameters such as solution concentration distribution, crystal dimensions, growth rate and velocity field were obtained by both experiment and numerical simulation. The influence of earth gravity during crystal growth process was analyzed. A reasonable theory model corresponding to the present experiment is advanced. The thickness of concentration boundary layer was investigated especially. The results from the experiment and numerical simulation match well.

A Continuation Method of Parameter Inversion for Non-Equilibrium Convection-Dispersion Equation

Based on the homotopy mapping, a globally convergent method of parameter inversion for non-equilibrium convection-dispersion equations (CDEs) is developed. Moreover, in order to further improve the computational efficiency of the algorithm, a properly smooth function, which is derived from the sigmoid function, is employed to update the homotopy parameter during iteration. Numerical results show the feature of global convergence and high performance of this method. In addition, even the measurement quantities are heavily contaminated by noises, and a good solution can be found.

Generalized Maugis-Dugdale Model of an Elastic Cylinder in Non-Slipping Adhesive Contact with a Stretched Substrate

We have recently developed a generalized JKR model for non-slipping adhesive contact between an elastic cylinder and a stretched substrate where both tangential and normal tractions are transmitted across the contact interface. Here we extend this model to a generalized Maugis-Dugdale model by adopting a Dugdale-type adhesive interaction law to eliminate the stress singularity near the edge of the contact zone. The non-slipping Maugis-Dugdale model is expected to have a broader range of validity in comparison with the non-slipping JKR model. The solution shares a number of common features with experimentally observed behaviors of cell reorientation on a cyclically stretched substrate.

Non-Slipping Jkr Model For Transversely Isotropic Materials

A generalized plane strain JKR model is established for non-slipping adhesive contact between an elastic transversely isotropic cylinder and a dissimilar elastic transversely isotropic half plane, in which a pulling force acts on the cylinder with the pulling direction at an angle inclined to the contact interface. Full-coupled solutions are obtained through the Griffith energy balance between elastic and surface energies. The analysis shows that, for a special case, i.e., the direction of pulling normal to the contact interface, the full-coupled solution can be approximated by a non-oscillatory one, in which the critical pull-off force, pull-off contact half-width and adhesion strength can be expressed explicitly. For the other cases, i.e., the direction of pulling inclined to the contact interface, tangential tractions have significant effects on the pull-off process, it should be described by an exact full-coupled solution. The elastic anisotropy leads to an orientation-dependent pull-off force and adhesion strength. This study could not only supply an exact solution to the generalized JKR model of transversely isotropic materials, but also suggest a reversible adhesion sensor designed by transversely isotropic materials, such as PZT or fiber-reinforced materials with parallel fibers. (c) 2007 Elsevier Ltd. All rights reserved.

On dispersion-controlled principles for non-oscillatory shock-capturing schemes

The role of dispersions in the numerical solutions of hydrodynamic equation systems has been realized for long time. It is only during the last two decades that extensive studies on the dispersion-controlled dissipative (DCD) schemes were reported. The studies have demonstrated that this kind of the schemes is distinct from conventional dissipation-based schemes in which the dispersion term of the modified equation is not considered in scheme construction to avoid nonphysical oscillation occurring in shock wave simulations. The principle of the dispersion controlled aims at removing nonphysical oscillations by making use of dispersion characteristics instead of adding artificial viscosity to dissipate the oscillation as the conventional schemes do. Research progresses on the dispersion controlled principles are reviewed in this paper, including the exploration of the role of dispersions in numerical simulations, the development of the dispersion-controlled principles, efforts devoted to high-order dispersion-controlled dissipative schemes, the extension to both the finite volume and the finite element methods, scheme verification and solution validation, and comments on several aspects of the schemes from author's viewpoint.

Concentration distribution in crystallization from solution under microgravity

Concentration distribution in crystallization from solution under microgravity is numerically studied. A quasi-steady state growth and dissolution in a 2D rectangular enclosure filled with sodium chlorate (NaClO3) aqueous solution, in which one wall is the growth surface of the crystal and the opposite one is the dissolution surface, is considered. The solute transport process at the growth surface is described by the diffusion-reaction theory with finite interface kinetics coefficient. The results show that the concentration at the growth surface is supersaturated and the supersaturation distribution is of non-uniformity, i.e. the supersaturation in a region facing an incoming flow is high. On the other hand, the non-uniformity of supersaturation at the growth surface is closely related to the gravity level even under microgravity, it exponentially increases as the thermal Rayleigh number on behalf of the gravity level rises.

The effect of the surface tension on non-propagating solitary waves

In this paper, the effect of the surface tension is considered carefully in the study of non-propagating solitary waves. The parameter plane of the surface tension and the fluid depth is divided into three regions; in two of them a breather soliton can be produced. In literature the parameters of breather solitons are all in one of the parameter regions. The new region reported here has been confirmed by our experiments. In the third region, the theoretical solution is a kink soliton, but a kind of the non-propagating solitary wave similar to the breather soliton was found in our experiments besides the kink soliton.

Numerical Solution for Polymer Transient Flows in a Circle Bounded Composite Formation

A new numerical model for transient flows of polymer solution in a circular bounded composite formation is presented in this paper. Typical curves of the wellbore transient pressure are yielded by FEM. The effects of non-Newtonian power-law index, mobility and boundary distance have been considered. It is found that for the mobility ratio larger than 1, which is favorable for the polymer flooding, the pressure derivative curve in log-log form rises up without any hollow. On the other hand, if the pressure derivative curve has a hollow and then is raised up, we say that the polymer flooding fails. Finally, the new model has been extended to more complicated boundary case.

Orientation-Dependent Adhesion Strength Of A Rigid Cylinder In Non-Slipping Contact With A Transversely Isotropic Half-Space

Recently, Chen and Gao [Chen, S., Gao, H., 2007. Bio-inspired mechanics of reversible adhesion: orientation-dependent adhesion strength for non-slipping adhesive contact with transversely isotropic elastic materials. J. Mech. Phys. solids 55, 1001-1015] studied the problem of a rigid cylinder in non-slipping adhesive contact with a transversely isotropic solid subjected to an inclined pulling force. An implicit assumption made in their study was that the contact region remains symmetric with respect to the center of the cylinder. This assumption is, however, not self-consistent because the resulting energy release rates at two contact edges, which are supposed to be identical, actually differ from each other. Here we revisit the original problem of Chen and Gao and derive the correct solution by removing this problematic assumption. The corrected solution provides a proper insight into the concept of orientation-dependent adhesion strength in anisotropic elastic solids. (c) 2008 Elsevier Ltd. All rights reserved.