CN107557085B - A kind of the kinematic viscosity test and prediction technique of biomass-based liquid fuel - Google Patents
A kind of the kinematic viscosity test and prediction technique of biomass-based liquid fuel Download PDFInfo
- Publication number
- CN107557085B CN107557085B CN201710789659.3A CN201710789659A CN107557085B CN 107557085 B CN107557085 B CN 107557085B CN 201710789659 A CN201710789659 A CN 201710789659A CN 107557085 B CN107557085 B CN 107557085B
- Authority
- CN
- China
- Prior art keywords
- test
- kinematic viscosity
- temperature
- fuel
- viscosity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Abstract
The invention discloses a kind of test of the kinematic viscosity of fuel combination and prediction techniques, this method calculates the test and prediction technique combined with numerical value using test, dynamic viscosity-temperature curve is tested based on rotational rheometer, and multiple numerical fitting is carried out using 3 points of calibration, lignocellulosic base levulic acid esters-diesel oil kinematic viscosity database critical issue is established in test and prediction lignocellulosic base levulic acid esters-diesel blended fuel kinematic viscosity, solution.The method of the present invention is greatly saved the testing time relative to conventional test methodologies, can the convenient and efficient kinematic viscosity predicted in wide temperature range.When tested viscosity is worth biggish liquid fuel, " air pocket " phenomenon bring measurement error can be evaded, effectively conducive to kinematic viscosity database is established.
Description
Technical field
The present invention relates to a kind of test of the kinematic viscosity of fuel combination and prediction technique more particularly to a kind of lignocellulosics
Base levulic acid esters-diesel blended fuel kinematic viscosity test and prediction technique.
Background technique
Petroleum-based energy rapidly depletes so that countries in the world accelerate the paces of demand clean energy resource.Biomass energy is six
It can uniquely collect, store in kind of renewable energy, transporting and the renewable energy of fixed carbon.It is produced based on lignocellulosic
Material liq fuel is an important research direction in Biomass Energy Utilization technology, and exploitation biomass liquid fuel technology is not
Chinese energy safety strategic importance can only be enhanced, and had to alleviating the withered consumption pressure of China's energy, reducing environmental emission pollution
Significance.
Lignocellulosic base liquid fuel make an addition to diesel fuel applications in vehicular engine be in Biomass Energy Utilization technology
The technology of one relative maturity, the key of the technological development are to carry out evaluation analysis to fuel combination.Conventional evaluation side at present
Method is to carry out direct test analysis according to engine testing rack, and this method is simply direct, but evaluation cycle is longer, is expended big
The time of amount and financial resources.It is a kind of feasible succinct approach that fuel is screened by the physicochemical property progress pre assessment of fuel.And
Kinematic viscosity is a particularly important parameter in physicochemical property, can influence the evaporation atomization and combustion stability of fuel.It establishes
Diesel blended fuel kinematic viscosity database is an important link for improving fuel physicochemical property database.
The conventional test methodologies of kinematic viscosity use GB 265-88, and this method test is relatively accurate.But to different temperatures
Under the conditions of kinematic viscosity test consuming time is long, the measuring accuracy of kinematic viscosity the accuracy of manufacture of capillary is relied on it is larger,
To cause the less reproducible of test;The biggish liquid of viscosity is also easy to produce " air pocket " phenomenon, causes test error larger.According to
The detailed database of kinematic viscosity that conventional method establishes fuel combination is almost the project for being difficult to complete.
Summary of the invention
The disadvantages of present invention is the consuming overlong time for solving conventional test methodologies, heavy fuel test inaccuracy, adopts
The test and prediction technique combined is calculated with numerical value with test, dynamic viscosity-temperature curve is tested based on rotational rheometer, and
Multiple numerical fitting, test and prediction lignocellulosic base levulic acid esters-diesel blended fuel are carried out using 3 points of calibration
Lignocellulosic base levulic acid esters-diesel oil kinematic viscosity database critical issue is established in kinematic viscosity, solution.
Therefore, present invention provide the technical scheme that a kind of kinematic viscosity test of fuel combination and prediction technique, the party
Method the following steps are included:
(1) kinematic viscosity for using rheometer test fuel combination chooses suitable flat plate mold diameter D, gap delta, angle speed
Spend r and range of temperature T1~T2DEG C, the test data of dynamic viscosity μ Yu temperature T are measured, half predicted motion viscosity number is obtained
According to η=ρ μ, density p is test scope (T in formula1~T2DEG C) the corresponding density of neutral temperature;
(2) three test calibration points are chosen in test scope, temperature test range are equally divided into three sections, in three temperature
It spends section and distinguishes one temperature spot T of any selection1、T2And T3, as calibration kinematic viscosity, calibration point uses the kinematic viscosity of test
GB 265-88 tests kinematic viscosity η1、η2And η3, obtain calibration point test value (T1,η1)、(T2,η2) and (T3,η3);
(3) by the test data of flow measurement, Numerical regression analysis is carried out to the kinematic viscosity of calculating, obtains Empirical Mode
Type f1, f1Regression analysis type is y=me-nx, m and n is constant in formula, it is desirable that the degree of fitting R of empirical model2>0.99;
(4) second is carried out to calibration point according to the regression analysis type for calculating the acquisition of kinematic viscosity data to be fitted, obtain
Demarcate fitting empirical formula f2, test temperature value range T ∈ (T1~T2℃);
(5) operation: f=(af is carried out according to acquisition prediction model1+bf2)/(a+b), a and b are respectively weight factor, to f1
And f2Weight distribution is carried out, a+b=1, b > a in this formula determine weight factor using factor analysis;
(6) finally determine that prediction model is f, T ∈ (T1~T2℃);
(7) according to prediction model f, T ∈ (T1~T2DEG C) successively calculate kinematic viscosity value under the conditions of each temperature.
The method of the present invention is to carry out numerical value calculating with kinematic viscosity of the numerical model to certain temperature section, and use test
Means logarithm model is modified, and can obtain the kinematic viscosity number for meeting certain precision again to reach and reduce a large number of experiments
According to the purpose in library.
Method of the present invention, the fuel combination are preferably biomass-based diesel fuel;It is more preferably wooden fibre
Tie up plain base levulic acid esters-diesel blended fuel.
Method of the present invention, in a specific embodiment, using rotational rheometer test mixing in step (1)
The kinematic viscosity of fuel, flat plate mold use diameter D for 60mm, gap delta 1mm, and angular speed r is 50r/min, test scope T1
~T2DEG C be 0 DEG C~60 DEG C, gap;In step (6), a=0.3, b=0.7 in formula.
Compared with prior art, the present invention is had the following advantages and beneficial effects:
The present invention is greatly saved the testing time relative to conventional test methodologies, can convenient and efficient prediction wide temperature model
Enclose interior kinematic viscosity.When tested viscosity is worth biggish liquid fuel, it can effectively evade the measurement of " air pocket " phenomenon bring and miss
Difference, conducive to kinematic viscosity database is established.
Detailed description of the invention
Fig. 1 is E5 kinematic viscosity-temperature prediction model figure in the embodiment of the present invention.
Fig. 2 is B5 kinematic viscosity-temperature prediction model figure in the embodiment of the present invention.
Specific embodiment
Come that the present invention is furture elucidated below by the detailed description of specific embodiment, but is not to limit of the invention
System, only illustrates.
Embodiment one
The fuel that test is prepared is ethyl levulinate-diesel oil and Butyl acetylpropanoate-diesel blended fuel, and is added micro-
N-butanol is measured as cosolvent, respectively E5 and B5.E5 at being grouped as are as follows: 4.5~5%vol ethyl levulinate, 0~
0.5%vol n-butanol, 95%vol0# diesel oil;B5 at being grouped as are as follows: 4.5~5%vol Butyl acetylpropanoate, 0~0.5%
Vol n-butanol, 95%vol0# diesel oil.
1. using in rotational rheometer DHR-1 (TA, the U.S.), flat plate mold uses (diameter 60mm, gap 1mm), angular speed
For 50r/min, 0 DEG C~60 DEG C of test scope.Dynamic viscosity test data is measured, half predicted motion viscosity data η=ρ is obtained
μ, density p is the corresponding fuel combination density of 30 DEG C of (neutral temperature of test scope) temperature spots, ρ in formulaE5=0.818g
cm3, ρB5=0.822gcm3, μ is the dynamic viscosity of rheometer test.
2. successively choosing the calibration point of kinematic viscosity, test scope is divided into three sections (0~20,21~40,41~60),
Three temperature section difference one temperature spots of any selection are as calibration point T1、T2And T3, scale of thermometer fixed point is using GB 265-88 survey
Test run kinetic viscosity obtains η1、η2And η3.E5 choose calibration point be respectively (11,3.63E-3), (28,2.54E-3), (49,
1.63E-3);The calibration point that B5 chooses is respectively (8,2.81E-3), (51,2.88E-3), (51,1.14E-3).
3. carrying out Numerical regression analysis by the test data of flow measurement to the kinematic viscosity of calculating, obtaining Empirical Mode
Type f1, as depicted in figs. 1 and 2.The f of E51Empirical model are as follows: y=4.372e-0.02x,R2=0.9929;The f of B51Empirical model
Are as follows: y=4.635e-0.024x, R2=0.9999.
4. being fitted according to the regression analysis type that kinematic viscosity data obtain is calculated to calibration point, according to f1It can must return
Analysis type is y=me-nx, m and n is constant in formula, carries out quadratic fit to calibration point according to this recursive type, is demarcated
Fitting empirical formula f2, as depicted in figs. 1 and 2.The f of E52Empirical model are as follows: y=4.499e-0.02x,R2=0.9995;The f of B52
Empirical model are as follows: y=4.319e-0.02x, R2=0.9931.
5. carrying out operation: f=(af according to prediction model is obtained1+bf2)/(a+b), a and b are respectively weight factor, to f1
And f2Weight distribution is carried out, a=0.3 in formula, b=0.7 are acquired using factor analysis.Prediction model according to operational formula E5
Are as follows: y=0.3*4.372e-0.02x+0.7*4.499e-0.02x;The prediction model of B5 are as follows: y=0.3*4.635e-0.024x+0.7*
4.319e-0.02x。
6.f is finally determining prediction model, T ∈ (0 DEG C~60 DEG C).
E5: η=4.4609e-0.02T, T ∈ (0 DEG C~60 DEG C);The prediction model of B5 are as follows: η=1.391e-0.024T+3.023e-0.02x, T ∈ (0 DEG C~60 DEG C).
7. calculating E5 corresponding kinematic viscosity value at different temperatures according to prediction model:
B5 corresponding kinematic viscosity value at different temperatures is calculated according to prediction model:
Claims (5)
1. a kind of kinematic viscosity of fuel combination is tested and prediction technique, method includes the following steps:
(1) kinematic viscosity for using rheometer test fuel combination, chooses suitable flat plate mold diameter D, gap delta, angular speed r and
Range of temperature T1~T2DEG C, measure the test data of dynamic viscosity μ Yu temperature T, obtain half predicted motion viscosity data η=
ρ μ, density p is test scope (T in formula1~T2DEG C) the corresponding density of neutral temperature;
(2) three test calibration points are chosen in test scope, temperature test range are equally divided into three sections, in three temperature sections
One temperature spot T of any selection respectively1、T2And T3, for the kinematic viscosity of test as calibration kinematic viscosity, calibration point uses GB
265-88 tests kinematic viscosity η1、η2And η3, obtain calibration point test value (T1,η1)、(T2,η2) and (T3,η3);
(3) by the test data of flow measurement, Numerical regression analysis is carried out to the kinematic viscosity of calculating, obtains empirical model f1,
f1Regression analysis type is y=me-nx, m and n is constant in formula, it is desirable that the degree of fitting R of empirical model2>0.99;
(4) second is carried out to calibration point according to the regression analysis type for calculating the acquisition of kinematic viscosity data to be fitted, demarcated
Fitting empirical formula f2, test temperature value range T ∈ (T1~T2℃);
(5) operation: f=(af is carried out according to acquisition prediction model1+bf2)/(a+b), a and b are respectively weight factor, to f1And f2
Carry out weight distribution, a+b=1, b > a in this formula;
(6) finally determine that prediction model is f, T ∈ (T1~T2℃);
(7) according to prediction model f, T ∈ (T1~T2DEG C) successively calculate kinematic viscosity value under the conditions of each temperature.
2. according to the method for claim 1, it is characterised in that: the fuel combination is biomass-based diesel fuel.
3. according to the method for claim 2, it is characterised in that: flat plate mold diameter D is 60mm in step (1), and gap delta is
1mm, angular speed r are 50r/min, test scope T1~T2DEG C be 0 DEG C~60 DEG C.
4. according to the method for claim 2, it is characterised in that: in step (6), a=0.3, b=0.7 in formula.
5. according to the method for claim 2, it is characterised in that: the fuel combination is lignocellulosic base levulic acid
Esters-diesel blended fuel.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710789659.3A CN107557085B (en) | 2017-09-04 | 2017-09-04 | A kind of the kinematic viscosity test and prediction technique of biomass-based liquid fuel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710789659.3A CN107557085B (en) | 2017-09-04 | 2017-09-04 | A kind of the kinematic viscosity test and prediction technique of biomass-based liquid fuel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107557085A CN107557085A (en) | 2018-01-09 |
CN107557085B true CN107557085B (en) | 2019-03-15 |
Family
ID=60978051
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710789659.3A Expired - Fee Related CN107557085B (en) | 2017-09-04 | 2017-09-04 | A kind of the kinematic viscosity test and prediction technique of biomass-based liquid fuel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107557085B (en) |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105547918B (en) * | 2015-12-04 | 2018-01-30 | 广西玉柴机器股份有限公司 | The assay method of machine oil kinematic viscosity under high temperature |
CN106018176B (en) * | 2016-07-07 | 2019-01-22 | 中国石油集团川庆钻探工程有限公司 | A kind of evaluation method of water-base drilling fluid rheological characteristic |
-
2017
- 2017-09-04 CN CN201710789659.3A patent/CN107557085B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN107557085A (en) | 2018-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Alaoui et al. | Liquid density of 1-butanol at pressures up to 140 MPa and from 293.15 K to 403.15 K | |
Haputhanthri et al. | Ammonia and gasoline fuel blends for spark ignited internal combustion engines | |
CN103675217B (en) | A kind of gas detection method and device | |
CN102445406B (en) | Method for measuring liquid phase diffusion coefficient, and apparatus thereof | |
Burke et al. | Distillation-based droplet modeling of non-ideal oxygenated gasoline blends: investigating the role of droplet evaporation on PM emissions | |
CN202522295U (en) | System for calibrating turbine flow meter under condition of large flow of low-temperature medium | |
CN102607993A (en) | Modularized vacuum system of saturated vapor pressure and solid density of tested substance and applications | |
Meng et al. | An experimental study on the potential usage of acetone as an oxygenate additive in PFI SI engines | |
CN105547916A (en) | System for detecting water content in crude oil of oil field united station | |
CN107557085B (en) | A kind of the kinematic viscosity test and prediction technique of biomass-based liquid fuel | |
Kaletnik et al. | Study on performance of compression engine operated by biodiesel fuel. | |
CN104502227A (en) | Method for measuring fiber density | |
Patil et al. | Experimental investigation to assess the efficacy of gasoline surrogates with engine testing | |
Kumbár et al. | Temperature dependence viscosity and density of different biodiesel blends | |
CN109884666A (en) | A kind of tropospheric delay correction method based on data assimilation | |
Anish et al. | Impact of reactivity controlled compression ignition (RCCI) mode engine operation in diesel engine powered with B20 blend of waste cooking oil biodiesel | |
Zachiotis et al. | Methodology to estimate road grade effects on consumption and emissions from a light commercial vehicle running on the WLTC cycle | |
CN105549127B (en) | A kind of the atomized rain strength meter and measurement method of hydraulic model | |
Han et al. | Influence of coupling action of oxygenated fuel and gas circuit oxygen on soot generation in a diesel engine | |
CN103033592A (en) | Measuring method and device for measuring detonation resistance property of fuel | |
CN101769823A (en) | Method for optimizing experimental fuel proportioning based on thermogravimetric analysis and gas-liquid dual-phase calculation | |
Tesfa et al. | Effect of biodiesel blends on the transient performance of compression ignition engines | |
CN113933345A (en) | Construction method of cohesive soil thermal conductivity coefficient evaluation model | |
Song | Test study of the water quantity loss of cooling tower on U-type liquidometer | |
Vojtisek-Lom et al. | Measurement of consumption rates of viscous biofuels |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190315 Termination date: 20200904 |