Calculation of subdominant solutions of linear differential equations.
Subdominant solutions of linear differential equations calculated by numerical method
NUMERICAL INTEGRATION↗
SEARCH · Engineering Papers
Results for “numerical differentiation”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Numerical Solution of Symmetric Positive Differential Equations
Finite difference method for symmetric positive linear differential equations, discussing methods to solve finite difference equations generated
Katsanis, T.↗
Numerical integration of second order differential equations
Performance characteristics of higher order approximations of Runge-Kutta type are analyzed, and performance predictors for time required on machine and for error size are developed. Technique is useful in evaluating system performance, analyzing material characteristics, and designing inertial guidance and nuclear instrumentation and materials.
Shanks, E. B.↗
On the numerical solution of elliptic partial differential equations by the method of lines.
Abstract missing.
Jones, D. J.↗
Utility of a finite element solution algorithm for initial-value problems
The Galerkin criterion within a finite element Weighted Residuals formulation is employed to establish an implicit solution algorithm for an initial-value partial differential equation. Numerical solutions of a transient parabolic and a hyperbolic equation, obtained using linear, quadratic and two cubic finite element basis functions, are employed to quantize accuracy and confirm and refine theoretical convergence rate estimates. The linear basis algorithm for the hyperbolic equation displays excellent accuracy on a coarse computational grid and a high-order convergence rate with discretization refinement. Good accuracy and a strong convergence rate in surface flux are determined for a nonhomogeneous Neumann boundary constraint applied to a parabolic equation. The results amply demonstrate the impact of the nondiagonal finite element initial-value matrix structure on solution accuracy and/or convergence rate.
Baker, A. J.↗
Theoretical regime diagrams for thermally driven flows in a beta-plane channel in the presence of variable gravity
The effect of a power law gravity field on baroclinic instability is examined, with a focus on the case of inverse fifth power gravity, since this is the power law produced when terrestrial gravity is simulated in spherical geometry by a dielectric force. Growth rates are obtained of unstable normal modes as a function of parameters of the problem by solving a second order differential equation numerically. It is concluded that over the range of parameter space explored, there is no significant change in the character of theoretical regime diagrams if the vertically averaged gravity is used as parameter.
Geisler, J. E.↗
Higher-Order Panel Method for Aerodynamic Flow Analysis
PANAIR uses high-order panel method to predict inviscid subsonic or supersonic flows about arbitrary configuration. Panel method solves linear partial differential equation numerically by approximating configuration surface with panels on which unknown "singularity strengths" are defined. PANAIR includes advanced software technology as well as advanced aerodynamic technology.
Erickson, L.↗
Boundary-layer effects in composite laminates: Free-edge stress singularities, part 6
A rigorous mathematical model was obtained for the boundary-layer free-edge stress singularity in angleplied and crossplied fiber composite laminates. The solution was obtained using a method consisting of complex-variable stress function potentials and eigenfunction expansions. The required order of the boundary-layer stress singularity is determined by solving the transcendental characteristic equation obtained from the homogeneous solution of the partial differential equations. Numerical results obtained show that the boundary-layer stress singularity depends only upon material elastic constants and fiber orientation of the adjacent plies. For angleplied and crossplied laminates the order of the singularity is weak in general.
Wanag, S. S.↗
Analytical calculation of a single jet in cross-flow and comparison with experiment
A series of progressively finer grid systems to calculate the single jet in cross-flow experimentally measured were employed. The experimental measurements provide a collection of velocities, turbulence intensities, and jet concentration profiles with measurements of the inlet field. The use of a series of progressively finer grid systems allows a differentiation between numerical errors and the hydrodynamic modeling assumptions embodied in the 3D combustor code.
Claus, R. W.↗
A penalty finite element algorithm for parabolic flow problems
The thin-layer simplification of the two-dimensional Navier-Stokes equations for steady viscous flow are developed using an order of magnitude analysis. A space marching finite element solution algorithm is developed, wherein the first order continuity effects are enforced as a penalty function differential constraint. Numerical results are presented to document accuracy and convergence features of the numerical solution algorithm.
Baker, A. J.↗
An explicit predictor-corrector solver with applications to Burgers' equation
Forward Euler's explicit, finite-difference formula of extrapolation, is used as a predictor and a convex formula as a corrector to integrate differential equations numerically. An application has been made to Burger's equation.
Dey, S. K.↗
Stability of limit cycles in frictionally damped and aerodynamically unstable rotor stages
This paper deals with the stability of limit cycles (Steady-State Oscillations) associated with the multi-degree-of-freedom model of a frictionally damped and aerodynamically unstable rotor stage. By using the first order averaging technique, a generalized criterion has been established to sort out those unstable limit cycles which govern the maximum transient amplitude beyond which the rotor stage becomes unstable.The stability of the remaining steady-state solutions is analyzed by linearizing the averaged system of differential equations. Numerical results are discussed for three-, four- and five-bladed disks.
Sinha, A.↗
Numerical integration of asymptotic solutions of ordinary differential equations
Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration.
Thurston, Gaylen A.↗
An automatic multigrid method for the solution of sparse linear systems
An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.
Shapira, Yair↗
On the numeric integration of dynamic attitude equations
We describe new types of numerical integration algorithms developed by the authors. The main aim of the algorithms is to numerically integrate differential equations which evolve on geometric objects, such as the rotation group. The algorithms provide iterates which lie on the prescribed geometric object, either exactly, or to some prescribed accuracy, independent of the order of the algorithm. This paper describes applications of these algorithms to the evolution of the attitude of a rigid body.
Crouch, P. E.↗
Nonlinear interaction between a pair of oblique modes in a supersonic mixing layer: Long-wave limit
The nonlinear interaction between a pair of symmetric, oblique, and spatial instability modes is studied in the long-wave limit using asymptotic methods. The base flow is taken to be a supersonic mixing layer whose Mach number is such that the corresponding vortex sheet is marginally stable according to Miles' criterion. It is shown that the amplitude of the mode obeys a nonlinear integro-differential equation. Numerical solutions of this equation show that, when the obliqueness angle is less than pi/4, the effect of the nonlinearity is to enhance the growth rate of the instability. The solution terminates in a singularity at a finite streamwise location. This result is reminiscent of that obtained in the vicinity of the neutral point by other authors in several different types of flows. On the other hand, when the obliqueness angle is more than pi/4, the streamwise development of the amplitude is characterized by a series of modulations. This arises from the fact that the nonlinear term in the amplitude equation may be either stabilizing or destabilizing, depending on the value of the streamwise coordinate. However, even in this case the amplitude of the disturbance increases, though not as rapidly as in the case for which the angle is less than pi/4. Quite generally then, the nonlinear interaction between two oblique modes in a supersonic mixing layer enhances the growth of the disturbance.
Balsa, Thomas F.↗
Sensitivity analysis for the dynamic response of viscoplastic shells of revolution
A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric response of viscoplastic shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. Spatial discretization is performed using finite elements, with discontinuous stress resultants across element interfaces. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) with an implicit constitutive update. The sensitivity coefficients are evaluated using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading, and a circular plate subjected to impulsive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to thickness, mass density, Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients are presented, along with spatial distributions of these quantities at selected times.
Kulkarni, Makarand↗
Sensitivity analysis of frictional contact response of axisymmetric composite structures
A computational procedure is presented for evaluating the sensitivity coefficients of the static frictional contact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of nodal displacements, and Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for the frictional contact of a composite spherical cap pressed against a rigid plate.
Karaoglan, I.↗