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.展开更多
How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linea...How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linearizing technique based on the nonlinear term to reduce the effect of the nonlinear terms.We decompose the nonlinear terms in the NAEs through a splitting parameter and then linearize the NAEs around the values at the previous step to a linear system.Through the maximal orthogonal projection concept,to minimize a merit function within a selected interval of splitting parameters,the optimal parameters can be quickly determined.In each step,a linear system is solved by the Gaussian elimination method,and the whole iteration procedure is convergent very fast.Several numerical tests show the high performance of the optimal split-linearization iterative method(OSLIM).展开更多
In this paper,the equal-norm multiple-scale Trefftz method combined with the implicit Lie-group scheme is applied to solve the two-dimensional nonlinear sloshing problem with baffles.When considering solving sloshing ...In this paper,the equal-norm multiple-scale Trefftz method combined with the implicit Lie-group scheme is applied to solve the two-dimensional nonlinear sloshing problem with baffles.When considering solving sloshing problems with baffles by using boundary integral methods,degenerate geometry and problems of numerical instability are inevitable.To avoid numerical instability,the multiple-scale characteristic lengths are introduced into T-complete basis functions to efficiently govern the high-order oscillation disturbance.Again,the numerical noise propagation at each time step is eliminated by the vector regularization method and the group-preserving scheme.A weighting factor of the group-preserving scheme is introduced into a linear system and then used in the initial and boundary value problems(IBVPs)at each time step.More importantly,the parameters of the algorithm,namely,the T-complete function,dissipation factor,and time step,can obtain a linear relationship.The boundary noise interference and energy conservation are successfully overcome,and the accuracy of the boundary value problem is also improved.Finally,benchmark cases are used to verify the correctness of the numerical algorithm.The numerical results show that this algorithm is efficient and stable for nonlinear two-dimensional sloshing problems with baffles.展开更多
The inverse problems of wave equation to recover unknown space-time dependent functions of wave speed and wave source are solved in this paper, without needing of initial conditions and no internal measurement of data...The inverse problems of wave equation to recover unknown space-time dependent functions of wave speed and wave source are solved in this paper, without needing of initial conditions and no internal measurement of data being required. After a homogenization technique, a sequence of spatial boundary functions at least the fourth-order polynomials are derived, which satisfy the homogeneous boundary conditions. The boundary functions and the zero element constitute a linear space, and then a new boundary functional is proved in the linear space, of which the energy is preserved for each dynamic energetic boundary function. The linear systems and iterative algorithms used to recover unknown wave speed and wave source functions with the dynamic energetic boundary functions as bases are developed, which converge fast at each time step. The input data are parsimonious, merely the measured boundary strains and the boundary values and slopes of unknown functions to be recovered. The accuracy and robustness of present methods are confirmed by comparing exact solutions with estimated results under large noises up to 20%.展开更多
基金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.
基金support provided by the Ministry of Science and Technology,Taiwan,ROC under Contract No.MOST 110-2221-E-019-044.
文摘How to accelerate the convergence speed and avoid computing the inversion of a Jacobian matrix is important in the solution of nonlinear algebraic equations(NAEs).This paper develops an approach with a splitting-linearizing technique based on the nonlinear term to reduce the effect of the nonlinear terms.We decompose the nonlinear terms in the NAEs through a splitting parameter and then linearize the NAEs around the values at the previous step to a linear system.Through the maximal orthogonal projection concept,to minimize a merit function within a selected interval of splitting parameters,the optimal parameters can be quickly determined.In each step,a linear system is solved by the Gaussian elimination method,and the whole iteration procedure is convergent very fast.Several numerical tests show the high performance of the optimal split-linearization iterative method(OSLIM).
基金The second author greatly appreciates the financial support provided by the Ministry of Science and Technology,Taiwan,ROC,under Contract No.MOST 108-2221-E-019-015.
文摘In this paper,the equal-norm multiple-scale Trefftz method combined with the implicit Lie-group scheme is applied to solve the two-dimensional nonlinear sloshing problem with baffles.When considering solving sloshing problems with baffles by using boundary integral methods,degenerate geometry and problems of numerical instability are inevitable.To avoid numerical instability,the multiple-scale characteristic lengths are introduced into T-complete basis functions to efficiently govern the high-order oscillation disturbance.Again,the numerical noise propagation at each time step is eliminated by the vector regularization method and the group-preserving scheme.A weighting factor of the group-preserving scheme is introduced into a linear system and then used in the initial and boundary value problems(IBVPs)at each time step.More importantly,the parameters of the algorithm,namely,the T-complete function,dissipation factor,and time step,can obtain a linear relationship.The boundary noise interference and energy conservation are successfully overcome,and the accuracy of the boundary value problem is also improved.Finally,benchmark cases are used to verify the correctness of the numerical algorithm.The numerical results show that this algorithm is efficient and stable for nonlinear two-dimensional sloshing problems with baffles.
文摘The inverse problems of wave equation to recover unknown space-time dependent functions of wave speed and wave source are solved in this paper, without needing of initial conditions and no internal measurement of data being required. After a homogenization technique, a sequence of spatial boundary functions at least the fourth-order polynomials are derived, which satisfy the homogeneous boundary conditions. The boundary functions and the zero element constitute a linear space, and then a new boundary functional is proved in the linear space, of which the energy is preserved for each dynamic energetic boundary function. The linear systems and iterative algorithms used to recover unknown wave speed and wave source functions with the dynamic energetic boundary functions as bases are developed, which converge fast at each time step. The input data are parsimonious, merely the measured boundary strains and the boundary values and slopes of unknown functions to be recovered. The accuracy and robustness of present methods are confirmed by comparing exact solutions with estimated results under large noises up to 20%.
