Showing papers on "Finite difference method published in 1974"
••
TL;DR: In this paper, an unconditionally stable second order accurate, implicit, finite difference method is described, where the coefficient matrix is tridiagonal and always diagonally dominant, and it is used to solve Burgers' equation.
492 citations
••
TL;DR: In this article, the field region is discretised by triangular Finite-elements and the solution to the field problem is obtained by minimizing the energy functional with respect to each of the vertex values of the vector potential.
Abstract: Analysis of the eddy-currentproblem in magnetic structures by the method of Finite-elements is presented. The linear diffusion equation representing the appropriate energy functional is described. The field region is discretised by triangular Finite-elements and the solution to the field problem is obtained by minimizing the energy functional with respect to each of the vertex values of the vector potential. Expressions for the magnetic field, electric field and eddy-current losses are presented. The method is applied to a few cases of engineering interest and compared with results of classical analysis and tests.
160 citations
••
153 citations
••
TL;DR: The finite element method, using smooth splines as basis functions, applied to the model problem $u_t = cu_x $ with periodic data generates a differential-difference equation whose phase error is closely estimated and compared with the phase error of both explicit and high order implicit centered differencing as mentioned in this paper.
Abstract: The finite element method, using smooth splines as basis functions, applied to the model problem $u_t = cu_x $ with periodic data generates a differential-difference equation whose phase error is closely estimated and compared with the phase error of both explicit and high order implicit centered differencing. We also compute and compare the minimum work required to obtain a fixed error for several fully discrete schemes.
99 citations
••
TL;DR: An adaptive finite difference method for first order nonlinear systems of ordinary differential equations subject to multipoint nonlinear boundary conditions is presented and is based on a discretization studied earlier by H. B. Keller.
Abstract: An adaptive finite difference method for first order nonlinear systems of ordinary differential equations subject to multipoint nonlinear boundary conditions is presented. The method is based on a discretization studied earlier by H. B. Keller. Variable order is provided through deferred corrections, while a built-in natural asymptotic estimator is used to automatically refine the mesh in order to achieve a required tolerance. Extensive numerical experimentation and a FORTRAN program
87 citations
••
TL;DR: In this article, a linear programming (LP) model for ground-water flow is proposed, which replaces the differential equations of groundwater flow by finite-difference approximations that include unknown sink/source terms.
Abstract: The method first replaces the differential equations of ground-water flow by finite-difference approximations that include unknown sink/source terms. The resulting system of algebraic linear equations has a rectangular matrix of coefficients. This system, together with linear inequalities relating sink/source terms, heads or both, and together with an objective function, forms a linear programming (LP) model. The method is applied to small-scale models of confined and unconfined saturated flow for steady-state and transient cases. The steady-state LP models are solved using available computer codes. For the transient confined model, the Crank-Nicolson scheme is used, and a single LP problem is solved covering all of the time steps. For the transient unconfined model, a predictor technique is used, and a LP problem is solved at each corrector step. The optimal solutions are consistent with the results of traditional analyses.
86 citations
••
TL;DR: In this paper, the orthogonal collocation method is used to obtain approximate solutions to the differential equations modeling chemical reactors, which is very often useful in engineering work, where valid approximations are accepted.
Abstract: The orthogonal collocation method is used to obtain approximate solutions to the differential equations modeling chemical reactors. There are two ways to view applications of the orthogonal collocation method. In one view it is a numerical method for which the convergence to the exact answer can be seen as the approximation is refined in successive calculations by using more collocation points, which are similar to grid points in a finite difference method. Another viewpoint considers only the first approximation, which can often be found analytically, and which gives valuable insight to the qualitative behavior of the solution. The answers, however, are of uncertain accuracy, so that the calculation must be refined to obtain useful numbers. However, with experience and appropriate caution, the first approximation is often sufficient and is easy to obtain. Thus it is very often useful in engineering work, where valid approximations are accepted. We present both viewpoints here: we use the first a...
73 citations
••
64 citations
••
TL;DR: In this paper, a detailed theoretical and experimental study of the problem for pipes of square and rectangular cross-sections is presented, where the main concern is the form that the secondary flows take and their influence on the flow rate through the pipe.
61 citations
••
01 Jul 1974TL;DR: In this paper, a finite difference method for the solution of the transonic flow about a harmonically oscillating wing is presented, where the flow is divided into steady and unsteady perturbation velocity potentials.
Abstract: A finite difference method for the solution of the transonic flow about a harmonically oscillating wing is presented. The partial differential equation for the unsteady transonic flow was linearized by dividing the flow into separate steady and unsteady perturbation velocity potentials and by assuming small amplitudes of harmonic oscillation. The resulting linear differential equation is of mixed type, being elliptic or hyperbolic whereever the steady flow equation is elliptic or hyperbolic. Central differences were used for all derivatives except at supersonic points where backward differencing was used for the streamwise direction. Detailed formulas and procedures are described in sufficient detail for programming on high speed computers. To test the method, the problem of the oscillating flap on a NACA 64A006 airfoil was programmed. The numerical procedure was found to be stable and convergent even in regions of local supersonic flow with shocks.
56 citations
••
TL;DR: In this paper, a new theory is presented for the approximation of second order elliptic Dirichlet boundary value problems by a well-known finite difference scheme on irregular networks which cannot be handled by the finite element theory.
Abstract: A new theory is presented for the approximation of second order elliptic Dirichlet boundary value problems by a well-known finite difference scheme on irregular networks which cannot be handled by the finite element theory. Adequate regularity conditions are imposed on the network, a finite difference scheme is set up, and the convergence of its solution towards the solution of the continuous boundary value problems is established, as well as a bound for the discretization error. A numerical example illustrates the theory.
••
TL;DR: In this paper, a subclass of general functional difference equations is derived which retain their stability under appropriate perturbations, and the results of Hale and Cruz are extended to include difference equations on the Banach space of p-th power integrable functions and essentially bounded functions.
••
TL;DR: In this article, the authors derived analytical solutions for a semi-ellipsoid and a pair of parallel semi-cylinders immersed in an otherwise homogeneous semi-infinite medium.
Abstract: Summary Refraction of heat flow across near surface boundaries between dissimilar media results in local variations in heat flow. The magnitude of the perturbation from the regional value is a function of position, conductivity ratio, and the configuration of media. New analytical solutions are given for a semi-ellipsoid and a pair of parallel semi-cylinders immersed in an otherwise homogeneous semi-infinite medium. An analytical solution is also derived for the case of a buried cylinder. For two-dimensional but arbitrary geometric configurations, a numerical scheme based on the finite element method is formulated to include the effects of conductivity contrast as well as relief, non-uniform surface temperature, and internal heat generation. The applicability of this scheme is demonstrated by comparing the numerical results with analytical solutions and results obtained from a finite difference method.
••
TL;DR: A numerical procedure that uses an explicit finite difference method to solve the wave equation is described and results for several specific empty cavity confocal unstable resonators are presented and compared to results obtained using an established calculation technique.
Abstract: A numerical procedure that uses an explicit finite difference method to solve the wave equation is described. This technique results in a propagation algorithm that can accurately propagate an arbitrary electric field through a uniform medium or a medium that is nonuniform, transversely flowing, saturable, and contains index inhomogeneities. By using the propagation algorithm to propagate an arbitrary field back and forth between two resonator mirrors, the three-dimensional transverse mode and the output beam characteristics for a laser resonator can be determined. The advantage of the finite difference method is that unlike integral techniques the computational accuracy and efficiency improve as the resonator Fresnel number increases. The computational techniques are explained, and results for several specific empty cavity confocal unstable resonators are presented and compared to results obtained using an established calculation technique. The application of the finite difference method to inhomogeneous laser media is described, and computational results for an existing CO2 gas dynamic laser are presented and compared to measured data. The medium kinetics and shock wave models used in the calculations are described.
••
TL;DR: The Galerkin method of approximation permits a high-order approximation in time as well as in space as mentioned in this paper, and the resulting approximate equations have been successfully solved by using a prismatic element with triangular cross section.
Abstract: Whereas considerable effort has been expended in generating approximations to the spatial derivatives encountered in porous media flow, the time derivative has received relatively little attention. In spite of the fact that sophisticated finite element formulations have been developed for the spatial derivatives, finite difference methods are generally applied to the time derivative. The Galerkin method of approximation permits a high-order approximation in time as well as in space. The resulting approximate equations have been successfully solved by using a prismatic element with triangular cross section. The time axis runs the length of the prism and is subdivided into elements that may be linear, quadratic, or cubic. Because this formulation requires in general the solution for several time levels simultaneously, there is a resulting increase in computer time required to solve the larger matrix. Numerical experiments indicate that the selection of an optimum numerical scheme is dependent not only on the particular problem considered but also on the sequence of time steps used.
••
TL;DR: In this paper, the approximate analytical solutions to the finite-geometry moving-boundary problem were presented and compared with those obtained by a finite-difference method based on a simple model of the system.
••
01 Jan 1974TL;DR: In this article, the relative efficiency of three types of spatial differencing used in setting up differential-difference equations for the 1-periodic parabolic problem ut = uxx was evaluated.
Abstract: 0. Summary.
Three subjects are considered in this paper. First, the notion of the resolving power of approximation methods, i.e. the number of intervals (or function values) per wavelength necessary to attain a preassigned error when approximating a given frequency, permits the evaluation of the relative efficiency of three types of spatial differencing used in setting up differential-difference equations for the 1-periodic parabolic problem ut = uxx. The three types of high-order spatial schemes considered are (explicit) centered differencing, the super-convergent smooth spline-Galerkin schemes discovered by Thomee, and the very high-order implicit schemes (Mehrstellenverfahren) which generalize Numerov's method. Seven time discretizations are introduced, namely Euler, backward-Euler, duFort-Frankel, 2nd order explicit, trapezoidal, Calahan-Zlamal, and 4th order Pade. The computational work necessary to solve each full discretization is minimized, for given error requirements, by balancing the number of intervals per wavelength against the number of time intervals per eth-life. The resulting data is used to compare the relative efficiency of these finite element and finite difference methods.
Secondly, some corresponding results for the hyperbolic problem ut = ux are briefly reviewed.
Finally, as Numerov-trapezoidal differencing turns out to be almost optimal for the heat equation, a tridiagonal implicit difference formula which extends Numerov differencing to the general 2nd order linear differential operator in one space dimension is presented. The technique used in deriving this scheme inspires certain difference analogs of some finite element schemes. It also leads to a curious modification of the diamond difference scheme for ut = ux. This “alternating kite” scheme has O(Δ2 + Δx3) truncation error with no more computational work than the diamond scheme itself.
••
TL;DR: In this article, Taylor's theorem is used to obtain the field values around the point '0' as: for the real part of the field value and similar equations for the imaginary part (a misplaced minus sign which appeared in the earlier paper has been correctly located above).
Abstract: It is apparent that the finite difference approximation to the Laplacian obtained by Williamson, Hewlett & Tammemagi (1974) differs from that given by Jones R: Pascoe (1971). The difference is related to the manner in which first derivatives are approximated in the two cases. Jones & Pascoe assunzed a central diflerence form for first derivatives from the outset, whereas the method of Williamson et al. implies a linear cornhination of backward and forward differences for first derivatives. It is instructive to examine the two approaches in greater detail. Following Jones & Pascoe, Taylor's theorem may bc used to obtain the field values around the point '0' as: for the real part of the field values and similar equations for the imaginary part (a misplaced minus sign which appeared in the earlier paper has been correctly located above). Combining the equations for f1 and f3 gives:
••
TL;DR: In this paper, a new quantum mechanical procedure, using the finite-difference boundary value method and an equation derived by Hazi and Taylor for the stabilization method, is proposed for finding the widths or lifetimes of resonance states in atomic collisions.
••
TL;DR: In this paper, the equations of motion for a high Prandtl number Boussinesq fluid in a square 2D cavity with side-wail heating and cooling and perfectly conducting end walls have been solved by means of Galerkin as well as ADI (alternating-direction-implicit) finite difference methods for Rayleigh numbers up to 8 × 10 6 and two angles of tilt.
••
TL;DR: An incremental finite element method of the large elastic-plasic deformation of non-axisymmetric metal diaphragms which are edge-clamped and bulged by hydrostatic pressure has been formulated in this paper.
Abstract: An incremental finite element method of the large elastic-plasic deformation of non-axisymmetric metal diaphragms which are edge-clamped and bulged by hydrostatic pressure has been formulated. The diaphragm is divided into a number of flat triangular elements and the behavior of a typical element is described in terms of the displacements of its nodes. The authors derive stiffness matrices, taking account of effects of shape change on the equilibrium equations, from a total differential of the equivalent nodal forces. In order to check the theory, numerical calculations were carried out for the bulging of a circular diaphragm under hydrostatic pressure. Theoretical results were in reasonable agreement with experiments for aluminium sheet and also with numerical solutions by means of a finite difference method.
••
TL;DR: In this paper, the problem of steady state free convection in an unconfined aquifer bounded by ocean on the sides with geothermal heating from below is investigated, and the governing nonlinear partial differential equations with nonlinear boundary conditions are approximated by a set of linear subproblems on the basis of perturbation method.
Abstract: The problem of steady state free convection in an unconfined aquifer bounded by ocean on the sides with geothermal heating from below is investigated in this paper. The governing nonlinear partial differential equations with nonlinear boundary conditions are approximated by a set of linear subproblems on the basis of perturbation method. The equations for the zero- and first-order approximations are of the elliptic type that can be solved numerically by the finite difference method. Numerical results, accurate to the first-order approximations, are obtained for temperature, pressure, and stream function as well as for the shape of the water table. The influence of the location and the size of the heat source as well as various parameters on heat transfer and fluid flow characteristics in a rectangular geothermal aquifer is discussed.
••
17 Jun 1974••
TL;DR: In this paper, the Colorado method for the solution of the non-linear form of Laplace's, Poisson's, and the diffusion partial differential equations is explained, and two methods of acceleration of convergence are described.
Abstract: The Colorado method for the solution of the non-linear form of Laplace's, Poisson's, and the diffusion partial differential equations is explained. Various boundary conditions can be satisfied. The transformation of the partial differential equation into a large set of finite difference equations is given. The discretization is based on a grid system consisting of two sets of orthogonal grid lines. The resulting meshes are nonuniform. Successive line overrelaxation method is used for the solution of the nonlinear equations in two steps. For the improvement of convergence, two methods of acceleration of convergence are described.
••
••
TL;DR: In this article, a finite difference method for the numerical solution of electromagnetic waveguide discontinuity problems is presented, where boundary relaxation is applied, using finite difference techniques in the nonuniform section of the waveguide and using a modal representation of the field in the uniform sections.
Abstract: Finite difference method for the numerical solution of electromagnetic waveguide discontinuity problems is presented. The method of boundary relaxation is applied, using finite difference techniques in the nonuniform section of the waveguide and using a modal representation of the field in the uniform sections of the waveguide. To illustrate the process some two-dimensional diffraction problems in an electromagnetic waveguide with rectangular cross section are solved.
••
TL;DR: In this paper, a finite-difference approximation to the three boundary value problems for Poisson's equation is given, with discretization errors of O(H^3) for the mixed boundary value problem, O(h^3 |ln(h)| for the Neumann problem, and O( h^4 ) for the Dirichlet problem respectively.
Abstract: In this thesis finite-difference approximations to the three boundary value problems for Poisson’s equation are given, with discretization errors of O(H^3) for the mixed boundary value problem, O(H^3 |ln(h)| for the Neumann problem and O(H^4)for the Dirichlet problem respectively . First an operator is constructed, which approximates the nor-ma1 derivative with a truncation errors of O(H^3). The derivation by which this result is obtained contains an improvement upon the one used for a similar operator by Bramble and Hubbard; it became thus possible to make their results valid under more general conditions. For points in a square net where the nine-point approximation to the Laplace operator cannot be used, because of their position near the boundary of the region under consideration, several approximations to the Laplace operator are given, dependent on the particular point configuration, which are all O(H^2). The above-mentioned operators are then used to formulate finite-difference problems with solutions approximating the corresponding continuous problems with the desired accuracy. These error bounds are an improvement of O(H) upon the most accurate approximations known for the Neumann and Robin problems, while the approximation given for the Dirichlet problem has as an advantage over a similar approximation given by Bramble and Hubbard, which is of the same order, that its res-ulting coefficient matrix is of positive type. All our approximations share this important property. For all three problems simple numerical examples are given. It is also pointed out that certain results of Bramble and Hubbard, with respect to the error in the numerical partial derivatives of the solutions, are valid for our approximations.
••
TL;DR: In this article, the effects of finite rates of reaction and ambipolar diffusion are considered, and the inner boundary conditions for both the electron energy and the electron continuity equations are obtained by investigating a sheath analysis.
Abstract: The compressible laminar and turbulent boundary-layer equations representing the fiow of cesium seeded argon plasma over successive segments of the cathode wall of an MHD channel were investigated and solved for both finite and intinite rates of reaction. The equations solved include the electron energy and the electron continuity equations as well as the usual global boundary-layer equations. The effects of finite rates of reaction and ambipolar diffusion are considered, and the inner boundary conditions for both the electron energy and the electron continuity equations are obtained by investigating a sheath analysis. The governing equations are solved using a finite-difference scheme without introducing a local similarity assumption. The finite rate results are considerably different from the infinite rate ones. The main reason is the ambipolar diffusion which depletes the electrons near the wall. No special problems have been found in the solution of the turbulent flow equations. This is due to the high mobility of electrons. (auth)
••
TL;DR: In this article, an alternative approach to the solution of the two-center electronic Schrodinger equation involving a finite-difference Newton-Raphson algorithm is described, where the usual separation in confocal elliptical coordinates is employed leading to two coupled one-dimensional differential equations with split boundary conditions.
••
TL;DR: In this article, the boundary conditions for large deflection of rectangular sandwich plates were derived using the principle of complementary energy, and the governing differential equations were transformed into systems of nonlinear algebraic equations using the finite difference method, and solved by successive iteration.