Engineering Papers⌕ Search

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.

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

At least 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Improved Discrete Ordinate Solutions in the Presence of an Anisotropically Reflecting Lower Boundary: Upgrades of the DISORT Computational Tool

A successor version 3 of DISORT (DISORT3) is presented with important upgrades that improve the accuracy, efficiency, and stability of the algorithm. Compared with version 2 (DISORT2 released in 2000) these upgrades include (a) a redesigned BRDF computation that improves both speed and accuracy, (b) a revised treatment of the single scattering correction, and (c) additional efficiency and stability upgrades for beam sources. In DISORT3 the BRDF computation is improved in the following three ways: (i) the Fourier decomposition is prepared "off-line", thus avoiding the repeated internal computations done in DISORT2; (ii) a large enough number of terms in the Fourier expansion of the BRDF is employed to guarantee accurate values of the expansion coefficients (default is 200 instead of 50 in DISORT2); (iii) in the post processing step the reflection of the direct attenuated beam from the lower boundary is included resulting in a more accurate single scattering correction. These improvements in the treatment of the BRDF have led to improved accuracy and a several-fold increase in speed. In addition, the stability of beam sources has been improved by removing a singularity occurring when the cosine of the incident beam angle is too close to the reciprocal of any of the eigenvalues. The efficiency for beam sources has been further improved from reducing by a factor of 2 (compared to DISORT2) the dimension of the linear system of equations that must be solved to obtain the particular solutions, and by replacing the LINPAK routines used in DISORT2 by LAPACK 3.5 in DISORT3. These beam source stability and efficiency upgrades bring enhanced stability and an additional 5-7% improvement in speed. Numerical results are provided to demonstrate and quantify the improvements in accuracy and efficiency of DISORT3 compared to DISORT2.

Lin, Z.↗
NASA NTRS1971

The problem of an elastic stiffener bonded to a half plane.

The contact problem of an elastic stiffener bonded to an elastic half plane with different mechanical properties is considered. The governing integral equation is reduced to an infinite system of linear algebraic equations. It is shown that, depending on the value of a parameter which is a function of the elastic constants and the thickness of the stiffener, the system is either regular or quasi-regular. A complete numerical example is given for which the strength of the stress singularity and the contact stresses are tabulated.

Erdogan, F.↗
NASA NTRS1978

Stress intensity factors in a cracked infinite elastic wedge loaded by a rigid punch

A plane elastic wedge-shaped solid was split through the application of a rigid punch. It was assumed that the coefficient of friction on the the contact area was constant, and the problem had a plane of symmetry with respect to loading and geometry, with the crack in the plane of symmetry. The problem was formulated in terms of a system of integral equations with the contact stress and the derivative of the crack surface displacement as the unknown functions. The solution was obtained for an internal crack and for an edge crack. The results include primarily the stress intensity factors at the crack tips, and the measure of the stress singularity at the wedge apex, and at the end points of the contact area.

Erdogan, F.↗
NASA NTRS1980

A general panel method for the analysis and design of arbitrary configurations in incompressible flows

A method for solving the linear integral equations of incompressible potential flow in three dimensions is presented. Both analysis (Neumann) and design (Dirichlet) boundary conditions are treated in a unified approach to the general flow problem. The method is an influence coefficient scheme which employs source and doublet panels as boundary surfaces. Curved panels possessing singularity strengths, which vary as polynomials are used, and all influence coefficients are derived in closed form. These and other features combine to produce an efficient scheme which is not only versatile but eminently suited to the practical realities of a user-oriented environment. A wide variety of numerical results demonstrating the method is presented.

Johnson, F. T.↗
NASA NTRS1985

Finite-volume scheme for transonic potential flow about airfoils and bodies in an arbitrarily-shaped channel

A conservative finite-volume difference scheme is developed for the potential equation to solve transonic flow about airfoils and bodies in an arbitrary channel. The scheme employs a mesh which is a nearly-conformal 'O' mesh about the airfoil and nearly orthogonal at the channel walls. The mesh extends to infinity upstream and downstream, where the mapping is singular. Special procedures are required to treat the singularities at infinity, including computation of the metrics near those points. Channels with exit areas different from inlet areas are solved; a body with a sting mount is an example of such a case.

South, J. C., Jr.↗
NASA NTRS1986

Finite-volume scheme for transonic potential flow about airfoils and bodies in an arbitrarily shaped channel

A conservative finite-volume difference scheme is developed for the potential equation to solve transonic flow about airfoils and bodies in an arbitrarily shaped channel. The scheme employs a mesh which is a nearly conformal O mesh about the airfoil and nearly orthogonal at the channel walls. The mesh extends to infinity upstream and downstream, where the mapping is singular. Special procedures are required to treat the singularities at infinity, including computation of the metrics near those points. Channels with exit areas different from inlet areas are solved; a body with a sting mount is an example of such a case.

South, Jerry C., Jr.↗
NASA NTRS1993

Nonlinear vortex trail dynamics. II - Analytic solutions

Spatially periodic large amplitude solutions of the von Karman model are obtained in the neighborhood of singularities. These singularities correspond to vortex clusters in the physical plane. The quasi-periodic and unbounded solutions found analytically confirm earlier numerical work and show qualitative agreement with experimental observations of large-scale phenomena of vortex trails. Separatrices or heteroclinic orbits were explicitly found for an integrable approximate equation, which indicate that the von Karman model itself supports chaotic solutions.

Lim, Chjan C.↗
NASA NTRS1989

Crack-face displacements for embedded elliptic and semi-elliptical surface cracks

Analytical expressions for the crack-face displacements of an embedded elliptic crack in infinite solid subjected to arbitrary tractions are obtained. The tractions on the crack faces are assumed to be expressed in a polynomial form. These displacements expressions complete the exact solution of Vijayakumar and Atluri, and Nishioki and Atluri. For the special case of an embedded crack in an infinite solid subjected to uniform pressure loading, the present displacements agree with those by Green and Sneddon. The displacement equations derived were used with the finite-element alternating method (FEAM) for the analysis of a semi-elliptic surface crack in a finite solid subjected to remote tensile loading. The maximum opening displacements obtained with FEAM are compared to those with the finite-element method with singularity elements. The maximum crack opening displacements by the two methods showed good agreement.

Raju, I. S.↗
NASA NTRS1990

Crack-face displacements for an embedded elliptic crack in an infinite solid subjected to arbitrary tractions

Analytical expression for the crack-face displacements of an embedded elliptic crack in infinite solid subjected to arbitrary traction are obtained. The tractions on the crack faces are assumed to be expressed in a polynomial form. These displacements expressions complete the exact solution of Vijayakumar and Atluri, and Nishioki and Atluri. For the special case of an embedded crack in an infinite solid subjected to uniform pressure loading, the present displacements agree with those by Green and Sneddon. The displacement equations derived were used with the finite-element alternating method (FEAM) for the analysis of a semi-elliptic surface crack in a finite solid subjected to remote tensile loading. The maximum opening displacements obtained with FEAM are compared to those with the finite-element method with singularity elements. The maximum crack opening displacements by the two methods showed good agreement.

Raju, I. S.↗
NASA NTRS1990

Scattering of impact wave by a crack in composite plate

The surface responses due to impact load on an infinite uniaxial graphite/epoxy plate with and without delamination cracks are investigated both in time and frequency domain by using a hybrid method combining the finite element discretization of the near-field with boundary integral representation of the field outside a contour enclosing completely the crack. This combined method leads to a set of linear unsymmetric complex matrix equations which are solved to obtain the response in the frequency domain by biconjugate gradient method. The time domain response is then obtained by using an FFT. In order to capture the time-domain characteristics accurately, high order finite elements have been used. Also, both the six node singular elements and eight node transition elements are used around the crack tips to model the crack-tip singularity. It is shown that from the numerical results for surface responses both depth and length of this crack can be identified.

Ju, T. H.↗
NASA NTRS1994

Tensor-GMRES method for large sparse systems of nonlinear equations

This paper introduces a tensor-Krylov method, the tensor-GMRES method, for large sparse systems of nonlinear equations. This method is a coupling of tensor model formation and solution techniques for nonlinear equations with Krylov subspace projection techniques for unsymmetric systems of linear equations. Traditional tensor methods for nonlinear equations are based on a quadratic model of the nonlinear function, a standard linear model augmented by a simple second order term. These methods are shown to be significantly more efficient than standard methods both on nonsingular problems and on problems where the Jacobian matrix at the solution is singular. A major disadvantage of the traditional tensor methods is that the solution of the tensor model requires the factorization of the Jacobian matrix, which may not be suitable for problems where the Jacobian matrix is large and has a 'bad' sparsity structure for an efficient factorization. We overcome this difficulty by forming and solving the tensor model using an extension of a Newton-GMRES scheme. Like traditional tensor methods, we show that the new tensor method has significant computational advantages over the analogous Newton counterpart. Consistent with Krylov subspace based methods, the new tensor method does not depend on the factorization of the Jacobian matrix. As a matter of fact, the Jacobian matrix is never needed explicitly.

Feng, Dan↗
NASA NTRS1979

Stability of the nakedness of Weyl singularities

The stability of the nakedness of the Weyl singularities against matter perturbations is investigated. Consideration is given to the effects of infalling test matter on the convergence of outgoing null rays. It is shown that the additional convergence induced by infalling test matter does not blow up sufficiently fast to reconverge diverging outgoing rays, at least in the equator, and that the nakedness seems to be stable in this limited sense.

Haugan, M. P.↗
NASA NTRS1992

A numerical method for solving systems of linear ordinary differential equations with rapidly oscillating solutions

The present numerical method for accurate and efficient solution of systems of linear equations proceeds by numerically developing a set of basis solutions characterized by slowly varying dependent variables. The solutions thus obtained are shown to have a computational overhead largely independent of the small size of the scale length which characterizes the solutions; in many cases, the technique obviates series solutions near singular points, and its known sources of error can be easily controlled without a substantial increase in computational time.

Bernstein, Ira B.↗
NASA NTRS2000

Higher Order Modeling In the BEM/FEM Hybrid Formulation

Hybrid formulations using low order curl-conforming bases to represent the total electric field within a finite element region and low order divergence-conforming bases to represent equivalent electric and magnetic currents on the boundary are well known. However, higher-order divergence and curl-conforming bases have been shown to provide significant benefits in convergence rates and accuracy when employed in strictly integral equation and strictly finite element formulations. In this paper, a hybrid electric field formulation employing higher order bases is presented. The paper addresses benefits and issues associated with using higher order divergence-and curl-conforming bases in the hybrid finite element/boundary element electric field formulation. The method of singularity subtraction may be used to compute the self terms of the boundary integral when the bases are of low order. But this method becomes laborious and requires great care when the divergence conforming bases are of higher order. In order to handle these singularities simply and accurately, a generalized Gaussian quadrature method is employed in which the expansion functions account for the singularity. In preliminary tests of the higher order hybrid formulation, the equivalent electric current induced by scattering of a plane wave from a square dielectric cylinder is examined. Accurate results are obtained using only a two-triangle mesh when the current basis is of order 4 or 5. Additional results are presented comparing the error obtained using higher order bases to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained with higher order bases are compared to those obtained with lower order bases for selected sample problems.

Fink, Patrick W.↗
NASA NTRS2003

Collective Excitations in InAs Well Intersubband Transitions

Intersubband transitions in semiconductor quantum well are studied using a density matrix theory that goes beyond the Hartree-Fock approximation by including the full second order electron-electron scattering terms in the polarization equation for the first time. Even though the spectral features remain qualitatively similar to the results obtained with dephasing rate approximation, significant quantitative changes result from such a more detailed treatment of dephasing. More specifically, we show how the interplay of the two fundamental collective excitations, the Fermi-edge singularity and the intersubband plasmon, leads to significant changes in lineshape as the electron density varies.

Li, Jian-Zhong↗
NASA NTRS1979

Solitary waves in stratified shear flows

An analysis is described of long, finite-amplitude internal waves in a stratified shear flow. Both regular and singular modes are considered with a nonlinear critical layer employed in the latter case. A three-layer model is used to develop the theory and closed-form expressions are obtained relating the phase speed to the Richardson number, the latter quantity being taken as O(1). The amplitude evolution equation is found to be either the Korteweg-de Vries equation or the Benjamin-Davis-Ono equation depending upon the distance of the more remote boundary from the edge of the shear layer.

Maslowe, S. A.↗
NASA NTRS1991

Three-dimensional kinematic reconnection of plasmoids

The kinematic reconnection model of Lau and Finn (1990) is applied in a theoretical investigation of three-dimensional plasmoid morphology, as seen in the solar corona and earth magnetotail. The derivation of the governing equations is outlined; long, short, and periodic plasmoid models are developed; and the evolution of the stable and unstable manifolds in these models is presented graphically. It is inferred that sheet currents and tangential discontinuities can form on surfaces topologically identical to those where Delta(phi) about equal to Delta(z) singularities occur in the kinematic reconnection model, and can be broadened in a similar way by nonideal effects. Two such surfaces exist in long plasmoids, and also in short plasmoids in the presence of finite resistivity; the intertwined surfaces characteristic of the periodic plasmoid form a fractal set but merge in the presence of finite resistivity, producing structures similar to those proposed by the sheet-current theory of Parker (1983).

Lau, Yun-Tung↗