To deal with the staircase approximation problem in the standard finite-difference time-domain(FDTD) simulation,the two-dimensional boundary condition equations(BCE) method is proposed in this paper.In the BCE met...To deal with the staircase approximation problem in the standard finite-difference time-domain(FDTD) simulation,the two-dimensional boundary condition equations(BCE) method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.展开更多
In this paper we give a priori estimates for the maximum modulus of generalizedsolulions of the quasilinear elliplic equations irith anisotropic growth condition.
This paper deals with the solution of a parametric equation with generalized boundary condiiton in transport theory. It gives the distribution of parameter (so called delta-eigenvalue [1]) with which the equation has ...This paper deals with the solution of a parametric equation with generalized boundary condiiton in transport theory. It gives the distribution of parameter (so called delta-eigenvalue [1]) with which the equation has non-zero solution. A necessary and sufficient condition for the existence of; he control critical eigenvalue delta0 is established.展开更多
A Fourier pseudospectral-finue difference scheme is proposed for three-dimensional vorticity. equalion with unilalerally periodic boundary. condilion. Thegeneralized stability. and conrergence are analyzed. The nunier...A Fourier pseudospectral-finue difference scheme is proposed for three-dimensional vorticity. equalion with unilalerally periodic boundary. condilion. Thegeneralized stability. and conrergence are analyzed. The nunierical results show. theadvantage of this method.展开更多
In this paper, some misunderstam igs concerning the necessary conditions for resonance for ordinary differential equations with turning point have been corrected, and a recursive process for finding the sequence of ne...In this paper, some misunderstam igs concerning the necessary conditions for resonance for ordinary differential equations with turning point have been corrected, and a recursive process for finding the sequence of necessary conditions for resonance has beenoffered.展开更多
Accurate boundary conditions of composite material plates with different holes are founded to settle boundary condition problems of complex holes by conformal mapping method upon the nonhomogeneous anisotropic elastic...Accurate boundary conditions of composite material plates with different holes are founded to settle boundary condition problems of complex holes by conformal mapping method upon the nonhomogeneous anisotropic elastic and complex function theory. And then the two stress functions required were founded on Cauchy integral by boundary conditions. The final stress distributions of opening structure and the analytical solution on composite material plate with rectangle hole and wing manholes were achieved. The influences on hole-edge stress concentration factors are discussed under different loads and fiber direction cases, and then contrast calculates are carried through FEM.展开更多
There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle th...There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle these kind of mixed boundary-value problems associated with the Laplace’s equation (or Helmholtz equation) arising in the study of waves propagating through solids or fluids. One of the widely used methods in wave structure interaction is the multipole expansion method. This expansion involves a general combination of a regular wave, a wave source, a wave dipole and a regular wave-free part. The wave-free part can be further expanded in terms of wave-free multipoles which are termed as wave-free potentials. These are singular solutions of Laplace’s equation or two-dimensional Helmholz equation. Construction of these wave-free potentials and multipoles are presented here in a systematic manner for a number of situations such as two-dimensional non-oblique and oblique waves, three dimensional waves in two-layer fluid with free surface condition with higher order partial derivative are considered. In particular, these are obtained taking into account of the effect of the presence of surface tension at the free surface and also in the presence of an ice-cover modelled as a thin elastic plate. Also for limiting case, it can be shown that the multipoles and wave-free potential functions go over to the single layer multipoles and wave-free potential.展开更多
With the help of skew-symmetric differential forms, the hidden properties of the mathematical physics equations that describe discrete quantum transitions and emergence the physical structures are investigated. It is ...With the help of skew-symmetric differential forms, the hidden properties of the mathematical physics equations that describe discrete quantum transitions and emergence the physical structures are investigated. It is shown that the mathematical physics equations possess a unique property. They can describe discrete quantum transitions, emergence of physical structures and occurrence observed formations. However, such a property possesses only equations on which no additional conditions, namely, the conditions of integrability, are imposed. The intergrability conditions are realized from the equations themselves. Just under realization of integrability conditions double solutions to the mathematical physics equations, which describe discrete transitions and so on, are obtained. The peculiarity consists in the fact that the integrability conditions do not directly follow from the mathematical physics equations;they are realized under the description of evolutionary process. The hidden properties of differential equations were discovered when studying the integrability of differential equations of mathematical physics that depends on the consistence between the derivatives in differential equations along different directions and on the consistence of equations in the set of equations. The results of this work were obtained with the help of skew-symmetric differential forms that possess a nontraditional mathematical apparatus such as nonidentical relations, degenerate transformations and the transition from nonintegrable manifolds to integrable structures. Such results show that mathematical physics equations can describe quantum processes.展开更多
The global boundness and existence are presented for the kind of the Rosseland equation with a general growth condition. A linearized map in a closed convex set is defined. The image set is precompact, and thus a fixe...The global boundness and existence are presented for the kind of the Rosseland equation with a general growth condition. A linearized map in a closed convex set is defined. The image set is precompact, and thus a fixed point exists. A multi-scale expansion method is used to obtain the homogenized equation. This equation satisfies a similar growth condition.展开更多
Debris flow simulations are useful for predicting the sediment supplied to watersheds from upstream areas. However, the topographic conditions upstream are more complicated than those downstream and the relationship b...Debris flow simulations are useful for predicting the sediment supplied to watersheds from upstream areas. However, the topographic conditions upstream are more complicated than those downstream and the relationship between the topographic conditions and debris flow initiation is not well understood. This study compared the use of several entrainment rate equations in numerical simulations of debris flows to examine the effect of topographic conditions on the flow. One-dimensional numerical simulations were performed based on the shallow water equations and three entrainment rate equations were tested. These entrainment rate equations were based on the same idea that erosion and the deposition of debris flows occur via the difference between the equilibrium and current conditions of debris flows, while they differed in the expression of the concentration, channel angle, and sediment amount. The comparison was performed using a straight channel with various channel angles and a channel with a periodically undulating surface. The three entrainment rate equations gave different amounts of channel bed degradation and hydrographs for a straight channel with a channel angle greater than 21° when water was supplied from upstream at a steady rate. The difference was caused by the expression of the entrainment rate equations. For channels with little undulation, the numerical simulations gave results almost identical to those for straight channels with the same channel angle. However, for channels with large undulations, the hydrographs differed from those for straight channels with the same channel angle when the channel angle was less than 21°. Rapid erosion occurred and the hydrograph showed a significant peak, especially in cases using the entrainment equation expressed by channel angle. This was caused by the effects of the steep undulating sections, since the effect increased with the magnitude of the undulation, suggesting that a debris flow in an upstream area develops differently according to the topographic conditions. These results also inferred that numerical simulations of debris flow can differ depending on the spatial resolution of the simulation domain, as the resolution determines the reproducibility of the undulations.展开更多
In this paper, we get the existence of a weak solution of the following inhomogeneous quasilinear elliptic equation with critical growth conditions: where N≥2, f(x,u)~|u|<sup>m-1</sup>e<sup>b|u|&...In this paper, we get the existence of a weak solution of the following inhomogeneous quasilinear elliptic equation with critical growth conditions: where N≥2, f(x,u)~|u|<sup>m-1</sup>e<sup>b|u|<sup>γ</sup></sup>at +∞, with γ=N/N-1, m≥1, b】0.展开更多
Using the Gleeble-1500 D simulator, the hot deformation behavior and dynamic recrystallization critical conditions of the 10%Ti C/Cu-Al2O3(volume fraction) composite were investigated by compression tests at the tempe...Using the Gleeble-1500 D simulator, the hot deformation behavior and dynamic recrystallization critical conditions of the 10%Ti C/Cu-Al2O3(volume fraction) composite were investigated by compression tests at the temperatures from 450 °C to 850 °C with the strain rates from 0.001 s-1 to 1 s-1. The results show that the softening mechanism of the dynamic recrystallization is a feature of high-temperature flow true stress-strain curves of the composite, and the peak stress increases with the decreasing deformation temperature or the increasing strain rate. The thermal deformation activation energy was calculated as 170.732 k J/mol and the constitutive equation was established. The inflection point in the lnθ-ε curve appears and the minimum value of-(lnθ)/ε-ε curve is presented when the critical state is attained for this composite. The critical strain increases with the increasing strain rate or the decreasing deformation temperature. There is linear relationship between critical strain and peak strain, i.e., εc=0.572εp. The predicting model of critical strain is described by the function of εc=1.062×10-2Z0.0826.展开更多
This paper describes the spectral method for numerically solving Zakharov equation with periodicboundary conditions. This method is spectral method for spatial variable and difference method fortime variable. We make ...This paper describes the spectral method for numerically solving Zakharov equation with periodicboundary conditions. This method is spectral method for spatial variable and difference method fortime variable. We make error estimation of approximate solution and prove the convergence of spectralmethod. We had given the convergence rate. Also, we prove the stability of approximate method forinitial values.展开更多
In this paper we are concerned with the regularity of solutions to the Navier-Stokes equations with the condition on the pressure on parts of the boundary where there is flow. For the steady Stokes problem a result si...In this paper we are concerned with the regularity of solutions to the Navier-Stokes equations with the condition on the pressure on parts of the boundary where there is flow. For the steady Stokes problem a result similar to L q-theory for the one with Dirichlet boundary condition is obtained. Using the result, for the steady Navier-Stokes equations we obtain regularity as the case of Dirichlet boundary conditions. Furthermore,for the time-dependent 2-D Navier-Stokes equations we prove uniqueness and existence of regular solutions,which is similar to J.M.Bernard's results[6]for the time-dependent 2-D Stokes equations.展开更多
In this paper, the concept of generalized ω-periodic solution is given for Riccati's equationy'=a(t)y^2+b(t)y+c(t)with perriodic coefficients, the relation between generalized ω-periodicsolutions and the cha...In this paper, the concept of generalized ω-periodic solution is given for Riccati's equationy'=a(t)y^2+b(t)y+c(t)with perriodic coefficients, the relation between generalized ω-periodicsolutions and the characteristic numbers of system x'_1=c(t)x_2, x'_2=-a(t)x_1-b(t)x_2 is indicated, andseveral necessary and sufficient conditions are given using the coefficients. Moreover, in the case of a(t)without zero, the relation between the number of continuous ω-periodic solutions of y'=a(t)y^2+b(t)y+c(t)+δand the parameter δ is given; thus the problem on the existence of continuous ω-periodic.solutions is basically solved.展开更多
This paper focuses on boundary stabilization of a one-dimensional wave equation with an unstable boundary condition,in which observations are subject to arbitrary fixed time delay.The observability inequality indicate...This paper focuses on boundary stabilization of a one-dimensional wave equation with an unstable boundary condition,in which observations are subject to arbitrary fixed time delay.The observability inequality indicates that the open-loop system is observable,based on which the observer and predictor are designed:The state of system is estimated with available observation and then predicted without observation.After that equivalently the authors transform the original system to the well-posed and exponentially stable system by backstepping method.The equivalent system together with the design of observer and predictor give the estimated output feedback.It is shown that the closed-loop system is exponentially stable.Numerical simulations are presented to illustrate the effect of the stabilizing controller.展开更多
To meet the requirements of fast and automatic computation of subsonic and transonic aerodynamics in aircraft conceptual design,a novel finite volume solver for full potential flows on adaptive Cartesian grids is deve...To meet the requirements of fast and automatic computation of subsonic and transonic aerodynamics in aircraft conceptual design,a novel finite volume solver for full potential flows on adaptive Cartesian grids is developed in this paper.Cartesian grids with geometric adaptation are firstly generated automatically with boundary cells processed by cell-cutting and cell-merging algorithms.The nonlinear full potential equation is discretized by a finite volume scheme on these Cartesian grids and iteratively solved in an implicit fashion with a generalized minimum residual(GMRES) algorithm.During computation,solution-based mesh adaptation is also applied so as to capture flow features more accurately.An improved ghost-cell method is proposed to implement the non-penetration wall boundary condition where the velocity-potential of a ghost cell is modified by an analytic method instead.According to the characteristics of the Cartesian grids,the Kutta condition is applied by specially computing the gradients on Kutta-faces without directly assigning the potential jump to cells adjacent wake faces,which can significantly improve the solution converging speed.The feasibility and accuracy of the proposed method are validated by several typical cases of sub/transonic flows around an ONERA M6 wing,a DLR-F4 wing-body,and an unconventional figuration of a blended wing body(BWB).The validation cases demonstrate a fast convergence with fully automatic grid treatment and computation,and the results suggest its capacity in application for aircraft conceptual design.展开更多
基金Project supported by the National Natural Science Foundation of China(Grant No.51025622)
文摘To deal with the staircase approximation problem in the standard finite-difference time-domain(FDTD) simulation,the two-dimensional boundary condition equations(BCE) method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.
文摘In this paper we give a priori estimates for the maximum modulus of generalizedsolulions of the quasilinear elliplic equations irith anisotropic growth condition.
基金Project supported by the National Natural Science Foundation of China.
文摘This paper deals with the solution of a parametric equation with generalized boundary condiiton in transport theory. It gives the distribution of parameter (so called delta-eigenvalue [1]) with which the equation has non-zero solution. A necessary and sufficient condition for the existence of; he control critical eigenvalue delta0 is established.
文摘A Fourier pseudospectral-finue difference scheme is proposed for three-dimensional vorticity. equalion with unilalerally periodic boundary. condilion. Thegeneralized stability. and conrergence are analyzed. The nunierical results show. theadvantage of this method.
基金Projects Supported by the Science Fund of the Chinese Academy of Sciences
文摘In this paper, some misunderstam igs concerning the necessary conditions for resonance for ordinary differential equations with turning point have been corrected, and a recursive process for finding the sequence of necessary conditions for resonance has beenoffered.
基金This project is supported by National Natural Science Foundation of China(No.50175031).
文摘Accurate boundary conditions of composite material plates with different holes are founded to settle boundary condition problems of complex holes by conformal mapping method upon the nonhomogeneous anisotropic elastic and complex function theory. And then the two stress functions required were founded on Cauchy integral by boundary conditions. The final stress distributions of opening structure and the analytical solution on composite material plate with rectangle hole and wing manholes were achieved. The influences on hole-edge stress concentration factors are discussed under different loads and fiber direction cases, and then contrast calculates are carried through FEM.
文摘There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle these kind of mixed boundary-value problems associated with the Laplace’s equation (or Helmholtz equation) arising in the study of waves propagating through solids or fluids. One of the widely used methods in wave structure interaction is the multipole expansion method. This expansion involves a general combination of a regular wave, a wave source, a wave dipole and a regular wave-free part. The wave-free part can be further expanded in terms of wave-free multipoles which are termed as wave-free potentials. These are singular solutions of Laplace’s equation or two-dimensional Helmholz equation. Construction of these wave-free potentials and multipoles are presented here in a systematic manner for a number of situations such as two-dimensional non-oblique and oblique waves, three dimensional waves in two-layer fluid with free surface condition with higher order partial derivative are considered. In particular, these are obtained taking into account of the effect of the presence of surface tension at the free surface and also in the presence of an ice-cover modelled as a thin elastic plate. Also for limiting case, it can be shown that the multipoles and wave-free potential functions go over to the single layer multipoles and wave-free potential.
文摘With the help of skew-symmetric differential forms, the hidden properties of the mathematical physics equations that describe discrete quantum transitions and emergence the physical structures are investigated. It is shown that the mathematical physics equations possess a unique property. They can describe discrete quantum transitions, emergence of physical structures and occurrence observed formations. However, such a property possesses only equations on which no additional conditions, namely, the conditions of integrability, are imposed. The intergrability conditions are realized from the equations themselves. Just under realization of integrability conditions double solutions to the mathematical physics equations, which describe discrete transitions and so on, are obtained. The peculiarity consists in the fact that the integrability conditions do not directly follow from the mathematical physics equations;they are realized under the description of evolutionary process. The hidden properties of differential equations were discovered when studying the integrability of differential equations of mathematical physics that depends on the consistence between the derivatives in differential equations along different directions and on the consistence of equations in the set of equations. The results of this work were obtained with the help of skew-symmetric differential forms that possess a nontraditional mathematical apparatus such as nonidentical relations, degenerate transformations and the transition from nonintegrable manifolds to integrable structures. Such results show that mathematical physics equations can describe quantum processes.
基金Supported by the National Basic Research Program of China(973 Program)(No.2012CB025904)the National Natural Science Foundation of China(No.90916027)
文摘The global boundness and existence are presented for the kind of the Rosseland equation with a general growth condition. A linearized map in a closed convex set is defined. The image set is precompact, and thus a fixed point exists. A multi-scale expansion method is used to obtain the homogenized equation. This equation satisfies a similar growth condition.
基金partially supported by Grant-in-Aid for Scientific Research 26292077, 2014, from the Ministry of Education, Science, Sports, and Culture of Japanby the River Fund in charge of the River Foundation, Japan
文摘Debris flow simulations are useful for predicting the sediment supplied to watersheds from upstream areas. However, the topographic conditions upstream are more complicated than those downstream and the relationship between the topographic conditions and debris flow initiation is not well understood. This study compared the use of several entrainment rate equations in numerical simulations of debris flows to examine the effect of topographic conditions on the flow. One-dimensional numerical simulations were performed based on the shallow water equations and three entrainment rate equations were tested. These entrainment rate equations were based on the same idea that erosion and the deposition of debris flows occur via the difference between the equilibrium and current conditions of debris flows, while they differed in the expression of the concentration, channel angle, and sediment amount. The comparison was performed using a straight channel with various channel angles and a channel with a periodically undulating surface. The three entrainment rate equations gave different amounts of channel bed degradation and hydrographs for a straight channel with a channel angle greater than 21° when water was supplied from upstream at a steady rate. The difference was caused by the expression of the entrainment rate equations. For channels with little undulation, the numerical simulations gave results almost identical to those for straight channels with the same channel angle. However, for channels with large undulations, the hydrographs differed from those for straight channels with the same channel angle when the channel angle was less than 21°. Rapid erosion occurred and the hydrograph showed a significant peak, especially in cases using the entrainment equation expressed by channel angle. This was caused by the effects of the steep undulating sections, since the effect increased with the magnitude of the undulation, suggesting that a debris flow in an upstream area develops differently according to the topographic conditions. These results also inferred that numerical simulations of debris flow can differ depending on the spatial resolution of the simulation domain, as the resolution determines the reproducibility of the undulations.
基金Supported by the Youth FoundationNatural Science Foundation, People's Republic of China.
文摘In this paper, we get the existence of a weak solution of the following inhomogeneous quasilinear elliptic equation with critical growth conditions: where N≥2, f(x,u)~|u|<sup>m-1</sup>e<sup>b|u|<sup>γ</sup></sup>at +∞, with γ=N/N-1, m≥1, b】0.
基金Project(51101052) supported by the National Natural Science Foundation of China
文摘Using the Gleeble-1500 D simulator, the hot deformation behavior and dynamic recrystallization critical conditions of the 10%Ti C/Cu-Al2O3(volume fraction) composite were investigated by compression tests at the temperatures from 450 °C to 850 °C with the strain rates from 0.001 s-1 to 1 s-1. The results show that the softening mechanism of the dynamic recrystallization is a feature of high-temperature flow true stress-strain curves of the composite, and the peak stress increases with the decreasing deformation temperature or the increasing strain rate. The thermal deformation activation energy was calculated as 170.732 k J/mol and the constitutive equation was established. The inflection point in the lnθ-ε curve appears and the minimum value of-(lnθ)/ε-ε curve is presented when the critical state is attained for this composite. The critical strain increases with the increasing strain rate or the decreasing deformation temperature. There is linear relationship between critical strain and peak strain, i.e., εc=0.572εp. The predicting model of critical strain is described by the function of εc=1.062×10-2Z0.0826.
基金Project supported by the Science Foundation of the Chinese Academy of Sciences
文摘This paper describes the spectral method for numerically solving Zakharov equation with periodicboundary conditions. This method is spectral method for spatial variable and difference method fortime variable. We make error estimation of approximate solution and prove the convergence of spectralmethod. We had given the convergence rate. Also, we prove the stability of approximate method forinitial values.
基金Supported by TWAS,UNESCO and AMSS in Chinese Academy of Sciences
文摘In this paper we are concerned with the regularity of solutions to the Navier-Stokes equations with the condition on the pressure on parts of the boundary where there is flow. For the steady Stokes problem a result similar to L q-theory for the one with Dirichlet boundary condition is obtained. Using the result, for the steady Navier-Stokes equations we obtain regularity as the case of Dirichlet boundary conditions. Furthermore,for the time-dependent 2-D Navier-Stokes equations we prove uniqueness and existence of regular solutions,which is similar to J.M.Bernard's results[6]for the time-dependent 2-D Stokes equations.
文摘In this paper, the concept of generalized ω-periodic solution is given for Riccati's equationy'=a(t)y^2+b(t)y+c(t)with perriodic coefficients, the relation between generalized ω-periodicsolutions and the characteristic numbers of system x'_1=c(t)x_2, x'_2=-a(t)x_1-b(t)x_2 is indicated, andseveral necessary and sufficient conditions are given using the coefficients. Moreover, in the case of a(t)without zero, the relation between the number of continuous ω-periodic solutions of y'=a(t)y^2+b(t)y+c(t)+δand the parameter δ is given; thus the problem on the existence of continuous ω-periodic.solutions is basically solved.
基金supported by the National Natural Science Foundation of China under Grant No.61203058the Training Program for Outstanding Young Teachers of North China University of Technology under Grant No.XN131+1 种基金the Construction Plan for Innovative Research Team of North China University of Technology under Grant No.XN129the Laboratory construction for Mathematics Network Teaching Platform of North China University of Technology under Grant No.XN041
文摘This paper focuses on boundary stabilization of a one-dimensional wave equation with an unstable boundary condition,in which observations are subject to arbitrary fixed time delay.The observability inequality indicates that the open-loop system is observable,based on which the observer and predictor are designed:The state of system is estimated with available observation and then predicted without observation.After that equivalently the authors transform the original system to the well-posed and exponentially stable system by backstepping method.The equivalent system together with the design of observer and predictor give the estimated output feedback.It is shown that the closed-loop system is exponentially stable.Numerical simulations are presented to illustrate the effect of the stabilizing controller.
基金co-supported by the National Natural Science Foundation of China(No.11672133)the Fundamental Research Funds for the Central UniversitiesThe support from the Priority Academic Program Development(PAPD)of Jiangsu Higher Education Institutions
文摘To meet the requirements of fast and automatic computation of subsonic and transonic aerodynamics in aircraft conceptual design,a novel finite volume solver for full potential flows on adaptive Cartesian grids is developed in this paper.Cartesian grids with geometric adaptation are firstly generated automatically with boundary cells processed by cell-cutting and cell-merging algorithms.The nonlinear full potential equation is discretized by a finite volume scheme on these Cartesian grids and iteratively solved in an implicit fashion with a generalized minimum residual(GMRES) algorithm.During computation,solution-based mesh adaptation is also applied so as to capture flow features more accurately.An improved ghost-cell method is proposed to implement the non-penetration wall boundary condition where the velocity-potential of a ghost cell is modified by an analytic method instead.According to the characteristics of the Cartesian grids,the Kutta condition is applied by specially computing the gradients on Kutta-faces without directly assigning the potential jump to cells adjacent wake faces,which can significantly improve the solution converging speed.The feasibility and accuracy of the proposed method are validated by several typical cases of sub/transonic flows around an ONERA M6 wing,a DLR-F4 wing-body,and an unconventional figuration of a blended wing body(BWB).The validation cases demonstrate a fast convergence with fully automatic grid treatment and computation,and the results suggest its capacity in application for aircraft conceptual design.