SEARCH · Engineering Papers

Results for “singular equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Features of sound propagation through and stability of a finite shear layer

The plane wave propagation, the stability, and the rectangular duct mode problems of a compressible, inviscid, linearly sheared, parallel, homogeneous flow are shown to be governed by Whittaker's equation. The exact solutions for the perturbation quantities are essentially the Whittaker M-functions where the nondimensional quantities have precise physical meanings. A number of known results are obtained as limiting cases of the exact solutions. For the compressible finite thickness shear layer it is shown that no resonances and no critical angles exist for all Mach numbers, frequencies, and shear layer velocity profile slopes except in the singular case of the vortex sheet.

Koutsoyannis, S. P.↗

DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

NASA NTRS1973

Classical eighth- and lower-order Runge-Kutta-Nystroem formulas with a new stepsize control procedure for special second-order differential equations

New Runge-Kutta-Nystrom formulas of the eighth, seventh, sixth, and fifth order are derived for the special second-order (vector) differential equation x = f (t,x). In contrast to Runge-Kutta-Nystrom formulas of an earlier NASA report, these formulas provide a stepsize control procedure based on the leading term of the local truncation error in x. This new procedure is more accurate than the earlier Runge-Kutta-Nystrom procedure (with stepsize control based on the leading term of the local truncation error in x) when integrating close to singularities. Two central orbits are presented as examples. For these orbits, the accuracy and speed of the formulas of this report are compared with those of Runge-Kutta-Nystrom and Runge-Kutta formulas of earlier NASA reports.

Fehlberg, E.↗

NASA NTRS1977

The geopotential in nonsingular orbital elements

Singularities are eliminated from the geopotential and its partial derivatives for zero eccentricity and inclination, and an expression for the geopotential expansion based entirely on nonsingular orbital elements is developed. The argument relies on the treatments of the geopotential function given by Izsak (1964), Allan (1965) and Kaula (1966); the geopotential expansion developed does not involve mixed variables, and therefore does not require the chain rule to formulate the Lagrange planetary equations. The need for recursion relations to be used in conjunction with the nonsingular version of the geopotential expansion is also mentioned.

Nacozy, P. E.↗

DOE OSTI2020

Hydrodynamics of two-dimensional compressible fluid with broken parity: Variational principle and free surface dynamics in the absence of dissipation

In this paper, we consider an isotropic compressible nondissipative fluid with broken parity subject to free surface boundary conditions in two spatial dimensions. The hydrodynamic equations describing the bulk dynamics of the fluid and the free surface boundary conditions depend explicitly on the parity-breaking nondissipative odd viscosity term. We construct an effective action which gives both bulk hydrodynamic equations and free surface boundary conditions. The free surface boundary conditions require an additional boundary term in the action which resembles a 1+1D chiral boson field coupled to the background geometry. We solve the linearized hydrodynamic equations for the deep water case and derive the dispersion of chiral surface waves. We show that in the long-wavelength limit the flow profile exhibits an oscillating vortical boundary layer near the free surface. The layer thickness is controlled by the ratio between the odd viscosity (ν o ) and the sound velocity (c s ), δ ~ ν o /c s . In the incompressible limit, c s → ∞, the vortical boundary layer becomes singular with the vorticity within the layer diverging as ω ~ c s . The boundary layer is formed by odd viscosity coupling the divergence of velocity ∇ • v to vorticity ∇ × v. It results in nontrivial chiral free surface dynamics even in the absence of external forces. The structure of the odd-viscosity-induced boundary layer is very different from the conventional free surface boundary layer associated with dissipative shear viscosity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1982

Power series expansions for the frequency and period of the limit cycle of the van der Pol equation

An equation reported by van der Pol (1926) in connection with relaxation-oscillations studies is considered. The equation contains the factor epsilon which can assume values in the range from zero to infinity. The period T(epsilon), or equivalently the frequency nu(epsilon) of the limit cycle has been studied. However, to date there has been little success in discovering the analytical structure of T(epsilon) as a function of epsilon. The present investigation has the objectives to present the Taylor series expansion of nu(epsilon), to locate the singularities which determine the radius of convergence of that expansion, to introduce a new damping variable in terms of which the expansion converges for all epsilon, to form a new expansion for the period T(epsilon) which improves the rate of convergence of the series, to attempt to 'complete' the series, and to compare the obtained results with the numerically determined values of T(epsilon) and with the asymptotic approximation valid for large epsilon.

Andersen, C. M.↗

NASA NTRS1992

Determination of stress intensity factors for interface cracks under mixed-mode loading

A simple technique was developed using conventional finite element analysis to determine stress intensity factors, K1 and K2, for interface cracks under mixed-mode loading. This technique involves the calculation of crack tip stresses using non-singular finite elements. These stresses are then combined and used in a linear regression procedure to calculate K1 and K2. The technique was demonstrated by calculating three different bimaterial combinations. For the normal loading case, the K's were within 2.6 percent of an exact solution. The normalized K's under shear loading were shown to be related to the normalized K's under normal loading. Based on these relations, a simple equation was derived for calculating K1 and K2 for mixed-mode loading from knowledge of the K's under normal loading. The equation was verified by computing the K's for a mixed-mode case with equal and normal shear loading. The correlation between exact and finite element solutions is within 3.7 percent. This study provides a simple procedure to compute K2/K1 ratio which has been used to characterize the stress state at the crack tip for various combinations of materials and loadings. Tests conducted over a range of K2/K1 ratios could be used to fully characterize interface fracture toughness.

Naik, Rajiv A.↗

NASA NTRS1993

Nested Krylov methods and preserving the orthogonality

Recently the GMRESR inner-outer iteraction scheme for the solution of linear systems of equations was proposed by Van der Vorst and Vuik. Similar methods have been proposed by Axelsson and Vassilevski and Saad (FGMRES). The outer iteration is GCR, which minimizes the residual over a given set of direction vectors. The inner iteration is GMRES, which at each step computes a new direction vector by approximately solving the residual equation. However, the optimality of the approximation over the space of outer search directions is ignored in the inner GMRES iteration. This leads to suboptimal corrections to the solution in the outer iteration, as components of the outer iteration directions may reenter in the inner iteration process. Therefore we propose to preserve the orthogonality relations of GCR in the inner GMRES iteration. This gives optimal corrections; however, it involves working with a singular, non-symmetric operator. We will discuss some important properties, and we will show by experiments that, in terms of matrix vector products, this modification (almost) always leads to better convergence. However, because we do more orthogonalizations, it does not always give an improved performance in CPU-time. Furthermore, we will discuss efficient implementations as well as the truncation possibilities of the outer GCR process. The experimental results indicate that for such methods it is advantageous to preserve the orthogonality in the inner iteration. Of course we can also use iteration schemes other than GMRES as the inner method; methods with short recurrences like GICGSTAB are of interest.

Desturler, Eric↗

NASA NTRS1992

A study of the control problem of the shoot side environment delivery system of a closed crop growth research chamber

The details of our initial study of the control problem of the crop shoot environment of a hypothetical closed crop growth research chamber (CGRC) are presented in this report. The configuration of the CGRC is hypothetical because neither a physical subject nor a design existed at the time the study began, a circumstance which is typical of large scale systems control studies. The basis of the control study is a mathematical model which was judged to adequately mimic the relevant dynamics of the system components considered necessary to provide acceptable realism in the representation. Control of pressure, temperature, and flow rate of the crop shoot environment, along with its oxygen, carbon dioxide, and water concentration is addressed. To account for mass exchange, the group of plants is represented in the model by a source of oxygen, a source of water vapor, and a sink for carbon dioxide. In terms of the thermal energy exchange, the group of plants is represented by a surface with an appropriate temperature. Most of the primitive equations about an experimental operating condition and a state variable representation which was extracted from the linearized equations are presented. Next, we present the results of a real Jordan decomposition and the repositioning of an undesirable eigenvalue via full state feedback. The state variable representation of the modeling system is of the nineteenth order and reflects the eleven control variables and eight system disturbances. Five real eigenvalues are very near zero, with one at zero, three having small magnitude positive values, and one having a small magnitude negative value. A Singular Value Decomposition analysis indicates that these non-zero eigenvalues are not results of numerical error.

Blackwell, C. C.↗

DOE OSTI2020

Numerical methods for nonlocal and fractional models

Partial differential equations (PDEs) are used with huge success to model phenomena across all scientific and engineering disciplines. However, across an equally wide swath, there exist situations in which PDEs fail to adequately model observed phenomena, or are not the best available model for that purpose. On the other hand, in many situations, nonlocal models that account for interaction occurring at a distance have been shown to more faithfully and effectively model observed phenomena that involve possible singularities and other anomalies. Here, we consider a generic nonlocal model, beginning with a short review of its definition, the properties of its solution, its mathematical analysis and of specific concrete examples. We then provide extensive discussions about numerical methods, including finite element, finite difference and spectral methods, for determining approximate solutions of the nonlocal models considered. In that discussion, we pay particular attention to a special class of nonlocal models that are the most widely studied in the literature, namely those involving fractional derivatives. The article ends with brief considerations of several modelling and algorithmic extensions, which serve to show the wide applicability of nonlocal modelling.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Annihilation-to-nothing: a quantum gravitational boundary condition for the Schwarzschild black hole

The interior of a static Schwarzschild metric can be written in terms of two functions, similar to some models of anisotropic cosmology. With a suitable choice of canonical variables, we solve the Wheeler-DeWitt equation (WDW) inside the horizon of a Schwarzschild black hole. By imposing classicality near the horizon, and requiring boundedness of the wave function, we get a rather generic solution of the WDW equation, whose steepest-descent solution, i.e., the ridge of the wave function, coincides nicely with the classical trajectory. However, there is an ambiguity in defining the arrow of time which leads to two possible interpretations—(i) if there is only one arrow of time, one can infer that the steepest-descent of the wave function follows the classical trajectory throughout: coming from the event horizon and going all the way down to the singularity, while (ii) if there are two different arrows of time in two separate regimes, it can be inferred that the steepest-descent of the wave function comes inwards from the event horizon in one region while it moves outwards from the singularity in the other region, and there exists an annihilation process of these two parts of the wave function inside the horizon. Additionally, adopting the second interpretation could shed light on the information loss paradox: as time goes on, probabilities for histories that include black holes and singularities decay to zero and eventually only trivial geometries dominate.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1976

The development of the Poincare-similar elements with eccentric anomaly as the independent variable

A new set of element differential equations for the perturbed two-body motion is derived. The elements are canonical and are similar to the classical canonical Poincare elements, which have time as the independent variable. The phase space is extended by introducing the total energy and time as canonically conjugated variables. The new independent variable is, to within an additive constant, the eccentric anomaly. These elements are compared to the Kustaanheimo-Stiefel (KS) element differential equations, which also have the eccentric anomaly as the independent variable. For several numerical examples, the accuracy and stability of the new set are equal to those of the KS solution. This comparable accuracy result can probably be attributed to the fact that both sets have the same time element and very similar energy elements. The new set has only 8 elements, compared to 10 elements for the KS set. Both sets are free from singularities due to vanishing eccentricity and inclination.

Bond, V. R.↗

NASA NTRS1979

Investigations on the hierarchy of reference frames in geodesy and geodynamics

Problems related to reference directions were investigated. Space and time variant angular parameters are illustrated in hierarchic structures or towers. Using least squares techniques, model towers of triads are presented which allow the formation of linear observation equations. Translational and rotational degrees of freedom (origin and orientation) are discussed along with and the notion of length and scale degrees of freedom. According to the notion of scale parallelism, scale factors with respect to a unit length are given. Three-dimensional geodesy was constructed from the set of three base vectors (gravity, earth-rotation and the ecliptic normal vector). Space and time variations are given with respect to a polar and singular value decomposition or in terms of changes in translation, rotation, deformation (shear, dilatation or angular and scale distortions).

Grafarend, E. W.↗

NASA NTRS1986

Taylor-Goertler instabilities of Tollmien-Schlichting waves and other flows governed by the interactive boundary-layer equations

The Taylor-Goertler vortex instability equations are formulated for steady and unsteady interacting boundary-layer flows. The effective Goertler number is shown to be a function of the wall shape in the boundary layer and the possibility of both steady and unsteady Taylor-Goertler modes exists. As an example the steady flow in a symmetrically constricted channel is considered and it is shown that unstable Goertler vortices exist before the boundary layers at the wall develop the Goldstein singularity discussed by Smith and Daniels (1981). As an example of an unsteady spatially varying basic state, it is considered the instability of high-frequency large-amplitude two- and three-dimensional Tollmien-Schlichting waves in a curved channel. It is shown that they are unstable in the first 'Stokes-layer stage' of the hierarchy of nonlinear states discussed by Smith and Burggraf (1985). This instability of Tollmien-Schlichting waves in an internal flow can occur in the presence of either convex or concave curvature. Some discussion of this instability in external flows is given.

Hall, Philip↗

NASA NTRS1987

Analytical caustic surfaces

This document discusses the determination of caustic surfaces in terms of rays, reflectors, and wavefronts. Analytical caustics are obtained as a family of lines, a set of points, and several types of equations for geometries encountered in optics and microwave applications. Standard methods of differential geometry are applied under different approaches: directly to reflector surfaces, and alternatively, to wavefronts, to obtain analytical caustics of two sheets or branches. Gauss/Seidel aberrations are introduced into the wavefront approach, forcing the retention of all three coefficients of both the first- and the second-fundamental forms of differential geometry. An existing method for obtaining caustic surfaces through exploitation of the singularities in flux density is examined, and several constant-intensity contour maps are developed using only the intrinsic Gaussian, mean, and normal curvatures of the reflector. Numerous references are provided for extending the material of the present document to the morphologies of caustics and their associated diffraction patterns.

Schmidt, R. F.↗

NASA NTRS1990

Elliptic grid generation with orthogonality and spacing control on an arbitrary number of boundaries

A procedure for the generation of two and quasi-three-dimensional grids with control of orthogonality and spacing with respect to any and/or all boundaries of the domain is described. The elliptic grid generation equations of Thompson are solved implicitly. Control of the grid behavior is achieved through the introduction of forcing functions terms in the manner of Steger and Sorenson or in a modification of the method of Hilgenstock. The forcing function terms are constructed on the boundaries and propagated into the domain using transfinite Lagrangian bivariate interpolation. An anisotropic transfinite stencil is introduced and is shown to produce excellent grid behavior particularly in the vicinity of corner singularities. Emphasis is placed on the generation of viscous grids and the method is shown to be suited for use in the generation of grids for internal as well as external flow geometries. A FORTRAN program named PISCES has been written to implement the algorithm. Examples of grids for internal and external flows are given that highlight the characteristics and behavior of the algorithm.

White, J. A.↗

NASA NTRS1990

Bifurcations and sensitivity in parametric nonlinear programming

The parametric nonlinear programming problem is that of determining the behavior of solution(s) as a parameter or vector of parameters alpha belonging to R(sup r) varies over a region of interest for the problem: Minimize over x the set f(x, alpha):h(x, alpha) = 0, g(x, alpha) is greater than or equal to 0, where f:R(sup (n+r)) approaches R, h:R(sup (n+r)) approaches R(sup q) and g:R(sup (n+r)) approaches R(sup p) are assumed to be at least twice continuously differentiable. Some of these parameters may be fixed but not known precisely and others may be varied to enhance the performance of the system. In both cases a fundamentally important problem in the investigation of global sensitivity of the system is to determine the stability boundaries of the regions in parameter space which define regions of qualitatively similar solutions. The objective is to explain how numerical continuation and bifurcation techniques can be used to investigate the parametric nonlinear programming problem in a global sense. Thus, first the problem is converted to a closed system of parameterized nonlinear equations whose solution set contains all local minimizers of the original problem. This system, which will be represented as F(z,alpha) = O, will include all Karush-Kuhn-Tucker and Fritz John points, both feasible and infeasible solutions, and relative minima, maxima, and saddle points of the problem. The local existence and uniqueness of a solution path (z(alpha), alpha) of this system as well as the solution type persist as long as a singularity in the Jacobian D(sub z)F(z,alpha) is not encountered. Thus the nonsingularity of this Jacobian is characterized in terms of conditions on the problem itself. Then, a class of efficient predictor-corrector continuation procedures for tracing solution paths of the system F(z,alpha) = O which are tailored specifically to the parametric programming problem are described. Finally, these procedures and the obtained information are illustrated within the context of design optimization.

Lundberg, Bruce N.↗

NASA NTRS2016

Diagonally Implicit Runge-Kutta Methods for Ordinary Differential Equations. A Review

A review of diagonally implicit Runge-Kutta (DIRK) methods applied to rst-order ordinary di erential equations (ODEs) is undertaken. The goal of this review is to summarize the characteristics, assess the potential, and then design several nearly optimal, general purpose, DIRK-type methods. Over 20 important aspects of DIRKtype methods are reviewed. A design study is then conducted on DIRK-type methods having from two to seven implicit stages. From this, 15 schemes are selected for general purpose application. Testing of the 15 chosen methods is done on three singular perturbation problems. Based on the review of method characteristics, these methods focus on having a stage order of two, sti accuracy, L-stability, high quality embedded and dense-output methods, small magnitudes of the algebraic stability matrix eigenvalues, small values of aii, and small or vanishing values of the internal stability function for large eigenvalues of the Jacobian. Among the 15 new methods, ESDIRK4(3)6L[2]SA is recommended as a good default method for solving sti problems at moderate error tolerances.

Kennedy, Christopher A.↗