CN113505483A - 一种考虑风机转动惯量的风速折算出力的精确方法 - Google Patents
一种考虑风机转动惯量的风速折算出力的精确方法 Download PDFInfo
- Publication number
- CN113505483A CN113505483A CN202110781808.8A CN202110781808A CN113505483A CN 113505483 A CN113505483 A CN 113505483A CN 202110781808 A CN202110781808 A CN 202110781808A CN 113505483 A CN113505483 A CN 113505483A
- Authority
- CN
- China
- Prior art keywords
- output
- wind speed
- fan
- wind
- speed
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/30—Wind power
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mechanical Engineering (AREA)
- Databases & Information Systems (AREA)
- Combustion & Propulsion (AREA)
- Chemical & Material Sciences (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Sustainable Energy (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Sustainable Development (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Wind Motors (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明涉及一种考虑风机转动惯量的风速折算出力的精确方法,包括以下步骤:获取水平轴风机的历史数据;风机的历史数据的预处理;采用功率曲线折算出未来的风机出力,得到从有效风速veq到出力p的具体折算方法;水平轴风机设计中,出力p和风轮转速n之间为分段的确定型函数关系,考虑机械惯性的更精确的出力折算,根据刚体定轴转动定律,忽略阻尼,建立风轮机械运动的时域微分方程;梯形法来求解风轮机械运动的时域微分方程,得到出力的折算值pj+1;分别采用具体模型求解梯形法的公式,得到pj+1的单项折算结果;采用偏最小二乘回归,优化组合成一个最终采用的出力值pj+1。
Description
技术领域
本发明属于风电技术领域,具体涉及一种计及风机机械惯性的“风速-出力”折算新方法。
背景技术
包括风能在内的新能源技术发展,不仅是碳中和的基本技术,而且对于我国的能源安全具有基础性的积极意义。自然界里风能的随机性、波动性和不可控性,使得高精度、高可靠性风电功率预测成为现代电网安全经济运行必须的关键基础信息。在我国,自然界里风的预报误差,是风电功率预测的第一误差来源;功率曲线(从风速到风力发电机组输出功率的折算)误差往往是第二误差来源。功率曲线误差可以引起百分之几的风电功率预测误差。
水平轴风力发电机的功率曲线,是风速和风力发电机组输出功率(以下简称“出力”)之间的函数关系,往往用图形或表格表示。一般由一周以上连续时间的风机实际运行的高频采样数据,考虑空气密度、风切变等多种因素的影响后,处理成10min时间分辨率的“风速-出力”函数关系。
目前的功率曲线,一般是将由多种因素引起的风速、出力的散点图,拟合成一条曲线,即“以线拟合带”。或者进一步加入可能的功率曲线变化的合理范围。散点图的分散性,大致可以分为人工和自然两类因素引起。如偏航角、桨距角的控制等人工因素,很难预先精确计及。由风机叶片的机械惯性等自然因素引起的分散性,反而可以表现出某种规律性。
功率曲线里的风速、出力,是一段时间内现场运行的一个可靠平均状态;但是,这个平均的风速、出力,并不利于及时反映实际运行中各种影响因素快速变化引起的偏差。所以,只采用一个风速数据点进行出力折算,无法有效考虑风速变化引起的出力变化。只有采用连续两点以上的风速数据点,才能反映风速变化率引起的出力偏移。
实际运行中风速快速变化,会由于风机叶片等的机械惯性,引起出力的时间延迟等。例如,我国杨茂等学者,发现了机械惯性引起的规律性变化:风速快速上升时,风机的实际出力往往比标准功率曲线折算出的出力小;反之亦然。已有考虑机械惯性的几项研究都提高了出力折算的准确性。总的来说,现有的考虑机械惯性的出力折算方法,通常都将历史数据进行离散性的近似分组处理,并采用各种先进的智能技术以得到更好的出力折算值。由于出力的历史数据受到多种因素的影响,因此仅依据历史数据,不利于更精细更可靠地依据叶片(叶轮)等的机械惯性来提高出力折算的准确性。
发明内容
本发明的目的是提供一种准确率高、能够反映风速快速变化作用的从风速到出力的折算新方法。在转速和出力为正相关的情况下,依据刚体定轴转动定律建立叶片(叶轮)的转动方程,并采用高精度、高稳定性的梯形法进行求解,实现风速快速变化时的出力折算。技术方案如下:
一种考虑风机转动惯量的风速折算出力的精确方法,包括以下步骤:
(2)风机的历史数据的预处理:对上述各种历史观测值进行有效性检验,数据修复预处理;
(4)水平轴风机设计中,出力p和风轮转速n之间为分段的确定型函数关系,在风速到达额定风速之前,当风速在一定范围变化时,风轮转速增大则引起出力增大;考虑机械惯性的更精确的出力折算,根据刚体定轴转动定律,忽略阻尼,建立风轮机械运动的时域微分方程:其中,J是叶片等的等效机械转动惯量,为风轮转速n对时间t的导数,为瞬时有效风速veq引起的动力转矩,TP为发电机出力折算的阻力转矩;
(6)对未来时间点j=1,2,…,将有效风速的veq,j+1、veq,j,以及由功率曲线折算的出力pPC,j,代入该梯形法公式,得到出力的折算值pj+1;
(7)分别采用三次多项式拟合、最小二乘支持向量机LSSVR、随机森林RF,以及深度学习里的长短期记忆网络LSTM、深度信念网络DBN共计5种具体模型求解梯形法的公式,得到pj+1的5种单项折算结果,即pCUBIC,j+1、pLSSVR,j+1、pRF,j+1、pLSTM,j+1、pDBN,j+1;
(8)采用偏最小二乘回归,5种单项折算的出力结果pCUBIC,j+1、pLSSVR,j+1、pRF,j+1、pLSTM,j+1、pDBN,j+1,优化组合成一个最终采用的出力值pj+1。
步骤(2)的预处理可以包括对历史风速、风向时间序列进行3点平滑去噪、双正交小波去噪。
附图说明
图1是本发明依据刚体定轴转动定律将未来风速风向折算为风机出力的流程图;
图2是用于验证本发明效果的历史风速、风轮转速和风机出力;
图3为实例采用的水平轴风机的历史风速和出力散点图;
图4是功率曲线散点图及其拟合函数(10m/s以下);
图5为实例风机的风轮转速与出力之间的确定性函数关系;
图6是功率曲线散点图及其拟合函数,(10r/min、0.15MW)、(18r/min、0.85MW)之间;
图7为实例中梯形法求解的三次多项式的回归次数a、b的具体变化;
图9是梯形法中三次多项式拟合折算的出力,和对应的出力历史值;
图10是梯形法中4种智能模型折算的出力,和对应的出力历史值;
图11是本发明的最终出力折算值。
具体实施方式
参见图1,本发明的采用精度高且稳定性好的梯形法来求解水平轴风机风速和风轮转速(发电机转速)之间刚体定轴转动定律方程,以实现风速精确折算为出力的方法,包括以下步骤:
本实例采用我国季风区某风场的2018年1月前后冬季风时期连续21天的实测数据,采样间隔为10秒钟。实测数据包括风速、风向、出力、转速、偏航角、桨距角和气温等。
以及未来3天的风速、风向等预报值。
用于验证本发明效果的历史出力、风轮转速、风速参见图2。这是一段较为快速变化的风速。
(2)风机的历史数据的预处理。
对上述各种历史观测值进行有效性检验,数据修复等常见的预处理。如对历史风速、风向等时间序列进行3点平滑去噪、双正交小波去噪。
(3)采用功率曲线折算出未来的风机出力。参照我国标准GB/T 33225-2016,并参考IEC 61400-12-1(2017-03),得到从有效风速veq到出力p的具体折算方法:这里分别为出力p和有效风速veq的10min平均值。在不同的统计周期(历史风速、出力等时间序列的时间分辨率)下,功率曲线p=f(veq)与10min下的变化不明显。将未来的风速代入即可得到由功率曲线折算出的未来出力。
本专利重点研究到达额定风速之前的出力折算,具体为0.15MW<p<0.85MW时的折算。
在本实例中,出力和有效风速的散点图请看图3。统计时间分别是10s和20s,以及10min(标准功率曲线的统计周期)。
本实例中,参照我国标准以及国际标准IEC 61400-12-1建议得到的功率曲线请看图4。风速在2.5m/s到10m/s之间的拟合函数为
图5为叶片(风轮)的转速与出力之间往往为分段函数关系。
图6是本实例选用的(10r/min、0.15MW)、(18r/min、0.85MW)之间的功率曲线细节。
在本实例中,以风机的额定功率1.6MW为基准值,在15min统计周期下,采用该标准功率曲线的出力折算误差的平均误差0.0000%,均方根误差1.5702%,折算误差绝对值的最大值2.4295%。可参见后面图10、图11。
(4)常见的水平轴风机设计中,出力p和风轮转速n之间为分段的确定型函数关系。在风速到达额定风速之前,当风速在一定范围变化时,风轮转速增大则引起出力增大。
如某些类型的风机,
这里Ak,k=1,2,3,4是常数。n的单位r/min。
考虑机械惯性的更精确的出力折算,可以分为两种途径:
1)将连续2点以上的未来有效风速veq,j+1、veq,j作为输入变量,输入到支持向量机、深度学习等智能模型,直接折算出未来的出力pj+1。尽管这是一类通用的方法,但当风速在分段函数特性转折点周围波动时,以及历史数据中的未能有效识别的不良数据等,会引起较为明显的误差。
2)在转速和出力为正相关的情况下,借助刚体定轴转动定律,从veq,j+1、veq,j等折算出pj+1。尽管这样会受到某些情况的限制,但是这是一种机理性的方法。可以有效区分分段函数特性转折点,并抵抗不良数据历史数据(不明显依赖历史数据)的不利影响。
在实际中,可以用下述方程很好地逼近上面的刚体定轴转动方程:
这里系数a、b是实数常数。与功率系数曲线CP类似,系数a、b随风速的变化有一定的变化。
a、b可以从历史观测值中获取。将上述方程离散化,即可得
对于符合发电机的分段特性的风速、出力和转速,采用连续多点观测值,即可回归出常数a、b。
假设当前时刻为(i+1),采用最近的连续S点历史观测值,回归常数a、b的具体输入数据格式如下:
图7为本实例中回归常数a、b的变化情况。
(5)之后,再采用精度高且稳定性好的梯形法来求解上面的微分方程。梯形法求解的公式为这里Δt为出力等变量的统计周期;系数a、b是实数常数。一般地,记时间点i=1,2,…。在由未来风速折算出力时,时间点记为j=1,2,…。
该迭代公式里有两个相互依存的未知量ni+1、pi+1。在水平轴风机设计时,它们之间设计成分段的显式函数关系。
上式可以进一步变形为:
与标准的功率曲线相比,上式的pi+1由veq,i+1、veq,i两个自变量引起,这就考虑了风速变化时风轮等机械惯性引起的出力偏移。
显然,在风电功率预测要求的15min统计周期下,由上述梯形法求解的定轴刚体转动定律方程效果并不会好。这是因为选用Δt=15min作为微分方程求解的步长太长了。实际上,选用Δt=20s是比较合适的步长(各种历史数据的统计周期)。
在本专利中,采用5种具体模型求解梯形法公式。即三次多项式拟合、最小二乘支持向量机LSSVR、随机森林RF,以及深度学习里的长短期记忆网络LSTM、深度信念网络DBN等,得到pj+1的5种单项折算结果,即pCUBIC,j+1、pLSSVR,j+1、pRF,j+1、pLSTM,j+1、pDBN,j+1。
第一:采用三次多项式拟合
在风机设计时,发电机转速和出力之间通常有确定的显示函数关系。发电机转速和风轮转速n之间有确定的关系(通常由齿轮箱传动)。因此风轮转速n可以用出力p表示出。类似地,也可以用出力p表示出。本专利采用3次多项式拟合它们,即:
n=A3·p3+A2·p2+A1·p+A0
这里Ak,Bk,k=1,2,3,4是常数。Ak,Bk可以从历史数据中拟合出,它们实质上是发电机的设计性能造成的分段显式函数;这些显式函数关系不随风速变化。因此,可以提前拟合出。
这样,在将未来有效风速的veq,j+1、veq,j,以及由功率曲线计算出出力pPC,j代入
图9是本实例中梯形法中三次多项式拟合折算的出力,和对应的出力历史值。历史数据采用20s统计周期。
以风机的额定功率1.6MW为基准值,在15min统计周期下,采用三次多项式拟合的出力折算误差的平均误差-0.4551%,均方根误差1.4976%,折算误差绝对值的最大值2.6735%。
第二:采用最小二乘支持向量机LSSVR、随机森林RF,以及深度学习里的长短期记忆网络LSTM、深度信念网络DBN求解pj+1的各个单项折算值时,先采用连续S个历史数据,以
Y1×S=(pj-S+1,…,pj-1,pj)
方式训练出智能模型F,即Y1×S=F(X3×S)。
之后将veq,j+1、veq,j,以及pPC,j代入模型F即可,即p*,j+1=F(veq,j+1,veq,j,pPC,j)。这里用*表示上面的4种智能模型。
图10是梯形法中4种智能模型折算的出力,和对应的出力历史值。
这4种智能模型折算的误差见表1。
表1实例中各种方法的出力折算误差(对额定功率的相对误差)
(7)采用偏最小二乘回归,将上面5种单项折算的出力结果pCUBIC,j+1、pLSSVR,j+1、pRF,j+1、pLSTM,j+1、pDBN,j+1,再优化组合成一个最终采用的出力值pj+1。
图11是本实例的最终出力折算值。
以风机的额定功率1.6MW为基准值,在15min统计周期下,采用误差PLSR组合后的最终折算误差:平均误差-0.4949%,均方根误差0.8584%,折算误差绝对值的最大值1.6138%。与标准的功率曲线方法相比,均方根误差降低了0.7118%,折算误差绝对值的最大值降低了0.8157%。
Claims (2)
1.一种考虑风机转动惯量的风速折算出力的精确方法,包括以下步骤:
(2)风机的历史数据的预处理:对上述各种历史观测值进行有效性检验,数据修复预处理;
(4)水平轴风机设计中,出力p和风轮转速n之间为分段的确定型函数关系,在风速到达额定风速之前,当风速在一定范围变化时,风轮转速增大则引起出力增大;考虑机械惯性的更精确的出力折算,根据刚体定轴转动定律,忽略阻尼,建立风轮机械运动的时域微分方程:其中,J是叶片等的等效机械转动惯量,为风轮转速n对时间t的导数,为瞬时有效风速veq引起的动力转矩,TP为发电机出力折算的阻力转矩;
(6)对未来时间点j=1,2,…,将有效风速的veq,j+1、veq,j,以及由功率曲线折算的出力pPC,j,代入该梯形法公式,得到出力的折算值pj+1;
(7)分别采用三次多项式拟合、最小二乘支持向量机LSSVR、随机森林RF,以及深度学习里的长短期记忆网络LSTM、深度信念网络DBN共计5种具体模型求解梯形法的公式,得到pj+1的5种单项折算结果,即pCUBIC,j+1、pLSSVR,j+1、pRF,j+1、pLSTM,j+1、pDBN,j+1;
(8)采用偏最小二乘回归,5种单项折算的出力结果pCUBIC,j+1、pLSSVR,j+1、pRF,j+1、pLSTM,j+1、pDBN,j+1,优化组合成一个最终采用的出力值pj+1。
2.根据权利要求1所述的方法,其特征在于,步骤(2)的预处理包括对历史风速、风向时间序列进行3点平滑去噪、双正交小波去噪。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110781808.8A CN113505483A (zh) | 2021-07-09 | 2021-07-09 | 一种考虑风机转动惯量的风速折算出力的精确方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110781808.8A CN113505483A (zh) | 2021-07-09 | 2021-07-09 | 一种考虑风机转动惯量的风速折算出力的精确方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113505483A true CN113505483A (zh) | 2021-10-15 |
Family
ID=78012216
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110781808.8A Pending CN113505483A (zh) | 2021-07-09 | 2021-07-09 | 一种考虑风机转动惯量的风速折算出力的精确方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113505483A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115964816A (zh) * | 2022-12-06 | 2023-04-14 | 曲阜师范大学 | 基于深度强化学习的大功率风电磁齿轮箱优化设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103268366A (zh) * | 2013-03-06 | 2013-08-28 | 辽宁省电力有限公司电力科学研究院 | 一种适用于分散式风电场的组合风电功率预测方法 |
CN110457800A (zh) * | 2019-07-30 | 2019-11-15 | 天津大学 | 考虑机械惯性的水平轴风机风速出力折算方法 |
-
2021
- 2021-07-09 CN CN202110781808.8A patent/CN113505483A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103268366A (zh) * | 2013-03-06 | 2013-08-28 | 辽宁省电力有限公司电力科学研究院 | 一种适用于分散式风电场的组合风电功率预测方法 |
CN110457800A (zh) * | 2019-07-30 | 2019-11-15 | 天津大学 | 考虑机械惯性的水平轴风机风速出力折算方法 |
Non-Patent Citations (3)
Title |
---|
刘丁一: ""提高风电功率预测效果的改进多步法和风速出力折算"", 《万方数据》 * |
李真真: ""风力发电机组推力曲线分散性的数据分析"", 《硕士电子期刊 工程科技Ⅱ辑》 * |
杨正瓴 等: ""基于季风和大气压分布的我国风电功率预测研究"", 《分布式能源》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115964816A (zh) * | 2022-12-06 | 2023-04-14 | 曲阜师范大学 | 基于深度强化学习的大功率风电磁齿轮箱优化设计方法 |
CN115964816B (zh) * | 2022-12-06 | 2024-02-09 | 曲阜师范大学 | 基于深度强化学习的大功率风电磁齿轮箱优化设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104156575B (zh) | 基于测风塔数据外推法的风电场理论功率计算方法 | |
CN109667727B (zh) | 一种基于功率曲线分析的风电机组偏航误差固有偏差辨识及补偿方法 | |
Fuglsang et al. | Optimization method for wind turbine rotors | |
AU2020399503B2 (en) | Method and apparatus for detecting yaw-to-wind abnormality, and device and storage medium thereof | |
DE112016005511T5 (de) | Verfahren zum Bewerten von Leistungscharakteristiken von Windkraftanlagen, Vorrichtung und Speichermedium | |
CN109063276B (zh) | 适用于长时域风速波动的风电场动态等值建模方法 | |
CN103887815A (zh) | 基于运行数据的风电场参数辨识及动态等值方法 | |
CN107944175B (zh) | 一种考虑风湍流强度的风机真实功率曲线获取方法 | |
CN105649875B (zh) | 风力发电机组的变桨控制方法及装置 | |
CN111322205B (zh) | 风电机组风向标零位误差辨识方法及校正方法 | |
WO2019165743A1 (zh) | 确定对风角度偏差及修正对风角度的方法、装置和系统 | |
CN109002650A (zh) | 一种风电机组功率曲线建模方法 | |
CN107002636A (zh) | 用于估计风速,包括计算针对叶片扭转调节的桨距角的方法 | |
CN103984986B (zh) | 实时校正的自学习arma模型风电功率超短期预测方法 | |
CN113935534A (zh) | 一种基于卷积transformer架构的光伏日前预测方法、系统及设备 | |
CN113505483A (zh) | 一种考虑风机转动惯量的风速折算出力的精确方法 | |
Porto et al. | Power performance of starting-improved and multi-bladed horizontal-axis small wind turbines | |
CN111720271B (zh) | 一种风电机组载荷在线预测的智能方法及风电机组 | |
Van Der Laan et al. | Faster wind farm AEP calculations with CFD using a generalized wind turbine model | |
Kumar et al. | Global sensitivity analysis of model uncertainty in aeroelastic wind turbine models | |
CN110457800B (zh) | 考虑机械惯性的水平轴风机风速出力折算方法 | |
CN117251995A (zh) | 基于可变遗忘因子最小二乘法的双馈风机惯量评估方法 | |
CN103903087A (zh) | 一种基于bp神经网络的汽动引风机全工况在线监测方法 | |
Li et al. | Wind power forecasting based on time series and neural network | |
CN116306236A (zh) | 一种基于gru-lstm-fc组合网络的风电场动态等值建模方法及系统 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211015 |
|
RJ01 | Rejection of invention patent application after publication |