Search results

The purpose of this paper is to develop a new transversal method of lines for one-dimensional reaction–diffusion equations that is conservative and provides piecewise–analytical solutions in space, analyze its truncation errors and linear stability, compare it with other finite-difference discretizations and assess the effects of the nonlinear diffusion coefficients, reaction rate terms and initial conditions on wave propagation and merging.

A conservative, transversal method of lines based on the discretization of time and piecewise analytical integration of the resulting two-point boundary-value problems subject to the continuity of the dependent variables and their fluxes at the control-volume boundaries, is presented. The method provides three-point finite difference expressions for the nodal values and continuous solutions in space, and its accuracy has been determined first analytically and then assessed in numerical experiments of reaction-diffusion problems, which exhibit interior and/or boundary layers.

The transversal method of lines presented here results in three-point finite difference equations for the nodal values, treats the diffusion terms implicitly and is unconditionally stable if the reaction terms are treated implicitly. The method is very accurate for problems with the interior and/or boundary layers. For a system of two nonlinearly-coupled, one-dimensional reaction–diffusion equations, the formation, propagation and merging of reactive fronts have been found to be strong function of the diffusion coefficients and reaction rates. For asymmetric ignition, it has been found that, after front merging, the temperature and concentration profiles are almost independent of the ignition conditions.

A new, conservative, transversal method of lines that treats the diffusion terms implicitly and provides piecewise exponential solutions in space without the need for interpolation is presented and applied to someone.

The purpose of this study is to address one-dimensional, unsteady heat conduction in a large plane wall exchanging heat convection with a nearby fluid under “small time” conditions.

The Transversal Method of Lines (TMOL) was used to reformulate the unsteady, one-dimensional heat conduction equation in the space coordinate and time into a transformed “quasi-steady”, one-dimensional heat conduction equation in the space coordinate housing the time as an embedded parameter. The resulting ordinary differential equation of second order with heat convection boundary conditions is solved analytically with the method of undetermined coefficients.

Semi-analytical TMOL dimensionless temperature profiles of compact form with/without regressed terms are obtained for the whole spectrum of Biot number (0 < Bi < ∞) in the “small time” sub-domain. In addition, a new “large time” sub-domain is redefined, that is, setting a smaller critical dimensionless time or critical Fourier number τ_cr = 0.18.

The computed dimensionless center, surface and mean temperature profiles in the large plane wall accounting for all Biot number (0 < Bi < ∞) in the “small time” sub-domain τ < τ_cr = 0.18 exhibit excellent quality while carrying reasonable relative errors for engineering applications. The exemplary level of accuracy indicates that the traditional evaluation of the center, surface and mean temperatures with the standard infinite series retaining a large number of terms is no longer necessary.

The purpose of this paper is to address a novel method for solving parabolic partial differential equations (PDEs) in general, wherein the heat conduction equation constitutes an important particular case. The new method, appropriately named the Improved Transversal Method of Lines (ITMOL), is inspired in the Transversal Method of Lines (TMOL), with strong insight from the method of separation of variables.

The essence of ITMOL revolves around an exponential variation of the dependent variable in the parabolic PDE for the evaluation of the time derivative. As will be demonstrated later, this key step is responsible for improving the accuracy of ITMOL over its predecessor TMOL. Throughout the paper, the theoretical properties of ITMOL, such as consistency, stability, convergence and accuracy are analyzed in depth. In addition, ITMOL has proven to be unconditionally stable in the Fourier sense.

In a case study, the 1-D heat conduction equation for a large plate with symmetric Dirichlet boundary conditions is transformed into a nonlinear ordinary differential equation by means of ITMOL. The numerical solution of the resulting differential equation is straightforward and brings forth a nearly zero truncation error over the entire time domain, which is practically nonexistent.

Accurate levels of the analytical/numerical solution of the 1-D heat conduction equation by ITMOL are easily established in the entire time domain.

The purpose of this study is to develop a new method of lines for one-dimensional (1D) advection-reaction-diffusion (ADR) equations that is conservative and provides piecewise analytical solutions in space, compare it with other finite-difference discretizations and assess the effects of advection and reaction on both 1D and two-dimensional (2D) problems.

A conservative method of lines based on the piecewise analytical integration of the two-point boundary value problems that result from the local solution of the advection-diffusion operator subject to the continuity of the dependent variables and their fluxes at the control volume boundaries is presented. The method results in nonlinear first-order, ordinary differential equations in time for the nodal values of the dependent variables at three adjacent grid points and triangular mass and source matrices, reduces to the well-known exponentially fitted techniques for constant coefficients and equally spaced grids and provides continuous solutions in space.

The conservative method of lines presented here results in three-point finite difference equations for the nodal values, implicitly treats the advection and diffusion terms and is unconditionally stable if the reaction terms are implicitly treated. The method is shown to be more accurate than other three-point, exponentially fitted methods for nonlinear problems with interior and/or boundary layers and/or source/reaction terms. The effects of linear advection in 1D reacting flow problems indicates that the wave front steepens as it approaches the downstream boundary, whereas its back corresponds to a translation of the initial conditions; for nonlinear advection, the wave front exhibits steepening but the wave back shows a linear dependence on space. For a system of two nonlinearly coupled, 2D ADR equations, it is shown that a counter-clockwise rotating vortical field stretches the spiral whose tip drifts about the center of the domain, whereas a clock-wise rotating one compresses the wave and thickens its arms.

A new, conservative method of lines that implicitly treats the advection and diffusion terms and provides piecewise-exponential solutions in space is presented and applied to some 1D and 2D advection reactions.

The purpose of this paper is to develop a new finite-volume method of lines for one-dimensional reaction-diffusion equations that provides piece-wise analytical solutions in space and is conservative, compare it with other finite-difference discretizations and assess the effects of the nonlinear diffusion coefficient on wave propagation.

A conservative, finite-volume method of lines based on piecewise integration of the diffusion operator that provides a globally continuous approximate solution and is second-order accurate is presented. Numerical experiments that assess the accuracy of the method and the time required to achieve steady state, and the effects of the nonlinear diffusion coefficients on wave propagation and boundary values are reported.

The finite-volume method of lines presented here involves the nodal values and their first-order time derivatives at three adjacent grid points, is linearly stable for a first-order accurate Euler’s backward discretization of the time derivative and has a smaller amplification factor than a second-order accurate three-point centered discretization of the second-order spatial derivative. For a system of two nonlinearly-coupled, one-dimensional reaction-diffusion equations, the amplitude, speed and separation of wave fronts are found to be strong functions of the dependence of the nonlinear diffusion coefficients on the concentration and temperature.

A new finite-volume method of lines for one-dimensional reaction-diffusion equations based on piecewise analytical integration of the diffusion operator and the continuity of the dependent variables and their fluxes at the cell boundaries is presented. The method may be used to study heat and mass transfer in layered media.

The purpose of this paper is to address unsteady heat conduction in two subsets of ordinary bodies. One subset consists of a large plane wall, a long cylinder and a sphere in one dimension. The other subset consists of a short cylinder and a large rectangular bar in two dimensions. The prevalent assumptions in the two subsets are: constant initial temperature, uniform surface heat flux and thermo-physical properties invariant with temperature. The engineering applications of the unsteady heat conduction deal with the determination of temperature–time histories in the two subsets using electric resistance heating, radiative heating and fire pool heating.

To this end, a novel numerical procedure named the enhanced method of discretization in time (EMDT) transforms the linear one-dimensional unsteady, heat conduction equations with non-homogeneous boundary conditions into equivalent nonlinear “quasi–steady” heat conduction equations having the time variable embedded as a time parameter. The equivalent nonlinear “quasi–steady” heat conduction equations are solved with a finite difference method.

Based on the numerical computations, it is demonstrated that the approximate temperature–time histories in the simple subset of ordinary bodies (large plane wall, long cylinder and sphere) exhibit a perfect matching over the entire time domain 0 < t < ∞ when compared against the rigorous exact temperature–time histories expressed by classical infinite series. Furthermore, using the method of superposition of solutions in the convoluted subset (short cylinder and large rectangular crossbar), the same level of agreement in the approximate temperature–time histories in the simple subset of ordinary bodies is evident.

The performance of the proposed EMDT coupled with a finite difference method is exhaustively assessed in the solution of the unsteady, one-dimensional heat conduction equations with prescribed surface heat flux for: a subset of one-dimensional bodies (plane wall, long cylinder and spheres) and a subset of two-dimensional bodies (short cylinder and large rectangular bar).

The purpose of this paper is to address the practically important problem of the load dependence of transverse vibrations for helical springs. At the beginning, the author develops the equations for transverse vibrations of the axially loaded helical springs. The method is based on the concept of an equivalent column. Second, the author reveals the effect of axial load on the fundamental frequency of transverse vibrations and derive the explicit formulas for this frequency. The fundamental natural frequency of the transverse vibrations of the spring depends on the variable length of the spring. The reduction of frequency with the load is demonstrated. Finally, when the frequency nullifies, the side buckling spring occurs.

Helical springs constitute an integral part of many mechanical systems. A coil spring is a special form of spatially curved column. The center of each cross-section is located on a helix. The helix is a curve that winds around with a constant slope of the surface of a cylinder. An exact stability analysis based on the theory of spatially curved bars is complicated and difficult for further applications. Hence, in most engineering applications a concept of an equivalent column is introduced. The spring is substituted for the simplification of the basic equations by an equivalent column. Such a column must account for compressibility of axis and shear effects. The transverse vibration is represented by a differential equation of fourth order in place and second order in time. The solution of the undamped model equation could be obtained by separation of variables. The fundamental natural frequency of the transverse vibrations depends on the current length of the spring. Natural frequency is the function of the deflection and slenderness ratio. Is the fundamental natural frequency of transverse oscillations nullifies, the lateral buckling of the spring with the natural form occurs. The mode shape corresponds to the buckling of the spring with moment-free, simply supported ends. The mode corresponds to the buckling of the spring with clamped ends. The author finds the critical spring compression.

Buckling refers to the loss of stability up to the sudden and violent failure of seed straight bars or beams under the action of pressure forces, whose line of action is the column axis. The known results for the buckling of axially overloaded coil springs were found using the static stability criterion. The author uses an alternative approach method for studying the stability of the spring. This method is based on dynamic equations. In this paper, the author derives the equations for transverse vibrations of the pressure-loaded coil springs. The fundamental natural frequency of the transverse vibrations of the column is proved to be the certain function of the axial force, as well as the variable length of the spring. Is the fundamental natural frequency of transverse oscillations turns to be to zero, is the lateral buckling of the spring occurs.

The spring is substituted for the simplification of the basic equations by an equivalent column. Such a column must account for compressibility of axis and shear effects. The more accurate model is based on the equations of motion of loaded helical Timoshenko beams. The dimensionless for beams of circular cross-section and the number of parameters governing the problem is reduced to four (helix angle, helix index, Poisson coefficient, and axial strain) is to be derived. Unfortunately, that for the spatial beam models only numerical results could be obtained.

The closed form analytical formulas for fundamental natural frequency of the transverse vibrations of the column as function of the axial force, as well as the variable length of the spring are derived. The practically important formulas for lateral buckling of the spring are obtained.

In this paper, the author derives the new equations for transverse vibrations of the pressure-loaded coil springs. The author demonstrates that the fundamental natural frequency of the transverse vibrations of the column is the function of the axial force. For study of the stability of the spring the author uses an alternative approach method. This method is based on dynamic equations. The new, original expressions for lateral buckling of the spring are also obtained.

This paper aims to present an analytical approach for the determination of helical gear tooth geometry and introduces the necessary parameters. Tooth geometry including tooth chamfer, involute curve, root fillet, helix as well as tooth microgeometry can be obtained using the presented approach.

The presented analytical approach involves deriving the equivalent equations at the transverse plane rather than the normal plane. Moreover, numerical evaluation of microgeometry modifications is presented for tooth profile, tooth lead and flank twist.

An analytical approach is presented and equations are derived and explained in detail for helical gear tooth geometry calculation, including tooth microgeometry. Method 1, which was presented by Lopez and Wheway (1986) for obtaining the root fillet, is examined and it is proven that it does not work accurately for helical gears, but rather it works perfectly in the case of spur gears. Changing the normal plane parameters in Method 1 to the transverse plane ones does not give correct results. Two alternative methods, namely, Methods 2 and 3, are developed in the current research for the calculation of the tooth root fillet of helical gears. The presented methods and also the numerical evaluation presented for microgeometry modification are examined against the geometry obtained from Windows LDP software. The results show very good agreement, and it is feasible to apply the approach using the presented equations.

In the gear design process, it is important to model the correct gear tooth geometry and deliver all related dimensions and calculations accurately. However, the determination of helical gear tooth geometry has not been presented adequately by equations to facilitate gear modelling. The detailed helical gear tooth root has been enveloped using software tools that can simulate the cutter motion. Deriving those equations, presented in this article, provides gear design engineers and researchers with the possibility to model helical gears and perform design calculations in a structured, applicable and accurate method.

Gives introductory remarks about chapter 1 of this group of 31 papers, from ISEF 1999 Proceedings, in the methodologies for field analysis, in the electromagnetic community. Observes that computer package implementation theory contributes to clarification. Discusses the areas covered by some of the papers ‐ such as artificial intelligence using fuzzy logic. Includes applications such as permanent magnets and looks at eddy current problems. States the finite element method is currently the most popular method used for field computation. Closes by pointing out the amalgam of topics.

This paper attempts to evaluate the transverse stresses that are generated within the interface between two layers of laminated composite and sandwich laminates by using Cℴ finite element formulation of higher‐order theories. These theories do not require the use of a fictitious shear correction coefficient which is usually associated with the first‐order Reissner‐Mindlin theory. The in‐plane stresses are evaluated by using constitutive relations. The transverse stresses are evaluated through the use of equilibrium equations. The integration of the equilibrium equations is attempted through forward and central direct finite difference techniques and a new approach, named as, an exact surface fitting method. Sixteen and nine‐noded quadrilateral Lagrangian elements are used. The numerical results obtained by the present approaches in general and the exact surface fitting method in particular, show excellent agreement with available elasticity solutions. New results for symmetric sandwich laminates are also presented for future comparisons.