Search results

The purpose of this paper is to propose the generalized integral transform technique (GITT) to the solution of convection-diffusion problems with nonlinear boundary conditions by employing the corresponding nonlinear eigenvalue problem in the construction of the expansion basis.

The original nonlinear boundary condition coefficients in the problem formulation are all incorporated into the adopted eigenvalue problem, which may be itself integral transformed through a representative linear auxiliary problem, yielding a nonlinear algebraic eigenvalue problem for the associated eigenvalues and eigenvectors, to be solved along with the transformed ordinary differential system. The nonlinear eigenvalues computation may also be accomplished by rewriting the corresponding transcendental equation as an ordinary differential system for the eigenvalues, which is then simultaneously solved with the transformed potentials.

An application on one-dimensional transient diffusion with nonlinear boundary condition coefficients is selected for illustrating some important computational aspects and the convergence behavior of the proposed eigenfunction expansions. For comparison purposes, an alternative solution with a linear eigenvalue problem basis is also presented and implemented.

This novel approach can be further extended to various classes of nonlinear convection-diffusion problems, either already solved by the GITT with a linear coefficients basis, or new challenging applications with more involved nonlinearities.

To introduce an efficient methodology for the computation of the periodic steady state solution of power systems with nonlinear and time‐varying components which combines a Newton method based on a numerical differentiation procedure to obtain a fast steady state solution in the time domain and parallel process techniques.

Nonlinear electric systems are represented by a set of differential equations, the conventional solution in the time domain is accelerated by a Newton method based on a numerical differentiation procedure for the convergence of state variables to the limit cycle and thus to the network periodic steady state solution. The efficiency of the solution is further enhanced with the application of parallel processing technology based on parallel virtual machine (PVM) and multi‐threading (MT).

The periodic steady state solution of nonlinear electric systems, even of large‐scale, can be efficiently obtained in the time domain with the application of Newton methods for the fast converge of state variables to the limit cycle. The efficiency of the computer solution can be dramatically enhanced with the application of parallel processing technology. The potential of the PVM and MT platforms is shown in the investigation. A comparison of advantages and disadvantages associated with each parallel processing platforms is given; a quantitative comparison between PVM and MT is provided.

The steady state solution of nonlinear electric systems can be efficiently obtained with a combination of Newton methods for the convergence acceleration to the limit cycle and parallel processing techniques.

The steady state solution of nonlinear electric systems using a Newton method based on a numerical differentiation procedure for the convergence acceleration to the limit cycle and parallel processing based on the PVM and MT platforms has not, to the authors' knowledge, reported before.

The purpose of this paper is to present a comparison of nonlinear and linear simulations of aircraft dynamics to determine the divergence of the linear solution from the nonlinear solution.

The general equations of motion of a transport aircraft are presented both in nonlinear and linear form. The nonlinear equations are solved by using the Runge Kutta method. Linear equations are solved numerically by using the Runge Kutta method and they are also solved exactly by using the Laplace transformation method. All of these solutions are obtained by using the body axis system. The results of the simulations are plotted for different control deflections.

Solution of linear equations by both methods gave the same results as expected. There are important differences in amplitude and frequency of oscillations which are obtained by using nonlinear and linear equations. These differences increase with growing input control deflection. Therefore, it is appropriate to prefer nonlinear approach to obtain more satisfactory results.

Accurate determination of the aerodynamic derivatives is important for the accuracy of the nonlinear solutions.

Many classical approaches use stability axis system for the solution of linear equations. However, in this paper transfer functions of the aircraft are redefined in the body axis system, because stability axes change with angle of attack and some of the stability derivatives need to be re‐evaluated for each angle of attack. Moreover, in addition to classical text book, linear equations are also solved by using the 4th order Runge Kutta medhod.

This paper aims to propose a new adaptive strategy for two-dimensional (2D) nonlinear finite element (FE) analysis of the minimal surface problem (MSP) based on the element energy projection (EEP) technique.

By linearizing nonlinear problems into a series of linear problems via the Newton method, the EEP technique, which is an effective and reliable point-wise super-convergent displacement recovery strategy for linear FE analysis, can be directly incorporated into the solution procedure. Accordingly, a posteriori error estimate in maximum norm was established and an adaptive 2D nonlinear FE strategy of h-version mesh refinement was developed.

Three classical known surfaces, including a singularity problem, were analysed. Moreover, an example whose analytic solution is unavailable was considered and a comparison was made between present results and those computed by the MATLAB PDE toolbox. The results show that the adaptively-generated meshes reflect the difficulties inherent in the problems and the proposed adaptive analysis can produce FE solutions satisfying the user-preset error tolerance in maximum norm with a fair adaptive convergence rate.

The EEP technique for linear FE analysis was extended to the nonlinear procedure of MSP and can be expected to apply to other 2D nonlinear problems. The employment of the maximum norm makes point-wisely error control on the sought surfaces possible and makes the proposed method distinguished from other adaptive FE analyses.

The purpose of this paper is to apply, for the first time, the authors' newly developed perturbation iteration method to heat transfer problems. The effectiveness of the new method in nonlinear heat transfer problems will be tested.

Nonlinear heat transfer problems are solved by perturbation iteration method. They are also solved by the well‐known technique variational iteration method in the literature.

It is found that perturbation iteration solutions converge faster to the numerical solutions. More accurate results can be achieved with this new method for nonlinear heat transfer problems.

A few iterations are actually sufficient. Further iterations need symbolic packages to calculate the solutions.

This new technique can practically be applied to many heat and flow problems.

The new perturbation iteration technique is successfully implemented to nonlinear heat transfer problems. Results show good agreement with the direct numerical simulations and the method performs better than the existing variational iteration method.

The purpose of this study is to develop a reliable finite element algorithm based on the transmission line method (TLM) to solve the nonlinear problem in electromagnetic field calculation.

In this paper, the TLM has been researched and applied to solve nonlinear iteration in FEM. LU decomposition method and the Jacobi preconditioned conjugate gradient method have been investigated to solve the equations in transmission line finite element method (FEM-TLM). The algorithms have been developed in C++ language. The algorithm is applied to analyze the magnetic field of a long straight current-carrying wire and a single-phase transformer.

FEM-TLM is more effective than traditional FEM in nonlinear iteration. The results of FEM-TLM have been compared and analyzed under different calculation scales. It is found that the LU decomposition method is more suitable for FEM-TLM because there is no need to repeatedly assemble the global coefficient matrix in the iterative solution process and it is not affected by the disturbance of the right-hand vector.

An effective algorithm is provided for solving nonlinear problems in the electromagnetic field, which can save a lot of computing costs. The efficiency of LU decomposition and CG method in FEM-TLM nonlinear iteration is investigated, which also makes a preliminary exploration for the research of FEM-TLM parallel algorithms.

In this study, a second-order standard wave equation extended to a two-dimensional viscous wave equation with timely differentiated advection-diffusion terms has been solved by differential quadrature methods (DQM) using a modification of cubic B-spline functions. Two numerical schemes are proposed and compared to achieve numerical approximations for the solutions of nonlinear viscous wave equations.

Two schemes are adopted to reduce the given system into two systems of nonlinear first-order partial differential equations (PDE). For each scheme, modified cubic B-spline (MCB)-DQM is used for calculating the spatial variables and their derivatives that reduces the system of PDEs into a system of nonlinear ODEs. The solutions of these systems of ODEs are determined by SSP-RK43 scheme. The CPU time is also calculated and compared. Matrix stability analysis has been performed for each scheme and both are found to be unconditionally stable. The results of both schemes have been extensively discussed and compared. The accuracy and reliability of the methods have been successfully tested on several examples.

A comparative study has been carried out for two different schemes. Results from both schemes are also compared with analytical solutions and the results available in literature. Experiments show that MCB-DQM with Scheme II yield more accurate and reliable results in solving viscous wave equations. But Scheme I is comparatively less expensive in terms of CPU time. For MCB-DQM, less depository requirements lead to less aggregation of approximation errors which in turn enhances the correctness and readiness of the numerical techniques. Approximate solutions to the two-dimensional nonlinear viscous wave equation have been found without linearizing the equation. Ease of implementation and low computation cost are the strengths of the method.

For the first time, a comparative study has been carried out for the solution of nonlinear viscous wave equation. Comparisons are done in terms of accuracy and CPU time. It is concluded that Scheme II is more suitable.

This paper presents a Newton‐Raphson algorithm for determining the Fourier spectrum of two‐periodic solutions for dynamic systems described by nonlinear ordinary differential equations. Assuming that two basic frequencies are known, the coefficients of a double Fourier series result from this algorithm. An application to the analysis of electromagnetic phenomena in electromechanical converters is described. In an example, of the steady‐state performances of current in a simple converter, the algorithm is tested with very good results.

This paper aims to apply an iterative numerical method to find the numerical solution of the nonlinear non-self-adjoint singular boundary value problems that arises in the theory of epitaxial growth.

The proposed problem has multiple solutions and it is singular too; so not every technique can capture all the solutions. This study proposes to use variational iterative numerical method and compute both the solutions. The computed solutions are very close to the exact result.

It turns out that the existence or nonexistence of numerical solutions fully depends on the value of a parameter. The authors show that numerical solutions exist for small positive values of this parameter. For large positive values of the parameter, they find nonexistence of solutions. They also observe existence of solutions for negative values of the parameter and determine the range of parameter values which separates existence and nonexistence of solutions. This parameter has a clear physical meaning, as it describes the rate at which new material is deposited onto the system. This fact allows interpreting the physical significance of the results.

The authors could capture both the solutions and got accurate estimation of the parameter. This method will be a great tool to handle such types of nonlinear non-self-adjoint equations that have multiple solutions in engineering and mathematical sciences. The results in this paper will have an impact on the understanding of theoretical models of epitaxial growth in near future.

The purpose of this paper is to propose a calibration method that can calibrate the relationships among the robot manipulator, the camera and the workspace.

The method uses a laser pointer rigidly mounted on the manipulator and projects the laser beam on the work plane. Nonlinear constraints governing the relationships of the geometrical parameters and measurement data are derived. The uniqueness of the solution is guaranteed when the camera is calibrated in advance. As a result, a decoupled multi‐stage closed‐form solution can be derived based on parallel line constraints, line/plane intersection and projective geometry. The closed‐form solution can be further refined by nonlinear optimization which considers all parameters simultaneously in the nonlinear model.

Computer simulations and experimental tests using actual data confirm the effectiveness of the proposed calibration method and illustrate its ability to work even when the eye cannot see the hand.

Only a laser pointer is required for this calibration method and this method can work without any manual measurement. In addition, this method can also be applied when the robot is not within the camera field of view.