CN1687773A - Method for accurate testing performance of pressure concentration-temperature of hydrogen storage material - Google Patents

Method for accurate testing performance of pressure concentration-temperature of hydrogen storage material Download PDF

Info

Publication number
CN1687773A
CN1687773A CN 200510046320 CN200510046320A CN1687773A CN 1687773 A CN1687773 A CN 1687773A CN 200510046320 CN200510046320 CN 200510046320 CN 200510046320 A CN200510046320 A CN 200510046320A CN 1687773 A CN1687773 A CN 1687773A
Authority
CN
China
Prior art keywords
hydrogen
pressure
temperature
storage material
air storage
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.)
Granted
Application number
CN 200510046320
Other languages
Chinese (zh)
Other versions
CN100371711C (en
Inventor
程宏辉
陈德敏
杨柯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Metal Research of CAS
Original Assignee
Institute of Metal Research of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Metal Research of CAS filed Critical Institute of Metal Research of CAS
Priority to CNB2005100463201A priority Critical patent/CN100371711C/en
Publication of CN1687773A publication Critical patent/CN1687773A/en
Application granted granted Critical
Publication of CN100371711C publication Critical patent/CN100371711C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Hydrogen, Water And Hydrids (AREA)

Abstract

This invention is a kind of method to measure the H2 storing material's PCT curved line. Especially it offers a kind of precise testing method which has wide use area to measure the H2 storing material's PCT curved line. This invention uses the volume method, that is to say, we will calculate the amount of H2 absorbed of the H2 storing material corresponding to different H2 pressure according to MBWR equation by continuously recording each part of pressure and temperature of H2 during the process of the H2 storing material absorbing or sending out H2, so as to get PCT curved line. After that we will realize the compression divisor's whole process of calculating and drawing picture by computer programming. As this invention has thought about the effect of the H2 storing room's temperature fluctuating, the PCT curved line will still be credible as though part of the room can not realize stable temperature. To sum up, this invention has good use value, making the instrument simple and easy to operate. At the same time we get higher degree of belief when we check the PCT curved line of carbon Nanometer pipe and V and so on which is the H2 storing material having high balance pressure.

Description

A kind of method for accurate testing performance of pressure concentration-temperature of hydrogen storage material
Technical field
The present invention provides the method for accurate testing of the PCT curve of the wide mensuration hydrogen storage material of a kind of scope of application especially for measuring a kind of method of hydrogen storage material Pressure-concentration-temperature (PCT) (PCT) curve.
Background technology
The method of mensuration PCT curve commonly used has volumetric method and mass method.
The capacity ratio juris is to extrapolate the variation of hydrogen storage material hydrogen content according to the pressure of the hydrogen storage material sample chamber of constant volume and air storage chamber.Filling hydrogen or putting the variation that hydrogen pressure changes and hydrogen concentration in the hydrogen storage material is calculated in accumulation in the METHOD FOR CONTINUOUS DETERMINATION constant volume obtains the PCT curve according to the state equation of hydrogen.When volumetric determination fills hydrogen process PCT curve, fill hydrogen according to the hydrogen amount of filling of estimation.Owing to fill the hydrogen amount by the Hydrogen Vapor Pressure adjustment sample of adjusting in the constant volume, test more convenient.And measure when putting the PCT curve of hydrogen, control hydrogen desorption capacity by the hydrogen pressure of adjusting in the constant volume of putting.The equipment investment of this method is less, generally all adopts this method to measure the PCT curve of hydrogen storage material.
Mass method is by the suction hydrogen of continuous recording hydrogen storage material or puts the hydrogen mass change and the variation of mutually deserved equalized pressure, thereby obtains the PCT curve.The advantage of mass method is to the direct measurement of being measured as of hydrogen or hydrogen desorption capacity, can be dynamically or measure the PCT curve statically, fill hydrogen and put the hydrogen process very convenient for mensuration.Yet this method equipment investment is huge, needs the little balance of high-resolution quartz oscillation.
Chinese patent application " a kind of pressure, density and temperature performance test method for hydrogen storing alloy, application number: 02109169.2 " reported method only is applicable to the measurement with low pressure or ultralow pressure hydrogen storage material PCT curve.And Chinese patent application " a kind of pressure, density and temperature performance test method for hydrogen storing alloy; application number: 02109169.2 " and " rare metal ", 1997,21 (1), the simple measuring of hydrogen-storage alloy PCT curve, reported method all adopts the perfect condition equation to measure, and this kind method precision is lower.When higher in equilibrium pressure especially, temperature was low, the state equation of actual hydrogen and perfect condition equation had bigger deviation, if the PCT curve that at this moment still adopts the perfect condition equation to measure hydrogen storage material will certainly cause very big error.
Summary of the invention
The object of the present invention is to provide the wide method for accurate testing performance of pressure concentration-temperature of hydrogen storage material of a kind of scope of application.
Technical scheme of the present invention is:
A kind of method of accurate mensuration hydrogen storage material PCT performance comprises the steps:
(1) takes by weighing the hydrogen storage material of quality m, put into the sample chamber;
(2) measure air storage chamber temperature, sample chamber temperature by temperature transmitter, pressure unit is measured the experimental pressure of air storage chamber and sample chamber, and the sample chamber keeps constant temperature;
(3) at first the sample chamber being placed temperature is to open valve under 100~120 ℃, whole test system is vacuumized continue 60~120 minutes; Stop to vacuumize, valve-off supposes that the hydrogen storage content of sample this moment is 0, the gaseous tension 0 of sample chamber;
(4) in air storage chamber, charge into a certain amount of hydrogen and make its pressure remain on 0.001-20MPa, measure and write down the temperature T c of air storage chamber this moment 1, pressure P c 1Open valve, the hydrogen in the hydrogen storage chamber enters sample chamber while hydrogen storage material by connecting tube and begins to inhale hydrogen, and the force value that shows when pressure unit did not change in 10 minutes, think that promptly hydrogen storage material suction hydrogen is saturated under this pressure, record system balancing pressure P this moment b 1And this moment air storage chamber temperature T b 1
(5) valve-off, charged pressure is the hydrogen of 1kPa~20MPa in air storage chamber, measures and write down the temperature T c of air storage chamber this moment 2, pressure P c 2, open valve, the hydrogen in the air storage chamber enters sample chamber while hydrogen storage material and begins to inhale hydrogen, and the force value that shows when pressure unit did not change in 10 minutes, thought that promptly hydrogen storage material suction hydrogen is saturated under this pressure, write down the equalized pressure record Pb of system's this moment 2And this moment air storage chamber temperature T b 2
(6) repeat above-mentioned steps 3,4 processes, it is saturated thoroughly to inhale hydrogen up to hydrogen storage material;
(7) wait to inhale hydrogen saturated after, after valve-off, air storage chamber are emitted the hydrogen of 1e-5~0.01mol, measure the temperature T c of air storage chamber this moment i, pressure P c i, open valve, the hydrogen of sample chamber enters air storage chamber hydrogen storage material releasing hydrogen gas simultaneously, and the force value that shows when pressure unit does not change in 10 minutes that to think promptly that under this pressure hydrogen storage material is put hydrogen abundant, record system balancing pressure P this moment b iAnd this moment air storage chamber temperature T b i
(8) repeat to put the hydrogen process repeatedly described in the above-mentioned steps 6, fully discharge, write down the temperature of last system balancing pressure and air storage chamber up to the hydrogen of hydrogen storage material;
(9) with the hydrogen under MBWR Equation for Calculating relevant temperature and pressure state compression factor Z (P, T), then according to formula n = PV RTZ Calculate the molar weight of balance front and back gaseous hydrogen, calculate its difference ordered series of numbers { Δ n}, this ordered series of numbers being obtained add up and ordered series of numbers cumsum{ Δ n}, is the y axle with equalized pressure Pb then, be that the x axle is drawn suction and put hydrogen PCT curve with hydrogen storage material hydrogen concentration cumsum{ Δ n} * 2/ (m/M); Wherein, R is a gas law constant, and T is an absolute temperature, and P is the pressure of hydrogen, and M is the mole fraction of hydrogen storage material;
(10) whole calculating and drawing course are finished by following program, and whole procedure comprises 5 functions: mbwr-comfactor function, mbwr-eos function, f function, hydron function, hydapi function; Wherein the mbwr-comfactor function is used to calculate the compressibility factor of hydrogen, and input variable is the pressure and temperature of hydrogen, and output variable is the compressibility factor of hydrogen; The mbwr-eos function is used to calculate the density of hydrogen, and input variable is the pressure and temperature of hydrogen, and output variable is the density of hydrogen, and used numerical computation method is a secant method; The f function is that the mbwr equation by above-mentioned hydrogen converts, and its input variable is density, the pressure and temperature of hydrogen, and output variable is a pressure, carries out interative computation for the mbwr-eos function call; The hydron function is used to calculate the molar weight of hydrogen, and its input variable is the temperature and pressure of volume of a container, gas, and output variable is the molar weight of hydrogen; The hydapi function is used to calculate the adding up and ordered series of numbers of mol ratio H/M increment of hydrogen atom and hydrogen storage material molecule, and hydrogen storage material hydrogen and hydrogen desorption capacity, and make the PCT curve map, its input variable is an air storage chamber pressure, sample chamber pressure, the air storage chamber temperature, the sample chamber temperature, the air storage chamber volume, the volume of connecting tube, the volume of sample chamber behind the removal sample volume, the molar weight of sample, its output variable is for inhaling adding up and ordered series of numbers of the mol ratio increment put hydrogen process hydrogen atom and hydrogen storage material molecule, hydrogen process equilibrium pressure ordered series of numbers is put in suction, maximum hydrogen of hydrogen storage material and maximum hydrogen desorption capacity.
Vacuumize in the described step 3, make vacuum tightness reach 10 -1-10 -2Pa.
Described step 8 repeats to put the hydrogen process 20-40 time described in the above-mentioned steps 6.
The invention has the beneficial effects as follows:
1, the present invention is owing to having considered the influence of air storage chamber temperature fluctuation to experimental result, so even the air storage chamber of equipment part can not realize constant temperature, PCT curve that its test obtains is still true and reliable in program.
2, the present invention has good practical value, and it makes equipment become simple and easy operating, has higher confidence level simultaneously when test such as CNT, vanadium etc. has the PCT curve of hydrogen storage material of high equilibrium pressure.
3, the present invention can obviously improve testing efficiency, and it is all finished whole data processing and drawing course automatically by computing machine, and the operator only needs related pressure and temperature data inserted and gets final product.
Description of drawings
Fig. 1 is the simple and easy synoptic diagram of the inventive method operative installations; Among the figure, 1 temperature transmitter; 2 pressure units; 3 air storage chambers; 4 valves; 5 connecting tubes; 6 sample chambers; 7 hydrogen storage materials; 8 temperature transmitters.
Fig. 2 is before valve 4 is opened among Fig. 1, the molar weight of hydrogen in the each several part.
Fig. 3 opens the back hydrogen storage material for valve 4 among Fig. 1 and inhales and put hydrogen and reached balance, the molar weight of hydrogen in the each several part.
The LaNi of Fig. 4 for using the inventive method to obtain 5Inhale 40 ℃ the first time and put hydrogen PC curve.
Fig. 5 a-d is the MATLAB program flow diagram.
Embodiment
The inventive method operative installations as shown in Figure 1, this device comprises temperature transmitter 1, pressure unit 2, air storage chamber 3, valve 4, connecting tube 5, sample chamber 6, hydrogen storage material 7, temperature transmitter 8.Temperature transmitter 1, pressure unit 2 are housed on the air storage chamber 3; Temperature transmitter 8 is housed on the sample chamber 6, hydrogen storage material 7 is housed in it; Air storage chamber 3 links to each other by connecting tube 5 with sample chamber 6, and connecting tube 5 is provided with valve 4.
The invention provides a kind of method of accurate mensuration hydrogen storage material PCT performance, comprise the steps:
1, takes by weighing the hydrogen storage material of certain mass m, put into the sample chamber;
2, measure air storage chamber temperature, sample chamber temperature by temperature transmitter 1, pressure unit 2 is measured the experimental pressure of air storage chamber 3 and sample chamber 6, and sample chamber 6 keeps constant temperature;
3, at first the sample chamber is placed under the high temperature (about 100 ℃) to open valve 4, by the direct connection vacuum pump whole test system is vacuumized and continue 60~120 minutes, make vacuum tightness reach 10 -1Pa.Stop to vacuumize, valve-off 4 supposes that the hydrogen storage content of sample this moment is 0, the gaseous tension 0 of sample chamber;
4, in air storage chamber 3, charge into the hydrogen (pressure is 0.01-20MPa) of certain pressure, measure and write down the temperature T c of air storage chamber 3 this moment 1, pressure P c 1Open valve 4, hydrogen in the hydrogen storage chamber 3 enters sample chamber 6 while hydrogen storage material 7 by connecting tube 5 and begins to inhale hydrogen, and the force value that shows when pressure unit did not have variation to think that promptly hydrogen storage material suction hydrogen is saturated under this pressure in 10 minutes, write down system balancing pressure P b this moment 1(this force value be the pressure of sample chamber also be the pressure of air storage chamber) and this moment air storage chamber temperature T b 1
5, valve-off 4, and (pressure is 1kPa~20MPa), measures the also temperature T c of record air storage chamber this moment to charge into the hydrogen of certain pressure in air storage chamber 2, pressure P c 2Open valve 4, hydrogen in the air storage chamber enters sample chamber while hydrogen storage material 7 and begins to inhale hydrogen, and the force value that shows when pressure unit did not have variation to think that promptly hydrogen storage material suction hydrogen is saturated under this pressure in 10 minutes, and the equalized pressure that writes down system's this moment writes down Pb 2(this force value be the pressure of sample chamber also be the pressure of air storage chamber) and this moment air storage chamber temperature T b 2
6, repeat above-mentioned steps 3,4 processes, it is saturated thoroughly to inhale hydrogen up to hydrogen storage material;
7, wait to inhale hydrogen saturated after, valve-off 4, air storage chamber are emitted a certain amount of hydrogen and (behind the 1e-5~0.01mol), are measured the temperature T c of air storage chamber this moment i, pressure P c i, open valve 4, the hydrogen of sample chamber enters air storage chamber hydrogen storage material releasing hydrogen gas simultaneously, and the force value that shows when pressure unit does not change in 10 minutes that to think promptly that under this pressure hydrogen storage material is put hydrogen abundant, record system balancing pressure P this moment b i(this force value be the pressure of sample chamber also be the pressure of air storage chamber) and this moment air storage chamber temperature T b i
8, the hydrogen process of putting described in the repeating step 6 repeatedly (present embodiment is 20 times) fully discharges up to the hydrogen of hydrogen storage material, writes down the temperature of last system balancing pressure and air storage chamber;
9, with the equation MBWR equation of the description hydrogen virtual condition the best of generally acknowledging at present, calculate corresponding to the hydrogen uniform temperature and pressure state under compress factor Z (P, T), then according to formula n = PV RTZ Calculate the molar weight of balance front and back gaseous hydrogen, calculate its difference ordered series of numbers { Δ n}, this ordered series of numbers obtained add up and ordered series of numbers cumsum{ Δ n}, be the y axle with equalized pressure Pb then, with hydrogen storage material hydrogen concentration cumsum{ Δ n} * 2/ (m/M) (M: be that the x axle is drawn to inhale and put hydrogen PCT curve the mole fraction of hydrogen storage material);
10, whole calculating and drawing course are finished by working out the MATLAB program voluntarily.
The present invention is owing to having considered the influence of air storage chamber temperature fluctuation to experimental result, so even the air storage chamber of equipment part can not realize constant temperature, PCT curve that its test obtains is still true and reliable in program.This feature has good practical value, and it makes equipment become simple and easy operating, has higher confidence level simultaneously when test such as CNT, vanadium etc. has the PCT curve of hydrogen storage material of high equilibrium pressure.
What the present invention adopted is volumetric method, i.e. suction hydrogen by the continuous recording hydrogen storage material or put the each several part Hydrogen Vapor Pressure of hydrogen process and temperature is extrapolated the hydrogen of corresponding different hydrogen pressure hydrogen storage material according to the MBWR equation obtains the PCT curve.Key problem in technology is the calculating that how realizes compressibility factor by computer programming.Very complicated equation of MBWR equation, its form is:
P=ρRT+ρ 2(N 0T+N 1T 1/2+N 2+N 3/T+N 4/T 2)
3(N 5T+N 6+N 7/T+N 8/T 2)
4(N 9T+N 10+N 11/T)+ρ 5(N 12)
6(N 13/T+N 14/T 2)+ρ 7(N 15/T)
8(N 16/T+N 17/T 2)+ρ 9(N 18/T 2)
3(N 19/T 2+N 20/T 3)exp(γρ 2)
5(N 21/T 2+N 22/T 4)exp(γρ 2)
7(N 23/T 2+N 24/T 3)exp(γρ 2)
9(N 25/T 2+N 26/T 4)exp(γρ 2)
11(N 27/T 2+N 28/T 3)exp(γρ 2)
13(N 29/T 2+N 30/T 3+N 31/T 4)exp(γρ 2)
N wherein 0~N 31, γ is for describing the parameter of hydrogen character, and ρ is a hydrogen density, and R is a gas law constant, and T is an absolute temperature, and P is the pressure of hydrogen, and what this equation was described is the relation of Hydrogen Vapor Pressure and its density and temperature.Can measure the pressure and the temperature of gas by pressure unit and temperature transmitter, realize that by programming computing machine calculates hydrogen density ρ automatically, the numerical computation method that is adopted is a secant method.Then, according to formula Z = P ρRT Calculate corresponding to this pressure, and the hydrogen compression factor Z under the temperature (P, T).Before can getting valve 4 and open by Fig. 2 in the air storage chamber hydrogen molar weight be n 1i, the hydrogen molar weight is n in the connecting tube 2i, hydrogen molar weight n in the sample chamber 3i, after can getting valve 4 and open by Fig. 3, when suction put the hydrogen process when reaching balance in the air storage chamber hydrogen molar weight be n 4i, the hydrogen molar weight is n in the connecting tube 5i, the hydrogen molar weight is n in the sample chamber 6i, then the sucking/placing hydrogen amount of hydrogen storage material is When
Figure A20051004632000083
The time, the suction hydrogen of hydrogen storage material, when The time, hydrogen storage material put hydrogen.Because in the actual measurement process, hydrogen storage material is inhaled and to be put hydrogen balance and charge and discharge time interval very short (less than one minute) of hydrogen process to air storage chamber next time, and ambient temperature fluctuates very little, so can draw tc I+1=tb i, the temperature the when meaning of this equation is this air storage chamber temperature that records with the balance that recorded last time equates (i be one more than or equal to 1 variable).In addition, because last hydrogen storage material inhales and puts hydrogen and reach balance, so can think that the gaseous tension of sample chamber remains unchanged before filling hydrogen, i.e. Py next time I+1=Pb i, the pressure the when meaning of this equation is this sample chamber pressure that records with the balance that recorded last time equates (i be one more than or equal to 1 variable).Through filling hydrogen and putting the hydrogen process repeatedly, relevant temperature and pressure data are handled to obtain hydrogen storage material sucking/placing hydrogen amount ordered series of numbers
Figure A20051004632000091
This ordered series of numbers is added up and processing obtains one and adds up and ordered series of numbers cumsum This be one about the adding up and ordered series of numbers of hydrogen molecule, and in the actual test in the PCT curve map its horizontal ordinate be the mol ratio H/M of hydrogen atom and hydrogen storage material molecule, so x i = cumsum { Δn i } × 2 n M . Y in addition i=Pb i, be horizontal ordinate with xi, yi is that the ordinate mapping just obtains the PCT curve map.
As shown in Figure 2, before valve 4 is opened among Fig. 1, the molar weight of hydrogen in the each several part; As shown in Figure 3, valve 4 is opened the back hydrogen storage material and is inhaled and to put hydrogen and reached balance, the molar weight of hydrogen in the each several part among Fig. 1.
MATLAB program flow diagram shown in Fig. 5 a-d, whole procedure comprise 5 functions: mbwr-comfactor function, mbwr-eos function, f function, hydron function, hydapi function.Wherein, shown in Fig. 5 c, the mbwr-comfactor function is used to calculate the compressibility factor of hydrogen, and input variable is the pressure and temperature of hydrogen, and output variable is the compressibility factor of hydrogen; Shown in Fig. 5 d, the mbwr-eos function is used to calculate the density of hydrogen, and input variable is the pressure and temperature of hydrogen, and output variable is the density of hydrogen, and used numerical computation method is a secant method; Shown in Fig. 5 d, the f function is that the mbwr equation by above-mentioned hydrogen converts, and its input variable is density, the pressure and temperature of hydrogen, and output variable is a pressure, and this function does not have physical significance, and purpose is to carry out interative computation for the mbwr-eos function call; Shown in Fig. 5 b, the hydron function is used to calculate the molar weight of hydrogen, and its input variable is a volume of a container, the temperature and pressure of gas, and output variable is the molar weight of hydrogen; Shown in Fig. 5 a, the hydapi function is used to calculate the adding up and ordered series of numbers of mol ratio H/M increment of hydrogen atom and hydrogen storage material molecule, and hydrogen storage material hydrogen and hydrogen desorption capacity, and makes the PCT curve map.Its input variable is the volume of air storage chamber pressure, sample chamber pressure, air storage chamber temperature, sample chamber temperature, air storage chamber volume, connecting tube, the volume of removing sample chamber behind the sample volume, the molar weight of sample, and its output variable is for inhaling adding up and ordered series of numbers, inhaling and put hydrogen process equilibrium pressure ordered series of numbers, hydrogen storage material maximum hydrogen and maximum hydrogen desorption capacity of the mol ratio increment of putting hydrogen process hydrogen atom and hydrogen storage material molecule.
It below is the MATLAB program of whole computation process.
This subprogram is used to calculate the compressibility factor of hydrogen:
function?z=mbwr_comfactor(P,T)
%The?32-term?equation?of?state;
%Based?on?the?Modified?Bennedict-webb-rubin?equation?of?state;
%P--kPa;
%T--k;
The compressibility factor that % might calculate under extreme environment is inaccurate;
% is only at conventional temperature 273.15~573.15k, and under conventional pressure 0.001~100Mpa, the compression that this program calculates is reliably because of %, and the PCT curve test environment of hydrogen storage material is contained in this scope.
The numerical computation method that % uses is secant method.
R=8.314472;
ro=mbwr_eos(P,T);
z=P./ro./R./T;
function?ro=mbwr_eos(P,T)
R=8.31447;
ro0=0.002;
ro1=P./R./T;
y=1;
while?ab(y)>1e-7
ro2=(ro0.*feval(@f,ro1,P,T)-ro1.*feval(@f,ro0,P,T))./...(feval(@f,ro1,P,T)-feval(@f,ro0,P,T));
ro0=ro1;
ro1=ro2;
y=feval(@f,ro1,P,T);
end
ro=ro1;
function?y=f(ro,P,T)
rr=0.00831434;
gm=-0.41000000E-02;
g0=.4675528393416E-04;g1=.4289274251454E-02;g2=.5164085596504E-01;
g3=.2961790279801E+00;g4=-.3027194968412E+01;g5=.1908100320379E-05;
g6=-.1339776859288E-03;g7=.3056473115421E-01;g8=.5161197159532E+01;
g9=.1999981550224E-07;g10=.2896367059356E-04;g11=-.2257803939041E-02;
g12=-.2287392761826E-06;g13=.2446261478645E-05;g14=-.1718181601119E-03;
g15=-.5465142603459E-07;g16=.4051941401315E-09;g17=.1157595123961E-06;
g18=-.1269162728389E-08;g19=-.4983023605519E+01;g20=-.1606676092098E+02;
g21=-.1926799185310E-01;g22=.9319894638928E+00;g23=-.3222596554434E-04;
g24=.1206839307669E-03;g25=-.3841588197470E-07;g26=-.4036157453608E-05;
g27=-.1250868123513E-10;g28=.1976107321888E-09;g29=.2411883474011E-13;
g30=-.4127551498251E-13;g31=.8917972883610E-12;
d0=ro;d1=ro.^2;d2=ro.^3;d3=ro.^4;d4=ro.^5;d5=ro.^6;d6=ro.^7;d7=ro.^8;d8=ro.^9;
d9=ro.^10;d10=ro.^11;d11=ro.^12;d12=ro.^13;
t0=sqrt(T);t1=T.^2;t2=T.^3;t3=T.^4;
f=exp(gm.*d1);
b0=d1.*T;b1=d1.*t0;b2=d1;
b3=d1./T;b4=d1./t1;b5=d2.*T;
b6=d2;b7=d2./T;b8=d2./t1;
b9=d3.*T;b10=d3;b11=d3./T;
b12=d4;b13=d5./T;b14=d5./t1;
b15=d6./T;b16=d7./T;b17=d7./t1;
b18=d8./t1;b19=d2.*f./t1;b20=d2.*f./t2;
b21=d4.*f./t1;b22=d4.*f./t3;b23=d6.*f./t1;
b24=d6.*f./t2;b25=d8.*f./t1;b26=d8.*f./t3;
b27=d10.*f./t1;b28=d10.*f./t2;b29=d12.*f./t1;
b30=d12.*f./t2;b31=d12.*f./t3;
p=(b0.*g0+b1.*g1+b2.*g2+b3.*g3+b4.*g4+b5.*g5+b6.*g6+b7.*g7+b8.*g8+b9.*g9...+b10.*g10+b11.*g11+b12.*g12+b13.*g13+b14.*g14+b15.*g15+b16.*g16+b17.*g17...+b18.*g18+b19.*g19+b20.*g20+b21.*g21+b22.*g22+b23.*g23+b24.*g24+b25.*g25...+b26.*g26+b27.*g27+b28.*g28+b29.*g29+b30.*g30+b31.*g31+rr.*d0.*T).*1000;
y=p-P;
This subprogram is used to calculate the molar weight of hydrogen:
function?n=hydron(P,t,V)
% is used to ask the true molar weight of hydrogen
%P--kPa;
%t--℃;
%V-ml;
%n--mol
When % was 0kPa when Hydrogen Vapor Pressure, its molar weight was 0
R=8.314472;
T=t+273.15;
J=find(P==0);
n=P.*V./1000./R./T./mbwr_comfactor(P,T);
n(J)=0;
This subprogram is used to calculate the adding up and ordered series of numbers of mol ratio H/M increment of hydrogen atom and hydrogen storage material molecule, and hydrogen storage material hydrogen and hydrogen desorption capacity, and makes the PCT curve map.
function[cumlsum,cumlsuma,cumlsumd,Pb,Pba,Pbd,maxa,maxd]=hydapi(Pc,Py,tc,ty,Vc,Vy1,Vy2,nm)%[cumlsum?cumlsuma?cumlsund?Pb?Pba?Pbd?maxa
%maxd]=hydapi(Pc,Py,tc,ty,Vc,Vy1,Vy2,nm);
The air storage chamber pressure that writes down in the %Pc-test process, unit are kPa;
The sample chamber pressure that writes down in the %Py-test process, unit are kPa;
The temperature of the air storage chamber that writes down in the %tc-test process, unit is ℃;
The temperature of the sample chamber of writing down in the %ty-test process, unit is ℃;
The volume of %Vc-air storage chamber, unit are ml;
The volume of %Vy1-connecting tube, unit are ml;
The volume of sample chamber behind the %Vy2-removal sample volume, unit is ml;
The molar weight of %nm-sample, unit is mol;
%cumlsum-comprises the adding up and ordered series of numbers of mol ratio H/M increment of the hydrogen atom inhaling hydrogen and put the hydrogen process and hydrogen storage material molecule;
%cumlsuma-only inhales the adding up and ordered series of numbers of mol ratio H/M increment of the hydrogen atom of hydrogen process and hydrogen storage material molecule;
%cumlsumd-is only put the adding up and ordered series of numbers of mol ratio H/M increment of the hydrogen atom of hydrogen process and hydrogen storage material molecule;
The %Pb-hydrogen storage material is inhaled and is put the equilibrium pressure ordered series of numbers of hydrogen process;
The %Pba-hydrogen storage material is inhaled the equilibrium pressure ordered series of numbers of hydrogen process;
The %Pbd-hydrogen storage material is put the equilibrium pressure ordered series of numbers of hydrogen process;
The maximum hydrogen of %maxa-hydrogen storage material;
The maximum hydrogen desorption capacity of %maxb-hydrogen storage material;
J=length (Py); % obtains the length of Py array;
% presses array according to equation Py (i+1)=Pb (i) equilibrium establishment;
Pb=zeros(J-1,1);
for?i=1:J-1
Pb(i)=Py(i+1);
end
% sets up the equilibrium pressure array of inhaling the hydrogen process;
[CI]=max(Pb);
Pba=zeros(I,1);
fori=1:I
Pba(i)=Pb(i);
end
% sets up the equilibrium pressure array of putting the hydrogen process;
Pbd=zeros(J-I,1);
for?i=1:J-I
Pbd(i)=Pb(J-i);
end
% sets up the array of the sample chamber pressure that is used for the use of actual computation process;
Py1=zeros(J-1,1);
for?i=1:J-1
Py1(i)=Py(i);
end
% sets up the array of the preceding air storage chamber temperature of molecular balance that is used for the actual computation process;
tc1=zeros(J-1,1);
for?i=1:J-1
tc1(i)=tc(i);
end
% sets up the array of air storage chamber temperature behind the molecular balance that is used for the actual computation process according to equation tc (i+1)=tb (i);
tb=zeros(J-1,1);
for?i=1:J-1
tb(i)=tc(i+1);
end
% calculates the sucking/placing hydrogen amount array of inhaling the hydrogen storage material of putting the hydrogen process;
N1=hydron (Pc, tc1, Vc); The molar weight of Gaseous Hydrogen molecule in the air storage chamber before the % reaction;
N2=hydron (Py1, tc1, Vy1); The molar weight of Gaseous Hydrogen molecule in the connecting tube before the % reaction;
N3=hydron (Py1, ty, Vy2); % reacts the molar weight of Gaseous Hydrogen molecule in the preceding sample chamber;
N4=hydron (Pb, tb, Vc); The molar weight of Gaseous Hydrogen molecule in the air storage chamber after the % balance;
N5=hydron (Pb, tb, Vy1); The molar weight of Gaseous Hydrogen molecule in the connecting tube after the % balance
N6=hydron (Pb, ty, Vy2); The molar weight of Gaseous Hydrogen molecule in the sample chamber after the % balance;
Detan=n1+n2+n3-n4-n5-n6; The % hydrogen storage material is inhaled and is put the hydrogen molar weight;
cumlsum=cumsum(detan).*2./nm;
cumlsuma=zeros(I,1);
for?i=1:I
cumlsuma(i)=cumlsum(i);
end
cumlsumd=zeros(J-I,1);
for?i=1:J-I
cumlsumd(i)=cumlsum(J-i);
end
maxa=cumlsum(I);
maxd=maxa-cumlsum(J-1);
plot(cumlsum,Pb);
Preparation hydrogen storage material LaNi 5Alloy button ingot is pressed the design mix preparation, and used metal purity (at%) is La 99.3%, and Ni 99.9%.Above-mentioned raw materials is cleaned with ultrasonator, press the design mix counterweight afterwards, the material for preparing is placed water cooled copper mould; under argon gas atmosphere protection, melting in electric arc furnaces, in fusion process with alloy turn-over remelting five times; carry out electromagnetic agitation simultaneously, to guarantee the homogeneity of alloy.After the melting alloy pig is enclosed in the vitreosil pipe, puts into heat-treatment furnace, be warming up to 1100 ℃ of insulations with stove and quench after 8 hours.On this integral test system, carry out LaNi then 5The hydrogen storage performance test of hydrogen storage material, used hydrogen purity 99.999%.The sample that 1~2 gram is ground is put into reactor, and sample vacuumizes (to 10 with mechanical pump in the time of 40 ℃ -2Pa) after 50~60 minutes, in device shown in Figure 1, test by the inventive method.
The LaNi of Fig. 4 for using the inventive method to obtain 5Inhale 40 ℃ the first time and put hydrogen pressure-composition (P-C) curve map.Maximum hydrogen is 6.2961, and maximum hydrogen desorption capacity is 5.841.As seen from the figure, the PCT curve map that utilizes the perfect condition equation to record obviously can not react truth.
Adopt the present invention to can be used as the hydrogen storage material performance testing equipment method of testing in (Chinese invention patent application, application number 200410050594.3, October 15 2004 applying date).

Claims (3)

1, a kind of method for accurate testing performance of pressure concentration-temperature of hydrogen storage material is characterized in that comprising the steps:
(1) takes by weighing the hydrogen storage material of quality m, put into the sample chamber;
(2) measure air storage chamber temperature, sample chamber temperature by temperature transmitter, pressure unit is measured the experimental pressure of air storage chamber and sample chamber, and the sample chamber keeps constant temperature;
(3) at first the sample chamber being placed temperature is to open valve under 100~120 ℃, whole test system is vacuumized continue 60~120 minutes; Stop to vacuumize, valve-off supposes that the hydrogen storage content of sample this moment is 0, the gaseous tension 0 of sample chamber;
(4) in air storage chamber, charge into a certain amount of hydrogen and make its pressure remain on 0.001-20MPa, measure and write down the temperature T c of air storage chamber this moment 1, pressure P c 1Open valve, the hydrogen in the hydrogen storage chamber enters sample chamber while hydrogen storage material by connecting tube and begins to inhale hydrogen, and the force value that shows when pressure unit did not change in 10 minutes, think that promptly hydrogen storage material suction hydrogen is saturated under this pressure, record system balancing pressure P this moment b 1And this moment air storage chamber temperature T b 1
(5) valve-off, charged pressure is the hydrogen of 1kPa~20MPa in air storage chamber, measures and write down the temperature T c of air storage chamber this moment 2, pressure P c 2, open valve, the hydrogen in the air storage chamber enters sample chamber while hydrogen storage material and begins to inhale hydrogen, and the force value that shows when pressure unit did not change in 10 minutes, thought that promptly hydrogen storage material suction hydrogen is saturated under this pressure, write down the equalized pressure record Pb of system's this moment 2And this moment air storage chamber temperature T b 2
(6) repeat above-mentioned steps 3,4 processes, it is saturated thoroughly to inhale hydrogen up to hydrogen storage material;
(7) wait to inhale hydrogen saturated after, after valve-off, air storage chamber are emitted the hydrogen of 1e-5~0.01mol, measure the temperature T c of air storage chamber this moment i, pressure P c i, open valve, the hydrogen of sample chamber enters air storage chamber hydrogen storage material releasing hydrogen gas simultaneously, and the force value that shows when pressure unit does not change in 10 minutes that to think promptly that under this pressure hydrogen storage material is put hydrogen abundant, record system balancing pressure P this moment b iAnd this moment air storage chamber temperature T b i
(8) repeat to put the hydrogen process repeatedly described in the above-mentioned steps 6, fully discharge, write down the temperature of last system balancing pressure and air storage chamber up to the hydrogen of hydrogen storage material;
(9) with the hydrogen under MBWR Equation for Calculating relevant temperature and pressure state compression factor Z (P, T), then according to formula n = PV RTZ Calculate the molar weight of balance front and back gaseous hydrogen, calculate its difference ordered series of numbers { Δ n}, this ordered series of numbers being obtained add up and ordered series of numbers cumsum{ Δ n}, is the y axle with equalized pressure Pb then, be that the x axle is drawn suction and put hydrogen PCT curve with hydrogen storage material hydrogen concentration cumsum{ Δ n} * 2/ (m/M); Wherein, R is a gas law constant, and T is an absolute temperature, and P is the pressure of hydrogen, and M is the mole fraction of hydrogen storage material;
(10) whole calculating and drawing course are finished by following program, and whole procedure comprises 5 functions: mbwr-comfactor function, mbwr-eos function, f function, hydron function, hydapi function; Wherein the mbwr-comfactor function is used to calculate the compressibility factor of hydrogen, and input variable is the pressure and temperature of hydrogen, and output variable is the compressibility factor of hydrogen; The mbwr-eos function is used to calculate the density of hydrogen, and input variable is the pressure and temperature of hydrogen, and output variable is the density of hydrogen, and used numerical computation method is a secant method; The f function is that the mbwr equation by above-mentioned hydrogen converts, and its input variable is density, the pressure and temperature of hydrogen, and output variable is a pressure, carries out interative computation for the mbwr-eos function call; The hydron function is used to calculate the molar weight of hydrogen, and its input variable is the temperature and pressure of volume of a container, gas, and output variable is the molar weight of hydrogen; The hydapi function is used to calculate the adding up and ordered series of numbers of mol ratio H/M increment of hydrogen atom and hydrogen storage material molecule, and hydrogen storage material hydrogen and hydrogen desorption capacity, and make the PCT curve map, its input variable is an air storage chamber pressure, sample chamber pressure, the air storage chamber temperature, the sample chamber temperature, the air storage chamber volume, the volume of connecting tube, the volume of sample chamber behind the removal sample volume, the molar weight of sample, its output variable is for inhaling adding up and ordered series of numbers of the mol ratio increment put hydrogen process hydrogen atom and hydrogen storage material molecule, hydrogen process equilibrium pressure ordered series of numbers is put in suction, maximum hydrogen of hydrogen storage material and maximum hydrogen desorption capacity.
2, according to the described method for accurate testing performance of pressure concentration-temperature of hydrogen storage material of claim 1, it is characterized in that: vacuumize in the described step 3, make vacuum tightness reach 10 -1-10 -2Pa.
3, according to the described method for accurate testing performance of pressure concentration-temperature of hydrogen storage material of claim 1, it is characterized in that: described step 8 repeats to put the hydrogen process 20-40 time described in the above-mentioned steps 6.
CNB2005100463201A 2005-04-27 2005-04-27 Method for accurate testing performance of pressure concentration-temperature of hydrogen storage material Expired - Fee Related CN100371711C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2005100463201A CN100371711C (en) 2005-04-27 2005-04-27 Method for accurate testing performance of pressure concentration-temperature of hydrogen storage material

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2005100463201A CN100371711C (en) 2005-04-27 2005-04-27 Method for accurate testing performance of pressure concentration-temperature of hydrogen storage material

Publications (2)

Publication Number Publication Date
CN1687773A true CN1687773A (en) 2005-10-26
CN100371711C CN100371711C (en) 2008-02-27

Family

ID=35305815

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2005100463201A Expired - Fee Related CN100371711C (en) 2005-04-27 2005-04-27 Method for accurate testing performance of pressure concentration-temperature of hydrogen storage material

Country Status (1)

Country Link
CN (1) CN100371711C (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477018B (en) * 2008-12-30 2011-04-20 中国科学院上海微系统与信息技术研究所 Full-automatic hydrogen storage material performance tester and its testing method
CN101876618B (en) * 2009-12-11 2012-08-22 北京有色金属研究总院 Measuring and correcting method for error of hydrogen storage capacity isasteric method system
CN103234860A (en) * 2013-04-28 2013-08-07 扬州大学 Carbon dioxide adsorption isotherm precision testing method
CN103439474A (en) * 2013-08-14 2013-12-11 扬州大学 Method for determining hydrogen absorption/desorption thermodynamic parameter of hydrogen storage alloy
CN105928855A (en) * 2016-04-20 2016-09-07 中国原子能科学研究院 Measuring apparatus and method for surface micropressure tritium absorption rate constant of material
CN106838610A (en) * 2017-03-02 2017-06-13 扬州大学 The activation of multichannel alloy-type hydrogen-storage tank, performance test and encapsulation integration apparatus and its technological process
CN107561134A (en) * 2017-09-20 2018-01-09 扬州大学 A kind of thin-film material inhales the method for testing for putting hydrogen PCT curves
CN107782487A (en) * 2017-09-23 2018-03-09 南京律智诚专利技术开发有限公司 A kind of pressure transmitter suitable for microwave hydrothermal reactor
CN110297066A (en) * 2019-07-15 2019-10-01 中国船舶重工集团公司第七一八研究所 A kind of VOCs concentration on-line measurement device
JP2020197248A (en) * 2019-06-04 2020-12-10 清水建設株式会社 Hydrogen storage state estimation device, hydrogen storage state estimation program and hydrogen storage state estimation method
CN113739067A (en) * 2021-09-02 2021-12-03 西北工业大学太仓长三角研究院 Dual-cavity-based hydrogen charging control device, hydrogen charging control method and equipment
CN114558028A (en) * 2022-03-18 2022-05-31 北京工业大学 Hydrogen molecule therapeutic agent for recovering brain function after TBI injury

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2940730B2 (en) * 1991-09-30 1999-08-25 三洋電機株式会社 Evaluation method of equilibrium characteristics of hydrogen storage alloy
JPH11101729A (en) * 1997-09-29 1999-04-13 Matsushita Electric Ind Co Ltd Apparatus for measuring reaction speed of hydrogen-storing alloy
CN1230668C (en) * 2002-02-09 2005-12-07 中国科学院金属研究所 Pressure, density and temperature performance test method for hydrogen storing alloy
CN1243241C (en) * 2003-12-10 2006-02-22 复旦大学 Measuring method for hydrogen storage quantity of nanometer hydrogen storage material
CN1580752A (en) * 2004-05-17 2005-02-16 华南理工大学 Method and device for synchronized determining hydrogen storage alloy structure and tissue change during charging-discharing circle
CN100487459C (en) * 2004-10-15 2009-05-13 中国科学院金属研究所 Hydrogen storage alloy performance integrated experimental apparatus

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477018B (en) * 2008-12-30 2011-04-20 中国科学院上海微系统与信息技术研究所 Full-automatic hydrogen storage material performance tester and its testing method
CN101876618B (en) * 2009-12-11 2012-08-22 北京有色金属研究总院 Measuring and correcting method for error of hydrogen storage capacity isasteric method system
CN103234860A (en) * 2013-04-28 2013-08-07 扬州大学 Carbon dioxide adsorption isotherm precision testing method
CN103234860B (en) * 2013-04-28 2015-04-01 扬州大学 Carbon dioxide adsorption isotherm precision testing method
CN103439474A (en) * 2013-08-14 2013-12-11 扬州大学 Method for determining hydrogen absorption/desorption thermodynamic parameter of hydrogen storage alloy
CN105928855A (en) * 2016-04-20 2016-09-07 中国原子能科学研究院 Measuring apparatus and method for surface micropressure tritium absorption rate constant of material
CN106838610A (en) * 2017-03-02 2017-06-13 扬州大学 The activation of multichannel alloy-type hydrogen-storage tank, performance test and encapsulation integration apparatus and its technological process
CN106838610B (en) * 2017-03-02 2018-11-20 扬州大学 The activation of multichannel alloy-type hydrogen-storage tank, performance test and encapsulation integration apparatus and its process flow
CN107561134B (en) * 2017-09-20 2019-04-05 扬州大学 A kind of thin-film material inhales the test method of hydrogen release PCT curve
CN107561134A (en) * 2017-09-20 2018-01-09 扬州大学 A kind of thin-film material inhales the method for testing for putting hydrogen PCT curves
CN107782487A (en) * 2017-09-23 2018-03-09 南京律智诚专利技术开发有限公司 A kind of pressure transmitter suitable for microwave hydrothermal reactor
JP2020197248A (en) * 2019-06-04 2020-12-10 清水建設株式会社 Hydrogen storage state estimation device, hydrogen storage state estimation program and hydrogen storage state estimation method
CN110297066A (en) * 2019-07-15 2019-10-01 中国船舶重工集团公司第七一八研究所 A kind of VOCs concentration on-line measurement device
CN110297066B (en) * 2019-07-15 2024-04-19 中国船舶重工集团公司第七一八研究所 VOCs concentration on-line measuring device
CN113739067A (en) * 2021-09-02 2021-12-03 西北工业大学太仓长三角研究院 Dual-cavity-based hydrogen charging control device, hydrogen charging control method and equipment
CN114558028A (en) * 2022-03-18 2022-05-31 北京工业大学 Hydrogen molecule therapeutic agent for recovering brain function after TBI injury

Also Published As

Publication number Publication date
CN100371711C (en) 2008-02-27

Similar Documents

Publication Publication Date Title
CN1687773A (en) Method for accurate testing performance of pressure concentration-temperature of hydrogen storage material
CN1214849C (en) Carbon dioxide absorption agent
CN1673693A (en) Precision calibrating method of container volume
CN1237337C (en) Isotope gas spectrometer instrument
CN1790616A (en) Substrate processing apparatus, program for performing operation and control method thereof, and computer readable storage medium storing the program
CN1235665C (en) Composite ceramic membrane oxygen separation method
CN1053200C (en) New fluoroelastomers comprising monomeric units deriving from a bis-olefin
CN101056906A (en) Modified polytetrafluoroethylen fine powder and modified polytetrafluoroethylene molded product
CN1035200C (en) Heat treatment method and device for metal materials
CN1017556B (en) Method of determining quality of steam for stimulating hydrocarbon production
CN1782728A (en) Battery power detection device
CN1108215A (en) Method for producing alpha-alumina powder
CN101074940A (en) Device for analysis gas in insulating oil and method for analysis gas in insulating oil
CN1934679A (en) Method and apparatus for forming silicon dots
CN1518350A (en) Frame data correction output device, correction device and method, display device
CN1898282A (en) Non melt processable polytetrafluoroethylene and fine powder thereof
CN1724986A (en) Method of monitoring pressure of a gas species and apparatus to do so
CN1517314A (en) Method for producing fibre-optical precast device and fibre-optical precast device and fibre-optical produced by the method
CN1269968C (en) Analytic method and use for reducing real-time fluorescent PCR instrument quantitative analytic system error
CN1210848A (en) Method and device for removing formic acid and acetic acid from water of glycol production
CN1652908A (en) Weighing equipment for concrete material
CN1701220A (en) Load sensor and its manufacturing method
CN1054778A (en) Novel type fluorubber based on vinylidene fluoride, R 1216 and tetrafluoroethylene
CN1187902C (en) Decoding device
CN1154238C (en) Interleave address generating device and interleave address generating method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080227

Termination date: 20110427