The Broadwell model with a subsonic boundary

In this paper, we focus on the time-asymptotic behavior of an initial boundary value problem (IBVP) for the Broadwell model with a subsonic physical boundary. By using the Green’s function for the initial problem established in [C.-Y. Lan, H.-E. Lin, S.-H. Yu, The Green’s functions for the Broadwell model in half space problem, Netw. Heterog. Media 1 (1) (2006)] and the weighted energy estimates, we construct the Green’s function for IBVP and show that the solution converges pointwise to the equilibrium state when the perturbations are sufficiently small.     

A Green’s function formulation of nonlocal finite-difference schemes for reaction-diffusion equations

Reaction-diffusion equations are commonly used in different science and engineering fields to describe spatial patterns arising from the interaction of chemical or biochemical reactions and diffusive transport mechanisms. The aim of this work is to show that a Green’s function formulation of reaction-diffusion PDEs is a suitable framework to derive FD schemes incorporating both O(h²) accuracy and nonlocal approximations in the whole domain (including boundary nodes). By doing so, the approach departs from a Green’s function formulation of the boundary-value problem to pose an approximation problem based on a domain decomposition. Within each subdomain, the corresponding integral equation is forced to have zero residual at given grid points. Different FD schemes are obtained depending on the numerical scheme used for computing the Green’s integral over each subdomain. Dirichlet and Neumann boundary conditions are considered, showing that the FD scheme based on the Green’s function formulation incorporates, in a natural way, the effects of boundary nodes in the discretization approximation.     

On the Green’s function for the Helmholtz operator in an impedance circular cylindrical waveguide

This paper addresses the problem of finding a series representation for the Green’s function of the Helmholtz operator in an infinite circular cylindrical waveguide with impedance boundary condition. Resorting to the Fourier transform, complex analysis techniques and the limiting absorption principle (when the undamped case is analyzed), a detailed deduction of the Green’s function is performed, generalizing the results available in the literature for the case of a complex impedance parameter. Procedures to obtain numerical values of the Green’s function are also developed in this article.     

Variational iteration method for elastodynamic Green’s functions

In this study, a new application of the variational iteration method is presented. We use this method to obtain approximations to 3D Green’s function for the dynamic system of anisotropic elasticity. The numerical results obtained from convolution of Green’s function and data of the Cauchy problem are compared with the exact solution for cubic media.     

2D Green’s functions for semi-infinite orthotropic thermoelastic plane

Based on the 2D general solutions of orthotropic thermoelastic material, the Green’s function for a steady point heat source in the interior of semi-infinite orthotropic thermoelastic plane is constructed by three newly introduced harmonic functions. All components of coupled field in semi-infinite thermoelastic plane are expressed in terms of elementary functions. Numerical results are given graphically by contours.     

The space fractional diffusion equation with Feller’s operator

This paper investigates the space fractional diffusion equation with fractional Feller’s operator. The Green’s function is obtained by using Fourier transform, and the analytical solutions of some space fractional diffusion equations with initial (or initial and boundary) condition are obtained in terms of Green’s function. In addition, numerical simulations are discussed. The results indicate that the effect range of skewness parameter θ has more effect on probability density than that of parameter α. The results also explain the property of the skewness and long tail in the asymmetry diffusion process.     

Using spectral element method for analyzing continuous beams and bridges subjected to a moving load

Spectral element method in frequency domain is employed to analyze continuous beams and bridges subjected to a moving load. The formulation is developed for an Euler beam under a moving load with an arbitrary amplitude and velocity. It is shown that the procedure is simplified for a moving load with a constant amplitude and velocity. Static Green’s function is used as a modifying function to improve the moment and shear force results. It is further shown that while modifying function is used in conjunction with spectral element method, fewer elements will be required to achieve proper results. The numerical examples show the accuracy of the method.     

Analytical modelling for a three-dimensional hydrofoil with winglets operating beneath a free surface

The current study focuses on establishing a theoretical lifting surface model for predicting the hydrodynamic loads acting on the three-dimensional hydrofoil with winglets, which is considerably influenced by the proximity to the free surface through finding the three-dimensional Green’s function for the planar and vertical horseshoe vortices operating below a free surface. The hydrofoil surface is decomposed into a finite number of elements along the span direction and the chord directions, each of which can then be represented by a horseshoe vortex. The linearized free surface boundary condition is applied to analyze the influence of the free surface on the hydrofoil as well as the winglets. The thickness problem is considered using the source distribution among the hydrofoil and winglets surfaces and the analytical Green’s function that satisfies the linearized free surface boundary condition is used. As a sample application, numerical examples were conducted to show the performance of the hydrodynamic characteristics for the hydrofoil with winglets as a function of the Froude number. It was concluded that there are significant efficiency benefits from using winglets inside the free surface proximity effect. These results are substantiated by the comparison with the available published data.     

New closed-form Green function and integral formula for a thermoelastic quadrant

In this study a new Green’s function and a new Green-type integral formula for a boundary value problem (BVP) in thermoelastostatics for a quadrant are derived in closed form. On the boundary semi-straight-lines twice mixed homogeneous mechanical boundary conditions (one boundary semi-straight-line is free of loadings and normal displacements and tangential stresses are prescribed on the other one) are prescribed. The thermoelastic displacements are subject by a heat source applied in the inner points of the quadrant and by mixed non-homogeneous boundary heat conditions (on one boundary semi-straight-line the temperature is prescribed and the heat flux is given on the other one). When thermoelastic Green’s function is derived the thermoelastic displacements are generated by an inner unit point heat source, described by δ Dirac’s function. All results are obtained in elementary functions that are formulated in a special theorem. A closed-form solution for a particular BVP of thermoelastostatics for a quadrant also is included. Using the proposed approach it is possible to extend the obtained for quadrant results to any other canonical Cartesian domain.     

Coupled boundary element method and finite difference method for the heat conduction in laser processing

This paper is presented as a way to model transient heat conduction in a 3-D axisymmetric case where large rates of heat fluxes are applied on the surfaces as done in the case of laser processing. This would result in large temperature gradients in a small area irradiated by the laser on the incident surface that could also reach melting and subsequent vaporization. BEM can handle large fluxes very easily and it also can be formulated if needed to incorporate the moving boundary problem in a unique manner while on the other hand FDM is a fast and efficient method. For these reasons a coupled BEM–FDM method is formulated to simulate the heat conduction process. In the BEM method linear elements for the boundary and quadratic elements for the domain were used. The integrals in BEM were integrated in time using the asymptotic expansion for the modified Bessel functions in the Green’s function. To further improve the accuracy, special techniques were employed in the spatial integration. As for the FDM formulation, a flux conservation scheme with a 4th order formula for the fluxes was used. The FDM and BEM were coupled at the interface by the temperature from the FDM formulation being imposed on the BEM and the flux from the BEM being utilized by the FDM elements near to the interface. To advance in time, the Crank–Nicholson scheme was used on the FDM directly and due to coupling indirectly on the BEM. The relative errors for the simulation of constant and variable flux cases demonstrate the successful nature of the numerical model.     

Special Green’s function boundary element approach for steady-state axisymmetric heat conduction across low and high conducting planar interfaces

The problem of determining the steady-state axisymmetric temperature distribution in a bimaterial with a planar interface is considered here. The interface is either low or high conducting. Special Green’s functions satisfying the thermal conditions on the interface are derived and employed to obtain boundary integral equations whose path of integration does not include the interface. Boundary element procedures that do not require the interface to be discretized into elements are proposed for solving the problem under consideration.     

Mathematical modelling and analytical solution for workpiece temperature in grinding

This paper deals with modelling the workpiece temperature field produced during the grinding process. The proposed model is given in terms of a two-dimensional boundary-value problem where the interdependence among the grinding wheel, the workpiece and the coolant is described by two variable functions in the boundary condition. An explicit integral form solution is constructed using the Laplace and Fourier transforms and the Green’s function method.     

Green’s function,temperature in a convectively cooled sphere with arbitrarily located spherical heat sources

Steady state heat conduction in a convectively cooled sphere having arbitrarily located spherical heat sources inside is treated with the method of Green’s function accompanied by a coordinate transform. Green’s function of the heat diffusion operator for a finite sphere with Robin boundary condition is obtained by spherical harmonics expansion. Verification of the analytical solution is exemplified in some generic cases related to the pebbles of South-African PBMR as of year 2000 with 268 MW thermal power. Analytical results for different sectors of the sphere (pebble) are compared with the results of computational fluid dynamics code FLUENT^™. This work is motivated through a modest effort to assess the stochastic effects of distribution and volumetric effects of fuel kernels within the pebbles of future-promising pebble bed reactors.     

A new criterion of optimization of the simple multipole coefficients in a modified Green’s function for the elastic two-dimensional case

The question of non-uniqueness in the integral formulation of an exterior boundary value problem in the elastic two-dimensional case has been resolved using the modified Green’s function technique. In this work, a new criterion of optimality based on the minimization of the norm of the kernel of the modified integral operator is established.     

A discrete fourth-order Lidstone problem with parameters

Various existence, multiplicity, and nonexistence results for nontrivial solutions to a nonlinear discrete fourth-order Lidstone boundary value problem with dependence on two parameters are given, using a symmetric Green’s function approach. An existence result is also given for a related semipositone problem, thus relaxing the usual assumption of nonnegativity on the nonlinear term.     

Strain gradient beam model for dynamics of microscale pipes conveying fluid

Based on the strain gradient theory, we present a microstructure-dependent Bernoulli–Euler model to analyze the vibration and stability of microscale pipes conveying fluid. The equation of motion and boundary conditions are derived using Hamilton’s principle. The proposed strain gradient beam model contains three material length scale parameters to capture the size effect. This new model may be reduced to the modified couple stress beam model when two of these three material length scale parameters vanish and may be reduced to the classical beam model in the absence of all the material length scale parameters. From the numerical calculations for micropipes with both ends positively supported, it is found that the natural frequency and the critical flow velocity are size-dependent. The results show that the microscale pipe displays remarkable size effect when its outside diameter becomes comparable to the material length scale parameter, while the size effect is almost diminishing as the diameter is far greater than the material length scale parameter. Moreover, the size effect predicted by the current strain gradient beam model is stronger than that predicted by the modified couple stress beam model, since two other material length scale parameters have been accounted for in the former.     

Effects of Young’s modulus on response of railway sleeper

As a main part of a railroad system, sleepers have important duty in conveying the load from rails to the ballast. The different situations in which the sleepers should function necessitate making them from different materials, such as various types of wood, reinforced concrete and even steel. In this work, the effects of Young’s modulus on response of railway sleeper are evaluated. As a main consideration, Winkler’s theorem is used to model the behavior of the elastic foundation. First, the response of a sleeper on a Winkler’s foundation is found. To evaluate the results of the closed form solution, a finite element model is used. Good agreement between the results of the closed form solution and the finite element model proves the validity of the results. In the next stage, the Young’s modulus is considered as a variable and the fundamental diagrams of the beam are plotted based on the variation of Young’s modulus.     

Three-dimensional dynamic Green’s functions in transversely isotropic tri-materials

An analytical derivation of the elastodynamic fundamental solutions for a transversely isotropic tri-material full-space is presented by means of a complete representation using two displacement potentials. The complete set of three-dimensional point-load, patch-load, and ring-load Green’s functions for stresses and displacements are given, for the first time, in the complex-plane line-integral representations. The formulation includes a complete set of transformed stress-potential and displacement-potential relations in the framework of Fourier expansions and Hankel integral transforms, that is useful in a variety of elastodynamic as well as elastostatic problems. For the numerical computation of the integrals, a robust and effective methodology is laid out. Selected numerical results for point-load and patch-load Green’s functions are presented to portray the dependence of the response on layering, the frequency of excitation, and type of loading.     

A study on the method of fundamental solutions using an image concept

In this paper, both analytical and semi-analytical solutions for Green’s functions are obtained by using the image method which can be seen as a special case of method of fundamental solutions (MFS). The image method is employed to solve the Green’s function for the annular, eccentric and half-plane Laplace problems. In addition, an analytical solution is derived for the fixed-free annular case. For the half-plane problem with a circular hole and an eccentric annulus, semi-analytical solutions are both obtained by using the image concept after determining the strengths of two frozen image points and a free constant by matching boundary conditions. It is found that two frozen images terminated at the two focuses in the bipolar coordinates for the problems with two circular boundaries. A boundary value problem of an eccentric annulus without sources is also considered. Error distribution is plotted after comparing with the analytical solution derived by Lebedev et al. using the bipolar coordinates. The optimal locations for the source distribution in the MFS are also examined by using the image concept. It is observed that we should locate singularities on the two focuses to obtain better results in the MFS. Besides, whether the free constant is required or not in the MFS is also studied. The results are compared well with the analytical solutions.     

The three-body Coulomb scattering problem in a discrete Hilbert-space basis representation

We propose modified Faddeev-Merkuriev integral equations for solving the 2→2, 3 quantum three-body Coulomb scattering problem. We show that the solution of these equations can be obtained using a discrete Hilbert-space basis and that the error in the scattering amplitudes due to truncating the basis can be made arbitrarily small. The Coulomb Green’s function is also confined to the two-body sector of the three-body configuration space by this truncation and can be constructed in the leading order using convolution integrals of two-body Green’s functions. To evaluate the convolution integral, we propose an integration contour that is applicable for all energies including bound-state energies and scattering energies below and above the three-body breakup threshold.