A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variable...A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.展开更多
Actuator faults can be critical in turbofan engines as they can lead to stall,surge,loss of thrust and failure of speed control.Thus,fault diagnosis of gas turbine actuators has attracted considerable attention,from b...Actuator faults can be critical in turbofan engines as they can lead to stall,surge,loss of thrust and failure of speed control.Thus,fault diagnosis of gas turbine actuators has attracted considerable attention,from both academia and industry.However,the extensive literature that exists on this topic does not address identifying the severity of actuator faults and focuses mainly on actuator fault detection and isolation.In addition,previous studies of actuator fault identification have not dealt with multiple concurrent faults in real time,especially when these are accompanied by sudden failures under dynamic conditions.This study develops component-level models for fault identification in four typical actuators used in high-bypass ratio turbofan engines under both dynamic and steady-state conditions and these are then integrated with the engine performance model developed by the authors.The research results reported here present a novel method of quantifying actuator faults using dynamic effect compensation.The maximum error for each actuator is less than0.06%and 0.07%,with average computational time of less than 0.0058 s and 0.0086 s for steady-state and transient cases,respectively.These results confirm that the proposed method can accurately and efficiently identify concurrent actuator fault for an engine operating under either transient or steady-state conditions,even in the case of a sudden malfunction.The research results emonstrate the potential benefit to emergency response capabilities by introducing this method of monitoring the health of aero engines.展开更多
The alternative working modes and flexible working states are the outstanding features of an adaptive cycle engine, with a proper control schedule design being the only way to exploit the performance of such an engine...The alternative working modes and flexible working states are the outstanding features of an adaptive cycle engine, with a proper control schedule design being the only way to exploit the performance of such an engine. However, unreasonable design in the control schedule causes not only performance deterioration but also serious aerodynamic stability problems. Thus, in this work,a hybrid optimization method that automatically chooses the working modes and identifies the optimal and smooth control schedules is proposed, by combining the differential evolution algorithm and the Latin hypercube sampling method. The control schedule architecture does not only optimize the engine steady-state performance under different working modes but also solves the control-schedule discontinuity problem, especially during mode transition. The optimal control schedules are continuous and almost monotonic, and hence are strongly suitable for a control system, and are designed for two different working conditions, i.e., supersonic and subsonic throttling,which proves that the proposed hybrid method applies to various working conditions. The evaluation demonstrates that the proposed control method optimizes the engine performance, the surge margin of the compression components, and the range of the thrust during throttling.展开更多
A Newton multigrid method is developed for one-dimensional （1D） and two- dimensional （2D） steady-state shallow water equations （SWEs） with topography and dry areas. The nonlinear system arising from the well-bal...A Newton multigrid method is developed for one-dimensional （1D） and two- dimensional （2D） steady-state shallow water equations （SWEs） with topography and dry areas. The nonlinear system arising from the well-balanced finite volume discretization of the steady-state SWEs is solved by the Newton method as the outer iteration and a geometric multigrid method with the block symmetric Gauss-Seidel smoother as the inner iteration. The proposed Newton multigrid method makes use of the local residual to regularize the Jacobian matrix of the Newton iteration, and can handle the steady- state problem with wet/dry transition. Several numerical experiments are conducted to demonstrate the efficiency, robustness, and well-balanced property of the proposed method. The relation between the convergence behavior of the Newton multigrid method and the distribution of the eigenvalues of the iteration matrix is detailedly discussed.展开更多
In this paper,we considered the improved element-free Galerkin(IEFG)method for solving 2D anisotropic steadystate heat conduction problems.The improved moving least-squares(IMLS)approximation is used to establish the ...In this paper,we considered the improved element-free Galerkin(IEFG)method for solving 2D anisotropic steadystate heat conduction problems.The improved moving least-squares(IMLS)approximation is used to establish the trial function,and the penalty method is applied to enforce the boundary conditions,thus the final discretized equations of the IEFG method for anisotropic steady-state heat conduction problems can be obtained by combining with the corresponding Galerkin weak form.The influences of node distribution,weight functions,scale parameters and penalty factors on the computational accuracy of the IEFG method are analyzed respectively,and these numerical solutions show that less computational resources are spent when using the IEFG method.展开更多
In this paper, we design high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (multi-resolution WENO) limiters to compute compressible steady-state pr...In this paper, we design high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (multi-resolution WENO) limiters to compute compressible steady-state problems on triangular meshes. A troubled cell indicator extended from structured meshes to unstructured meshes is constructed to identify triangular cells in which the application of the limiting procedures is required. In such troubled cells, the multi-resolution WENO limiting methods are used to the hierarchical L^(2) projection polynomial sequence of the DG solution. Through using the RKDG methods with multi-resolution WENO limiters, the optimal high-order accuracy can be gradually reduced to first-order in the triangular troubled cells, so that the shock wave oscillations can be well suppressed. In steady-state simulations on triangular meshes, the numerical residual converges to near machine zero. The proposed spatial reconstruction methods enhance the robustness of classical DG methods on triangular meshes. The good results of these RKDG methods with multi-resolution WENO limiters are verified by a series of two-dimensional steady-state problems.展开更多
In this paper,we consider solving the topology optimization for steady-state incompressibleNavier-Stokes problems via a new topology optimization method called parameterized level set method,which can maintain a relat...In this paper,we consider solving the topology optimization for steady-state incompressibleNavier-Stokes problems via a new topology optimization method called parameterized level set method,which can maintain a relatively smooth level set function with a local optimality condition.The objective of topology optimization is tond an optimal conguration of theuid and solid materials that minimizes power dissipation under a prescribeduid volume fraction constraint.An articial friction force is added to the Navier-Stokes equations to apply the no-slip boundary condition.Although a great deal of work has been carried out for topology optimization ofuidow in recent years,there are few researches on the topology optimization ofuidow with physical body forces.To simulate theuidow in reality,the constant body force(e.g.,gravity)is considered in this paper.Several 2D numerical examples are presented to discuss the relationships between the proposed method with Reynolds number and initial design,and demonstrate the feasibility and superiority of the proposed method in dealing with unstructuredmesh problems.Three 3D numerical examples demonstrate the proposedmethod is feasible in three-dimensional.展开更多
An analysis technique of steady state and stability for closed-loop PWM DC/DC switching converters is presented. Using this method, the closed-loop switching converter is transformed into an open-loop system. By means...An analysis technique of steady state and stability for closed-loop PWM DC/DC switching converters is presented. Using this method, the closed-loop switching converter is transformed into an open-loop system. By means of the fact that in steady state, the two boundary values are equal in one switching period. The exponential matrix is evaluated by precise time-domain-integration method, and then the related curve between feedback duty cycle and the input one is obtained. Not only can the steady-state duty cycle be found from the curve, but also the stability and stable domain of the system. Compared with other methods, it features with simplicity and less calculation, and fit for numerical simulation and analysis for closed-loop switching converters. The simulation results of examples indicate the correctness of the presented method.展开更多
Based on the indirect Trefftz approach, a wave number method （WNM） is proposed to deal with three-dimensional steady-state acoustic problems. In the WNM, the dynamic pressure response variable is approximated by a s...Based on the indirect Trefftz approach, a wave number method （WNM） is proposed to deal with three-dimensional steady-state acoustic problems. In the WNM, the dynamic pressure response variable is approximated by a set of wave functions, which exactly satisfy the Helmholtz equation. The set of wave functions comprise the exact solutions of the homogeneous part of the governing equations and some particular solution functions. The unknown coefficients of the wave functions can be obtained by enforcing the pressure approximation to satisfy the boundary conditions. Compared with the boundary element method （BEM）, the WNM have a smaller system matrix, and is applicable to the radiation problems since the wave functions are independent of the domain size. A 3D acoustic cavity is exemplified to show the properties of the method. The results show that the wave number method is more efficient than the BEM, and it is fairly accurate.展开更多
NT build 443, or profile fitting method, is often used to measure the diffusion coefficient of chloride in concrete. However, this method is quite laborious and needs special equipment. Colorimetric method is a quick ...NT build 443, or profile fitting method, is often used to measure the diffusion coefficient of chloride in concrete. However, this method is quite laborious and needs special equipment. Colorimetric method is a quick and simple method to measure the penetration depth of chloride by spraying 0.1 mol/L silver nitrate solution. The objective of this work is to study the possibility of the use of colorimetric method in the calculation of non-steady-state diffusion coefficient. Twelve concrete mixtures with different supplementary cementitious materials and water-to-cement ratios of 0.35, 0.48 and 0.6 were used for study. According to NT build 443, the concrete specimens were immersed in 165 g/L NaC1 （2.8 mol/L） solution for 42 d. Both water-soluble （convert to free chloride） chloride and acid-soluble chloride at different layers of specimens were measured. The results show that the mean value of free chloride concentration at the color change boundary Cd was 0.306 mol/L. The surface free chloride concentration cs was obtained by profile fitting method, which was 40% lower than the chloride concentration of exposure solution after an immersion period of 42 d. Chloride diffusion coefficients obtained by the colorimetric method was not well correlated with those obtained by profile fitting method.展开更多
This study presents a numerical analysis of the steady-state solution for transient magnetohydrodynamic(MHD)dissipative and radiative fluid flow,incorporating an inducedmagnetic field(IMF)and considering a relatively ...This study presents a numerical analysis of the steady-state solution for transient magnetohydrodynamic(MHD)dissipative and radiative fluid flow,incorporating an inducedmagnetic field(IMF)and considering a relatively high concentration of foreign mass(accounting for Soret and Dufour effects)over a vertically oriented semi-infinite plate.The governing equations were normalized using boundary layer(BL)approximations.The resulting nonlinear system of partial differential equations(PDEs)was discretized and solved using an efficient explicit finite difference method(FDM).Numerical simulations were conducted using MATLAB R2015a,and the developed numerical code was verified through comparison with another code written in FORTRAN 6.6a.To ensure the reliability of the results,both mesh refinement and steady-state time validation tests were performed.Furthermore,a comparison with existing published studies was made to confirm the accuracy of the findings.The dimensionless equations revealed the impacts of several key parameters.The IMF initially intensifies near the plate before gradually diminishing as the magnetic parameter increases.For the range 0≤y≤1.8(where y is the horizontal direction),the IMF decreases with a rise in the magnetic Prandtl number;however,for 1.8≤y≤7(approximately),the magnetic field begins to increase.Beyond this,the profile of the magnetic field becomes somewhat irregular through the remaining part of the BL.展开更多
The isogeometric boundary element technique(IGABEM)is presented in this study for steady-state inhomogeneous heat conduction analysis.The physical unknowns in the boundary integral formulations of the governing equati...The isogeometric boundary element technique(IGABEM)is presented in this study for steady-state inhomogeneous heat conduction analysis.The physical unknowns in the boundary integral formulations of the governing equations are discretized using non-uniform rational B-spline(NURBS)basis functions,which are utilized to build the geometry of the structures.To speed up the assessment of NURBS basis functions,the Bezier extraction´approach is used.To solve the extra domain integrals,we use a radial integration approach.The numerical examples show the potential of IGABEM for dimension reduction and smooth integration of CAD and numerical analysis.展开更多
Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed. It lead...Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed. It leads to the high quality locM error bounds in the problem of the direct-solution steady-state dynamic analysis with a frequency-domain finite element, which involves the enrichments with plural variable basis functions. The solution of the steady-state dynamic procedure calculates the harmonic response directly in terms of the physical degrees of freedom in the model, which uses the mass, damping, and stiffness matrices of the system. A three-dimensional finite element example is carried out to illustrate the computational procedures.展开更多
A mathematical model of an amperometric biosensor with the substrate inhibition for steady-state condition is discussed. The model is based on the system of non-stationary diffusion equation containing a non-linear te...A mathematical model of an amperometric biosensor with the substrate inhibition for steady-state condition is discussed. The model is based on the system of non-stationary diffusion equation containing a non-linear term related to non-Michaelis–Menten kinetics of the enzymatic reaction. This paper presents the analytical expression of concentrations and current for all values of parameters φ2s φ2s α and β . Here the Adomian decomposition method (ADM) is used to find the analytical expressions for substrate, product concentration and current. A comparison of the analytical approximation and numerical simulation is also presented. A good agreement between theoretical predictions and numerical results is observed.展开更多
A mathematical modelling of diagnostic biosensors system at three basic types of enzyme kinetics is discussed in the presence of diffusion. Enzyme kinetics is adopted to be first order, Michaelis-Menten and ping-pong ...A mathematical modelling of diagnostic biosensors system at three basic types of enzyme kinetics is discussed in the presence of diffusion. Enzyme kinetics is adopted to be first order, Michaelis-Menten and ping-pong mechanism. In this paper, approximate analytical solutions are obtained for the non-linear equations under steady-state conditions by using the new Homotopy perturbation method. Simple and closed forms of analytical expressions for concentrations of substrate, product and co-substrate and corresponding current response have been derived for all possible values of parameters. Furthermore, the numerical simulation of the problem is also reported here by using Matlab program. Good agreement between analytical and numerical results is noted.展开更多
In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to ...In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.展开更多
The low frequency oscillation is a serious threat to security and stability of a power grid.How to locate the disturbance source accurately is an important issue to low frequency oscillation disposal.Existing methods ...The low frequency oscillation is a serious threat to security and stability of a power grid.How to locate the disturbance source accurately is an important issue to low frequency oscillation disposal.Existing methods have poor adaptability to the low frequency oscillation with time-varying steady-state points because of the limitations in the location criterion derivation.A disturbance source location method on a low frequency oscillation with good generality is presented in the paper.Firstly,the reasons why the steady-state points are time-varying on a low frequency oscillation are analyzed.Then,based on the energy function construction form,the branch transmission energy is decomposed into state energy,reciprocating energy and dissipation energy by mathematical derivation.The flow direction of the dissipation energy shows the source and destination of the disturbance energy,and the specific location of a disturbance source can be identified according to its flow direction.Meanwhile,to meet the needs of energy calculation,a recognition method on the electrical quantities steady-state points is also presented by using the cubic spline interpolation.Simulation results show the correctness of the derivation and analysis on energy structure in the paper,and the disturbance source can be located accurately according to the dissipation energy.展开更多
This paper proposes a steady-state errors correction(SSEC)method for eliminating measurement errors.This method is based on the detections of error signal E(s)and output C(s)which generate an expected output R(s).In c...This paper proposes a steady-state errors correction(SSEC)method for eliminating measurement errors.This method is based on the detections of error signal E(s)and output C(s)which generate an expected output R(s).In comparison with the conventional solutions which are based on detecting the expected output R(s)and output C(s)to obtain error signal E(s),the measurement errors are eliminated even the error might be at a significant level.Moreover,it is possible that the individual debugging by regulating the coefficient K for every member of the multiple objectives achieves the optimization of the open loop gain.Therefore,this simple method can be applied to the weak coupling and multiple objectives system,which is usually controlled by complex controller.The principle of eliminating measurement errors is derived analytically,and the advantages comparing with the conventional solutions are depicted.Based on the SSEC method analysis,an application of this method for an active power filter(APF)is investigated and the effectiveness and viability of the scheme are demonstrated through the simulation and experimental verifications.展开更多
Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery...Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery.We rst derive the projection formulas for a vector onto the feasible sets.The centralized circumcentered-reection method is designed to solve the convex feasibility problem.Some numerical experiments demonstrate the feasibility and e ectiveness of the proposed algorithm,showing superior performance compared to conventional alternating projection methods.展开更多
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
基金supported by the National Special Fund for Major Research Instrument Development(2011YQ140145)111 Project (B07009)+1 种基金the National Natural Science Foundation of China(11002013)Defense Industrial Technology Development Program(A2120110001 and B2120110011)
文摘A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.
基金support by the National Natural Science Foundation of China(Grant No.52402520)。
文摘Actuator faults can be critical in turbofan engines as they can lead to stall,surge,loss of thrust and failure of speed control.Thus,fault diagnosis of gas turbine actuators has attracted considerable attention,from both academia and industry.However,the extensive literature that exists on this topic does not address identifying the severity of actuator faults and focuses mainly on actuator fault detection and isolation.In addition,previous studies of actuator fault identification have not dealt with multiple concurrent faults in real time,especially when these are accompanied by sudden failures under dynamic conditions.This study develops component-level models for fault identification in four typical actuators used in high-bypass ratio turbofan engines under both dynamic and steady-state conditions and these are then integrated with the engine performance model developed by the authors.The research results reported here present a novel method of quantifying actuator faults using dynamic effect compensation.The maximum error for each actuator is less than0.06%and 0.07%,with average computational time of less than 0.0058 s and 0.0086 s for steady-state and transient cases,respectively.These results confirm that the proposed method can accurately and efficiently identify concurrent actuator fault for an engine operating under either transient or steady-state conditions,even in the case of a sudden malfunction.The research results emonstrate the potential benefit to emergency response capabilities by introducing this method of monitoring the health of aero engines.
基金funded by National Nature Science Foundation of China(Nos.51776010 and 91860205)supported by the Academic Excellence Foundation of BUAA for PhD Students,China。
文摘The alternative working modes and flexible working states are the outstanding features of an adaptive cycle engine, with a proper control schedule design being the only way to exploit the performance of such an engine. However, unreasonable design in the control schedule causes not only performance deterioration but also serious aerodynamic stability problems. Thus, in this work,a hybrid optimization method that automatically chooses the working modes and identifies the optimal and smooth control schedules is proposed, by combining the differential evolution algorithm and the Latin hypercube sampling method. The control schedule architecture does not only optimize the engine steady-state performance under different working modes but also solves the control-schedule discontinuity problem, especially during mode transition. The optimal control schedules are continuous and almost monotonic, and hence are strongly suitable for a control system, and are designed for two different working conditions, i.e., supersonic and subsonic throttling,which proves that the proposed hybrid method applies to various working conditions. The evaluation demonstrates that the proposed control method optimizes the engine performance, the surge margin of the compression components, and the range of the thrust during throttling.
基金Project supported by the National Natural Science Foundation of China(Nos.91330205and 11421101)the National Key Research and Development Program of China(No.2016YFB0200603)
文摘A Newton multigrid method is developed for one-dimensional （1D） and two- dimensional （2D） steady-state shallow water equations （SWEs） with topography and dry areas. The nonlinear system arising from the well-balanced finite volume discretization of the steady-state SWEs is solved by the Newton method as the outer iteration and a geometric multigrid method with the block symmetric Gauss-Seidel smoother as the inner iteration. The proposed Newton multigrid method makes use of the local residual to regularize the Jacobian matrix of the Newton iteration, and can handle the steady- state problem with wet/dry transition. Several numerical experiments are conducted to demonstrate the efficiency, robustness, and well-balanced property of the proposed method. The relation between the convergence behavior of the Newton multigrid method and the distribution of the eigenvalues of the iteration matrix is detailedly discussed.
基金supported by Natural Science Foundation of Shanxi Province(Grant No.20210302124388).
文摘In this paper,we considered the improved element-free Galerkin(IEFG)method for solving 2D anisotropic steadystate heat conduction problems.The improved moving least-squares(IMLS)approximation is used to establish the trial function,and the penalty method is applied to enforce the boundary conditions,thus the final discretized equations of the IEFG method for anisotropic steady-state heat conduction problems can be obtained by combining with the corresponding Galerkin weak form.The influences of node distribution,weight functions,scale parameters and penalty factors on the computational accuracy of the IEFG method are analyzed respectively,and these numerical solutions show that less computational resources are spent when using the IEFG method.
基金supported by the NSFC Grant No.11872210 and Grant No.MCMS-I-0120G01Chi-Wang Shu:Research is supported by the AFOSR Grant FA9550-20-1-0055 and the NSF Grant DMS-2010107Jianxian Qiu:Research is supported by the NSFC Grant No.12071392.
文摘In this paper, we design high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (multi-resolution WENO) limiters to compute compressible steady-state problems on triangular meshes. A troubled cell indicator extended from structured meshes to unstructured meshes is constructed to identify triangular cells in which the application of the limiting procedures is required. In such troubled cells, the multi-resolution WENO limiting methods are used to the hierarchical L^(2) projection polynomial sequence of the DG solution. Through using the RKDG methods with multi-resolution WENO limiters, the optimal high-order accuracy can be gradually reduced to first-order in the triangular troubled cells, so that the shock wave oscillations can be well suppressed. In steady-state simulations on triangular meshes, the numerical residual converges to near machine zero. The proposed spatial reconstruction methods enhance the robustness of classical DG methods on triangular meshes. The good results of these RKDG methods with multi-resolution WENO limiters are verified by a series of two-dimensional steady-state problems.
基金supported by the National Natural Science Foundation of China (Grant No.12072114)the National Key Research and Development Plan (Grant No.2020YFB1709401)the Guangdong Provincial Key Laboratory of Modern Civil Engineering Technology (2021B1212040003).
文摘In this paper,we consider solving the topology optimization for steady-state incompressibleNavier-Stokes problems via a new topology optimization method called parameterized level set method,which can maintain a relatively smooth level set function with a local optimality condition.The objective of topology optimization is tond an optimal conguration of theuid and solid materials that minimizes power dissipation under a prescribeduid volume fraction constraint.An articial friction force is added to the Navier-Stokes equations to apply the no-slip boundary condition.Although a great deal of work has been carried out for topology optimization ofuidow in recent years,there are few researches on the topology optimization ofuidow with physical body forces.To simulate theuidow in reality,the constant body force(e.g.,gravity)is considered in this paper.Several 2D numerical examples are presented to discuss the relationships between the proposed method with Reynolds number and initial design,and demonstrate the feasibility and superiority of the proposed method in dealing with unstructuredmesh problems.Three 3D numerical examples demonstrate the proposedmethod is feasible in three-dimensional.
文摘An analysis technique of steady state and stability for closed-loop PWM DC/DC switching converters is presented. Using this method, the closed-loop switching converter is transformed into an open-loop system. By means of the fact that in steady state, the two boundary values are equal in one switching period. The exponential matrix is evaluated by precise time-domain-integration method, and then the related curve between feedback duty cycle and the input one is obtained. Not only can the steady-state duty cycle be found from the curve, but also the stability and stable domain of the system. Compared with other methods, it features with simplicity and less calculation, and fit for numerical simulation and analysis for closed-loop switching converters. The simulation results of examples indicate the correctness of the presented method.
文摘Based on the indirect Trefftz approach, a wave number method （WNM） is proposed to deal with three-dimensional steady-state acoustic problems. In the WNM, the dynamic pressure response variable is approximated by a set of wave functions, which exactly satisfy the Helmholtz equation. The set of wave functions comprise the exact solutions of the homogeneous part of the governing equations and some particular solution functions. The unknown coefficients of the wave functions can be obtained by enforcing the pressure approximation to satisfy the boundary conditions. Compared with the boundary element method （BEM）, the WNM have a smaller system matrix, and is applicable to the radiation problems since the wave functions are independent of the domain size. A 3D acoustic cavity is exemplified to show the properties of the method. The results show that the wave number method is more efficient than the BEM, and it is fairly accurate.
基金Project(2010QZZD018) supported by Central South Universtiy,China,and Ghent University,Belgium
文摘NT build 443, or profile fitting method, is often used to measure the diffusion coefficient of chloride in concrete. However, this method is quite laborious and needs special equipment. Colorimetric method is a quick and simple method to measure the penetration depth of chloride by spraying 0.1 mol/L silver nitrate solution. The objective of this work is to study the possibility of the use of colorimetric method in the calculation of non-steady-state diffusion coefficient. Twelve concrete mixtures with different supplementary cementitious materials and water-to-cement ratios of 0.35, 0.48 and 0.6 were used for study. According to NT build 443, the concrete specimens were immersed in 165 g/L NaC1 （2.8 mol/L） solution for 42 d. Both water-soluble （convert to free chloride） chloride and acid-soluble chloride at different layers of specimens were measured. The results show that the mean value of free chloride concentration at the color change boundary Cd was 0.306 mol/L. The surface free chloride concentration cs was obtained by profile fitting method, which was 40% lower than the chloride concentration of exposure solution after an immersion period of 42 d. Chloride diffusion coefficients obtained by the colorimetric method was not well correlated with those obtained by profile fitting method.
基金supported by the NST Fellowship under the Ministry of Science and Technology,Government of the People’s Republic of Bangladesh(Session:2019–2020,merit number:334,serial number:714,physical science).
文摘This study presents a numerical analysis of the steady-state solution for transient magnetohydrodynamic(MHD)dissipative and radiative fluid flow,incorporating an inducedmagnetic field(IMF)and considering a relatively high concentration of foreign mass(accounting for Soret and Dufour effects)over a vertically oriented semi-infinite plate.The governing equations were normalized using boundary layer(BL)approximations.The resulting nonlinear system of partial differential equations(PDEs)was discretized and solved using an efficient explicit finite difference method(FDM).Numerical simulations were conducted using MATLAB R2015a,and the developed numerical code was verified through comparison with another code written in FORTRAN 6.6a.To ensure the reliability of the results,both mesh refinement and steady-state time validation tests were performed.Furthermore,a comparison with existing published studies was made to confirm the accuracy of the findings.The dimensionless equations revealed the impacts of several key parameters.The IMF initially intensifies near the plate before gradually diminishing as the magnetic parameter increases.For the range 0≤y≤1.8(where y is the horizontal direction),the IMF decreases with a rise in the magnetic Prandtl number;however,for 1.8≤y≤7(approximately),the magnetic field begins to increase.Beyond this,the profile of the magnetic field becomes somewhat irregular through the remaining part of the BL.
基金supported by Key Scientific Research Projects of Universities and Key Scientific and Technological Projects in Henan Province,which numbers are 21A440015,22A570007 and 212102310601,respectively.
文摘The isogeometric boundary element technique(IGABEM)is presented in this study for steady-state inhomogeneous heat conduction analysis.The physical unknowns in the boundary integral formulations of the governing equations are discretized using non-uniform rational B-spline(NURBS)basis functions,which are utilized to build the geometry of the structures.To speed up the assessment of NURBS basis functions,the Bezier extraction´approach is used.To solve the extra domain integrals,we use a radial integration approach.The numerical examples show the potential of IGABEM for dimension reduction and smooth integration of CAD and numerical analysis.
基金Project supported by the National Natural Science Foundation of China (No. 10876100)
文摘Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed. It leads to the high quality locM error bounds in the problem of the direct-solution steady-state dynamic analysis with a frequency-domain finite element, which involves the enrichments with plural variable basis functions. The solution of the steady-state dynamic procedure calculates the harmonic response directly in terms of the physical degrees of freedom in the model, which uses the mass, damping, and stiffness matrices of the system. A three-dimensional finite element example is carried out to illustrate the computational procedures.
文摘A mathematical model of an amperometric biosensor with the substrate inhibition for steady-state condition is discussed. The model is based on the system of non-stationary diffusion equation containing a non-linear term related to non-Michaelis–Menten kinetics of the enzymatic reaction. This paper presents the analytical expression of concentrations and current for all values of parameters φ2s φ2s α and β . Here the Adomian decomposition method (ADM) is used to find the analytical expressions for substrate, product concentration and current. A comparison of the analytical approximation and numerical simulation is also presented. A good agreement between theoretical predictions and numerical results is observed.
文摘A mathematical modelling of diagnostic biosensors system at three basic types of enzyme kinetics is discussed in the presence of diffusion. Enzyme kinetics is adopted to be first order, Michaelis-Menten and ping-pong mechanism. In this paper, approximate analytical solutions are obtained for the non-linear equations under steady-state conditions by using the new Homotopy perturbation method. Simple and closed forms of analytical expressions for concentrations of substrate, product and co-substrate and corresponding current response have been derived for all possible values of parameters. Furthermore, the numerical simulation of the problem is also reported here by using Matlab program. Good agreement between analytical and numerical results is noted.
基金Supported in part by Natural Science Foundation of Guangxi(2023GXNSFAA026246)in part by the Central Government's Guide to Local Science and Technology Development Fund(GuikeZY23055044)in part by the National Natural Science Foundation of China(62363003)。
文摘In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.
基金This work was supported in part by National Natural key R&D Program of China(2016YFB0900100).
文摘The low frequency oscillation is a serious threat to security and stability of a power grid.How to locate the disturbance source accurately is an important issue to low frequency oscillation disposal.Existing methods have poor adaptability to the low frequency oscillation with time-varying steady-state points because of the limitations in the location criterion derivation.A disturbance source location method on a low frequency oscillation with good generality is presented in the paper.Firstly,the reasons why the steady-state points are time-varying on a low frequency oscillation are analyzed.Then,based on the energy function construction form,the branch transmission energy is decomposed into state energy,reciprocating energy and dissipation energy by mathematical derivation.The flow direction of the dissipation energy shows the source and destination of the disturbance energy,and the specific location of a disturbance source can be identified according to its flow direction.Meanwhile,to meet the needs of energy calculation,a recognition method on the electrical quantities steady-state points is also presented by using the cubic spline interpolation.Simulation results show the correctness of the derivation and analysis on energy structure in the paper,and the disturbance source can be located accurately according to the dissipation energy.
基金National Natural Science Foundation of China(No.61273172)
文摘This paper proposes a steady-state errors correction(SSEC)method for eliminating measurement errors.This method is based on the detections of error signal E(s)and output C(s)which generate an expected output R(s).In comparison with the conventional solutions which are based on detecting the expected output R(s)and output C(s)to obtain error signal E(s),the measurement errors are eliminated even the error might be at a significant level.Moreover,it is possible that the individual debugging by regulating the coefficient K for every member of the multiple objectives achieves the optimization of the open loop gain.Therefore,this simple method can be applied to the weak coupling and multiple objectives system,which is usually controlled by complex controller.The principle of eliminating measurement errors is derived analytically,and the advantages comparing with the conventional solutions are depicted.Based on the SSEC method analysis,an application of this method for an active power filter(APF)is investigated and the effectiveness and viability of the scheme are demonstrated through the simulation and experimental verifications.
基金Supported by the Natural Science Foundation of Guangxi Province(Grant Nos.2023GXNSFAA026067,2024GXN SFAA010521)the National Natural Science Foundation of China(Nos.12361079,12201149,12261026).
文摘Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery.We rst derive the projection formulas for a vector onto the feasible sets.The centralized circumcentered-reection method is designed to solve the convex feasibility problem.Some numerical experiments demonstrate the feasibility and e ectiveness of the proposed algorithm,showing superior performance compared to conventional alternating projection methods.
