Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is a...Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is analyzed. The dynamic model uses Euler-Bernoulli beam theory for cantilever beam, in which the driving force (or driving moment) of the electromagnetic system is taken into account, and the contact force between moving contact and stationary contact is simulated by the Kelvin-Voigt vis-coelastic...展开更多
Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with ...Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with theoretical method, the finite difference method has been verified to be feasible by a case study. It is found that under seismic loading, the dynamic response of anchorage system is synchronously fluctuated with the seismic vibration. The change of displacement amplitude of material points is slight, and comparatively speaking, the displacement amplitude of the outside point is a little larger than that of the inside point, which shows amplification effect of surface. While the axial force amplitude transforms considerably from the inside to the outside. It increases first and reaches the peak value in the intersection between the anchoring section and free section, then decreases slowly in the free section. When considering damping effect of anchorage system, the finite difference method can reflect the time attenuation characteristic better, and the calculating result would be safer and more reasonable than the dynamic steady-state theoretical method. What is more, the finite difference method can be applied to the dynamic response analysis of harmonic and seismic random vibration for all kinds of anchor, and hence has a broad application prospect.展开更多
For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study prop...For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.展开更多
Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the ...Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference to展开更多
For the three-dimensional convection-dominated problem of dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Fract...For the three-dimensional convection-dominated problem of dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Fractional steps techniques are needed to convert a multi-dimensional problem into a series of successive one-dimensional problems. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators, and the theory of prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources and predicting the consequences of seawater intrusion and protection projects.展开更多
Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been...Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.展开更多
The coupled nonlocal nonlinear Schrödinger equations with variable coefficients are researched using the nonstandard Hirota bilinear method.The two-soliton and double-hump one-soliton solutions for the equations ...The coupled nonlocal nonlinear Schrödinger equations with variable coefficients are researched using the nonstandard Hirota bilinear method.The two-soliton and double-hump one-soliton solutions for the equations are first obtained.By assigning different functions to the variable coefficients,we obtain V-shaped,Y-shaped,wave-type,exponential solitons,and so on.Next,we reveal the influence of the real and imaginary parts of the wave numbers on the double-hump structure based on the soliton solutions.Finally,by setting different wave numbers,we can change the distance and transmission direction of the solitons to analyze their dynamic behavior during collisions.This study establishes a theoretical framework for controlling the dynamics of optical fiber in nonlocal nonlinear systems.展开更多
The value of form factor k at different drafts is important in predicting full-scale total resistance and speed for different types of ships. In the ITTC community, most organizations predict form factor k using a low...The value of form factor k at different drafts is important in predicting full-scale total resistance and speed for different types of ships. In the ITTC community, most organizations predict form factor k using a low-speed model test. However, this method is problematic for ships with bulbous bows and transom. In this article, a Computational Fluid Dynamics(CFD)-based method is introduced to obtain k for different type of ships at different drafts, and a comparison is made between the CFD method and the model test. The results show that the CFD method produces reasonable k values. A grid generating method and turbulence model are briefly discussed in the context of obtaining a consistent k using CFD.展开更多
This paper presents a numerical analysis of the dynamic transient behaviors of undersea cables.In this numerical study,the governing equations based on Euler-Bernoulli beam theory are adopted,and they can satisfy many...This paper presents a numerical analysis of the dynamic transient behaviors of undersea cables.In this numerical study,the governing equations based on Euler-Bernoulli beam theory are adopted,and they can satisfy many applications no matter what the magnitude of the cable tension is.The nonlinear coupled equations are solved by a popular central finite difference method,and the numerical results of transient behaviors are presented when several kinds of surrounding conditions,such as different towing speeds of surface vessel,different currents and waves with various frequencies and amplitudes,are exerted.Then a detailed comparison of the results,including the upper end tension and cable shape in time-domain,is made under the above external excitations,and finally the possible reasons for these are further explained.展开更多
A mathematical model has been developed to describe the dynamic heat transfer in the clothing microclimate under transient wear conditions. This model is solved numerically by the implicit finite difference method. If...A mathematical model has been developed to describe the dynamic heat transfer in the clothing microclimate under transient wear conditions. This model is solved numerically by the implicit finite difference method. If the physical activity and ambient conditions are specified, the model can predict the thermoregulatory response of the body. Experimental measurements with garments made of fibers with different levels of hygroscopicity are compared with predictions by the model. There is good agreement between prediction and experiment for the temperature of the clothing microclimate.展开更多
This paper presents a numerical investigation into the dynamics of marine cables which are extensively used in offshore industry. In this numerical study, the Euler-Bernoulli beam model is adopted to develop the gover...This paper presents a numerical investigation into the dynamics of marine cables which are extensively used in offshore industry. In this numerical study, the Euler-Bernoulli beam model is adopted to develop the governing equations of the cable. Bending stiffness is considered to cope with the low tension problem in local area of towing cable, and thus a more accurate solution with the consideration of the axial elongation can be given.The derived strongly-coupled and nonlinear governing equations are solved by a second-order accurate, implicit,and large time step stable central finite difference method. The quadratically convergent Newton-Raphson iteration method is applied to solving the discrete nonlinear algebraic equations. Then a towed array sonar system(TASS)problem is studied. The numerical solutions agree reasonably well with the experimental data and the simulated results of the references. The specified program of the present paper shows great robustness with high efficiency.展开更多
In this paper.the equations of motion of axisymmetrically laminated cylindrical orthotropic spherical shells are derived.Theeffects of transverse shear deformation and rotatory inertia are considered.On this basis,th...In this paper.the equations of motion of axisymmetrically laminated cylindrical orthotropic spherical shells are derived.Theeffects of transverse shear deformation and rotatory inertia are considered.On this basis,the dynamic response of spherical shells under axisymmetric dynamic load is calculated using the finite difference method The effects of material parameters.structural parameters and transverse shear dgformation are discussed.展开更多
We used simulated data to investigate both the small and large sample properties of the within-groups (WG) estimator and the first difference generalized method of moments (FD-GMM) estimator of a dynamic panel data (D...We used simulated data to investigate both the small and large sample properties of the within-groups (WG) estimator and the first difference generalized method of moments (FD-GMM) estimator of a dynamic panel data (DPD) model. The magnitude of WG and FD-GMM estimates are almost the same for square panels. WG estimator performs best for long panels such as those with time dimension as large as 50. The advantage of FD-GMM estimator however, is observed on panels that are long and wide, say with time dimension at least 25 and cross-section dimension size of at least 30. For small-sized panels, the two methods failed since their optimality was established in the context of asymptotic theory. We developed parametric bootstrap versions of WG and FD-GMM estimators. Simulation study indicates the advantages of the bootstrap methods under small sample cases on the assumption that variances of the individual effects and the disturbances are of similar magnitude. The boostrapped WG and FD-GMM estimators are optimal for small samples.展开更多
For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-d...For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as implicit-explicit difference scheme, calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L 2 norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources. Keywords: combinatorial system, multilayer dynamics of fluids in porous media, two-class upwind finite difference fractional steps method, convergence, numerical simulation of energy sources.展开更多
The generalized finite difference method (GFDM) used for irregular grids is first introduced into the numerical study of thelevel set equation, which is coupled with the theory of detonation shock dynamics (DSD) descr...The generalized finite difference method (GFDM) used for irregular grids is first introduced into the numerical study of thelevel set equation, which is coupled with the theory of detonation shock dynamics (DSD) describing the propagation of thedetonation shock front. The numerical results of a rate-stick problem, a converging channel problem and an arc channel prob-lem for specified boundaries show that GFDM is effective on solving the level set equation in the irregular geometrical domain.The arrival time and the normal velocity distribution of the detonation shock front of these problems can then be obtainedconveniently with this method. The numerical results also confirm that when there is a curvature effect, the theory of DSDmust be considered for the propagation of detonation shock surface, while classic Huygens construction is not suitable anymore.展开更多
For nonlinear coupled system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward, trod...For nonlinear coupled system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward, trod two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L2 norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources.展开更多
Data, including the spatial data and the non-spatial data, are the basis of all digital scientific engineering projects, such as the digital earth and the digital nation, the digital mine. The spatial data have the ch...Data, including the spatial data and the non-spatial data, are the basis of all digital scientific engineering projects, such as the digital earth and the digital nation, the digital mine. The spatial data have the characteristics of many sources, multi-dimension, multi-type, many time states and different accuracy. The spatial data firstly must be processed before using these data. The parameter estimation model to process the data is commonly the more complex nonlinear model including random parameters and non-random parameters. So a generalized nonlinear dynamic least squares method to process these data is put forward. According to the special structure of the generalized nonlinear dynamic least squares problem and the solution to the first order, a new solving model and a corresponding method to process the problem are put forward. The complex problem can be divided into two sub-problems so that the number of the unknown parameters is reduced largely. Therefore it reduces the computing difficulty and load.展开更多
We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is ...We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.展开更多
文摘Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is analyzed. The dynamic model uses Euler-Bernoulli beam theory for cantilever beam, in which the driving force (or driving moment) of the electromagnetic system is taken into account, and the contact force between moving contact and stationary contact is simulated by the Kelvin-Voigt vis-coelastic...
基金Projects(51308273,41372307,41272326) supported by the National Natural Science Foundation of ChinaProjects(2010(A)06-b) supported by Science and Technology Fund of Yunan Provincial Communication Department,China
文摘Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with theoretical method, the finite difference method has been verified to be feasible by a case study. It is found that under seismic loading, the dynamic response of anchorage system is synchronously fluctuated with the seismic vibration. The change of displacement amplitude of material points is slight, and comparatively speaking, the displacement amplitude of the outside point is a little larger than that of the inside point, which shows amplification effect of surface. While the axial force amplitude transforms considerably from the inside to the outside. It increases first and reaches the peak value in the intersection between the anchoring section and free section, then decreases slowly in the free section. When considering damping effect of anchorage system, the finite difference method can reflect the time attenuation characteristic better, and the calculating result would be safer and more reasonable than the dynamic steady-state theoretical method. What is more, the finite difference method can be applied to the dynamic response analysis of harmonic and seismic random vibration for all kinds of anchor, and hence has a broad application prospect.
基金National Natural Science Foundation of China under Grant Nos.51978213 and 51778190the National Key Research and Development Program of China under Grant Nos.2017YFC0703605 and 2016YFC0701106。
文摘For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.
文摘Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference to
基金Project supported by the Major State Basic Research Program of China (No.G1999032803)the National Tackling Key Problems Program (No.20050200069)the National Natural Science Foundation of China (Nos.10372052, 10271066)the Doctoral Foundation of Ministry of Education of China (No.20030422047).
文摘For the three-dimensional convection-dominated problem of dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Fractional steps techniques are needed to convert a multi-dimensional problem into a series of successive one-dimensional problems. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators, and the theory of prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources and predicting the consequences of seawater intrusion and protection projects.
文摘Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.
基金supported by the National Key R&D Program of China(Grant No.2022YFA1604200)the National Natural Science Foundation of China(Grant No.12261131495)Institute of Systems Science,Beijing Wuzi University(Grant No.BWUISS21).
文摘The coupled nonlocal nonlinear Schrödinger equations with variable coefficients are researched using the nonstandard Hirota bilinear method.The two-soliton and double-hump one-soliton solutions for the equations are first obtained.By assigning different functions to the variable coefficients,we obtain V-shaped,Y-shaped,wave-type,exponential solitons,and so on.Next,we reveal the influence of the real and imaginary parts of the wave numbers on the double-hump structure based on the soliton solutions.Finally,by setting different wave numbers,we can change the distance and transmission direction of the solitons to analyze their dynamic behavior during collisions.This study establishes a theoretical framework for controlling the dynamics of optical fiber in nonlocal nonlinear systems.
基金Supported by Ministry of Industry and Information(No.K24097)
文摘The value of form factor k at different drafts is important in predicting full-scale total resistance and speed for different types of ships. In the ITTC community, most organizations predict form factor k using a low-speed model test. However, this method is problematic for ships with bulbous bows and transom. In this article, a Computational Fluid Dynamics(CFD)-based method is introduced to obtain k for different type of ships at different drafts, and a comparison is made between the CFD method and the model test. The results show that the CFD method produces reasonable k values. A grid generating method and turbulence model are briefly discussed in the context of obtaining a consistent k using CFD.
基金the National High Technology Research and Development Program (863) of China(No. 2008AA092301-1)the National Natural Science Foundation of China(No. 50909061)the Ph.D. Programs Foundation of Ministry of Education of China(No. 20070248103)
文摘This paper presents a numerical analysis of the dynamic transient behaviors of undersea cables.In this numerical study,the governing equations based on Euler-Bernoulli beam theory are adopted,and they can satisfy many applications no matter what the magnitude of the cable tension is.The nonlinear coupled equations are solved by a popular central finite difference method,and the numerical results of transient behaviors are presented when several kinds of surrounding conditions,such as different towing speeds of surface vessel,different currents and waves with various frequencies and amplitudes,are exerted.Then a detailed comparison of the results,including the upper end tension and cable shape in time-domain,is made under the above external excitations,and finally the possible reasons for these are further explained.
文摘A mathematical model has been developed to describe the dynamic heat transfer in the clothing microclimate under transient wear conditions. This model is solved numerically by the implicit finite difference method. If the physical activity and ambient conditions are specified, the model can predict the thermoregulatory response of the body. Experimental measurements with garments made of fibers with different levels of hygroscopicity are compared with predictions by the model. There is good agreement between prediction and experiment for the temperature of the clothing microclimate.
基金the National Science and Technology Major Project(No.2011ZX05027-004)the National Natural Science Foundation of China(No.51279107)
文摘This paper presents a numerical investigation into the dynamics of marine cables which are extensively used in offshore industry. In this numerical study, the Euler-Bernoulli beam model is adopted to develop the governing equations of the cable. Bending stiffness is considered to cope with the low tension problem in local area of towing cable, and thus a more accurate solution with the consideration of the axial elongation can be given.The derived strongly-coupled and nonlinear governing equations are solved by a second-order accurate, implicit,and large time step stable central finite difference method. The quadratically convergent Newton-Raphson iteration method is applied to solving the discrete nonlinear algebraic equations. Then a towed array sonar system(TASS)problem is studied. The numerical solutions agree reasonably well with the experimental data and the simulated results of the references. The specified program of the present paper shows great robustness with high efficiency.
文摘In this paper.the equations of motion of axisymmetrically laminated cylindrical orthotropic spherical shells are derived.Theeffects of transverse shear deformation and rotatory inertia are considered.On this basis,the dynamic response of spherical shells under axisymmetric dynamic load is calculated using the finite difference method The effects of material parameters.structural parameters and transverse shear dgformation are discussed.
文摘We used simulated data to investigate both the small and large sample properties of the within-groups (WG) estimator and the first difference generalized method of moments (FD-GMM) estimator of a dynamic panel data (DPD) model. The magnitude of WG and FD-GMM estimates are almost the same for square panels. WG estimator performs best for long panels such as those with time dimension as large as 50. The advantage of FD-GMM estimator however, is observed on panels that are long and wide, say with time dimension at least 25 and cross-section dimension size of at least 30. For small-sized panels, the two methods failed since their optimality was established in the context of asymptotic theory. We developed parametric bootstrap versions of WG and FD-GMM estimators. Simulation study indicates the advantages of the bootstrap methods under small sample cases on the assumption that variances of the individual effects and the disturbances are of similar magnitude. The boostrapped WG and FD-GMM estimators are optimal for small samples.
基金This work was partly supported by the National Natural Science Foundation of China （11501581）, the Project by Central South Uni- versity （502042032）, and the China Postdoctoral Science Foundation （2015M570683）.
基金This work was supported by the Major State Basic Research Program of China(Grant No. 1990328) the National Tackling Key Problem Program, the National Natural Science Foundation of China (Grant Nos. 19871051 and 19972039) the Doctorate Foundation
文摘For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as implicit-explicit difference scheme, calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L 2 norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources. Keywords: combinatorial system, multilayer dynamics of fluids in porous media, two-class upwind finite difference fractional steps method, convergence, numerical simulation of energy sources.
基金supported by the National Natural Science Foundation of China (Grant No. 11002029)
文摘The generalized finite difference method (GFDM) used for irregular grids is first introduced into the numerical study of thelevel set equation, which is coupled with the theory of detonation shock dynamics (DSD) describing the propagation of thedetonation shock front. The numerical results of a rate-stick problem, a converging channel problem and an arc channel prob-lem for specified boundaries show that GFDM is effective on solving the level set equation in the irregular geometrical domain.The arrival time and the normal velocity distribution of the detonation shock front of these problems can then be obtainedconveniently with this method. The numerical results also confirm that when there is a curvature effect, the theory of DSDmust be considered for the propagation of detonation shock surface, while classic Huygens construction is not suitable anymore.
文摘For nonlinear coupled system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward, trod two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L2 norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources.
基金Project (40174003) supported by the National Natural Science Foundation of China
文摘Data, including the spatial data and the non-spatial data, are the basis of all digital scientific engineering projects, such as the digital earth and the digital nation, the digital mine. The spatial data have the characteristics of many sources, multi-dimension, multi-type, many time states and different accuracy. The spatial data firstly must be processed before using these data. The parameter estimation model to process the data is commonly the more complex nonlinear model including random parameters and non-random parameters. So a generalized nonlinear dynamic least squares method to process these data is put forward. According to the special structure of the generalized nonlinear dynamic least squares problem and the solution to the first order, a new solving model and a corresponding method to process the problem are put forward. The complex problem can be divided into two sub-problems so that the number of the unknown parameters is reduced largely. Therefore it reduces the computing difficulty and load.
基金supported by the National Natural Science Foundation of China(Nos.12132001 and 52192632)。
文摘We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.
