For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of calculus of variations, multiplicative commutation r...For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of 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 approximate solution.展开更多
After discretization by the finite volume method,the numerical solution of fractional diffusion equations leads to a linear system with the Toeplitz-like structure.The theoretical analysis gives sufficient conditions ...After discretization by the finite volume method,the numerical solution of fractional diffusion equations leads to a linear system with the Toeplitz-like structure.The theoretical analysis gives sufficient conditions to guarantee the positive-definite property of the discretized matrix.Moreover,we develop a class of positive-definite operator splitting iteration methods for the numerical solution of fractional diffusion equations,which is unconditionally convergent for any positive constant.Meanwhile,the iteration methods introduce a new preconditioner for Krylov subspace methods.Numerical experiments verify the convergence of the positive-definite operator splitting iteration methods and show the efficiency of the proposed preconditioner,compared with the existing approaches.展开更多
In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor depos...In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor deposition) processes, means the flow of species to a gas-phase, which are influenced by an electric field. Such a field we can model by wave equations. The main contributions are to improve the standard discretization schemes of each part of the coupling equation. So we discuss an improvement with implicit Runge- Kutta methods instead of the Yee’s algorithm. Further we balance the solver method between the Maxwell and Transport equation.展开更多
The operator-splitting methods for the mathematic model of one kind of oin reactions for the problem of groundwater are considered.Optimal error estimates in L 2 and H 1 norm are obtained for the approximation solut...The operator-splitting methods for the mathematic model of one kind of oin reactions for the problem of groundwater are considered.Optimal error estimates in L 2 and H 1 norm are obtained for the approximation solution.展开更多
The Operator Splitting method is applied to differential equations occurring as mathematical models in financial models. This paper provides various operator splitting methods to obtain an effective and accurate solut...The Operator Splitting method is applied to differential equations occurring as mathematical models in financial models. This paper provides various operator splitting methods to obtain an effective and accurate solution to the Black-Scholes equation with appropriate boundary conditions for a European option pricing problem. Finally brief comparisons of option prices are given by different models.展开更多
Let H;, H;, H;be real Hilbert spaces, let A : H;→ H;, B : H;→ H;be two bounded linear operators. The split equality common fixed point problem（SECFP） in the infinite-dimensional Hilbert spaces introduced by Moudaf...Let H;, H;, H;be real Hilbert spaces, let A : H;→ H;, B : H;→ H;be two bounded linear operators. The split equality common fixed point problem（SECFP） in the infinite-dimensional Hilbert spaces introduced by Moudafi（Alternating CQ-algorithm for convex feasibility and split fixed-point problems. Journal of Nonlinear and Convex Analysis）is to find x ∈ F（U）, y ∈ F（T） such that Ax = By,（1）where U : H;→ H;and T : H;→ H;are two nonlinear operators with nonempty fixed point sets F（U） = {x ∈ H;: Ux = x} and F（T） = {x ∈ H;: Tx = x}. Note that,by taking B = I and H;= H;in（1）, we recover the split fixed point problem originally introduced in Censor and Segal. Recently, Moudafi introduced alternating CQ-algorithms and simultaneous iterative algorithms with weak convergence for the SECFP（1） of firmly quasi-nonexpansive operators. In this paper, we introduce two viscosity iterative algorithms for the SECFP（1） governed by the general class of quasi-nonexpansive operators. We prove the strong convergence of algorithms. Our results improve and extend previously discussed related problems and algorithms.展开更多
This article deals with the design of energy efficient water utilization systems allowing operation split. Practical features such as operating flexibility and capital cost have made the number of sub operations an im...This article deals with the design of energy efficient water utilization systems allowing operation split. Practical features such as operating flexibility and capital cost have made the number of sub operations an important parameter of the problem. By treating the direct and indirect heat transfers separately, target freshwater and energy consumption as well as the operation split conditions are first obtained. Subsequently, a mixed integer non-linear programming （MINLP） model is established for the design of water network and the heat exchanger network （HEN）. The proposed systematic approach is limited to a single contaminant. Example from literature is used to illustrate the applicability of the approach.展开更多
Operability problem of dividing wall column (DWC) raised by vapor split was investigated by numerically analyzing four cases defined by different compositions of a three-component mixture. DWCs were firstly designed f...Operability problem of dividing wall column (DWC) raised by vapor split was investigated by numerically analyzing four cases defined by different compositions of a three-component mixture. DWCs were firstly designed for each case by optimizing the vapor split to the two sides of the dividing wall, and then their feasibilities and total annual costs in operation were evaluated against different vapor split ratios. The analysis on the operability of the DWC for four cases was made based on two scenarios: (1) vapor split is shifted by the vapor resistance difference between the column sections in the two sides of the dividing wall and (2) the feed composition is changed. It was demonstrated that the positioning of the dividing wall and the decision on the vapor split may affect significantly the operability of a DWC.展开更多
We propose two schemes for splitting single- and two-qubit states by using four-particle genuine entangled state as the quantum channel. After the sender performs Bell-basis （or three-partite GHZ- basis） measurement...We propose two schemes for splitting single- and two-qubit states by using four-particle genuine entangled state as the quantum channel. After the sender performs Bell-basis （or three-partite GHZ- basis） measurements on her particles, and the cooperators operate single-particle measurements on their particles, the state receiver can reconstruct the original state of the sender by applying the appropriate unitary operation. In particular, in the scheme for splitting two-qubit state, the receiver needs to introduce an auxiliary particle and carries out a C-NOT operation.展开更多
Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by t...Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.展开更多
Based on non-maximally entangled four-particle cluster states, we propose a new hierarchical information splitting protocol to probabilistically realize the quantum state sharing of an arbitrary unknown two-qubit stat...Based on non-maximally entangled four-particle cluster states, we propose a new hierarchical information splitting protocol to probabilistically realize the quantum state sharing of an arbitrary unknown two-qubit state. In this scheme, the sender transmits the two-qubit secret state to three agents who are divided into two grades with two Bell-state measurements,and broadcasts the measurement results via a classical channel. One agent is in the upper grade and two agents are in the lower grade. The agent in the upper grade only needs to cooperate with one of the other two agents to recover the secret state but both of the agents in the lower grade need help from all of the agents. Every agent who wants to recover the secret state needs to introduce two ancillary qubits and performs a positive operator-valued measurement（POVM） instead of the usual projective measurement. Moreover, due to the symmetry of the cluster state, we extend this protocol to multiparty agents.展开更多
The simplified four-qubit cluster state （i.e., （｜0000〉 ＋ ｜0011〉 ＋ ｜1100〉 -｜1111〉）/2） is explored for splitting an arbitrary single-qubit quantum information （QI）. Various feasible distributions of the ...The simplified four-qubit cluster state （i.e., （｜0000〉 ＋ ｜0011〉 ＋ ｜1100〉 -｜1111〉）/2） is explored for splitting an arbitrary single-qubit quantum information （QI）. Various feasible distributions of the four qubits among the Q,I sender and receivers for tri-splitting or hi-splitting are found out. For the distribution representations the corresponding splitting schemes and their LOCCs （local operation and classical communication） are presented amply while others are mentioned concisely.展开更多
A differentially weighted operator splitting Monte Carlo （DWOSMC） method is developed to solve com- plex aerosol dynamic processes by coupling the differentially weighted Monte Carlo method and the operator splittin...A differentially weighted operator splitting Monte Carlo （DWOSMC） method is developed to solve com- plex aerosol dynamic processes by coupling the differentially weighted Monte Carlo method and the operator splitting technique. This method is validated by analytical solutions and a sectional method in different aerosol dynamic processes. It is first validated in coagulation and condensation processes, and nucleation and coagulation processes, and then validated through simultaneous nucleation, coagulation, and condensation processes. The results show that the DWOSMC method is a computationally efficient and quantitatively accurate method for simulating complex aerosol dynamic processes.展开更多
An operator-splitting algorithm for the three-dimensional convection-diffusion equa- tion is presented.The flow region is discretized into tetrahedronal elements which are fixed in time. The transport equation is spli...An operator-splitting algorithm for the three-dimensional convection-diffusion equa- tion is presented.The flow region is discretized into tetrahedronal elements which are fixed in time. The transport equation is split into two successive initial value problems:a pure convection problem and a pure diffusion problem.For the pure convection problem,solutions are found by the method of characteristiCS.The solution algorithm involves tracing the characteristic lines backwards in time from a vertex of an element to an interior point.A cubic polynomial is used to interpolate the concentration and its derivatives within each element.For the diffusion problem,an explicit finite element algorithm is employed.Numerical examples are given which agree well with the analytical solutions.展开更多
This paper presents an Operator -Splitting Method (OSM) for the solution of the universal Reynolds equation. Jakoobsson-Floberg-Olsson (JFO) pressure conditions were incorporated for the study of cavitation in a liqui...This paper presents an Operator -Splitting Method (OSM) for the solution of the universal Reynolds equation. Jakoobsson-Floberg-Olsson (JFO) pressure conditions were incorporated for the study of cavitation in a liquid-lubricated journal bearings. Shear flow component of the oil film was first solved by a modified upwind finite difference method. The solution of the pressure gradient flow component was completed by the Galerkin finite element method. Present OSM solutions for a slider bearing are in agreement with Elord's results. OSM was then applied to herringbone grooved journal bearing in this work. The film pressure, cavitation areas, load capacity and attitude angle were obtained with JFO pressure conditions. The calculated load capacities are in agreement with Hirs's experimental data. A comparison of the present results and those predicted by the Reynolds pressure conditions shows some differences. The numerical results indicate that the load capacity and the critical mass of journal (linear stability indicator) are higher, and the attitude angle is lower than those predicted by Reynolds pressure conditions in cases of high eccentricities.展开更多
For separating some specific four component mixtures into four products, the four-product dividing wall column(FPDWC) with two partition walls can provide the same utility consumption with the extended Petlyuk configu...For separating some specific four component mixtures into four products, the four-product dividing wall column(FPDWC) with two partition walls can provide the same utility consumption with the extended Petlyuk configuration, although with structure simplicity. However, the reluctance to implement this kind of four product dividing wall column industrially also consists in the two uncontrollable vapor splits associated with it. The vapor split ratios are set at the design stage and might not be the optimal value for changed feed composition, thus minimum energy consumption could not be ensured. In the present work, a sequential iterative optimization approach was initially employed to determine the parameters of cost-effective FPDWC. Then the effect of maintaining the vapor split ratios at their nominal value on the energy penalty was investigated for the FPDWC with two partition walls, in case of feed composition disturbance. The result shows that no more than + 2% above the optimal energy requirements could be ensured for 20% feed composition disturbances, which is encouraging for industrial implementation.展开更多
The Douglas–Peaceman–Rachford–Varga operator splitting methods are a class ofefficient methods for finding a zero of the sum of two maximal monotone operatorsin a real Hilbert space;however, they are sometimes diff...The Douglas–Peaceman–Rachford–Varga operator splitting methods are a class ofefficient methods for finding a zero of the sum of two maximal monotone operatorsin a real Hilbert space;however, they are sometimes difficult or even impossible tosolve the subproblems exactly. In this paper, we suggest an inexact version in whichsome relative error criterion is discussed. The corresponding convergence propertiesare established, and some preliminary numerical experiments are reported to illustrateits efficiency.展开更多
This paper is concerned with efficient numerical methods for the advectiondiffusion equation in a heterogeneous porous medium containing fractures.A dimensionally reduced fracture model is considered,in which the frac...This paper is concerned with efficient numerical methods for the advectiondiffusion equation in a heterogeneous porous medium containing fractures.A dimensionally reduced fracture model is considered,in which the fracture is represented as an interface between subdomains and is assumed to have larger permeability than the surrounding area.We develop three global-in-time domain decomposition methods coupled with operator splitting for the reduced fracture model,where the advection and the diffusion are treated separately by different numerical schemes and with different time steps.Importantly,smaller time steps can be used in the fracture-interface than in the subdomains.The first two methods are based on the physical transmission conditions,while the third one is based on the optimized Schwarz waveform relaxation approach with Ventcel-Robin transmission conditions.A discrete space-time interface system is formulated for each method and is solved iteratively and globally in time.Numerical results for two-dimensional problems with various P′eclet numbers and different types of fracture are presented to illustrate and compare the convergence and accuracy in time of the proposed methods with local time stepping.展开更多
In this article, we discuss a numerical method for the computation of the minimal and maximal solutions of a steady scalar Eikonal equation. This method relies on a penalty treatment of the nonlinearity, a biharmonic ...In this article, we discuss a numerical method for the computation of the minimal and maximal solutions of a steady scalar Eikonal equation. This method relies on a penalty treatment of the nonlinearity, a biharmonic regularization of the resulting variational problem, and the time discretization by operator-splitting of an initial value problem associated with the Euler-Lagrange equations of the regularized variational problem. A low-order finite element discretization is advocated since it is well-suited to the low regularity of the solutions. Numerical experiments show that the method sketched above can capture efficiently the extremal solutions of various two-dimensional test problems and that it has also the ability of handling easily domains with curved boundaries.展开更多
基金the Major State Basic Research Program of China(19990328)NNSF of China(19871051,19972039) the Doctorate Foundation of the State Education Commission
文摘For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of 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 approximate solution.
基金This work was supported by the National Natural Science Foundation of China(No.11971354)The author Yi-Shu Du acknowledges the financial support from the China Scholarship Council(File No.201906260146).
文摘After discretization by the finite volume method,the numerical solution of fractional diffusion equations leads to a linear system with the Toeplitz-like structure.The theoretical analysis gives sufficient conditions to guarantee the positive-definite property of the discretized matrix.Moreover,we develop a class of positive-definite operator splitting iteration methods for the numerical solution of fractional diffusion equations,which is unconditionally convergent for any positive constant.Meanwhile,the iteration methods introduce a new preconditioner for Krylov subspace methods.Numerical experiments verify the convergence of the positive-definite operator splitting iteration methods and show the efficiency of the proposed preconditioner,compared with the existing approaches.
文摘In this article a new approach is considered for implementing operator splitting methods for transport problems, influenced by electric fields. Our motivation came to model PE-CVD (plasma-enhanced chemical vapor deposition) processes, means the flow of species to a gas-phase, which are influenced by an electric field. Such a field we can model by wave equations. The main contributions are to improve the standard discretization schemes of each part of the coupling equation. So we discuss an improvement with implicit Runge- Kutta methods instead of the Yee’s algorithm. Further we balance the solver method between the Maxwell and Transport equation.
基金the National Sciences Foundation of China and the Doctorial Program of Higher Edua-tion
文摘The operator-splitting methods for the mathematic model of one kind of oin reactions for the problem of groundwater are considered.Optimal error estimates in L 2 and H 1 norm are obtained for the approximation solution.
文摘The Operator Splitting method is applied to differential equations occurring as mathematical models in financial models. This paper provides various operator splitting methods to obtain an effective and accurate solution to the Black-Scholes equation with appropriate boundary conditions for a European option pricing problem. Finally brief comparisons of option prices are given by different models.
基金supported by National Natural Science Foundation of China(61503385)Fundamental Research Funds for the Central Universities of China(3122016L002)
文摘Let H;, H;, H;be real Hilbert spaces, let A : H;→ H;, B : H;→ H;be two bounded linear operators. The split equality common fixed point problem（SECFP） in the infinite-dimensional Hilbert spaces introduced by Moudafi（Alternating CQ-algorithm for convex feasibility and split fixed-point problems. Journal of Nonlinear and Convex Analysis）is to find x ∈ F（U）, y ∈ F（T） such that Ax = By,（1）where U : H;→ H;and T : H;→ H;are two nonlinear operators with nonempty fixed point sets F（U） = {x ∈ H;: Ux = x} and F（T） = {x ∈ H;: Tx = x}. Note that,by taking B = I and H;= H;in（1）, we recover the split fixed point problem originally introduced in Censor and Segal. Recently, Moudafi introduced alternating CQ-algorithms and simultaneous iterative algorithms with weak convergence for the SECFP（1） of firmly quasi-nonexpansive operators. In this paper, we introduce two viscosity iterative algorithms for the SECFP（1） governed by the general class of quasi-nonexpansive operators. We prove the strong convergence of algorithms. Our results improve and extend previously discussed related problems and algorithms.
基金Supported by the Major Project of National Natural Science Foundation of China （No.20409205） and National High Technology Research and Development Program of China （No.G20070040）.
文摘This article deals with the design of energy efficient water utilization systems allowing operation split. Practical features such as operating flexibility and capital cost have made the number of sub operations an important parameter of the problem. By treating the direct and indirect heat transfers separately, target freshwater and energy consumption as well as the operation split conditions are first obtained. Subsequently, a mixed integer non-linear programming （MINLP） model is established for the design of water network and the heat exchanger network （HEN）. The proposed systematic approach is limited to a single contaminant. Example from literature is used to illustrate the applicability of the approach.
基金Supported by the State Key Fundamental Research Program(2012CB720500)
文摘Operability problem of dividing wall column (DWC) raised by vapor split was investigated by numerically analyzing four cases defined by different compositions of a three-component mixture. DWCs were firstly designed for each case by optimizing the vapor split to the two sides of the dividing wall, and then their feasibilities and total annual costs in operation were evaluated against different vapor split ratios. The analysis on the operability of the DWC for four cases was made based on two scenarios: (1) vapor split is shifted by the vapor resistance difference between the column sections in the two sides of the dividing wall and (2) the feed composition is changed. It was demonstrated that the positioning of the dividing wall and the decision on the vapor split may affect significantly the operability of a DWC.
文摘We propose two schemes for splitting single- and two-qubit states by using four-particle genuine entangled state as the quantum channel. After the sender performs Bell-basis （or three-partite GHZ- basis） measurements on her particles, and the cooperators operate single-particle measurements on their particles, the state receiver can reconstruct the original state of the sender by applying the appropriate unitary operation. In particular, in the scheme for splitting two-qubit state, the receiver needs to introduce an auxiliary particle and carries out a C-NOT operation.
文摘Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.
基金Project supported by the National Natural Science Foundation of China(Grant No.61671087)
文摘Based on non-maximally entangled four-particle cluster states, we propose a new hierarchical information splitting protocol to probabilistically realize the quantum state sharing of an arbitrary unknown two-qubit state. In this scheme, the sender transmits the two-qubit secret state to three agents who are divided into two grades with two Bell-state measurements,and broadcasts the measurement results via a classical channel. One agent is in the upper grade and two agents are in the lower grade. The agent in the upper grade only needs to cooperate with one of the other two agents to recover the secret state but both of the agents in the lower grade need help from all of the agents. Every agent who wants to recover the secret state needs to introduce two ancillary qubits and performs a positive operator-valued measurement（POVM） instead of the usual projective measurement. Moreover, due to the symmetry of the cluster state, we extend this protocol to multiparty agents.
基金Supported by the Program for New Century Excellent Talents at the University of China under Grant No.NCET-06-0554the National Natural Science Foundation of China under Grant Nos.10975001,60677001,10747146,and 10874122+3 种基金the Science-technology Fund of Anhui Province for Outstanding Youth under Grant No.06042087the Key Fund of the Ministry of Education of China under Grant No.206063 the General Fund of the Educational Committee of Anhui Province under Grant No.2006KJ260Bthe Natural Science Foundation of Guangdong Province under Grant Nos.06300345 and 7007806
文摘The simplified four-qubit cluster state （i.e., （｜0000〉 ＋ ｜0011〉 ＋ ｜1100〉 -｜1111〉）/2） is explored for splitting an arbitrary single-qubit quantum information （QI）. Various feasible distributions of the four qubits among the Q,I sender and receivers for tri-splitting or hi-splitting are found out. For the distribution representations the corresponding splitting schemes and their LOCCs （local operation and classical communication） are presented amply while others are mentioned concisely.
文摘A differentially weighted operator splitting Monte Carlo （DWOSMC） method is developed to solve com- plex aerosol dynamic processes by coupling the differentially weighted Monte Carlo method and the operator splitting technique. This method is validated by analytical solutions and a sectional method in different aerosol dynamic processes. It is first validated in coagulation and condensation processes, and nucleation and coagulation processes, and then validated through simultaneous nucleation, coagulation, and condensation processes. The results show that the DWOSMC method is a computationally efficient and quantitatively accurate method for simulating complex aerosol dynamic processes.
文摘An operator-splitting algorithm for the three-dimensional convection-diffusion equa- tion is presented.The flow region is discretized into tetrahedronal elements which are fixed in time. The transport equation is split into two successive initial value problems:a pure convection problem and a pure diffusion problem.For the pure convection problem,solutions are found by the method of characteristiCS.The solution algorithm involves tracing the characteristic lines backwards in time from a vertex of an element to an interior point.A cubic polynomial is used to interpolate the concentration and its derivatives within each element.For the diffusion problem,an explicit finite element algorithm is employed.Numerical examples are given which agree well with the analytical solutions.
文摘This paper presents an Operator -Splitting Method (OSM) for the solution of the universal Reynolds equation. Jakoobsson-Floberg-Olsson (JFO) pressure conditions were incorporated for the study of cavitation in a liquid-lubricated journal bearings. Shear flow component of the oil film was first solved by a modified upwind finite difference method. The solution of the pressure gradient flow component was completed by the Galerkin finite element method. Present OSM solutions for a slider bearing are in agreement with Elord's results. OSM was then applied to herringbone grooved journal bearing in this work. The film pressure, cavitation areas, load capacity and attitude angle were obtained with JFO pressure conditions. The calculated load capacities are in agreement with Hirs's experimental data. A comparison of the present results and those predicted by the Reynolds pressure conditions shows some differences. The numerical results indicate that the load capacity and the critical mass of journal (linear stability indicator) are higher, and the attitude angle is lower than those predicted by Reynolds pressure conditions in cases of high eccentricities.
基金Supported by Open Research Project of State Key Laboratory of Chemical Engineering(SKL-ChE-16B06)Yangtze Scholars and Innovative Research Team in Chinese University(IRT-17R81)
文摘For separating some specific four component mixtures into four products, the four-product dividing wall column(FPDWC) with two partition walls can provide the same utility consumption with the extended Petlyuk configuration, although with structure simplicity. However, the reluctance to implement this kind of four product dividing wall column industrially also consists in the two uncontrollable vapor splits associated with it. The vapor split ratios are set at the design stage and might not be the optimal value for changed feed composition, thus minimum energy consumption could not be ensured. In the present work, a sequential iterative optimization approach was initially employed to determine the parameters of cost-effective FPDWC. Then the effect of maintaining the vapor split ratios at their nominal value on the energy penalty was investigated for the FPDWC with two partition walls, in case of feed composition disturbance. The result shows that no more than + 2% above the optimal energy requirements could be ensured for 20% feed composition disturbances, which is encouraging for industrial implementation.
基金This work was partially supported by the National Natural Science Foundations of China(Nos.11471102 and 11701150)the Key Basic Research Foundation of the Higher Education Institutions of Henan Province(No.16A110012).
文摘The Douglas–Peaceman–Rachford–Varga operator splitting methods are a class ofefficient methods for finding a zero of the sum of two maximal monotone operatorsin a real Hilbert space;however, they are sometimes difficult or even impossible tosolve the subproblems exactly. In this paper, we suggest an inexact version in whichsome relative error criterion is discussed. The corresponding convergence propertiesare established, and some preliminary numerical experiments are reported to illustrateits efficiency.
基金partially supported by the US National Science Foundation under grant number DMS-1912626.
文摘This paper is concerned with efficient numerical methods for the advectiondiffusion equation in a heterogeneous porous medium containing fractures.A dimensionally reduced fracture model is considered,in which the fracture is represented as an interface between subdomains and is assumed to have larger permeability than the surrounding area.We develop three global-in-time domain decomposition methods coupled with operator splitting for the reduced fracture model,where the advection and the diffusion are treated separately by different numerical schemes and with different time steps.Importantly,smaller time steps can be used in the fracture-interface than in the subdomains.The first two methods are based on the physical transmission conditions,while the third one is based on the optimized Schwarz waveform relaxation approach with Ventcel-Robin transmission conditions.A discrete space-time interface system is formulated for each method and is solved iteratively and globally in time.Numerical results for two-dimensional problems with various P′eclet numbers and different types of fracture are presented to illustrate and compare the convergence and accuracy in time of the proposed methods with local time stepping.
基金supported by the National Science Foundation(No.DMS-0913982)
文摘In this article, we discuss a numerical method for the computation of the minimal and maximal solutions of a steady scalar Eikonal equation. This method relies on a penalty treatment of the nonlinearity, a biharmonic regularization of the resulting variational problem, and the time discretization by operator-splitting of an initial value problem associated with the Euler-Lagrange equations of the regularized variational problem. A low-order finite element discretization is advocated since it is well-suited to the low regularity of the solutions. Numerical experiments show that the method sketched above can capture efficiently the extremal solutions of various two-dimensional test problems and that it has also the ability of handling easily domains with curved boundaries.
