Energy storage systems(ESSs)operate as independent market participants and collaborate with photovoltaic(PV)generation units to enhance the flexible power supply capabilities of PV units.However,the dynamic variations...Energy storage systems(ESSs)operate as independent market participants and collaborate with photovoltaic(PV)generation units to enhance the flexible power supply capabilities of PV units.However,the dynamic variations in the profitability of ESSs in the electricity market are yet to be fully understood.This study introduces a dual-timescale dynamics model that integrates a spot market clearing(SMC)model into a system dynamics(SD)model to investigate the profit-aware capacity growth of ESSs and compares the profitability of independent energy storage systems(IESSs)with that of an ESS integrated within a PV(PV-ESS).Furthermore,this study aims to ascertain the optimal allocation of the PV-ESS.First,SD and SMC models were set up.Second,the SMC model simulated on an hourly timescale was incorporated into the SD model as a subsystem,a dual-timescale model was constructed.Finally,a development simulation and profitability analysis was conducted from 2022 to 2040 to reveal the dynamic optimal range of PV-ESS allocation.Additionally,negative electricity prices were considered during clearing processes.The simulation results revealed differences in profitability and capacity growth between IESS and PV-ESS,helping grid investors and policymakers to determine the boundaries of ESSs and dynamic optimal allocation of PV-ESSs.展开更多
An accurate and novel small-signal equivalent circuit model for GaN high-electron-mobility transistors(HEMTs)is proposed,which considers a dual-field-plate(FP)made up of a gate-FP and a source-FP.The equivalent circui...An accurate and novel small-signal equivalent circuit model for GaN high-electron-mobility transistors(HEMTs)is proposed,which considers a dual-field-plate(FP)made up of a gate-FP and a source-FP.The equivalent circuit of the overall model is composed of parasitic elements,intrinsic transistors,gate-FP,and source-FP networks.The equivalent circuit of the gate-FP is identical to that of the intrinsic transistor.In order to simplify the complexity of the model,a series combination of a resistor and a capacitor is employed to represent the source-FP.The analytical extraction procedure of the model parameters is presented based on the proposed equivalent circuit.The verification is carried out on a 4×250μm GaN HEMT device with a gate-FP and a source-FP in a 0.45μm technology.Compared with the classic model,the proposed novel small-signal model shows closer agreement with measured S-parameters in the range of 1.0 to 18.0 GHz.展开更多
We estimate tree heights using polarimetric interferometric synthetic aperture radar(PolInSAR)data constructed by the dual-polarization(dual-pol)SAR data and random volume over the ground(RVoG)model.Considering the Se...We estimate tree heights using polarimetric interferometric synthetic aperture radar(PolInSAR)data constructed by the dual-polarization(dual-pol)SAR data and random volume over the ground(RVoG)model.Considering the Sentinel-1 SAR dual-pol(SVV,vertically transmitted and vertically received and SVH,vertically transmitted and horizontally received)configuration,one notes that S_(HH),the horizontally transmitted and horizontally received scattering element,is unavailable.The S_(HH)data were constructed using the SVH data,and polarimetric SAR(PolSAR)data were obtained.The proposed approach was first verified in simulation with satisfactory results.It was next applied to construct PolInSAR data by a pair of dual-pol Sentinel-1A data at Duke Forest,North Carolina,USA.According to local observations and forest descriptions,the range of estimated tree heights was overall reasonable.Comparing the heights with the ICESat-2 tree heights at 23 sampling locations,relative errors of 5 points were within±30%.Errors of 8 points ranged from 30%to 40%,but errors of the remaining 10 points were>40%.The results should be encouraged as error reduction is possible.For instance,the construction of PolSAR data should not be limited to using SVH,and a combination of SVH and SVV should be explored.Also,an ensemble of tree heights derived from multiple PolInSAR data can be considered since tree heights do not vary much with time frame in months or one season.展开更多
D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration a...D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration are used to solve seismic wave equation, which means difference scheme for space domain and analytic integration for time domain. Both the principle and algorithm of this method are introduced in the paper. Based on the theory, the numerical examples prove that this hybrid method can lead to higher accuracy than the traditional finite difference method and the solution is very close to the exact one. Also the seismic modeling examples show the good performance of this method even in the case of complex surface conditions and complicated structures.展开更多
Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transfo...Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transformation to curvilinear coordinates using body-fitted grids have been proposed, e.g., stand staggered grid(SSG) with interpolation, nonstaggered grid, rotated staggered grid(RSG), and fully staggered. The FD based on the RSG is somewhat superior to others because it satisfies the spatial distribution of the wave equation without additional memory and computational requirements; furthermore, it is simpler to implement. We use the RSG FD method to transform the firstorder stress–velocity equation in the curvilinear coordinates system and introduce the highprecision adaptive, unilateral mimetic finite-difference(UMFD) method to process the freeboundary conditions of an irregular surface. The numerical results suggest that the precision of the solution is higher than that of the vacuum formalism. When the minimum wavelength is low, UMFD avoids the surface wave dispersion. We compare FD methods based on RSG, SEM, and nonstaggered grid and infer that all simulation results are consistent but the computational efficiency of the RSG FD method is higher than the rest.展开更多
In the dual risk model, we consider the optimal dividend and capital injection problem, which involves a random time horizon and a ruin penalty. Both fixed and proportional costs from the transactions of capital injec...In the dual risk model, we consider the optimal dividend and capital injection problem, which involves a random time horizon and a ruin penalty. Both fixed and proportional costs from the transactions of capital injection are considered. The objective is to maximize the total value of the expected discounted dividends, and the penalized discounted both capital injections and ruin penalty during the horizon, which is described by the minimum of the time of ruin and an exponential random variable. The explicit solutions for optimal strategy and value function are obtained, when the income jumps follow a hyper-exponential distribution.Besides, some numerical examples are presented to illustrate our results.展开更多
Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with ...Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.展开更多
In this paper, a new lattice hydrodynamic model based on Nagatani's model INagatani T 1998 Physica A 261 5991 is presented by introducing the flow difference effect. The stability condition for the new model is obtai...In this paper, a new lattice hydrodynamic model based on Nagatani's model INagatani T 1998 Physica A 261 5991 is presented by introducing the flow difference effect. The stability condition for the new model is obtained by using the linear stability theory. The result shows that considering the flow difference effect leads to stabilization of the system compared with the original lattice hydrodynamic model. The jamming transitions among the freely moving phase, the coexisting phase, and the uniform congested phase are studied by nonlinear analysis. The modified KdV equation near the critical point is derived to describe the traffic jam, and kink -antikink soliton solutions related to the traffic density waves are obtained. The simulation results are consistent with the theoretical analysis for the new model.展开更多
As one of the basic inventory cost models, the (Q, τ)inventory cost model of dual suppliers with random procurement lead time is mostly formulated by using the concepts of "effective lead time" and "lead time de...As one of the basic inventory cost models, the (Q, τ)inventory cost model of dual suppliers with random procurement lead time is mostly formulated by using the concepts of "effective lead time" and "lead time demand", which may lead to an imprecise inventory cost. Through the real-time statistic of the inventory quantities, this paper considers the precise (Q, τ) inventory cost model of dual supplier procurement by using an infinitesimal dividing method. The traditional modeling method of the inventory cost for dual supplier procurement includes complex procedures. To reduce the complexity effectively, the presented method investigates the statistics properties in real-time of the inventory quantities with the application of the infinitesimal dividing method. It is proved that the optimal holding and shortage costs of dual supplier procurement are less than those of single supplier procurement respectively. With the assumption that both suppliers have the same distribution of lead times, the convexity of the cost function per unit time is proved. So the optimal solution can be easily obtained by applying the classical convex optimization methods. The numerical examples are given to verify the main conclusions.展开更多
A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordin...A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.展开更多
The influence of initial groove angle on strain rate inside and outside groove of Ti6Al4V alloy was investigated.Based on the evolution of strain rate inside and outside groove,the effect of strain rate difference on ...The influence of initial groove angle on strain rate inside and outside groove of Ti6Al4V alloy was investigated.Based on the evolution of strain rate inside and outside groove,the effect of strain rate difference on the evolution of normal stress and effective stress inside and outside groove was also analyzed.The results show that when linear loading path changes from uniaxial tension to equi-biaxial tension,the initial groove angle plays a weaker role in the evolution of strain rate in the M-K model.Due to the constraint of force equilibrium between inside and outside groove,the strain rate difference makes the normal stress inside groove firstly decrease and then increase during calculation,which makes the prediction algorithm of forming limit convergent at elevated temperature.The decrease of normal stress inside groove is mainly caused by high temperature softening effect and the rotation of groove,while the increase of normal stress inside groove is mainly due to strain rate hardening effect.展开更多
To realize high accurate control of relative position and attitude between two spacecrafts, the coupling between position and attitude must be fully considered and a more precise model should be established. This pape...To realize high accurate control of relative position and attitude between two spacecrafts, the coupling between position and attitude must be fully considered and a more precise model should be established. This paper breaks the traditional divide and conquer idea, and uses a mathematical tool, namely dual quaternion to establish the integrated 6 degree-of-freedom(6-DOF) model of relative position and attitude, which describes the coupled relative motion in a compact and efficient form and needs less information of the target. Considering the complex operation rules and the unclarity of the current relative motion model in dual quaternion, necessary mathematical foundations are given at first, followed by clear and detailed modeling process and analysis. Finally a generalized proportion-derivative(PD) controller law is designed. The simulation results show that based on the integrated model established by dual quaternion, this control law can achieve a high control accuracy of relative motion.展开更多
In this paper, we consider the dual risk model in which periodic taxation are paid according to a loss-carry-forward system and dividends are paid under a threshold strategy. We give an analytical approach to derive t...In this paper, we consider the dual risk model in which periodic taxation are paid according to a loss-carry-forward system and dividends are paid under a threshold strategy. We give an analytical approach to derive the expression of gδ(u) (i.e. the Laplace transform of the first upper exit time). We discuss the expected discounted tax payments for this model and obtain its corresponding integro-differential equations. Finally, for Erlang (2) inter-innovation distribution, closedform expressions for the expected discounted tax payments are given.展开更多
In this paper, we consider the dual of the generalized Erlang (n) risk model under a threshold dividend strategy. We derive an integrodifferential equation satisfied by the expectation of the discounted dividends unti...In this paper, we consider the dual of the generalized Erlang (n) risk model under a threshold dividend strategy. We derive an integrodifferential equation satisfied by the expectation of the discounted dividends until ruin. The case when profits follow an exponential distribution is solved.展开更多
Simulation of fluid flow in the fractured porous media is very important and challenging.Researchers have developed some models for fractured porous media.With the development of related research in recent years,the p...Simulation of fluid flow in the fractured porous media is very important and challenging.Researchers have developed some models for fractured porous media.With the development of related research in recent years,the prospect of embedded discrete fracture model(EDFM)is more and more bright.However,since the size of the fractures in the actual reservoir varies greatly,a very fine grid should be used which leads to a huge burden to the computing resources.To address this challenge,in the present paper,an upscaling based model is proposed.In this model,the flow in large-scale fractures is directly described by the EDFM while that in the small-scale fractures is upscaled through local simulation by EDFM.The EDFM is used to simulate the large-and small-scale fractures independently two times,so the new model is called dual embedded discrete fracture model(D-EDFM).In this paper,the detailed implementation process of D-EDFM is introduced and,through test cases,it is found the proposed model is a feasible method to simulate the flow in fractured porous media.展开更多
Surface subsidence induced by underground mining is a typical serious geohazard.Numerical approaches such as the discrete element method(DEM)and finite difference method(FDM)have been widely used to model and analyze ...Surface subsidence induced by underground mining is a typical serious geohazard.Numerical approaches such as the discrete element method(DEM)and finite difference method(FDM)have been widely used to model and analyze mining-induced surface subsidence.However,the DEM is typically computationally expensive,and is not capable of analyzing large-scale problems,while the mesh distortion may occur in the FDM modeling of largely deformed surface subsidence.To address the above problems,this paper presents a geometrically and locally adaptive remeshing method for the FDM modeling of largely deformed surface subsidence induced by underground mining.The essential ideas behind the proposed method are as follows:(i)Geometrical features of elements(i.e.the mesh quality),rather than the calculation errors,are employed as the indicator for determining whether to conduct the remeshing;and(ii)Distorted meshes with multiple attributes,rather than those with only a single attribute,are locally regenerated.In the proposed method,the distorted meshes are first adaptively determined based on the mesh quality,and then removed from the original mesh model.The tetrahedral mesh in the distorted area is first regenerated,and then the physical field variables of old mesh are transferred to the new mesh.The numerical calculation process recovers when finishing the regeneration and transformation.To verify the effectiveness of the proposed method,the surface deformation of the Yanqianshan iron mine,Liaoning Province,China,is numerically investigated by utilizing the proposed method,and compared with the numerical results of the DEM modeling.Moreover,the proposed method is applied to predicting the surface subsidence in Anjialing No.1 Underground Mine,Shanxi Province,China.展开更多
In this work,we developed a simple and direct circuit model with a dual two-diode model that can be solved by a SPICE numerical simulation to comprehensively describe the monolithic perovskite/crystalline silicon(PVS/...In this work,we developed a simple and direct circuit model with a dual two-diode model that can be solved by a SPICE numerical simulation to comprehensively describe the monolithic perovskite/crystalline silicon(PVS/c-Si)tandem solar cells.We are able to reveal the effects of different efficiency-loss mechanisms based on the illuminated current density-voltage(J-V),semi-log dark J-V,and local ideality factor(m-V)curves.The effects of the individual efficiency-loss mechanism on the tandem cell’s efficiency are discussed,including the exp(V/VT)and exp(V/2VT)recombination,the whole cell’s and subcell’s shunts,and the Ohmic-contact or Schottky-contact of the intermediate junction.We can also fit a practical J-V curve and find a specific group of parameters by the trial-and-error method.Although the fitted parameters are not a unique solution,they are valuable clues for identifying the efficiency loss with the aid of the cell’s structure and experimental processes.This method can also serve as an open platform for analyzing other tandem solar cells by substituting the corresponding circuit models.In summary,we developed a simple and effective methodology to diagnose the efficiency-loss source of a monolithic PVS/c-Si tandem cell,which is helpful to researchers who wish to adopt the proper approaches to improve their solar cells.展开更多
基金supported by National Natural Science Foundation of China(U2066209)。
文摘Energy storage systems(ESSs)operate as independent market participants and collaborate with photovoltaic(PV)generation units to enhance the flexible power supply capabilities of PV units.However,the dynamic variations in the profitability of ESSs in the electricity market are yet to be fully understood.This study introduces a dual-timescale dynamics model that integrates a spot market clearing(SMC)model into a system dynamics(SD)model to investigate the profit-aware capacity growth of ESSs and compares the profitability of independent energy storage systems(IESSs)with that of an ESS integrated within a PV(PV-ESS).Furthermore,this study aims to ascertain the optimal allocation of the PV-ESS.First,SD and SMC models were set up.Second,the SMC model simulated on an hourly timescale was incorporated into the SD model as a subsystem,a dual-timescale model was constructed.Finally,a development simulation and profitability analysis was conducted from 2022 to 2040 to reveal the dynamic optimal range of PV-ESS allocation.Additionally,negative electricity prices were considered during clearing processes.The simulation results revealed differences in profitability and capacity growth between IESS and PV-ESS,helping grid investors and policymakers to determine the boundaries of ESSs and dynamic optimal allocation of PV-ESSs.
文摘An accurate and novel small-signal equivalent circuit model for GaN high-electron-mobility transistors(HEMTs)is proposed,which considers a dual-field-plate(FP)made up of a gate-FP and a source-FP.The equivalent circuit of the overall model is composed of parasitic elements,intrinsic transistors,gate-FP,and source-FP networks.The equivalent circuit of the gate-FP is identical to that of the intrinsic transistor.In order to simplify the complexity of the model,a series combination of a resistor and a capacitor is employed to represent the source-FP.The analytical extraction procedure of the model parameters is presented based on the proposed equivalent circuit.The verification is carried out on a 4×250μm GaN HEMT device with a gate-FP and a source-FP in a 0.45μm technology.Compared with the classic model,the proposed novel small-signal model shows closer agreement with measured S-parameters in the range of 1.0 to 18.0 GHz.
文摘We estimate tree heights using polarimetric interferometric synthetic aperture radar(PolInSAR)data constructed by the dual-polarization(dual-pol)SAR data and random volume over the ground(RVoG)model.Considering the Sentinel-1 SAR dual-pol(SVV,vertically transmitted and vertically received and SVH,vertically transmitted and horizontally received)configuration,one notes that S_(HH),the horizontally transmitted and horizontally received scattering element,is unavailable.The S_(HH)data were constructed using the SVH data,and polarimetric SAR(PolSAR)data were obtained.The proposed approach was first verified in simulation with satisfactory results.It was next applied to construct PolInSAR data by a pair of dual-pol Sentinel-1A data at Duke Forest,North Carolina,USA.According to local observations and forest descriptions,the range of estimated tree heights was overall reasonable.Comparing the heights with the ICESat-2 tree heights at 23 sampling locations,relative errors of 5 points were within±30%.Errors of 8 points ranged from 30%to 40%,but errors of the remaining 10 points were>40%.The results should be encouraged as error reduction is possible.For instance,the construction of PolSAR data should not be limited to using SVH,and a combination of SVH and SVV should be explored.Also,an ensemble of tree heights derived from multiple PolInSAR data can be considered since tree heights do not vary much with time frame in months or one season.
基金This project is sponsored by the Specialized Prophasic Basic Research of the"973"Programme,contract No:2001cca02300
文摘D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration are used to solve seismic wave equation, which means difference scheme for space domain and analytic integration for time domain. Both the principle and algorithm of this method are introduced in the paper. Based on the theory, the numerical examples prove that this hybrid method can lead to higher accuracy than the traditional finite difference method and the solution is very close to the exact one. Also the seismic modeling examples show the good performance of this method even in the case of complex surface conditions and complicated structures.
基金supported by the National Nature Science Foundation of China(Nos.41504102 and 41604037)National Science and Technology Major Project(No.2016ZX05015-006)Yangtze University Youth Found(No.2015cqn32)
文摘Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transformation to curvilinear coordinates using body-fitted grids have been proposed, e.g., stand staggered grid(SSG) with interpolation, nonstaggered grid, rotated staggered grid(RSG), and fully staggered. The FD based on the RSG is somewhat superior to others because it satisfies the spatial distribution of the wave equation without additional memory and computational requirements; furthermore, it is simpler to implement. We use the RSG FD method to transform the firstorder stress–velocity equation in the curvilinear coordinates system and introduce the highprecision adaptive, unilateral mimetic finite-difference(UMFD) method to process the freeboundary conditions of an irregular surface. The numerical results suggest that the precision of the solution is higher than that of the vacuum formalism. When the minimum wavelength is low, UMFD avoids the surface wave dispersion. We compare FD methods based on RSG, SEM, and nonstaggered grid and infer that all simulation results are consistent but the computational efficiency of the RSG FD method is higher than the rest.
基金Supported by the National Natural Science Foundation of China(11231005)Promotive research fund for excellent young and middle-aged scientists of Shandong Province(BS2014SF006)the Natural Science Foundation of the Jiangsu Higher Education Institutions of China(15KJB110009)
文摘In the dual risk model, we consider the optimal dividend and capital injection problem, which involves a random time horizon and a ruin penalty. Both fixed and proportional costs from the transactions of capital injection are considered. The objective is to maximize the total value of the expected discounted dividends, and the penalized discounted both capital injections and ruin penalty during the horizon, which is described by the minimum of the time of ruin and an exponential random variable. The explicit solutions for optimal strategy and value function are obtained, when the income jumps follow a hyper-exponential distribution.Besides, some numerical examples are presented to illustrate our results.
基金supported by China Scholarship Council and partially by the National "863" Program of China under contract No. 2007AA06Z218.
文摘Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.
基金Project supported by the National Basic Research Program of China (Grant No. G2006CB705500)the National Natural Science Foundation of China (Grant Nos. 70501004,70701004 and 70631001)Program for New Century Excellent Talents in University(Grant No. NCET-07-0057)
文摘In this paper, a new lattice hydrodynamic model based on Nagatani's model INagatani T 1998 Physica A 261 5991 is presented by introducing the flow difference effect. The stability condition for the new model is obtained by using the linear stability theory. The result shows that considering the flow difference effect leads to stabilization of the system compared with the original lattice hydrodynamic model. The jamming transitions among the freely moving phase, the coexisting phase, and the uniform congested phase are studied by nonlinear analysis. The modified KdV equation near the critical point is derived to describe the traffic jam, and kink -antikink soliton solutions related to the traffic density waves are obtained. The simulation results are consistent with the theoretical analysis for the new model.
基金supported by the National High Technology Research and Development Program of China(863 Program)(2007AA04Z102)the National Natural Science Foundation of China(6087407160574077).
文摘As one of the basic inventory cost models, the (Q, τ)inventory cost model of dual suppliers with random procurement lead time is mostly formulated by using the concepts of "effective lead time" and "lead time demand", which may lead to an imprecise inventory cost. Through the real-time statistic of the inventory quantities, this paper considers the precise (Q, τ) inventory cost model of dual supplier procurement by using an infinitesimal dividing method. The traditional modeling method of the inventory cost for dual supplier procurement includes complex procedures. To reduce the complexity effectively, the presented method investigates the statistics properties in real-time of the inventory quantities with the application of the infinitesimal dividing method. It is proved that the optimal holding and shortage costs of dual supplier procurement are less than those of single supplier procurement respectively. With the assumption that both suppliers have the same distribution of lead times, the convexity of the cost function per unit time is proved. So the optimal solution can be easily obtained by applying the classical convex optimization methods. The numerical examples are given to verify the main conclusions.
基金The National Natural Science Foundation of China (No10471023)
文摘A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.
基金Project(51775023)supported by the National Natural Science Foundation of ChinaProject(YWF-18-BJ-J-75)supported by the Fundamental Research Funds for the Central Universities,China
文摘The influence of initial groove angle on strain rate inside and outside groove of Ti6Al4V alloy was investigated.Based on the evolution of strain rate inside and outside groove,the effect of strain rate difference on the evolution of normal stress and effective stress inside and outside groove was also analyzed.The results show that when linear loading path changes from uniaxial tension to equi-biaxial tension,the initial groove angle plays a weaker role in the evolution of strain rate in the M-K model.Due to the constraint of force equilibrium between inside and outside groove,the strain rate difference makes the normal stress inside groove firstly decrease and then increase during calculation,which makes the prediction algorithm of forming limit convergent at elevated temperature.The decrease of normal stress inside groove is mainly caused by high temperature softening effect and the rotation of groove,while the increase of normal stress inside groove is mainly due to strain rate hardening effect.
基金supported by the National Natural Science Foundation of China(6107412761427809)
文摘To realize high accurate control of relative position and attitude between two spacecrafts, the coupling between position and attitude must be fully considered and a more precise model should be established. This paper breaks the traditional divide and conquer idea, and uses a mathematical tool, namely dual quaternion to establish the integrated 6 degree-of-freedom(6-DOF) model of relative position and attitude, which describes the coupled relative motion in a compact and efficient form and needs less information of the target. Considering the complex operation rules and the unclarity of the current relative motion model in dual quaternion, necessary mathematical foundations are given at first, followed by clear and detailed modeling process and analysis. Finally a generalized proportion-derivative(PD) controller law is designed. The simulation results show that based on the integrated model established by dual quaternion, this control law can achieve a high control accuracy of relative motion.
文摘In this paper, we consider the dual risk model in which periodic taxation are paid according to a loss-carry-forward system and dividends are paid under a threshold strategy. We give an analytical approach to derive the expression of gδ(u) (i.e. the Laplace transform of the first upper exit time). We discuss the expected discounted tax payments for this model and obtain its corresponding integro-differential equations. Finally, for Erlang (2) inter-innovation distribution, closedform expressions for the expected discounted tax payments are given.
文摘In this paper, we consider the dual of the generalized Erlang (n) risk model under a threshold dividend strategy. We derive an integrodifferential equation satisfied by the expectation of the discounted dividends until ruin. The case when profits follow an exponential distribution is solved.
基金supported by National Natural Science Foundations of China(Grant Nos.51706021,51936001 and 51804033)the Beijing Youth Talent Support Program(Grant No.CIT&TCD201804037)+1 种基金Joint Project of the Beijing Natural Science Foundation and the Beijing Municipal Education Commission(Grant No.KZ201810017023)the Great Wall Scholar program(Grant No.CIT&TCD20180313).
文摘Simulation of fluid flow in the fractured porous media is very important and challenging.Researchers have developed some models for fractured porous media.With the development of related research in recent years,the prospect of embedded discrete fracture model(EDFM)is more and more bright.However,since the size of the fractures in the actual reservoir varies greatly,a very fine grid should be used which leads to a huge burden to the computing resources.To address this challenge,in the present paper,an upscaling based model is proposed.In this model,the flow in large-scale fractures is directly described by the EDFM while that in the small-scale fractures is upscaled through local simulation by EDFM.The EDFM is used to simulate the large-and small-scale fractures independently two times,so the new model is called dual embedded discrete fracture model(D-EDFM).In this paper,the detailed implementation process of D-EDFM is introduced and,through test cases,it is found the proposed model is a feasible method to simulate the flow in fractured porous media.
基金supported by the National Natural Science Foundation of China(Grant Nos.11602235 and 41772326)the Fundamental Research Funds for the Central Universities of China(Grant No.2652018091)。
文摘Surface subsidence induced by underground mining is a typical serious geohazard.Numerical approaches such as the discrete element method(DEM)and finite difference method(FDM)have been widely used to model and analyze mining-induced surface subsidence.However,the DEM is typically computationally expensive,and is not capable of analyzing large-scale problems,while the mesh distortion may occur in the FDM modeling of largely deformed surface subsidence.To address the above problems,this paper presents a geometrically and locally adaptive remeshing method for the FDM modeling of largely deformed surface subsidence induced by underground mining.The essential ideas behind the proposed method are as follows:(i)Geometrical features of elements(i.e.the mesh quality),rather than the calculation errors,are employed as the indicator for determining whether to conduct the remeshing;and(ii)Distorted meshes with multiple attributes,rather than those with only a single attribute,are locally regenerated.In the proposed method,the distorted meshes are first adaptively determined based on the mesh quality,and then removed from the original mesh model.The tetrahedral mesh in the distorted area is first regenerated,and then the physical field variables of old mesh are transferred to the new mesh.The numerical calculation process recovers when finishing the regeneration and transformation.To verify the effectiveness of the proposed method,the surface deformation of the Yanqianshan iron mine,Liaoning Province,China,is numerically investigated by utilizing the proposed method,and compared with the numerical results of the DEM modeling.Moreover,the proposed method is applied to predicting the surface subsidence in Anjialing No.1 Underground Mine,Shanxi Province,China.
基金This work was supported by Zhejiang Energy Group(znkj-2018-118)Key Research and Development Program of Zhejiang Province(2021C01006)+5 种基金Key Project of Zhejiang Province(2021C04009)Science and technology projects in Liaoning Province 2021(2021JH1/10400104)Ningbo“Innovation 2025”Major Project(2020Z098)National Key R&D Program of China(2018YFB1500403)National Natural Science Foundation of China(61974178,61874177,62004199)Youth Innovation Promotion Association(2018333).
文摘In this work,we developed a simple and direct circuit model with a dual two-diode model that can be solved by a SPICE numerical simulation to comprehensively describe the monolithic perovskite/crystalline silicon(PVS/c-Si)tandem solar cells.We are able to reveal the effects of different efficiency-loss mechanisms based on the illuminated current density-voltage(J-V),semi-log dark J-V,and local ideality factor(m-V)curves.The effects of the individual efficiency-loss mechanism on the tandem cell’s efficiency are discussed,including the exp(V/VT)and exp(V/2VT)recombination,the whole cell’s and subcell’s shunts,and the Ohmic-contact or Schottky-contact of the intermediate junction.We can also fit a practical J-V curve and find a specific group of parameters by the trial-and-error method.Although the fitted parameters are not a unique solution,they are valuable clues for identifying the efficiency loss with the aid of the cell’s structure and experimental processes.This method can also serve as an open platform for analyzing other tandem solar cells by substituting the corresponding circuit models.In summary,we developed a simple and effective methodology to diagnose the efficiency-loss source of a monolithic PVS/c-Si tandem cell,which is helpful to researchers who wish to adopt the proper approaches to improve their solar cells.