This article aims to reveal the transient thermal characteristics of the solid domain in a nozzle flowmeter when measuring fluids of varying temperatures.Based on finite element method,the transient numerical calculat...This article aims to reveal the transient thermal characteristics of the solid domain in a nozzle flowmeter when measuring fluids of varying temperatures.Based on finite element method,the transient numerical calculation of the thermal characteristics of each component of the nozzle flowmeter has been conducted.The research shows that:as the fluid passes through the flowmeter,the high heat flux area inside the nozzle flowmeter gradually transfer from the center of the nozzle to the inlet and outlet,as well as the pressure tapping points upstream and downstream;High thermal stress zones are present near the upstream and downstream pressure tapping points and inlet and outlet area;High thermal deformation zones occur near the upstream and downstream pressure tapping points and eight slot nozzle.展开更多
With advantages of strong drive capability,nested-loop secondary linear machine(NLS-LM)has great potentiality in linear metro.For its secondary structure with multiple loops,it is difficult to calculate the electromag...With advantages of strong drive capability,nested-loop secondary linear machine(NLS-LM)has great potentiality in linear metro.For its secondary structure with multiple loops,it is difficult to calculate the electromagnetic thrust of NLS-LM reasonably.Hence,in this paper,one thrust calculation method is proposed considering variable loop inductance and transient loop current.Firstly,to establish the secondary winding function,the modeling domain is confined to a limited range,and the equivalent loop span is employed by analyzing the coupling relationship between primary and secondary.Then,in order to obtain the secondary flux density,the transient secondary current is solved based on the loop impedance and induced voltage.Finally,the electromagnetic thrust can be calculated reasonably by the given primary current sheet and the calculated secondary flux density.Comprehensive simulations and experiments have demonstrated the effectiveness of the proposed method.展开更多
A new method of 3D transient eddy current field calculation is proposed. The Maxwell equations with time component elimination (METCE) are derived under the assumption of magnetic quasi static approximation, especia...A new method of 3D transient eddy current field calculation is proposed. The Maxwell equations with time component elimination (METCE) are derived under the assumption of magnetic quasi static approximation, especially for the sample of low conductivity. Based on METCE, we deduce a more efficient reconstruction algorithm of a 3D transient eddy current field. The computational burden is greatly reduced through the new algorithm, and the computational efficiency is improved. This new algorithm decompounds the space-time variables into two individual variables. The idea is to solve the spatial vector component firstly, and then multiply it by the corresponded time component. The iterative methods based on METCE are introduced to recover the distribution of conductivity in magneto-acoustic tomography. The reconstructed images of conductivity are consistent with the original distribution, which validate the new method.展开更多
From basic equations of gas-liquid, solid-liquid, solid-gas two-phase flow,the calculating method on flow transients of two-phase flow is developed by means of characteristicmethod. As one example, a gas-liquid flow t...From basic equations of gas-liquid, solid-liquid, solid-gas two-phase flow,the calculating method on flow transients of two-phase flow is developed by means of characteristicmethod. As one example, a gas-liquid flow transient is calculated and it agrees well with theexperimental result. It is shown that the method is satisfactory for engineering demand.展开更多
Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given. Through...Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given. Through the experiments of one dimension transient flow of HPAM (Hydrolytic Polyacrylamide) foam in the artificial sandstone core, the HPAM foam generation and coalescence coefficient of the mathematical model mentioned above are determined. The profiles of the liquid phase saturation, the pressure drop and the number density of one dimension transient flow of HPAM foam with the dimensionless time in artificial sandstone core are numerically calculated and analyzed by using the numerical calculation method.展开更多
To perform nuclear reactor simulations in a more realistic manner,the coupling scheme between neutronics and thermal-hydraulics was implemented in the HNET program for both steady-state and transient conditions.For si...To perform nuclear reactor simulations in a more realistic manner,the coupling scheme between neutronics and thermal-hydraulics was implemented in the HNET program for both steady-state and transient conditions.For simplicity,efficiency,and robustness,the matrixfree Newton/Krylov(MFNK)method was applied to the steady-state coupling calculation.In addition,the optimal perturbation size was adopted to further improve the convergence behavior of the MFNK.For the transient coupling simulation,the operator splitting method with a staggered time mesh was utilized to balance the computational cost and accuracy.Finally,VERA Problem 6 with power and boron perturbation and the NEACRP transient benchmark were simulated for analysis.The numerical results show that the MFNK method can outperform Picard iteration in terms of both efficiency and robustness for a wide range of problems.Furthermore,the reasonable agreement between the simulation results and the reference results for the NEACRP transient benchmark verifies the capability of predicting the behavior of the nuclear reactor.展开更多
In the present paper the analytical formulas for calculating the equivalent deterministic transients for multivariable cross correlated random processes are developed.The formulas permit the determination of the Root-...In the present paper the analytical formulas for calculating the equivalent deterministic transients for multivariable cross correlated random processes are developed.The formulas permit the determination of the Root-mean-square of the responses of a linear time-invariant system to stationary multiple random inputs in the time domain.The method is applicable in the study of flight of airplanes in atmospheric turbulence and is also useful for general engineering applications of stochastc processes control.展开更多
针对油浸式电力变压器瞬态温升计算效率过低的问题,该文提出本征正交分解-αATS(proper orthogonal decomposition-adaptive time stepping based onαfactor,POD-αATS)降阶自适应变步长瞬态计算方法。首先,推导变压器绕组瞬态温升计...针对油浸式电力变压器瞬态温升计算效率过低的问题,该文提出本征正交分解-αATS(proper orthogonal decomposition-adaptive time stepping based onαfactor,POD-αATS)降阶自适应变步长瞬态计算方法。首先,推导变压器绕组瞬态温升计算的有限元离散方程;其次,采用POD降阶算法改善传统瞬态计算中存在的条件数过大及方程阶数过高的问题;同时对于瞬态计算中的时间步长选择问题,提出适用于非线性问题的αATS变步长策略;然后,为验证方法的有效性,基于110 kV油浸式电力变压器绕组的基本结构建立二维八分区数值计算模型,同时将计算结果与基于110 kV绕组的温升实验结果进行对比。数值计算及实验结果表明,所提算法与全阶定步长算法在流场和温度场中的精度几乎相同,且流场计算效率提升约45倍,温度场计算效率提升约38倍,计算速度得到显著提高。这一点在温升实验中同样得到验证,说明该文所提算法的准确性、高效性及一定的工程实用性。展开更多
Since the fault dynamic of droop-controlled inverter is different from synchronous generators (SGs), protection devices may become invalid, and the fault overcurrent may damage power electronic devices and threaten th...Since the fault dynamic of droop-controlled inverter is different from synchronous generators (SGs), protection devices may become invalid, and the fault overcurrent may damage power electronic devices and threaten the safety of the microgrid. Therefore, it is imperative to conduct a comprehensive fault analysis of the inverter to guide the design of protection schemes. However, due to the complexity of droop control strategy, existing literatures have simplified asymmetric fault analysis of droop-controlled inverters to varying degrees. Therefore, accurate fault analysis of a droop-controlled inverter is needed. In this paper, by analyzing the control system, an accurate fault model is established. Based on this, a calculation method for instantaneous asymmetrical fault current is proposed. In addition, the current components and current characteristics are analyzed. It was determined that fault currents are affected by control loops, fault types, fault distance and nonlinear limiters. In particular, the influences of limiters on the fault model, fault current calculation and fault current characteristics were analyzed. Through detailed analysis, it was found that dynamics of the control loop cannot be ignored, the fault type and fault distance determine fault current level, and part of the limiters will totally change the fault current trend. Finally, calculation and experimental results verify the correctness of the proposed method.展开更多
基于电缆外部精确电磁场分析方法,提出了一种拓展传输线理论,能够计及广义大地参数和电磁暂态仿真的影响计算电缆输入阻抗。采用电磁暂态仿真程序EMTP-RV,利用拓展传输线理论对典型高压交流电缆开展输入阻抗和金属护套电压特性研究,研...基于电缆外部精确电磁场分析方法,提出了一种拓展传输线理论,能够计及广义大地参数和电磁暂态仿真的影响计算电缆输入阻抗。采用电磁暂态仿真程序EMTP-RV,利用拓展传输线理论对典型高压交流电缆开展输入阻抗和金属护套电压特性研究,研究当合闸于空载电缆时的金属护套暂态电压特性;通过对正序激励下的电缆输入阻抗模值-频率响应进行研究,分析由三相金属护套交叉互联产生的多种传播模态。采用拓展传输线理论计算的结果与使用基于矩量法的改进数值电磁场分析方法(method of moments-surface admittanceoperator,MoM-SO)的计算结果基本拟合,与现有电磁暂态仿真软件中(ATP-EMTP、PSCAD/EMTDC和EMTP-RV)的CableConstants子程序计算结果出现显著性差异。当合闸于空载电缆线路时,受大地电阻率影响,采用Cable Constants子程序计算出的电缆金属护套最大暂态电压与采用拓展传输线理论所得出的结果相比,计算偏差位于1.8%至135%之间。考虑输入阻抗计算和暂态仿真的关键性差异,研究内容对今后开展高压交流电缆参数计算、模态分析、电磁暂态仿真、过电压研究具有参考意义。展开更多
In order to further achieve the balance between the calculation accuracy and efficiency of the transient analysis of the aero-engine disc cavity system,an Optimized Time-adaptive Aerother-mal Coupling calculation(OTAC...In order to further achieve the balance between the calculation accuracy and efficiency of the transient analysis of the aero-engine disc cavity system,an Optimized Time-adaptive Aerother-mal Coupling calculation(OTAC)method has been proposed.It combines one-dimensional tran-sient calculation of air system,Conventional Sequence Staggered(CSS)method,Time-adaptive Aerothermal Coupling calculation(TAC)method and differential evolution optimization algorithm to obtain an efficient and high-precision aerothermal coupling calculation method of air system.Considering both the heat conduction in the solid domain and the flow in the fluid domain as unsteady states in the OTAC,the interaction of fluid-solid information within a single coupling time step size was implemented based on the CSS method.Furthermore,the coupling time step size was automatically adjusted with the number of iterations by using the Proportional-Integral-Deri vative(PID)controller.Results show that when compared with the traditional loosely coupling method with a fixed time step size,the computational accuracy and efficiency of the OTAC method are improved by 8.9%and 30%,respectively.Compared with the tight coupling calculation,the OTAC method can achieve a speedup of 1 to 2 orders of magnitude,while the calculation error is maintained within 6.1%.展开更多
基金supported by the“Pioneer”and“Leading Goose”R&D Program of Zhejiang (Grant No.2022C03170)Science and Technology Project of Quzhou (Grant No.2022K98)Science and Technology Project of Zhejiang (Grant No.LGC21E050001).
文摘This article aims to reveal the transient thermal characteristics of the solid domain in a nozzle flowmeter when measuring fluids of varying temperatures.Based on finite element method,the transient numerical calculation of the thermal characteristics of each component of the nozzle flowmeter has been conducted.The research shows that:as the fluid passes through the flowmeter,the high heat flux area inside the nozzle flowmeter gradually transfer from the center of the nozzle to the inlet and outlet,as well as the pressure tapping points upstream and downstream;High thermal stress zones are present near the upstream and downstream pressure tapping points and inlet and outlet area;High thermal deformation zones occur near the upstream and downstream pressure tapping points and eight slot nozzle.
基金supported in part by the National Natural Science Foundation of China under Grants 52277050the Shenzhen International Collaboration under Grant GJHZ20210705142539007。
文摘With advantages of strong drive capability,nested-loop secondary linear machine(NLS-LM)has great potentiality in linear metro.For its secondary structure with multiple loops,it is difficult to calculate the electromagnetic thrust of NLS-LM reasonably.Hence,in this paper,one thrust calculation method is proposed considering variable loop inductance and transient loop current.Firstly,to establish the secondary winding function,the modeling domain is confined to a limited range,and the equivalent loop span is employed by analyzing the coupling relationship between primary and secondary.Then,in order to obtain the secondary flux density,the transient secondary current is solved based on the loop impedance and induced voltage.Finally,the electromagnetic thrust can be calculated reasonably by the given primary current sheet and the calculated secondary flux density.Comprehensive simulations and experiments have demonstrated the effectiveness of the proposed method.
基金Supported by the National Natural Science Foundation of China under Grant Nos 51137004,61427806 and 51577184the Equipment Development Project of Chinese Academy of Sciences under Grant No YZ201507
文摘A new method of 3D transient eddy current field calculation is proposed. The Maxwell equations with time component elimination (METCE) are derived under the assumption of magnetic quasi static approximation, especially for the sample of low conductivity. Based on METCE, we deduce a more efficient reconstruction algorithm of a 3D transient eddy current field. The computational burden is greatly reduced through the new algorithm, and the computational efficiency is improved. This new algorithm decompounds the space-time variables into two individual variables. The idea is to solve the spatial vector component firstly, and then multiply it by the corresponded time component. The iterative methods based on METCE are introduced to recover the distribution of conductivity in magneto-acoustic tomography. The reconstructed images of conductivity are consistent with the original distribution, which validate the new method.
文摘From basic equations of gas-liquid, solid-liquid, solid-gas two-phase flow,the calculating method on flow transients of two-phase flow is developed by means of characteristicmethod. As one example, a gas-liquid flow transient is calculated and it agrees well with theexperimental result. It is shown that the method is satisfactory for engineering demand.
文摘Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given. Through the experiments of one dimension transient flow of HPAM (Hydrolytic Polyacrylamide) foam in the artificial sandstone core, the HPAM foam generation and coalescence coefficient of the mathematical model mentioned above are determined. The profiles of the liquid phase saturation, the pressure drop and the number density of one dimension transient flow of HPAM foam with the dimensionless time in artificial sandstone core are numerically calculated and analyzed by using the numerical calculation method.
基金supported by the China Postdoctoral Science Foundation(No.2021M703045)the National Natural Science Foundation of China(No.12075067)the National Key R&D Program of China(No.2018YFE0180900).
文摘To perform nuclear reactor simulations in a more realistic manner,the coupling scheme between neutronics and thermal-hydraulics was implemented in the HNET program for both steady-state and transient conditions.For simplicity,efficiency,and robustness,the matrixfree Newton/Krylov(MFNK)method was applied to the steady-state coupling calculation.In addition,the optimal perturbation size was adopted to further improve the convergence behavior of the MFNK.For the transient coupling simulation,the operator splitting method with a staggered time mesh was utilized to balance the computational cost and accuracy.Finally,VERA Problem 6 with power and boron perturbation and the NEACRP transient benchmark were simulated for analysis.The numerical results show that the MFNK method can outperform Picard iteration in terms of both efficiency and robustness for a wide range of problems.Furthermore,the reasonable agreement between the simulation results and the reference results for the NEACRP transient benchmark verifies the capability of predicting the behavior of the nuclear reactor.
文摘In the present paper the analytical formulas for calculating the equivalent deterministic transients for multivariable cross correlated random processes are developed.The formulas permit the determination of the Root-mean-square of the responses of a linear time-invariant system to stationary multiple random inputs in the time domain.The method is applicable in the study of flight of airplanes in atmospheric turbulence and is also useful for general engineering applications of stochastc processes control.
文摘针对油浸式电力变压器瞬态温升计算效率过低的问题,该文提出本征正交分解-αATS(proper orthogonal decomposition-adaptive time stepping based onαfactor,POD-αATS)降阶自适应变步长瞬态计算方法。首先,推导变压器绕组瞬态温升计算的有限元离散方程;其次,采用POD降阶算法改善传统瞬态计算中存在的条件数过大及方程阶数过高的问题;同时对于瞬态计算中的时间步长选择问题,提出适用于非线性问题的αATS变步长策略;然后,为验证方法的有效性,基于110 kV油浸式电力变压器绕组的基本结构建立二维八分区数值计算模型,同时将计算结果与基于110 kV绕组的温升实验结果进行对比。数值计算及实验结果表明,所提算法与全阶定步长算法在流场和温度场中的精度几乎相同,且流场计算效率提升约45倍,温度场计算效率提升约38倍,计算速度得到显著提高。这一点在温升实验中同样得到验证,说明该文所提算法的准确性、高效性及一定的工程实用性。
基金supported by National Natural Science Foundation of China under Grant 51977066。
文摘Since the fault dynamic of droop-controlled inverter is different from synchronous generators (SGs), protection devices may become invalid, and the fault overcurrent may damage power electronic devices and threaten the safety of the microgrid. Therefore, it is imperative to conduct a comprehensive fault analysis of the inverter to guide the design of protection schemes. However, due to the complexity of droop control strategy, existing literatures have simplified asymmetric fault analysis of droop-controlled inverters to varying degrees. Therefore, accurate fault analysis of a droop-controlled inverter is needed. In this paper, by analyzing the control system, an accurate fault model is established. Based on this, a calculation method for instantaneous asymmetrical fault current is proposed. In addition, the current components and current characteristics are analyzed. It was determined that fault currents are affected by control loops, fault types, fault distance and nonlinear limiters. In particular, the influences of limiters on the fault model, fault current calculation and fault current characteristics were analyzed. Through detailed analysis, it was found that dynamics of the control loop cannot be ignored, the fault type and fault distance determine fault current level, and part of the limiters will totally change the fault current trend. Finally, calculation and experimental results verify the correctness of the proposed method.
文摘基于电缆外部精确电磁场分析方法,提出了一种拓展传输线理论,能够计及广义大地参数和电磁暂态仿真的影响计算电缆输入阻抗。采用电磁暂态仿真程序EMTP-RV,利用拓展传输线理论对典型高压交流电缆开展输入阻抗和金属护套电压特性研究,研究当合闸于空载电缆时的金属护套暂态电压特性;通过对正序激励下的电缆输入阻抗模值-频率响应进行研究,分析由三相金属护套交叉互联产生的多种传播模态。采用拓展传输线理论计算的结果与使用基于矩量法的改进数值电磁场分析方法(method of moments-surface admittanceoperator,MoM-SO)的计算结果基本拟合,与现有电磁暂态仿真软件中(ATP-EMTP、PSCAD/EMTDC和EMTP-RV)的CableConstants子程序计算结果出现显著性差异。当合闸于空载电缆线路时,受大地电阻率影响,采用Cable Constants子程序计算出的电缆金属护套最大暂态电压与采用拓展传输线理论所得出的结果相比,计算偏差位于1.8%至135%之间。考虑输入阻抗计算和暂态仿真的关键性差异,研究内容对今后开展高压交流电缆参数计算、模态分析、电磁暂态仿真、过电压研究具有参考意义。
基金support of the National Natural Science Foundation of China (No.52007002)the Science Center for Gas Turbine Project,China (No.P2022-A-II-007-001)the Fundamental Research Funds for the Central Universities,China (No.NS2023010).
文摘In order to further achieve the balance between the calculation accuracy and efficiency of the transient analysis of the aero-engine disc cavity system,an Optimized Time-adaptive Aerother-mal Coupling calculation(OTAC)method has been proposed.It combines one-dimensional tran-sient calculation of air system,Conventional Sequence Staggered(CSS)method,Time-adaptive Aerothermal Coupling calculation(TAC)method and differential evolution optimization algorithm to obtain an efficient and high-precision aerothermal coupling calculation method of air system.Considering both the heat conduction in the solid domain and the flow in the fluid domain as unsteady states in the OTAC,the interaction of fluid-solid information within a single coupling time step size was implemented based on the CSS method.Furthermore,the coupling time step size was automatically adjusted with the number of iterations by using the Proportional-Integral-Deri vative(PID)controller.Results show that when compared with the traditional loosely coupling method with a fixed time step size,the computational accuracy and efficiency of the OTAC method are improved by 8.9%and 30%,respectively.Compared with the tight coupling calculation,the OTAC method can achieve a speedup of 1 to 2 orders of magnitude,while the calculation error is maintained within 6.1%.