In this paper,the topological space(PF_(MP)(X),T) based on prime MP-filters of a lattice FI-algebra X is constructed firstly and we proved that it is a compact T_0-space if X with condition(P).Secondly,we restricted T...In this paper,the topological space(PF_(MP)(X),T) based on prime MP-filters of a lattice FI-algebra X is constructed firstly and we proved that it is a compact T_0-space if X with condition(P).Secondly,we restricted T to the set of all maximal MP-filters MF_(MP)(X) of X and concluded that(PF_(MP)(X),T |_(PF_(MP)(X)) )is a compact T_2 space if X with conditions(P) and(S).展开更多
The purpose of this paper is to establish a formula of higher derivative byFaa di Bruno formula, and apply it to some known results to get some identities involvingcomplete Bell polynomials.
In the present paper a numerical method is developed to approximate the solution of two-dimensional Nonlinear Schrodinger equation in the presence of a sin- gular potential. The method leads to generalized Lyapunov-Sy...In the present paper a numerical method is developed to approximate the solution of two-dimensional Nonlinear Schrodinger equation in the presence of a sin- gular potential. The method leads to generalized Lyapunov-Sylvester algebraic opera- tors that are shown to be invertible using original topological and differential calculus issued methods. The numerical scheme is proved to be consistent, convergent and sta- ble using the Lyapunov criterion, lax equivalence theorem and the properties of the generalized Lyapunov-Sylvester operators.展开更多
This investigation describes the nanofluid flow in a non-Darcy porous medium between two stretching and rotating disks. A nanofluid comprises of nanoparticles of silver and copper. Water is used as a base fluid. Heat ...This investigation describes the nanofluid flow in a non-Darcy porous medium between two stretching and rotating disks. A nanofluid comprises of nanoparticles of silver and copper. Water is used as a base fluid. Heat is being transferred with thermal radiation and the Joule heating. A system of ordinary differential equations is obtained by appropriate transformations. Convergent series solutions are obtained. Effects of various parameters are analyzed for the velocity and temperature. Numerical values of the skin friction coefficient and the Nusselt number are tabulated and examined. It can be seen that the radial velocity is affected in the same manner with both porous and local inertial parameters. A skin friction coefficient depicts the same impact on both disks for both nanofluids with larger stretching parameters.展开更多
We develop a new moving-water equilibria preserving partial relaxation(PR)scheme for the two-dimensional(2-D)Saint-Venant systemof shallowwater equations.The new scheme is a 2-D generalization of the one-dimensional(1...We develop a new moving-water equilibria preserving partial relaxation(PR)scheme for the two-dimensional(2-D)Saint-Venant systemof shallowwater equations.The new scheme is a 2-D generalization of the one-dimensional(1-D)PR scheme recently proposed in[X.Liu,X.Chen,S.Jin,A.Kurganov,andH.Yu,SIAMJ.Sci.Comput.,42(2020),pp.A2206–A2229].Our scheme is based on the PR approximation,which is designed in two steps.First,the geometric source terms are incorporated into the discharge fluxes,which results in a hyperbolic system with global fluxes.Second,the discharge equations are relaxed so that the nonlinearity is moved into the stiff right-hand side of the four added auxiliary equation.The obtained PR system is then numerically integrated using a semi-discrete hybrid upwind/central-upwind finitevolume method combined with an efficient semi-implicit ODE solver.The new 2-D PR scheme inherits the main advantages of the 1-D PR scheme:(i)no special treatment of the geometric source terms is required,(ii)no nonlinear(cubic)equations should be solved to obtain the point values of the water depth out of the reconstructed equilibriumvariables.The performance of the proposed PR scheme is illustrated on a number of numerical examples,in which we demonstrate that the PR scheme not only capable of exactly preserving quasi 1-D moving-water steady states and accurately capturing their small perturbations,but can also handle genuinely 2-D steady states and their small perturbations in a non-oscillatory manner.展开更多
In rarefied gas flows,the spatial grid size could vary by several orders of magnitude in a single flow configuration(e.g.,inside the Knudsen layer it is at the order of mean free path of gas molecules,while in the bul...In rarefied gas flows,the spatial grid size could vary by several orders of magnitude in a single flow configuration(e.g.,inside the Knudsen layer it is at the order of mean free path of gas molecules,while in the bulk region it is at a much larger hydrodynamic scale).Therefore,efficient implicit numerical method is urgently needed for time-dependent problems.However,the integro-differential nature of gas kinetic equations poses a grand challenge,as the gain part of the collision operator is non-invertible.Hence an iterative solver is required in each time step,which usually takes a lot of iterations in the(near)continuum flow regime where the Knudsen number is small;worse still,the solution does not asymptotically preserve the fluid dynamic limit when the spatial cell size is not refined enough.Based on the general synthetic iteration scheme for steady-state solution of the Boltzmann equation,we propose two numerical schemes to push the multiscale simulation of unsteady rarefied gas flows to a new boundary,that is,the numerical solution not only converges within dozens of iterations in each time step,but also asymptotically preserves the Navier-Stokes-Fourier limit in the continuum flow regime,when the spatial grid is coarse,and the time step is large(e.g.,in simulating the extreme slow decay of two-dimensional Taylor vortex,the time step is even at the order of vortex decay time).The properties of fast convergence and asymptotic preserving of the proposed schemes are not only rigorously proven by the Fourier stability analysis for simplified gas kinetic models,but also demonstrated by several numerical examples for the gas kinetic models and the Boltzmann equation.展开更多
For the planar Z2-equivariant cubic systems having twoelementary focuses, the characterization of a bi-center problem and shortened expressions of the first six Liapunov constants are completely discussed. The necessa...For the planar Z2-equivariant cubic systems having twoelementary focuses, the characterization of a bi-center problem and shortened expressions of the first six Liapunov constants are completely discussed. The necessary and sufficient conditions for the existence of the bi-center are obtained. All possible first integrals are given. Under small Z2-equivariant cubic perturbations, the conclusion that there exist at most 12 small-amplitude limit cycles with the scheme (6 II 6) is proved.展开更多
In recent years,physical informed neural networks(PINNs)have been shown to be a powerful tool for solving PDEs empirically.However,numerical analysis of PINNs is still missing.In this paper,we prove the convergence ra...In recent years,physical informed neural networks(PINNs)have been shown to be a powerful tool for solving PDEs empirically.However,numerical analysis of PINNs is still missing.In this paper,we prove the convergence rate to PINNs for the second order elliptic equations with Dirichlet boundary condition,by establishing the upper bounds on the number of training samples,depth and width of the deep neural networks to achieve desired accuracy.The error of PINNs is decomposed into approximation error and statistical error,where the approximation error is given in C2 norm with ReLU^(3)networks(deep network with activation function max{0,x^(3)})and the statistical error is estimated by Rademacher complexity.We derive the bound on the Rademacher complexity of the non-Lipschitz composition of gradient norm with ReLU^(3)network,which is of immense independent interest.展开更多
In this paper,a new earthquake location method based on the waveform inversion is proposed.As is known to all,the waveform misfit function under the L2 measure is suffering from the cycle skipping problem.This leads t...In this paper,a new earthquake location method based on the waveform inversion is proposed.As is known to all,the waveform misfit function under the L2 measure is suffering from the cycle skipping problem.This leads to a very small convergence domain of the conventional waveform based earthquake location methods.In present study,by introducing and solving two simple sub-optimization problems,we greatly expand the convergence domain of the waveform based earthquake location method.According to a large number of numerical experiments,the new method expands the range of convergence by several tens of times.This allows us to locate the earthquake accurately even from some relatively bad initial values.展开更多
This paper is concerned with numerical solutions of the LDG method for 1D wave equations.Superconvergence and energy conserving properties have been studied.We first study the superconvergence phenomenon for linear pr...This paper is concerned with numerical solutions of the LDG method for 1D wave equations.Superconvergence and energy conserving properties have been studied.We first study the superconvergence phenomenon for linear problems when alternating fluxes are used.We prove that,under some proper initial discretization,the numerical trace of the LDG approximation at nodes,as well as the cell average,converge with an order 2k+1.In addition,we establish k+2-th order and k+1-th order superconvergence rates for the function value error and the derivative error at Radau points,respectively.As a byproduct,we prove that the LDG solution is superconvergent with an order k+2 towards the Radau projection of the exact solution.Numerical experiments demonstrate that in most cases,our error estimates are optimal,i.e.,the error bounds are sharp.In the second part,we propose a fully discrete numerical scheme that conserves the discrete energy.Due to the energy conserving property,after long time integration,our method still stays accurate when applied to nonlinear Klein-Gordon and Sine-Gordon equations.展开更多
In this paper,the author considers a class of complete noncompact Riemannian manifoldswhich satisfy certain conditions on Ricci curvature and volume comparison. It is shown thatany harmonic map with finite energy from...In this paper,the author considers a class of complete noncompact Riemannian manifoldswhich satisfy certain conditions on Ricci curvature and volume comparison. It is shown thatany harmonic map with finite energy from such a manifold M into a normal geodesic ball inanother manifold N must be asymptotically constant at the infinity of each large end of M. Arelated existence theorem for harmonic maps is established.展开更多
Understanding the near boundary acoustic oscillation of microbubbles is critical for the effective design of ultrasonic biomedical devices and surface cleaning technologies.Accordingly,this study investigates the thre...Understanding the near boundary acoustic oscillation of microbubbles is critical for the effective design of ultrasonic biomedical devices and surface cleaning technologies.Accordingly,this study investigates the three-dimensional microbubble oscillation between two curved rigid plates experiencing a planar acoustic field using boundary integral method(BIM).The numerical model is validated via comparison with the nonlinear oscillation of the bubble governed by the modified Rayleigh-Plesset equation and with the axisymmetric model for an acoustic microbubble in infinite fluid domain.Then,the influence of the wave direction and horizontal standoff distance(h)on the bubble dynamics(including jet velocity,jet direction,centroid movement,total energy,and Kelvin impulse)were evaluated.It was concluded that the jet velocity,the maximum radius and the total energy of the bubble are not significantly influenced by the wave direction,while the jet direction and the high-pressure region depend strongly on it.More importantly,it was found that the jet velocity and the high-pressure region around the jet in acoustic bubble are drastically larger than their counterparts in the gas bubble.展开更多
A new method of the reproducing kernel Hilbert space is applied to a twodimensional parabolic inverse source problem with the final overdetermination. The exact and approximate solutions are both obtained in a reprodu...A new method of the reproducing kernel Hilbert space is applied to a twodimensional parabolic inverse source problem with the final overdetermination. The exact and approximate solutions are both obtained in a reproducing kernel space. The approximate solution and its partial derivatives are proved to converge to the exact solution and its partial derivatives, respectively. A technique is proposed to improve some existing methods. Numerical results show that the method is of high precision, and confirm the robustness of our method for reconstructing source parameter.展开更多
Time-dependent fractional partial differential equations typically require huge amounts of memory and computational time,especially for long-time integration,which taxes computational resources heavily for high-dimens...Time-dependent fractional partial differential equations typically require huge amounts of memory and computational time,especially for long-time integration,which taxes computational resources heavily for high-dimensional problems.Here,we first analyze existing numerical methods of sum-of-exponentials for approximating the kernel function in constant-order fractional operators,and identify the current pitfalls of such methods.In order to overcome the pitfalls,an improved sum-of-exponentials is developed and verified.We also present several sumof-exponentials for the approximation of the kernel function in variable-order fractional operators.Subsequently,based on the sum-of-exponentials,we propose a unified framework for fast time-stepping methods for fractional integral and derivative operators of constant and variable orders.We test the fast method based on several benchmark problems,including fractional initial value problems,the time-fractional Allen-Cahn equation in two and three spatial dimensions,and the Schr¨odinger equation with nonreflecting boundary conditions,demonstrating the efficiency and robustness of the proposed method.The results show that the present fast method significantly reduces the storage and computational cost especially for long-time integration problems.展开更多
How to obtain an effective design is a major concern of scientific research. This topic always involves high-dimensional inputs with limited resources. The foldover is a quick and useful technique in construction of f...How to obtain an effective design is a major concern of scientific research. This topic always involves high-dimensional inputs with limited resources. The foldover is a quick and useful technique in construction of fractional designs, which typically releases aliased factors or interactions.This paper takes the wrap-around L_2-discrepancy as the optimality measure to assess the optimal three-level combined designs. New and efficient analytical expressions and lower bounds of the wraparound L_2-discrepancy for three-level combined designs are obtained. The new lower bound is useful and sharper than the existing lower bound. Using the new analytical expression and lower bound as the benchmarks, the authors may implement an effective algorithm for constructing optimal three-level combined designs.展开更多
In this paper,we propose a novel PDE-based model for the multi-phase segmentation problem by using a complex version of Cahn-Hilliard equations.Specifically,we modify the original complex system of Cahn-Hilliard equat...In this paper,we propose a novel PDE-based model for the multi-phase segmentation problem by using a complex version of Cahn-Hilliard equations.Specifically,we modify the original complex system of Cahn-Hilliard equations by adding the mean curvature term and the fitting term to the evolution of its real part,which helps to render a piecewisely constant function at the steady state.By applying the K-means method to this function,one could achieve the desired multiphase segmentation.To solve the proposed system of equations,a semi-implicit finite difference scheme is employed.Numerical experiments are presented to demonstrate the feasibility of the proposed model and compare our model with other related ones.展开更多
基金Supported by the NSF of China(10371106,60774073)
文摘In this paper,the topological space(PF_(MP)(X),T) based on prime MP-filters of a lattice FI-algebra X is constructed firstly and we proved that it is a compact T_0-space if X with condition(P).Secondly,we restricted T to the set of all maximal MP-filters MF_(MP)(X) of X and concluded that(PF_(MP)(X),T |_(PF_(MP)(X)) )is a compact T_2 space if X with conditions(P) and(S).
基金Supported by the National Natural Science Foundation of China(Grant No.11601543,No.11601216,11701257)Supported by the NSF of Henan Province under Grant(No.172102410069)+1 种基金Supported by the NSF of Education Bureau of Henan Province under Grant(No.16B110009,18A110025)Supported by the Youth Foundation of Luoyang Normal university under Grant(No.2013-QNJJ-001)
文摘The purpose of this paper is to establish a formula of higher derivative byFaa di Bruno formula, and apply it to some known results to get some identities involvingcomplete Bell polynomials.
文摘In the present paper a numerical method is developed to approximate the solution of two-dimensional Nonlinear Schrodinger equation in the presence of a sin- gular potential. The method leads to generalized Lyapunov-Sylvester algebraic opera- tors that are shown to be invertible using original topological and differential calculus issued methods. The numerical scheme is proved to be consistent, convergent and sta- ble using the Lyapunov criterion, lax equivalence theorem and the properties of the generalized Lyapunov-Sylvester operators.
文摘This investigation describes the nanofluid flow in a non-Darcy porous medium between two stretching and rotating disks. A nanofluid comprises of nanoparticles of silver and copper. Water is used as a base fluid. Heat is being transferred with thermal radiation and the Joule heating. A system of ordinary differential equations is obtained by appropriate transformations. Convergent series solutions are obtained. Effects of various parameters are analyzed for the velocity and temperature. Numerical values of the skin friction coefficient and the Nusselt number are tabulated and examined. It can be seen that the radial velocity is affected in the same manner with both porous and local inertial parameters. A skin friction coefficient depicts the same impact on both disks for both nanofluids with larger stretching parameters.
基金The work of A.Kurganov was supported in part by NSFC grants 12111530004 and 12171226by the fund of the Guangdong Provincial Key Laboratory of Computational Science and Material Design(No.2019B030301001).
文摘We develop a new moving-water equilibria preserving partial relaxation(PR)scheme for the two-dimensional(2-D)Saint-Venant systemof shallowwater equations.The new scheme is a 2-D generalization of the one-dimensional(1-D)PR scheme recently proposed in[X.Liu,X.Chen,S.Jin,A.Kurganov,andH.Yu,SIAMJ.Sci.Comput.,42(2020),pp.A2206–A2229].Our scheme is based on the PR approximation,which is designed in two steps.First,the geometric source terms are incorporated into the discharge fluxes,which results in a hyperbolic system with global fluxes.Second,the discharge equations are relaxed so that the nonlinearity is moved into the stiff right-hand side of the four added auxiliary equation.The obtained PR system is then numerically integrated using a semi-discrete hybrid upwind/central-upwind finitevolume method combined with an efficient semi-implicit ODE solver.The new 2-D PR scheme inherits the main advantages of the 1-D PR scheme:(i)no special treatment of the geometric source terms is required,(ii)no nonlinear(cubic)equations should be solved to obtain the point values of the water depth out of the reconstructed equilibriumvariables.The performance of the proposed PR scheme is illustrated on a number of numerical examples,in which we demonstrate that the PR scheme not only capable of exactly preserving quasi 1-D moving-water steady states and accurately capturing their small perturbations,but can also handle genuinely 2-D steady states and their small perturbations in a non-oscillatory manner.
基金supported by the National Natural Science Foundation of China(12172162)the Guangdong-Hong Kong-Macao Joint Laboratory for Data-Driven Fluid Mechanics and Engineering Applications in China(2020B1212030001).
文摘In rarefied gas flows,the spatial grid size could vary by several orders of magnitude in a single flow configuration(e.g.,inside the Knudsen layer it is at the order of mean free path of gas molecules,while in the bulk region it is at a much larger hydrodynamic scale).Therefore,efficient implicit numerical method is urgently needed for time-dependent problems.However,the integro-differential nature of gas kinetic equations poses a grand challenge,as the gain part of the collision operator is non-invertible.Hence an iterative solver is required in each time step,which usually takes a lot of iterations in the(near)continuum flow regime where the Knudsen number is small;worse still,the solution does not asymptotically preserve the fluid dynamic limit when the spatial cell size is not refined enough.Based on the general synthetic iteration scheme for steady-state solution of the Boltzmann equation,we propose two numerical schemes to push the multiscale simulation of unsteady rarefied gas flows to a new boundary,that is,the numerical solution not only converges within dozens of iterations in each time step,but also asymptotically preserves the Navier-Stokes-Fourier limit in the continuum flow regime,when the spatial grid is coarse,and the time step is large(e.g.,in simulating the extreme slow decay of two-dimensional Taylor vortex,the time step is even at the order of vortex decay time).The properties of fast convergence and asymptotic preserving of the proposed schemes are not only rigorously proven by the Fourier stability analysis for simplified gas kinetic models,but also demonstrated by several numerical examples for the gas kinetic models and the Boltzmann equation.
基金Supported by National Natural Science Foundation of China (Grant Nos. 10831003 and 10771196)
文摘For the planar Z2-equivariant cubic systems having twoelementary focuses, the characterization of a bi-center problem and shortened expressions of the first six Liapunov constants are completely discussed. The necessary and sufficient conditions for the existence of the bi-center are obtained. All possible first integrals are given. Under small Z2-equivariant cubic perturbations, the conclusion that there exist at most 12 small-amplitude limit cycles with the scheme (6 II 6) is proved.
基金supported by the National Key Research and Development Program of China(No.2020YFA0714200)the National Science Foundation of China(No.12125103,No.12071362,No.11971468,No.11871474,No.11871385)+1 种基金the Natural Science Foundation of Hubei Province(No.2021AAA010,No.2019CFA007)the Fundamental Research Funds for the Central Universities.
文摘In recent years,physical informed neural networks(PINNs)have been shown to be a powerful tool for solving PDEs empirically.However,numerical analysis of PINNs is still missing.In this paper,we prove the convergence rate to PINNs for the second order elliptic equations with Dirichlet boundary condition,by establishing the upper bounds on the number of training samples,depth and width of the deep neural networks to achieve desired accuracy.The error of PINNs is decomposed into approximation error and statistical error,where the approximation error is given in C2 norm with ReLU^(3)networks(deep network with activation function max{0,x^(3)})and the statistical error is estimated by Rademacher complexity.We derive the bound on the Rademacher complexity of the non-Lipschitz composition of gradient norm with ReLU^(3)network,which is of immense independent interest.
基金This work was supported by the National Nature Science Foundation of China(Grant Nos.41230210,41390452)Hao Wu was also partially supported by the National Nature Science Foundation of China(Grant Nos.11101236,91330203)and SRF for ROCS,SEM.The authors are grateful to Prof.Shi Jin for his helpful suggestions and discussions that greatly improve the presentation.Hao Wu would like to thank Prof.Ping Tong for his valuable comments.The authors would also like to thank the referees for their valuable suggestions which helped to improve the content and presentation of this paper.
文摘In this paper,a new earthquake location method based on the waveform inversion is proposed.As is known to all,the waveform misfit function under the L2 measure is suffering from the cycle skipping problem.This leads to a very small convergence domain of the conventional waveform based earthquake location methods.In present study,by introducing and solving two simple sub-optimization problems,we greatly expand the convergence domain of the waveform based earthquake location method.According to a large number of numerical experiments,the new method expands the range of convergence by several tens of times.This allows us to locate the earthquake accurately even from some relatively bad initial values.
基金This work is supported in part by the National Natural Science Foundation of China(NSFC)under grants Nos.11201161,11471031,11501026,91430216,U1530401China Postdoctoral Science Foundation under grant Nos.2015M570026,2016T90027the US National Science Foundation(NSF)through grant DMS-1419040。
文摘This paper is concerned with numerical solutions of the LDG method for 1D wave equations.Superconvergence and energy conserving properties have been studied.We first study the superconvergence phenomenon for linear problems when alternating fluxes are used.We prove that,under some proper initial discretization,the numerical trace of the LDG approximation at nodes,as well as the cell average,converge with an order 2k+1.In addition,we establish k+2-th order and k+1-th order superconvergence rates for the function value error and the derivative error at Radau points,respectively.As a byproduct,we prove that the LDG solution is superconvergent with an order k+2 towards the Radau projection of the exact solution.Numerical experiments demonstrate that in most cases,our error estimates are optimal,i.e.,the error bounds are sharp.In the second part,we propose a fully discrete numerical scheme that conserves the discrete energy.Due to the energy conserving property,after long time integration,our method still stays accurate when applied to nonlinear Klein-Gordon and Sine-Gordon equations.
文摘In this paper,the author considers a class of complete noncompact Riemannian manifoldswhich satisfy certain conditions on Ricci curvature and volume comparison. It is shown thatany harmonic map with finite energy from such a manifold M into a normal geodesic ball inanother manifold N must be asymptotically constant at the infinity of each large end of M. Arelated existence theorem for harmonic maps is established.
文摘Understanding the near boundary acoustic oscillation of microbubbles is critical for the effective design of ultrasonic biomedical devices and surface cleaning technologies.Accordingly,this study investigates the three-dimensional microbubble oscillation between two curved rigid plates experiencing a planar acoustic field using boundary integral method(BIM).The numerical model is validated via comparison with the nonlinear oscillation of the bubble governed by the modified Rayleigh-Plesset equation and with the axisymmetric model for an acoustic microbubble in infinite fluid domain.Then,the influence of the wave direction and horizontal standoff distance(h)on the bubble dynamics(including jet velocity,jet direction,centroid movement,total energy,and Kelvin impulse)were evaluated.It was concluded that the jet velocity,the maximum radius and the total energy of the bubble are not significantly influenced by the wave direction,while the jet direction and the high-pressure region depend strongly on it.More importantly,it was found that the jet velocity and the high-pressure region around the jet in acoustic bubble are drastically larger than their counterparts in the gas bubble.
基金supported by the National Natural Science Foundation of China(No.91230119)
文摘A new method of the reproducing kernel Hilbert space is applied to a twodimensional parabolic inverse source problem with the final overdetermination. The exact and approximate solutions are both obtained in a reproducing kernel space. The approximate solution and its partial derivatives are proved to converge to the exact solution and its partial derivatives, respectively. A technique is proposed to improve some existing methods. Numerical results show that the method is of high precision, and confirm the robustness of our method for reconstructing source parameter.
基金supported by the NSF of China(Nos.12171283,12071301,12120101001)the National Key R&D Program of China(2021YFA1000202)+2 种基金the startup fund from Shandong University(No.11140082063130)the Shanghai Municipal Science and Technology Commission(No.20JC1412500)the science challenge project(No.TZ2018001).
文摘Time-dependent fractional partial differential equations typically require huge amounts of memory and computational time,especially for long-time integration,which taxes computational resources heavily for high-dimensional problems.Here,we first analyze existing numerical methods of sum-of-exponentials for approximating the kernel function in constant-order fractional operators,and identify the current pitfalls of such methods.In order to overcome the pitfalls,an improved sum-of-exponentials is developed and verified.We also present several sumof-exponentials for the approximation of the kernel function in variable-order fractional operators.Subsequently,based on the sum-of-exponentials,we propose a unified framework for fast time-stepping methods for fractional integral and derivative operators of constant and variable orders.We test the fast method based on several benchmark problems,including fractional initial value problems,the time-fractional Allen-Cahn equation in two and three spatial dimensions,and the Schr¨odinger equation with nonreflecting boundary conditions,demonstrating the efficiency and robustness of the proposed method.The results show that the present fast method significantly reduces the storage and computational cost especially for long-time integration problems.
基金supported by the National Natural Science Foundation of China under Grant Nos.11271147,11471135,11471136the UIC Grant R201409+1 种基金the Zhuhai Premier Discipline Grantthe Self-Determined Research Funds of CCNU from the Colleges Basic Research and Operation of MOE under Grant Nos.CCNU14A05041,CCNU16A02012
文摘How to obtain an effective design is a major concern of scientific research. This topic always involves high-dimensional inputs with limited resources. The foldover is a quick and useful technique in construction of fractional designs, which typically releases aliased factors or interactions.This paper takes the wrap-around L_2-discrepancy as the optimality measure to assess the optimal three-level combined designs. New and efficient analytical expressions and lower bounds of the wraparound L_2-discrepancy for three-level combined designs are obtained. The new lower bound is useful and sharper than the existing lower bound. Using the new analytical expression and lower bound as the benchmarks, the authors may implement an effective algorithm for constructing optimal three-level combined designs.
基金supported by the NSFC(Projects No.12025104,11871298,81930119).
文摘In this paper,we propose a novel PDE-based model for the multi-phase segmentation problem by using a complex version of Cahn-Hilliard equations.Specifically,we modify the original complex system of Cahn-Hilliard equations by adding the mean curvature term and the fitting term to the evolution of its real part,which helps to render a piecewisely constant function at the steady state.By applying the K-means method to this function,one could achieve the desired multiphase segmentation.To solve the proposed system of equations,a semi-implicit finite difference scheme is employed.Numerical experiments are presented to demonstrate the feasibility of the proposed model and compare our model with other related ones.