To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection techniq...To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).展开更多
This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
This paper discusses a fictitious domain method for the linear Dirichlet problem and its applications to the generalized Stokes problem. This method treats Dirichlet boundary condit ion via a Lagrange multiplier tec...This paper discusses a fictitious domain method for the linear Dirichlet problem and its applications to the generalized Stokes problem. This method treats Dirichlet boundary condit ion via a Lagrange multiplier technique and is well suited to the no-slip bound ary condition in viscous flow problems. In order to improve the accuracy of solu tions, meshes are refined according to the a posteriori error estimate. The mini -element discretization is applied to solve the generalized Stokes problem. Fin ally, some numerical results to validate this method are presented for partial d ifferential equations with Dirichlet boundary condition.展开更多
The dynamic response of pile in layered soil is theoretically investigated when considering the transverse inertia effect.Firstly, the fictitious soil-pile model is employed to simulate the dynamic interaction between...The dynamic response of pile in layered soil is theoretically investigated when considering the transverse inertia effect.Firstly, the fictitious soil-pile model is employed to simulate the dynamic interaction between the pile and the soil layers beneath pile toe. The dynamic interactions of adjacent soil layers along the vertical direction are simplified as distributed Voigt models.Meanwhile, the pile and fictitious soil-pile are assumed to be viscoelastic Rayleigh-Love rods, and both the radial and vertical displacement continuity conditions at the soil-pile interface are taken into consideration. On this basis, the analytical solution for dynamic response at the pile head is derived in the frequency domain and the corresponding quasi-analytical solution in the time domain is then obtained by means of the convolution theorem. Following this, the accuracy and parameter value of the hypothetical boundaries for soil-layer interfaces are discussed. Comparisons with published solution and measured data are carried out to verify the rationality of the present solution. Parametric analyses are further conducted by using the present solution to investigate the relationships between the transverse inertia effects and soil-pile parameters.展开更多
The potential field determined based on the fictitious compress recovery approach is influenced by the errors contained in the boundary (the Earth's surface or the surface corresponding to the satellite altitude) v...The potential field determined based on the fictitious compress recovery approach is influenced by the errors contained in the boundary (the Earth's surface or the surface corresponding to the satellite altitude) values. Given the boundary value with definite accuracy, the accuracy of the field determined based on the fictitious compress recovery approach is estimated, and it is theoretically shown that the determined field has the same accuracy level as the given boundary value.展开更多
The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully disc...The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.展开更多
The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which...The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which can be solved by traditional numerical methods and an infinite domain with a crack which can be solved by Muskhelishvili’s solutions.However,this alternating method cannot be directly applied to the edge crack problems since partial crack surface of Muskhelishvili’s solutions is located outside the computational domain.In this paper,an improved alternating method,the spline fictitious boundary element alternating method(SFBEAM),based on infinite domain with the combination of spline fictitious boundary element method(SFBEM)and Muskhelishvili’s solutions is proposed to solve the edge crack problems.Since the SFBEM and Muskhelishvili’s solutions are obtained in the framework of infinite domain,no special treatment is needed for solving the problem of edge cracks.Different mixed boundary conditions edge crack problems with varies of computational parameters are given to certify the high precision,efficiency and applicability of the proposed method compared with other alternating methods and extend finite element method.展开更多
In modern computer games, "bots" - intelligent realistic agents play a prominent role in the popularity of a game in the market. Typically, bots are modeled using finite-state machine and then programmed via simple ...In modern computer games, "bots" - intelligent realistic agents play a prominent role in the popularity of a game in the market. Typically, bots are modeled using finite-state machine and then programmed via simple conditional statements which are hard-coded in bots logic. Since these bots have become quite predictable to an experienced games' player, a player might lose interest in the game. We propose the use of a game theoretic based learning rule called fictitious play for improving behavior of these computer game bots which will make them less predictable and hence, more a enjoyable game.展开更多
This paper investigates the stability and convergence of some knowndifference schemes for the numerical solution to heat conduction equation withderivative boundary conditions by the fictitious domain method.The discr...This paper investigates the stability and convergence of some knowndifference schemes for the numerical solution to heat conduction equation withderivative boundary conditions by the fictitious domain method.The discrete vari-ables at the false mesh points are firstly eliminated from the difference schemes andthe local truncation errors are then analyzed in detail.The stability and convergenceof the schemes are proved by energy method.An improvement is proposed to obtainbetter schemes over the original ones.Several numerical examples and comparisonswith other schemes are presented.展开更多
This paper investigates some known difference schemes for the numerical solution to parabolic differential equation with derivative boundary conditions by the fictitious domain method.The stability and convergence in...This paper investigates some known difference schemes for the numerical solution to parabolic differential equation with derivative boundary conditions by the fictitious domain method.The stability and convergence in L ∞ are proven.展开更多
Abstract In this article we briefly studied the fictitious domain methods for steady and nonsteady Navier-Stokes equations based on Penalty factorεon the extended domain. The convergence u<sup>?</sup>→u ...Abstract In this article we briefly studied the fictitious domain methods for steady and nonsteady Navier-Stokes equations based on Penalty factorεon the extended domain. The convergence u<sup>?</sup>→u in H<sub>0</sub><sup>1</sup>(Ω)<sup>d</sup> and L<sup>2</sup>(Ω)<sup>d</sup>(d=2,3) is given as well as p<sup>?</sup>→p in L<sub>0</sub><sup>2</sup>(Ω).展开更多
In this paper, an efficient multigrid fictitious boundary method (MFBM) coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NAC...In this paper, an efficient multigrid fictitious boundary method (MFBM) coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.展开更多
This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the...This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as "fictitious" fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure' s rigid body shape and behavior, a rigid body constraint is enforced on the "fictitious" fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain (DLM/ FD) method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.展开更多
To apply the fictitious domain method and conduct numericalexperiments, a boundary value problem for an ordinary differential equation is considered. The results of numerical calculations for different valuesof the it...To apply the fictitious domain method and conduct numericalexperiments, a boundary value problem for an ordinary differential equation is considered. The results of numerical calculations for different valuesof the iterative parameter τ and the small parameter ε are presented. Astudy of the auxiliary problem of the fictitious domain method for NavierStokes equations with continuation into a fictitious subdomain by highercoefficients with a small parameter is carried out. A generalized solutionof the auxiliary problem of the fictitious domain method with continuationby higher coefficients with a small parameter is determined. After all theabove mathematical studies, a computational algorithm has been developedfor the numerical solution of the problem. Two methods were used to solvethe problem numerically. The first variant is the fictitious domain methodassociated with the modification of nonlinear terms in a fictitious subdomain.The model problem shows the effectiveness of using such a modification. Theproposed version of the method is used to solve two problems at once that arisewhile numerically solving systems of Navier-Stokes equations: the problem ofa curved boundary of an arbitrary domain and the problem of absence of aboundary condition for pressure in physical formulation of the internal flowproblem. The main advantage of this method is its universality in developmentof computer programs. The second method used calculation on a uniform gridinside the area. When numerically implementing the solution on a uniformgrid inside the domain, using this method it’s possible to accurately take intoaccount the boundaries of the curved domain and ensure the accuracy of thevalue of the function at the boundaries of the domain. Methodical calculationswere carried out, the results of numerical calculations were obtained. Whenconducting numerical experiments in both cases, quantitative and qualitativeindicators of numerical results coincide.展开更多
The fictitious compress recovery approach is introduced, which could be applied to the establishment of the Rungerarup theorem, the determination of the Bjerhammar's fictitious gravity anomaly, the solution of the "...The fictitious compress recovery approach is introduced, which could be applied to the establishment of the Rungerarup theorem, the determination of the Bjerhammar's fictitious gravity anomaly, the solution of the "downward con- tinuation" problem of the gravity field, the confirmation of the convergence of the spherical harmonic expansion series of the Earth's potential field, and the gravity field determination in three cases: gravitational potential case, gravitation case, and gravitational gradient case. Several tests using simulation experiments show that the fictitious compress recovery approach shows promise in physical geodesy applications.展开更多
Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the ...Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the fictitious crack model in fracture mechanics of concrete. The new analytical element can be implemented into FEM program systems to solve fictitious crack propagation problems for concrete cracked plates with arbitrary shapes and loads. Numerical results indicate that the method is more efficient and accurate than ordinary finite element method.展开更多
We propose a novel on-chip platform for controlling and manipulating cold atoms precisely and coherently. The scheme is achieved by producing optically induced fictitious magnetic traps(OFMTs) with 790 nm(for -(87...We propose a novel on-chip platform for controlling and manipulating cold atoms precisely and coherently. The scheme is achieved by producing optically induced fictitious magnetic traps(OFMTs) with 790 nm(for -(87)Rb) circularly polarized laser beams and state-dependent potentials simultaneously for two internal atomic states with microwave coplanar waveguides. We carry out numerical calculations and simulations for controlled collisional interactions between OFMTs and addressable single atoms' manipulation on our designed hybrid atom chips. The results show that our proposed platform is feasible and flexible, which has wide applications including collisional dynamics investigation, entanglement generation,and scalable quantum gates implementation.展开更多
As The Three-Body Problem translated by Ken Liu went popular in America, Chinese science fiction has got spotlightfrom all over the world. This article studies the translation of fictitious elements in The Three-Body ...As The Three-Body Problem translated by Ken Liu went popular in America, Chinese science fiction has got spotlightfrom all over the world. This article studies the translation of fictitious elements in The Three-Body Problem. Based on the specificmethods Ken Liu employs, general requirements for translators are proposed and strategies of translating fictitious elements areelicited. It intends to bring forth enlightening ideas in translating fictitious elements of science fiction and shed lights on this areaof study which few have set foot on.展开更多
A fictitious soil pile(FSP)model is developed to simulate the behavior of pipe piles with soil plugs undergoing high-strain dynamic impact loading.The developed model simulates the base soil with a fictitious hollow p...A fictitious soil pile(FSP)model is developed to simulate the behavior of pipe piles with soil plugs undergoing high-strain dynamic impact loading.The developed model simulates the base soil with a fictitious hollow pile fully filled with a soil plug extending at a cone angle from the pile toe to the bedrock.The friction on the outside and inside of the pile walls is distinguished using different shaft models,and the propagation of stress waves in the base soil and soil plug is considered.The motions of the pile—soil system are solved by discretizing them into spring-mass model based on the finite difference method.Comparisons of the predictions of the proposed model and conventional numerical models,as well as measurements for pipe piles in field tests subjected to impact loading,validate the accuracy of the proposed model.A parametric analysis is conducted to illustrate the influence of the model parameters on the pile dynamic response.Finally,the effective length of the FSP is proposed to approximate the affected soil zone below the pipe pile toe,and some guidance is provided for the selection of the model parameters.展开更多
To investigate the interaction of the bolt-reinforced rock and the surface support,an analytical model of the convergence-confinement type is proposed,considering the sequential installation of the fully grouted rockb...To investigate the interaction of the bolt-reinforced rock and the surface support,an analytical model of the convergence-confinement type is proposed,considering the sequential installation of the fully grouted rockbolts and the surface support.The rock mass is assumed to be elastic-brittle-plastic material,obeying the linear Mohr-Coulomb criterion or the non-linear Hoek-Brown criterion.According to the strain states of the tunnel wall at bolt and surface support installation and the relative magnitude between the bolt length and the plastic depth during the whole process,six cases are categorized upon solving the problem.Each case is divided into three stages due to the different effects of the active rockbolts and the passive surface support.The fictitious pressure is introduced to quantify the threedimensional(3D)effect of the tunnel face,and thus,the actual physical location along the tunnel axis of the analytical section can be considered.By using the bolt-rock strain compatibility and the rocksurface support displacement compatibility conditions,the solutions of longitudinal tunnel displacement and the reaction pressure of surface support along the tunnel axis are obtained.The proposed analytical solutions are validated by a series of 3D numerical simulations.Extensive parametric studies are conducted to examine the effect of the typical parameters of rockbolts and surface support on the tunnel displacement and the reaction pressure of the surface support under different rock conditions.The results show that the rockbolts are more effective in controlling the tunnel displacement than the surface support,which should be installed as soon as possible with a suitable length.For tunnels excavated in weak rocks or with restricted displacement control requirements,the surface support should also be installed or closed timely with a certain stiffness.The proposed method provides a convenient alternative approach for the optimization of rockbolts and surface support at the preliminary stage of tunnel design.展开更多
基金supported by the the National Science and Technology Council(Grant Number:NSTC 112-2221-E239-022).
文摘To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
文摘This paper discusses a fictitious domain method for the linear Dirichlet problem and its applications to the generalized Stokes problem. This method treats Dirichlet boundary condit ion via a Lagrange multiplier technique and is well suited to the no-slip bound ary condition in viscous flow problems. In order to improve the accuracy of solu tions, meshes are refined according to the a posteriori error estimate. The mini -element discretization is applied to solve the generalized Stokes problem. Fin ally, some numerical results to validate this method are presented for partial d ifferential equations with Dirichlet boundary condition.
基金Projects(51378464,51309207)supported by the National Natural Science Foundation of China
文摘The dynamic response of pile in layered soil is theoretically investigated when considering the transverse inertia effect.Firstly, the fictitious soil-pile model is employed to simulate the dynamic interaction between the pile and the soil layers beneath pile toe. The dynamic interactions of adjacent soil layers along the vertical direction are simplified as distributed Voigt models.Meanwhile, the pile and fictitious soil-pile are assumed to be viscoelastic Rayleigh-Love rods, and both the radial and vertical displacement continuity conditions at the soil-pile interface are taken into consideration. On this basis, the analytical solution for dynamic response at the pile head is derived in the frequency domain and the corresponding quasi-analytical solution in the time domain is then obtained by means of the convolution theorem. Following this, the accuracy and parameter value of the hypothetical boundaries for soil-layer interfaces are discussed. Comparisons with published solution and measured data are carried out to verify the rationality of the present solution. Parametric analyses are further conducted by using the present solution to investigate the relationships between the transverse inertia effects and soil-pile parameters.
基金Funded by the National Natural Science Foundation of China (No.40574004, No.40374004, No.40174004).
文摘The potential field determined based on the fictitious compress recovery approach is influenced by the errors contained in the boundary (the Earth's surface or the surface corresponding to the satellite altitude) values. Given the boundary value with definite accuracy, the accuracy of the field determined based on the fictitious compress recovery approach is estimated, and it is theoretically shown that the determined field has the same accuracy level as the given boundary value.
基金P.Sun was supported by NSF Grant DMS-1418806C.S.Zhang was partially supported by the National Key Research and Development Program of China(Grant No.2016YFB0201304)+1 种基金the Major Research Plan of National Natural Science Foundation of China(Grant Nos.91430215,91530323)the Key Research Program of Frontier Sciences of CAS.
文摘The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.
基金supported by the National Natural Science Foundation of China(51078150)the National Natural Science Foundation of China(11602087)+1 种基金the State Key Laboratory of Subtropical Building Science,South China University of Technology(2017ZB32)National Undergraduate Innovative and Entrepreneurial Training Program(201810561180).
文摘The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which can be solved by traditional numerical methods and an infinite domain with a crack which can be solved by Muskhelishvili’s solutions.However,this alternating method cannot be directly applied to the edge crack problems since partial crack surface of Muskhelishvili’s solutions is located outside the computational domain.In this paper,an improved alternating method,the spline fictitious boundary element alternating method(SFBEAM),based on infinite domain with the combination of spline fictitious boundary element method(SFBEM)and Muskhelishvili’s solutions is proposed to solve the edge crack problems.Since the SFBEM and Muskhelishvili’s solutions are obtained in the framework of infinite domain,no special treatment is needed for solving the problem of edge cracks.Different mixed boundary conditions edge crack problems with varies of computational parameters are given to certify the high precision,efficiency and applicability of the proposed method compared with other alternating methods and extend finite element method.
文摘In modern computer games, "bots" - intelligent realistic agents play a prominent role in the popularity of a game in the market. Typically, bots are modeled using finite-state machine and then programmed via simple conditional statements which are hard-coded in bots logic. Since these bots have become quite predictable to an experienced games' player, a player might lose interest in the game. We propose the use of a game theoretic based learning rule called fictitious play for improving behavior of these computer game bots which will make them less predictable and hence, more a enjoyable game.
文摘This paper investigates the stability and convergence of some knowndifference schemes for the numerical solution to heat conduction equation withderivative boundary conditions by the fictitious domain method.The discrete vari-ables at the false mesh points are firstly eliminated from the difference schemes andthe local truncation errors are then analyzed in detail.The stability and convergenceof the schemes are proved by energy method.An improvement is proposed to obtainbetter schemes over the original ones.Several numerical examples and comparisonswith other schemes are presented.
文摘This paper investigates some known difference schemes for the numerical solution to parabolic differential equation with derivative boundary conditions by the fictitious domain method.The stability and convergence in L ∞ are proven.
文摘Abstract In this article we briefly studied the fictitious domain methods for steady and nonsteady Navier-Stokes equations based on Penalty factorεon the extended domain. The convergence u<sup>?</sup>→u in H<sub>0</sub><sup>1</sup>(Ω)<sup>d</sup> and L<sup>2</sup>(Ω)<sup>d</sup>(d=2,3) is given as well as p<sup>?</sup>→p in L<sub>0</sub><sup>2</sup>(Ω).
基金Supported by National 863 Plan Project of Ministry of Science and Technology of China under Grant No. 2006AA09Z354National Natural Science Foundation of China under Grant No. 10672101.
文摘In this paper, an efficient multigrid fictitious boundary method (MFBM) coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.
基金This study is supported by the National Natural Science Foundation of China (Grant No50579046) the Science Foundation of Tianjin Municipal Commission of Science and Technology (Grant No043114711)
文摘This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as "fictitious" fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure' s rigid body shape and behavior, a rigid body constraint is enforced on the "fictitious" fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain (DLM/ FD) method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.
基金This research is funded by the Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan(Grant No.AP09058430)。
文摘To apply the fictitious domain method and conduct numericalexperiments, a boundary value problem for an ordinary differential equation is considered. The results of numerical calculations for different valuesof the iterative parameter τ and the small parameter ε are presented. Astudy of the auxiliary problem of the fictitious domain method for NavierStokes equations with continuation into a fictitious subdomain by highercoefficients with a small parameter is carried out. A generalized solutionof the auxiliary problem of the fictitious domain method with continuationby higher coefficients with a small parameter is determined. After all theabove mathematical studies, a computational algorithm has been developedfor the numerical solution of the problem. Two methods were used to solvethe problem numerically. The first variant is the fictitious domain methodassociated with the modification of nonlinear terms in a fictitious subdomain.The model problem shows the effectiveness of using such a modification. Theproposed version of the method is used to solve two problems at once that arisewhile numerically solving systems of Navier-Stokes equations: the problem ofa curved boundary of an arbitrary domain and the problem of absence of aboundary condition for pressure in physical formulation of the internal flowproblem. The main advantage of this method is its universality in developmentof computer programs. The second method used calculation on a uniform gridinside the area. When numerically implementing the solution on a uniformgrid inside the domain, using this method it’s possible to accurately take intoaccount the boundaries of the curved domain and ensure the accuracy of thevalue of the function at the boundaries of the domain. Methodical calculationswere carried out, the results of numerical calculations were obtained. Whenconducting numerical experiments in both cases, quantitative and qualitativeindicators of numerical results coincide.
基金Supported bythe National Natural Science Foundation of China (No.40637034, No. 40574004), the National 863 Program of China (No. 2006AA12Z211).
文摘The fictitious compress recovery approach is introduced, which could be applied to the establishment of the Rungerarup theorem, the determination of the Bjerhammar's fictitious gravity anomaly, the solution of the "downward con- tinuation" problem of the gravity field, the confirmation of the convergence of the spherical harmonic expansion series of the Earth's potential field, and the gravity field determination in three cases: gravitational potential case, gravitation case, and gravitational gradient case. Several tests using simulation experiments show that the fictitious compress recovery approach shows promise in physical geodesy applications.
文摘Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the fictitious crack model in fracture mechanics of concrete. The new analytical element can be implemented into FEM program systems to solve fictitious crack propagation problems for concrete cracked plates with arbitrary shapes and loads. Numerical results indicate that the method is more efficient and accurate than ordinary finite element method.
基金Project supported by the National Key Research and Development Program of China(Grant No.2016YFA0302800)the National Natural Science Foundation of China(Grant No.11674361)
文摘We propose a novel on-chip platform for controlling and manipulating cold atoms precisely and coherently. The scheme is achieved by producing optically induced fictitious magnetic traps(OFMTs) with 790 nm(for -(87)Rb) circularly polarized laser beams and state-dependent potentials simultaneously for two internal atomic states with microwave coplanar waveguides. We carry out numerical calculations and simulations for controlled collisional interactions between OFMTs and addressable single atoms' manipulation on our designed hybrid atom chips. The results show that our proposed platform is feasible and flexible, which has wide applications including collisional dynamics investigation, entanglement generation,and scalable quantum gates implementation.
文摘As The Three-Body Problem translated by Ken Liu went popular in America, Chinese science fiction has got spotlightfrom all over the world. This article studies the translation of fictitious elements in The Three-Body Problem. Based on the specificmethods Ken Liu employs, general requirements for translators are proposed and strategies of translating fictitious elements areelicited. It intends to bring forth enlightening ideas in translating fictitious elements of science fiction and shed lights on this areaof study which few have set foot on.
基金This work was supported by the Key Project of Natural Science Foundation of Zhejiang Province(No.LXZ22E080001)National Natural Science Foundation of China(Grant Nos.51779217,52178358,and 52108349)China Scholarship Council(No.202006320262).We gratefully acknowledge this support.
文摘A fictitious soil pile(FSP)model is developed to simulate the behavior of pipe piles with soil plugs undergoing high-strain dynamic impact loading.The developed model simulates the base soil with a fictitious hollow pile fully filled with a soil plug extending at a cone angle from the pile toe to the bedrock.The friction on the outside and inside of the pile walls is distinguished using different shaft models,and the propagation of stress waves in the base soil and soil plug is considered.The motions of the pile—soil system are solved by discretizing them into spring-mass model based on the finite difference method.Comparisons of the predictions of the proposed model and conventional numerical models,as well as measurements for pipe piles in field tests subjected to impact loading,validate the accuracy of the proposed model.A parametric analysis is conducted to illustrate the influence of the model parameters on the pile dynamic response.Finally,the effective length of the FSP is proposed to approximate the affected soil zone below the pipe pile toe,and some guidance is provided for the selection of the model parameters.
基金funding support from the Fundamental Research Funds for the Central Universities(Grant No.2023JBZY024)the National Natural Science Foundation of China(Grant Nos.52208382 and 52278387).
文摘To investigate the interaction of the bolt-reinforced rock and the surface support,an analytical model of the convergence-confinement type is proposed,considering the sequential installation of the fully grouted rockbolts and the surface support.The rock mass is assumed to be elastic-brittle-plastic material,obeying the linear Mohr-Coulomb criterion or the non-linear Hoek-Brown criterion.According to the strain states of the tunnel wall at bolt and surface support installation and the relative magnitude between the bolt length and the plastic depth during the whole process,six cases are categorized upon solving the problem.Each case is divided into three stages due to the different effects of the active rockbolts and the passive surface support.The fictitious pressure is introduced to quantify the threedimensional(3D)effect of the tunnel face,and thus,the actual physical location along the tunnel axis of the analytical section can be considered.By using the bolt-rock strain compatibility and the rocksurface support displacement compatibility conditions,the solutions of longitudinal tunnel displacement and the reaction pressure of surface support along the tunnel axis are obtained.The proposed analytical solutions are validated by a series of 3D numerical simulations.Extensive parametric studies are conducted to examine the effect of the typical parameters of rockbolts and surface support on the tunnel displacement and the reaction pressure of the surface support under different rock conditions.The results show that the rockbolts are more effective in controlling the tunnel displacement than the surface support,which should be installed as soon as possible with a suitable length.For tunnels excavated in weak rocks or with restricted displacement control requirements,the surface support should also be installed or closed timely with a certain stiffness.The proposed method provides a convenient alternative approach for the optimization of rockbolts and surface support at the preliminary stage of tunnel design.