The numerical methods of Fourier eigen transform FET and its inversion are discussed and applied to the boundary element method for elastodynamics. The program for solving elastodynamic problems with the boundary elem...The numerical methods of Fourier eigen transform FET and its inversion are discussed and applied to the boundary element method for elastodynamics. The program for solving elastodynamic problems with the boundary element method is developed and some examples are given. From the numerical results of the examples, we know the method can increase the computing speed 5 similar to 10 times and the accuracy is guaranteed.展开更多
For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The...For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The method can be used to remove the transverse vibration due to the accumulation of computational errors. A combined boundary condition of sponge layer and Sommerfeld radiation condition is used to remove the reflected waves from the computing domain. By following the water particle on the water surface, the third order Stokes wave transform is simulated by the numerical wave flume technique. The computed results are in good agreement with theoretical ones.展开更多
The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic diff...The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.展开更多
This paper focuses on the finite element method in the complex frequency domain(CFD-FEM)for the transient electric field.First,the initial value and boundary value problem of the transient electric field under the ele...This paper focuses on the finite element method in the complex frequency domain(CFD-FEM)for the transient electric field.First,the initial value and boundary value problem of the transient electric field under the electroquasistatic field in the complex frequency domain is given.In addition,the finite element equation and the constrained electric field equation on the boundary are derived.Secondly,the indirect algorithm of the numerical inverse Laplace transform is introduced.Based on it,the calculation procedures of the CFD-FEM are illustrated in detail.Thirdly,the step response,zero-state response under the positive periodic square waveform(PPSW)voltage,and the zero-input response by the CFD-FEM with direct algorithm and indirect algorithm are compared.Finally,the reason for the numerical oscillations of the zero-state response under the PPSW voltage is analyzed,and the method to reduce oscillations is proposed.The results show that the numerical accuracy of the indirect algorithm of the CFD-FEM is more than an order of magnitude higher than that of the direct algorithm when calculating the step response of the transient electric field.The proposed method can significantly reduce the numerical oscillations of the zero-state response under the PPSW voltage.The proposed method is helpful for the calculation of the transient electric field,especially in the case of frequency-dependent parameters.展开更多
In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier ser...In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.展开更多
The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology ...The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology in this work is suitable for modeling diffusion problems with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The underpinning of the modeling approach is to decompose the global initial-boundary value problem into a steady-state component and a transient component. The steady-state component is governed by the Laplace PDE and is modeled using the Complex Variable Boundary Element Method. The transient component is governed by the linear diffusion PDE and is modeled by a linear combination of basis functions that are the products of a two-dimensional Fourier sine series and an exponential function. The global approximation function is the sum of the approximate solutions from the two components. The boundary conditions of the steady-state problem are specified to match the boundary conditions from the global problem so that the CVBEM approximation function satisfies the global boundary conditions. Consequently, the boundary conditions of the transient problem are specified to be continuously zero. The initial condition of the transient component is specified as the difference between the initial condition of the global initial-boundary value problem and the CVBEM approximation of the steady-state solution. Therefore, when the approximate solutions from the two components are summed, the resulting global approximation function approximately satisfies the global initial condition. In this work, it will be demonstrated that the coupled global approximation function satisfies the governing diffusion PDE. Lastly, a procedure for developing streamlines at arbitrary model time is discussed.展开更多
This paper presents three boundary meshless methods for solving problems of steady-state and transient heat conduction in nonlinear functionally graded materials(FGMs).The three methods are,respectively,the method of ...This paper presents three boundary meshless methods for solving problems of steady-state and transient heat conduction in nonlinear functionally graded materials(FGMs).The three methods are,respectively,the method of fundamental solution(MFS),the boundary knot method(BKM),and the collocation Trefftz method(CTM)in conjunction with Kirchhoff transformation and various variable transformations.In the analysis,Laplace transform technique is employed to handle the time variable in transient heat conduction problem and the Stehfest numerical Laplace inversion is applied to retrieve the corresponding time-dependent solutions.The proposed MFS,BKM and CTM are mathematically simple,easyto-programming,meshless,highly accurate and integration-free.Three numerical examples of steady state and transient heat conduction in nonlinear FGMs are considered,and the results are compared with those from meshless local boundary integral equation method(LBIEM)and analytical solutions to demonstrate the effi-ciency of the present schemes.展开更多
A hybrid finite element-Laplace transform method is implemented to analyze the time domain electromagnetic scattering induced by a 2-D overfilled cavity embedded in the infinite ground plane.The algorithm divides the ...A hybrid finite element-Laplace transform method is implemented to analyze the time domain electromagnetic scattering induced by a 2-D overfilled cavity embedded in the infinite ground plane.The algorithm divides the whole scattering domain into two,interior and exterior,sub-domains.In the interior sub-domain which covers the cavity,the problem is solved via the finite element method.The problem is solved analytically in the exterior sub-domain which slightly overlaps the interior subdomain and extends to the rest of the upper half plane.The use of the Laplace transform leads to an analytical link condition between the overlapping sub-domains.The analytical link guides the selection of the overlapping zone and eliminates the need to use the conventional Schwartz iteration.This dramatically improves the efficiency for solving transient scattering problems.Numerical solutions are tested favorably against analytical ones for a canonical geometry.The perfect link over the artificial boundary between the finite element approximation in the interior and analytical solution in the exterior further indicates the reliability of the method.An error analysis is also performed.展开更多
The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an ...The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an orthogonal function system mixed a weight function is defined. Next the improved moving least-square approximation is discussed in detail. The improved method has higher computational efficiency and precision than the old method, and cannot form an ill-conditioned equation system. A boundary element-free method (BEFM) for elastodynamics problems is presented by combining the boundary integral equation method for elastodynamics and the improved moving least-square approximation. The boundary element-free method is a meshless method of boundary integral equation and is a direct numerical method compared with others, in which the basic unknowns are the real solutions of the nodal variables and the boundary conditions can be applied easily. The boundary element-free method has a higher computational efficiency and precision. In addition, the numerical procedure of the boundary element-free method for elastodynamics problems is presented in this paper. Finally, some numerical examples are given.展开更多
The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the sin...The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the singular element. In the computer program BEMTDY the Koizumi′s numerical inversion method is used and three examples of the 3 D vibrated foundation under harmonic forces and the influence with both adjacent foundations are studied.展开更多
The basic relations in linear isotropic photoviscoelasticity have been theoretically discussed in de- tails.A new routine has been found to obtain the time-dependent principal stress without the measurement of isoclin...The basic relations in linear isotropic photoviscoelasticity have been theoretically discussed in de- tails.A new routine has been found to obtain the time-dependent principal stress without the measurement of isoclinics.As a test of our method,examples are given at the end of this paper.展开更多
Based on the integral equation transformed from three dimensional Laplace equation and by the adoption of the division manner of sub-region boundary element method, the numerical computations of the velocity potential...Based on the integral equation transformed from three dimensional Laplace equation and by the adoption of the division manner of sub-region boundary element method, the numerical computations of the velocity potential of each sub-region are given considering the continuity conditions of potential and normal derivatives at the interface of sub-regions, Therefore, computation of wave deformation in offshore flow field is realized. The present numerical model provides a good solution for the application of boundary element method to the calculation of wave deformation in large areas.展开更多
基金This project is supported by National Natural Science Foundation of China
文摘The numerical methods of Fourier eigen transform FET and its inversion are discussed and applied to the boundary element method for elastodynamics. The program for solving elastodynamic problems with the boundary element method is developed and some examples are given. From the numerical results of the examples, we know the method can increase the computing speed 5 similar to 10 times and the accuracy is guaranteed.
基金National Natural Science Foundation of China(No.49876026)
文摘For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The method can be used to remove the transverse vibration due to the accumulation of computational errors. A combined boundary condition of sponge layer and Sommerfeld radiation condition is used to remove the reflected waves from the computing domain. By following the water particle on the water surface, the third order Stokes wave transform is simulated by the numerical wave flume technique. The computed results are in good agreement with theoretical ones.
基金supported by the National Natural Science Foundation of China
文摘The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.
基金supported by the National Natural Science Foundation of China(No.52077073).
文摘This paper focuses on the finite element method in the complex frequency domain(CFD-FEM)for the transient electric field.First,the initial value and boundary value problem of the transient electric field under the electroquasistatic field in the complex frequency domain is given.In addition,the finite element equation and the constrained electric field equation on the boundary are derived.Secondly,the indirect algorithm of the numerical inverse Laplace transform is introduced.Based on it,the calculation procedures of the CFD-FEM are illustrated in detail.Thirdly,the step response,zero-state response under the positive periodic square waveform(PPSW)voltage,and the zero-input response by the CFD-FEM with direct algorithm and indirect algorithm are compared.Finally,the reason for the numerical oscillations of the zero-state response under the PPSW voltage is analyzed,and the method to reduce oscillations is proposed.The results show that the numerical accuracy of the indirect algorithm of the CFD-FEM is more than an order of magnitude higher than that of the direct algorithm when calculating the step response of the transient electric field.The proposed method can significantly reduce the numerical oscillations of the zero-state response under the PPSW voltage.The proposed method is helpful for the calculation of the transient electric field,especially in the case of frequency-dependent parameters.
文摘In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.
文摘The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology in this work is suitable for modeling diffusion problems with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The underpinning of the modeling approach is to decompose the global initial-boundary value problem into a steady-state component and a transient component. The steady-state component is governed by the Laplace PDE and is modeled using the Complex Variable Boundary Element Method. The transient component is governed by the linear diffusion PDE and is modeled by a linear combination of basis functions that are the products of a two-dimensional Fourier sine series and an exponential function. The global approximation function is the sum of the approximate solutions from the two components. The boundary conditions of the steady-state problem are specified to match the boundary conditions from the global problem so that the CVBEM approximation function satisfies the global boundary conditions. Consequently, the boundary conditions of the transient problem are specified to be continuously zero. The initial condition of the transient component is specified as the difference between the initial condition of the global initial-boundary value problem and the CVBEM approximation of the steady-state solution. Therefore, when the approximate solutions from the two components are summed, the resulting global approximation function approximately satisfies the global initial condition. In this work, it will be demonstrated that the coupled global approximation function satisfies the governing diffusion PDE. Lastly, a procedure for developing streamlines at arbitrary model time is discussed.
文摘This paper presents three boundary meshless methods for solving problems of steady-state and transient heat conduction in nonlinear functionally graded materials(FGMs).The three methods are,respectively,the method of fundamental solution(MFS),the boundary knot method(BKM),and the collocation Trefftz method(CTM)in conjunction with Kirchhoff transformation and various variable transformations.In the analysis,Laplace transform technique is employed to handle the time variable in transient heat conduction problem and the Stehfest numerical Laplace inversion is applied to retrieve the corresponding time-dependent solutions.The proposed MFS,BKM and CTM are mathematically simple,easyto-programming,meshless,highly accurate and integration-free.Three numerical examples of steady state and transient heat conduction in nonlinear FGMs are considered,and the results are compared with those from meshless local boundary integral equation method(LBIEM)and analytical solutions to demonstrate the effi-ciency of the present schemes.
基金This workwas supported in part by the Air Force Office of Scientific Research.
文摘A hybrid finite element-Laplace transform method is implemented to analyze the time domain electromagnetic scattering induced by a 2-D overfilled cavity embedded in the infinite ground plane.The algorithm divides the whole scattering domain into two,interior and exterior,sub-domains.In the interior sub-domain which covers the cavity,the problem is solved via the finite element method.The problem is solved analytically in the exterior sub-domain which slightly overlaps the interior subdomain and extends to the rest of the upper half plane.The use of the Laplace transform leads to an analytical link condition between the overlapping sub-domains.The analytical link guides the selection of the overlapping zone and eliminates the need to use the conventional Schwartz iteration.This dramatically improves the efficiency for solving transient scattering problems.Numerical solutions are tested favorably against analytical ones for a canonical geometry.The perfect link over the artificial boundary between the finite element approximation in the interior and analytical solution in the exterior further indicates the reliability of the method.An error analysis is also performed.
基金supported by the National Natural Science Foundation of China(Grant No.10571118)the Shanghai Leading Academic Discipline Project(Grant No.Y0103).
文摘The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an orthogonal function system mixed a weight function is defined. Next the improved moving least-square approximation is discussed in detail. The improved method has higher computational efficiency and precision than the old method, and cannot form an ill-conditioned equation system. A boundary element-free method (BEFM) for elastodynamics problems is presented by combining the boundary integral equation method for elastodynamics and the improved moving least-square approximation. The boundary element-free method is a meshless method of boundary integral equation and is a direct numerical method compared with others, in which the basic unknowns are the real solutions of the nodal variables and the boundary conditions can be applied easily. The boundary element-free method has a higher computational efficiency and precision. In addition, the numerical procedure of the boundary element-free method for elastodynamics problems is presented in this paper. Finally, some numerical examples are given.
文摘The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the singular element. In the computer program BEMTDY the Koizumi′s numerical inversion method is used and three examples of the 3 D vibrated foundation under harmonic forces and the influence with both adjacent foundations are studied.
文摘The basic relations in linear isotropic photoviscoelasticity have been theoretically discussed in de- tails.A new routine has been found to obtain the time-dependent principal stress without the measurement of isoclinics.As a test of our method,examples are given at the end of this paper.
基金The present research work was financially supported by the National Natural Science Foundation of China (Grant No. 49876026) by Hongkong Research Grants Council (No. 49910161985)
文摘Based on the integral equation transformed from three dimensional Laplace equation and by the adoption of the division manner of sub-region boundary element method, the numerical computations of the velocity potential of each sub-region are given considering the continuity conditions of potential and normal derivatives at the interface of sub-regions, Therefore, computation of wave deformation in offshore flow field is realized. The present numerical model provides a good solution for the application of boundary element method to the calculation of wave deformation in large areas.
