Search NASA⌕ Search

SEARCH · Search NASA

Results for “SINGULAR INTEGRAL EQUATION”

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 253 records · Page 14

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.↗

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.↗

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.↗

Numerical solution methods for viscoelastic orthotropic materials

Numerical solution methods for viscoelastic orthotropic materials, specifically fiber reinforced composite materials, are examined. The methods include classical lamination theory using time increments, direction solution of the Volterra Integral, Zienkiewicz's linear Prony series method, and a new method called Nonlinear Differential Equation Method (NDEM) which uses a nonlinear Prony series. The criteria used for comparison of the various methods include the stability of the solution technique, time step size stability, computer solution time length, and computer memory storage. The Volterra Integral allowed the implementation of higher order solution techniques but had difficulties solving singular and weakly singular compliance function. The Zienkiewicz solution technique, which requires the viscoelastic response to be modeled by a Prony series, works well for linear viscoelastic isotropic materials and small time steps. The new method, NDEM, uses a modified Prony series which allows nonlinear stress effects to be included and can be used with orthotropic nonlinear viscoelastic materials. The NDEM technique is shown to be accurate and stable for both linear and nonlinear conditions with minimal computer time.

Gramoll, K. C.↗

Radiative transfer in the electrode region of an arc heater

Radiation from hot gases in the electrode region of an arc heater are discussed. Major original contributions are made: (1) an analytical expression for the radiative heat flux, for an axisymmetric geometry and a grey reflecting wall, is derived, (2) the expression is rendered analytically differentiable/integrable by using appropriate techniques to remove the singularities, and (3), finally a two band model is devised to compute the radiative heat load at the walls in NASA Ames' jet facility. The advantage of using these techniques is that the time taken for the numerical integration is considerably reduced and computation of a 39 x 39 axisymmetric geometry becomes manageable. Fluid equations, representative of a high current arc heater, are numerically solved for an axisymmetric geometry. The computed average heat load of 400 Watts/sq cm using the two band model agrees well with the experimental estimates made at NASA Ames' facility.

Durgapal, P.↗

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.↗

Lunar and solar perturbations on the orbit of a geosynchronous satellite

The luni-solar effects cause a large amplitude, long-period perturbation of the orbital plane. Canonical differential equations associated with this motion contain a singular point, and are expanded about this point to third order. A solution is given that is valid for long times and is not restricted to small inclinations. Higher order terms are investigated. Expressions are given for the inclination (I) and node as a function of time. Comparisons with a numerically integrated solution show a disagreement in inclination of only .03 deg after 11 years.

Graf, O. F., Jr.↗

Finite elements: Theory and application

Recent advances in FEM techniques and applications are discussed in reviews and reports presented at the ICASE/LaRC workshop held in Hampton, VA in July 1986. Topics addressed include FEM approaches for partial differential equations, mixed FEMs, singular FEMs, FEMs for hyperbolic systems, iterative methods for elliptic finite-element equations on general meshes, mathematical aspects of FEMS for incompressible viscous flows, and gradient weighted moving finite elements in two dimensions. Consideration is given to adaptive flux-corrected FEM transport techniques for CFD, mixed and singular finite elements and the field BEM, p and h-p versions of the FEM, transient analysis methods in computational dynamics, and FEMs for integrated flow/thermal/structural analysis.

Dwoyer, D. L.↗

Elastic interactions of a fatigue crack with a micro-defect by the mixed boundary integral equation method

In this paper, the mixed boundary integral equation method is developed to study the elastic interactions of a fatigue crack and a micro-defect such as a void, a rigid inclusion or a transformation inclusion. The method of pseudo-tractions is employed to study the effect of a transformation inclusion. An enriched element which incorporates the mixed-mode stress intensity factors is applied to characterize the singularity at a moving crack tip. In order to evaluate the accuracy of the numerical procedure, the analysis of a crack emanating from a circular hole in a finite plate is performed and the results are compared with the available numerical solution. The effects of various micro-defects on the crack path and fatigue life are investigated. The results agree with the experimental observations.

Lua, Yuan J.↗

Stress-intensity factors for internal surface cracks in cylindrical pressure vessels

The paper presents stress-intensity factors for a wide range of semi-elliptical surface cracks on the inside of pressurized cylinders. The ratio of crack depth to crack length ranged from 0.2 to 1; the ratio of crack depth to wall thickness ranged from 0.2 to 0.8; and the ratio of wall thickness to vessel radius was 0.1 to to 0.25. The stress-intensity factors were calculated by a three-dimensional finite-element method using singularity elements along the crack front and linear-strain elements elsewhere. An equation for the stress-intensity factors was obtained which applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was also made between the results and other analyses of internal surface cracks in cylinders. The results from a boundary-integral equation method were in agreement and those from another finite-element method were in fair agreement (+ or - 8 percent) with the results.

Newman, J. C., Jr.↗

A singularity free analytical solution of artificial satellite motion with drag

The connection between the existing Delaunay-Similar and Poincare-Similar satellite theories in the true anomaly version is outlined for the J(2) perturbation and the new drag approach. An overall description of the concept of the approach is given while the necessary expansions and the procedure to arrive at the computer program for the canonical forces is delineated. The procedure for the analytical integration of these developed equations is described. In addition, some numerical results are given. The computer program for the algebraic multiplication of the Fourier series which creates the FORTRAN coding in an automatic manner is described and documented.

Scheifele, G.↗

Role of pressure diffusion in non-homogeneous shear flows

A non-local model is presented for approximating the pressure diffusion in calculations of turbulent free shear and boundary layer flows. It is based on the solution of an elliptic relaxation equation which enables local diffusion sources to be distributed over lengths of the order of the integral scale. The pressure diffusion model was implemented in a boundary layer code within the framework of turbulence models based on both the kappa-epsilon-(bar)upsilon(exp 2) system of equations and the full Reynolds stress equations. Model computations were performed for mixing layers and boundary layer flows. In each case, the pressure diffusion model enabled the well-known free-stream edge singularity problem to be eliminated. There was little effect on near-wall properties. Computed results agreed very well with experimental and DNS data for the mean flow velocity, the turbulent kinetic energy, and the skin-friction coefficient.

Demuren, A. O.↗

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.↗

A robust direct-integration method for rotorcraft maneuver and periodic response

The Newmark-Beta method and the Newton-Raphson iteration scheme are combined to develop a direct-integration method for evaluating the maneuver and periodic-response expressions for rotorcraft. The method requires the generation of Jacobians and includes higher derivatives in the formulation of the geometric stiffness matrix to enhance the convergence of the system. The method leads to effective convergence with nonlinear structural dynamics and aerodynamic terms. Singularities in the matrices can be addressed with the method as they arise from a Lagrange multiplier approach for coupling equations with nonlinear constraints. The method is also shown to be general enough to handle singularities from quasisteady control-system models. The method is shown to be more general and robust than the similar 2GCHAS method for analyzing rotorcraft dynamics.

Panda, Brahmananda↗

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.↗

Energy transfer and constrained simulations in isotropic turbulence

The defining characteristic of turbulent flows is their ability to dissipate energy, even in the limit of zero viscosity. The Euler equations, if constrained in such a way that the velocity derivatives remain bounded, conserve energy. But when they arise as the limit of the Navier-Stokes (NS) equations, when the Reynolds number goes to infinity, there is persuasive empirical evidence that the gradients become singular as just the right function of Re for the dissipation to remain non-zero and to approach a well defined limit. It is generally believed that this limiting value of the dissipation is a property of the Euler equations themselves, independent of the particular dissipative mechanism involved, and that it can be normalized with the large scale properties of the turbulent flow (e.g. the kinetic energy per unit volume u'(exp 2)/2, and the integral scale L) without reference to the Reynolds number or to other dissipative quantities. This is usually taken to imply that the low wave number end of the energy spectrum, far from the dissipative range, is also independent of the particular mechanism chosen to dispose of the energy transfer. In the following sections, we present some numerical experiments on the effect of substituting different dissipation models into the truncated Euler equations. We will see that the effect is mainly felt in the 'near dissipation' range of the energy spectrum, but that this range can be quite wide in some cases, contaminating a substantial range of wave numbers. In the process, we will develop a 'practical' approximation to the subgrid energy transfer in isotropic turbulence, and we will gain insight into the structure of the nonlinear interactions among turbulent scales of comparable size, and into the nature of energy backscatter. Some considerations on future research directions are offered at the end.

Jimenez, Javier↗

Stress-intensity factors for internal surface cracks in cylindrical pressure vessels

The stress intensity factors were calculated by a three dimensional finite element method. The finite element models employed singularity elements along the crack front and linear strain elements elsewhere. The models had about 6500 degrees of freedom. The stress intensity factors were evaluated from a nodal force method. An equation for the stress intensity factors was obtained form the results of the present analysis. The equation applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was made between the present results and other analyses of internal surface cracks in cylinders. The results from a boundary integral equation method were in agreement (+ or - 2 percent) and those from another finite element were in fair agreement (+ or - 8 percent) with the present results.

Newman, J. C., Jr.↗

Differential Flatness and Cooperative Tracking in the Lorenz System

In this paper the control of the Lorenz system for both stabilization and tracking problems is studied via feedback linearization and differential flatness. By using the Rayleigh number as the control, only variable physically tunable, a barrier in the controllability of the system is incidentally imposed. This is reflected in the appearance of a singularity in the state transformation. Composite controllers that overcome this difficulty are designed and evaluated. The transition through the manifold defined by such a singularity is achieved by inducing a chaotic response within a boundary layer that contains it. Outside this region, a conventional feedback nonlinear control is applied. In this fashion, the authority of the control is enlarged to the whole. state space and the need for high control efforts is mitigated. In addition, the differential parametrization of the problem is used to track nonlinear functions of one state variable (single tracking) as well as several state variables (cooperative tracking). Control tasks that lead to integrable and non-integrable differential equations for the nominal flat output in steady-state are considered. In particular, a novel numerical strategy to deal with the non-integrable case is proposed. Numerical results validate very well the control design.

Crespo, Luis G.↗