A Fourth-Order Unstructured NURBS-Enhanced Finite Volume WENO Scheme for Steady Euler Equations in Curved Geometries

In Li and Ren(Int.J.Numer.Methods Fluids 70:742–763,2012),a high-order k-exact WENO finite volume scheme based on secondary reconstructions was proposed to solve the two-dimensional time-dependent Euler equations in a polygonal domain,in which the high-order numerical accuracy and the oscillations-free property can be achieved.In this paper,the method is extended to solve steady state problems imposed in a curved physical domain.The numerical framework consists of a Newton type finite volume method to linearize the nonlinear governing equations,and a geometrical multigrid method to solve the derived linear system.To achieve high-order non-oscillatory numerical solutions,the classical k-exact reconstruction with k=3 and the efficient secondary reconstructions are used to perform the WENO reconstruction for the conservative variables.The non-uniform rational B-splines(NURBS)curve is used to provide an exact or a high-order representation of the curved wall boundary.Furthermore,an enlarged reconstruction patch is constructed for every element of mesh to significantly improve the convergence to steady state.A variety of numerical examples are presented to show the effectiveness and robustness of the proposed method.

CALCULATIONS OF THE PHASE EQUILIBRIA AND CRITICAL CURVES WITH THE EQUATION OF STATE FOR BINARY MIXTURES

A rational equation of state of the perturbation type with a repulsion and attraction term has been applied to reproduce critical curves of six different binary systems up to high temperatures and pressures. A square well potential for intermolecular interaction is used. With pairwise combination rules for these potentials three adjustable parameters are needed. The experimental critical point and phase equilibrium data are compared with the values predicted using the equation of state. Good agreement is obtained for the analysis of the critical pressure composition data and molar volumes.

Furnace Temperature Curve Optimization Model Based on Differential Evolution Algorithm

When soldering electronic components onto circuit boards,the temperature curves of the reflow ovens across different zones and the conveyor belt speed significantly influence the product quality.This study focuses on optimizing the furnace temperature curve under varying settings of reflow oven zone temperatures and conveyor belt speeds.To address this,the research sequentially develops a heat transfer model for reflow soldering,an optimization model for reflow furnace conditions using the differential evolution algorithm,and an evaluation and decision model combining the differential evolution algorithm with the Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)method.This approach aims to determine the optimal furnace temperature curve,zone temperatures of the reflow oven,and the conveyor belt speed.

ANALYTIC INVARIANT CURVES OF A NONLINEAR SECOND ORDER DIFFERENCE EQUATION

This article studies the existence of analytic invariant curves for a nonlinear second order difference equation which was modeled from macroeconomics of the business cycle. The author not only discusses the case of the eigenvalue off the unit circle S^1 and the case on S^1 with the Diophantine condition but also considers the case of the eigenvalue at a root of the unity, which obviously violates the Diophantine condition.

Improved Nonlinear Equation Method for Numerical Prediction of Jominy End-Quench Curves

Without considering the effects of alloying interaction on the Jominy end-quench curves, the prediction resuits obtained by YU Bai-hai＇s nonlinear equation method for multi-alloying steels were different from those experimental ones reported in literature. Some alloying elements have marked influence on Jominy end-quench curves of steels. An improved mathematical model for simulating the Jominy end-quench curves is proposed by introducing a parameter named alloying interactions equivalent （Le）. With the improved model, the Jominy end-quench curves of steels so obtained agree very well with the experimental ones.

CONTINUOUS DEPENDENCE OF THE SOLUTIONS OF IMPULSIVE DIFFERENTIAL EQUATIONS ON THE INITIAL CONDITIONS AND BARRIER CURVES

The basic objects of investigation in this article are nonlinear impulsive dif- ferential equations. The impulsive moments coincide with the moments of meeting of the integral curve and some of the so-called barrier curves. For such type of equations, suf- ficient conditions are found under which the solutions are continuously dependent on the perturbations with respect to the initial conditions and barrier curves. The results are applied to a mathematical model of population dynamics.

Modified non-rectangular hyperbola equation with plant height for photosynthetic light-response curves of Potentilla anserina and Elymus nutans at various growth phases in the Heihe River Basin, Northwest China

The non-rectangular hyperbola(NRH)equation is the most popular method that plots the photosynthetic light-response(PLR)curve and helps to identify plant photosynthetic capability.However,the PLR curve can't be plotted well by the NRH equation at different plant growth phases due to the variations of plant development.Recently,plant physiological parameters have been considered into the NRH equation to establish the modified NRH equation,but plant height(H),an important parameter in plant growth phases,is not taken into account.In this study,H was incorporated into the NRH equation to establish the modified NRH equation,which could be used to estimate photosynthetic capability of herbage at different growth phases.To explore photosynthetic capability of herbage,we selected the dominant herbage species Potentilla anserina L.and Elymus nutans Griseb.in the Heihe River Basin,Northwest China as the research materials.Totally,twenty-four PLR curves and H at different growth phases were measured during the growing season in 2016.Results showed that the maximum net photosynthetic rate and the initial slope of PLR curve linearly increased with H.The modified NRH equation,which is established by introducing H and an H-based adjustment factor into the NRH equation,described better the PLR curves of P.anserina and E.nutans than the original ones.The results may provide an effective method to estimate the net primary productivity of grasslands in the study area.

Joule-Thomson inversion curve prediction by using equation of state

In this paper five equations of state are tested for checking their ability to predict the Joule-Thomson inversion curve.These five equations of state are:Mohsennia-Modarres-Mansoori(MMM),Ji-Lemp(JL),modified Soave-Redlich-Kwang(SRK)equation of state by Graboski(MSRK1),modified SRK equation of state by Peneloux and Rauzy(MSRK2),and modified Peng-Robinson (PR)equation of state by Rauzy(PRmr).The investigated equations of state give good prediction of the low-temperature branch of the inversion curve,except for MMM equation of state.The high-temperature branch and the peak of the inversion curve have been observed,in general,to be sensitive to the applied equation of state.The values of the maximum inversion temperature and maximum inversion pressure are calculated for each component used in this work.

BIFURCATIONS OF INVARIANT CURVES OF A DIFFERENCE EQUATION

Bifurcation of the invariant curves of a difference equation is studied. The system defined by the difference equation is integrable, so the study of the invariant curves of the difference system can become the study of topological classification of the planar phase portraits defined by a planar Hamiltonian system. By strict qualitative analysis, the classification of the invariant curves in parameter space can be obtained.

Investigations on if-E Curve in Cyclic Derivative Chronopotentiometry(Ⅵ)——Theoretical Equations of if-E Curves in the Case of Quasi-reversible and Irreversible Electrode Reactions

The theory of if-E curve in cyclic derivative chronopotentiometry is presented. Theoretical equations of if-E curves in the case of quasi-reversible and irreversible electrode reactions are deduced respectively.

Riemann-Hilbert problem for first order complex equations of mixed type with degenerate curve

This paper considers the Riemann-Hilbert problem for linear mixed（elliptichyperbolic） complex equations of first order with degenerate curve in a simply connected domain. We first give the representation theorem and uniqueness of solutions for such boundary value problem. Then by using the methods of successive iteration and parameter extension, the existence of solutions for this problem is proved.

Dynamic equations for curved submerged floating tunnel

In virtue of reference Cartesian coordinates, geometrical relations of spatial curved structure are presented in orthogonal curvilinear coordinates. Dynamic equations for helical girder are derived by Hamilton principle. These equations indicate that four generalized displacements are coupled with each other. When spatial structure degenerates into planar curvilinear structure, two generalized displacements in two perpendicular planes are coupled with each other. Dynamic equations for arbitrary curvilinear structure may be obtained by the method used in this paper.

Computer program of nonlinear, curved regression for ‘probacent’-probability equation in biomedicine

On the basis of experimental observations on animals, applications to clinical data on patients and theoretical statistical reasoning, the author developed a com-puter-assisted general mathematical model of the ‘probacent’-probability equation, Equation (1) and death rate (mortality probability) equation, Equation (2) derivable from Equation (1) that may be applica-ble as a general approximation method to make use-ful predictions of probable outcomes in a variety of biomedical phenomena [1-4]. Equations (1) and (2) contain a constant, γ and c, respectively. In the pre-vious studies, the author used the least maximum- difference principle to determine these constants that were expected to best fit reported data, minimizing the deviation. In this study, the author uses the method of computer-assisted least sum of squares to determine the constants, γ and c in constructing the ‘probacent’-related formulas best fitting the NCHS- reported data on survival probabilities and death rates in the US total adult population for 2001. The results of this study reveal that the method of com-puter-assisted mathematical analysis with the least sum of squares seems to be simple, more accurate, convenient and preferable than the previously used least maximum-difference principle, and better fit-ting the NCHS-reported data on survival probabili-ties and death rates in the US total adult population. The computer program of curved regression for the ‘probacent’-probability and death rate equations may be helpful in research in biomedicine.

Schrdinger Equation of a Particle on a Rotating Curved Surface

We derive the Schr6dinger equation of a particle constrained to move on a rotating curved surface S. Using the thin-layer quantization scheme to confine the particle on S, and with a proper choice of gauge transformation for the wave function, we obtain the well-known geometric potentiM Vg and an additive Coriolis-induced geometric potential in the co-rotationM curvilinear coordinates. This novel effective potential, which is included in the surface Schr6dinger equation and is coupled with the mean curvature of S, contains an imaginary part in the general case which gives rise to a non-Hermitian surface Hamiltonian. We find that the non-Hermitian term vanishes when S is a minimal surface or a revolution surface which is axially symmetric around the rolling axis.

SPEED UP RATIONAL POINT SCALAR MULTIPLICATIONS ON ELLIPTIC CURVES BY FROBENIUS EQUATIONS

Let q be a power of a prime and φ be the Frobenius endomorphism on E（Fqk）, then q = tφ - φ^2. Applying this equation, a new algorithm to compute rational point scalar multiplications on elliptic curves by finding a suitable small positive integer s such that q^s can be represented as some very sparse φ-polynomial is proposed. If a Normal Basis （NB） or Optimal Normal Basis （ONB） is applied and the precomputations are considered free, our algorithm will cost, on average, about 55% to 80% less than binary method, and about 42% to 74% less than φ-ary method. For some elliptic curves, our algorithm is also taster than Mǖller＇s algorithm. In addition, an effective algorithm is provided for finding such integer s.

Description of Evaporation Curve for Liquid Metals by the Generalized Van-der-Waals-Berthelot Equation. Part II

In the present work it is shown that the generalized Van-der-Waals-Berthelot equation describes the evaporation curves （saturation curves） for alkali metals with good accuracy. This result is obtained on the basis of the calculations performed by the author for thermodynamic parameters of the saturation curves described by the generalized Van-der-Waals-Berthelot equation.

Finite-difference time-domain modeling of curved material interfaces by using boundary condition equations method

To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.

Characterization of reservoir properties and pore structure based on micro-resistivity imaging logging: Porosity spectrum, permeability spectrum, and equivalent capillary pressure curve

According to the capillary theory,an equivalent capillary model of micro-resistivity imaging logging was built.On this basis,the theoretical models of porosity spectrum(Ф_(i)),permeability spectrum(K_(i))and equivalent capillary pressure curve(pe)were established to reflect the reservoir heterogeneity.To promote the application of the theoretical models,the Archie's equation was introduced to establish a general model for quantitatively characterizing bi,K,and pei.Compared with the existing models,it is shown that:(1)the existing porosity spectrum model is the same as the general equation of gi;(2)the Ki model can display the permeability spectrum as compared with Purcell's permeability model;(3)the per model is constructed on a theoretical basis and avoids the limitations of existing models that are built only based on the component of porosity spectrum,as compared with the empirical model of capillary pressure curve.The application in the Permian Maokou Formation of Well TsX in the Central Sichuan paleo-uplift shows that the Ф_(i),K_(i),and p_(ci) models can be effectively applied to the identification of reservoir types,calculation of reservoir properties and pore structure parameters,and evaluation of reservoir heterogeneity.

Establishment and application of logging saturation interpretation equation in vuggy reservoirs

Vuggy reservoirs are the most common, albeit important heterogeneous carbonate reservoirs in China. However, saturation calculations using logging data are not well developed, whereas Archie method is more common. In this study, electrical conduction in a vuggy reservoir is theoretically analyzed to establish a new saturation equation for vuggy reservoirs. We found that vugs have a greater effect on saturation than resistivity, which causes inflection in the rock-electricity curve. Using single-variable exPeriments, we evaluated the effects of rug size, vug number, and vug distribution on the rock-electricity relation. Based on the general saturation model, a saturation equation for vuggy reservoirs is derived, and the physical significance of the equation parameters is discussed based on the seepage-electricity similarity. The equation parameters depend on the pore structure, and vugs and matrix pore size distribution. Furthermore, a method for calculating the equation parameters is proposed, which uses nuclear magnetic resonance （NMR） data to calculate the capillary pressure curve. Field application of the proposed equation and parameter derivation method shows good match between calculated and experimental results, with an average absolute error of 5.8%.

Soil freezing process and different expressions for the soil-freezing characteristic curve

The soil-freezing characteristic curve(SFCC),which represents the relationship between unfrozen water content and subfreezing temperature(or suction at ice-water interface)in a freezing soil,can be used for understanding the transportation of heat,water,and solute in frozen soils.In this paper,the soil freezing process and the similarity between the SFCC of saturated frozen soil and soil-water characteristic curve(SWCC)of unfrozen unsaturated soil are reviewed.Based on similar characteristics between SWCC and SFCC,a conceptual SFCC is drawn for illustrating the main features of soil freezing and thawing processes.Various SFCC expressions from the literature are summarized.Four widely used expressions(i.e.,power relationship,exponential relationship,van Genuchten 1980 equation and Fredlund and Xing 1994 equation)are evaluated using published experimental data on four different soils(i.e.,sandy loam,silt,clay,and saline silt).Results show that the exponential relationship and van Genuchten(1980)equation are more suitable for sandy soils.The simple power relationship can be used to reasonably best-fit the SFCC for soils with different particle sizes;however,it exhibits limitations when fitting the saline silt data.The Fredlund and Xing(1994)equation is suitable for fitting the SFCCs for all soils studied in this paper.