A Numerical Study of Nonaxisymmetric Stokesian Flow in a Circular Tube

1981 ◽  
Vol 48 (3) ◽  
pp. 459-464
Author(s):  
J. Strigberger ◽  
A. Plotkin

A numerical study of the nonaxisymmetric Stokesian flow of a Newtonian fluid in a rigid circular tube of fixed radius has been performed. The analysis presented here is an integral part of the problem of modeling the flow of blood near the ostia of the intercostal arteries of rabbits in order to study a possible factor in the initiation of atherosclerosis. The method of lines is used to reduce the mathematical problem to one of solving a system of first-order ordinary differential equations along lines parallel to the tube axis. Solutions are obtained analytically using matrix eigenvalue techniques for the first two Fourier components of the flow and the accuracy of the numerical method is verified by suitable comparison with the results of independent computations.

Author(s):  
Antonio Campo

The present study on engineering education addresses the Method of Lines and its variant the Numerical Method of Lines as a reliable avenue for the numerical analysis of one-dimensional unsteady heat conduction in walls, cylinders, and spheres involving surface convection interaction with a nearby fluid. The Method of Lines transforms the one-dimensional unsteady heat conduction equation in the spatial and time variables x, t into an adjoint system of first-order ordinary differential equations in the time variable t. Subsequently, the adjoint system of first-order ordinary differential equations is channeled through the Numerical Method of Lines and the powerful fourth-order Runge–Kutta algorithm. The numerical solution of the adjoint system of first-order ordinary differential equations can be carried out by heat transfer students employing appropriate routines embedded in the computer codes Maple, Mathematica, Matlab, and Polymath. For comparison, the baseline solutions used are the exact, analytical temperature distributions that are available in the heat conduction literature.


1972 ◽  
Vol 39 (3) ◽  
pp. 689-695 ◽  
Author(s):  
W. W. Recker

The two-dimensional equations of magnetoelastodynamics are considered as a symmetric hyperbolic system of linear first-order partial-differential equations in three independent variables. The characteristic properties of the system are determined and a numerical method for obtaining the solution to mixed initial and boundary-value problems in plane magnetoelastodynamics is presented. Results on the von Neumann necessary condition are presented. Application of the method to a problem which has a known solution provides further numerical evidence of the convergence and stability of the method.


Sign in / Sign up

Export Citation Format

Share Document