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 541 records · Page 30

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Optimal Observations for Variational Data Assimilation

An important aspect of Ocean state estimation is the design of an observing system that allows the efficient study of climate aspects in the ocean. A solution of the design problem is presented here in terms of optimal observations that emerge as nondimensionalized singular vectors of the modified data resolution matrix. The actual computation is feasible only for scalar quantities in the limit of large observational errors. In the framework of a lo resolution North Atlantic primitive equation model it is demonstrated that such optimal observations when applied to determining the strength of the volume and heat transport across the Greenland-Scotland ridge, perform significantly better than traditional section data. On seasonal to inter-annual time-scales optimal observations are located primarily along the continental shelf and information about heat-transport, wind stress and stratification is being communicated via boundary waves and advective processes. On time-scales of about a month, sea surface height observations appear to be more efficient in reconstructing the cross-ridge heat transport than hydrographic observations. Optimal observations also provide a tool for understanding how the ocean state is effected by anomalies of integral quantities such as meridional heat transport.

Koehl, Armin↗

NASA NTRS1975

Steady subsonic flow around finite-thickness wings

The general method for analyzing steady subsonic potential aerodynamic flow around a lifting body having arbitrary shape is presented. By using the Green function method, an integral representation for the potential is obtained. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface of the body. Hence if the point P approaches the surface of the body, the representation reduces to an integral equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings. However, numerical results obtained for a rectangular wing in subsonic flow show that these problems do not appear even for thickness ratio tau = .001. Comparison with existing results shows that the proposed method is at least as fast and accurate as the lifting surface theories.

Kuo, C. C.↗

NASA NTRS1978

A singularity free analytical solution of artificial satellite motion with drag

An analytical satellite theory based on the regular, canonical Poincare-Similar (PS phi) elements is described along with an accurate density model which can be implemented into the drag theory. A computationally efficient manner in which to expand the equations of motion into a fourier series is discussed.

Mueller, A.↗

NASA NTRS1980

Stresses in a three-dimensional unidirectional composite containing broken fibers

An approximate solution is developed for the determination of the interlaminar normal and shear stresses in the vicinity of a crack in a three dimensional composite containing unidirectional linearly elastic fibers in an infinite linearly elastic matrix. In order to reduce the complexity of the formulation, certain assumptions are made as to the physically significant stresses to be retained. These simplifications reduce the partial differential equations of elasticity to differential-difference equations which are tractable using Fourier transform techniques. This 'material modeling' approach is in contrast with solutions developed by considering each lamina as a homogeneous, orthotropic layer. The resulting solution does not contain the classical singular stress field for the fibers and the influence of broken fibers on unbroken fibers is felt by a change in stress concentration factors. The matrix stresses however, are unbounded as the fiber spacing vanishes and an equivalent fiber-matrix geometry is proposed which gives the correct singular behavior. The numerical solution is considered in detail and several specific examples are presented. The potential for damaged or debonded zones to be generated by an embedded crack is discussed, and stress concentration factors for fibers near the crack are given. Detailed comparisons are made between the present solution, the analogous two-dimensional problem, and corresponding shear-lag models.

Goree, J. G.↗

NASA NTRS1985

Asymptotic matching by the symbolic manipulator MACSYMA

The delegation of the labor of calculating higher-order terms in singular perturbation (SP) expansions to a computer by the use of MACSYMA is considered. The method of matched asymptotic expansions is studied in detail for two model SP problems: a model resembling the boundary layer equation with a small parameter multiplying the highest derivatives; and a turning-point problem. It is shown that MACSYMA has successfully performed the higher-order matching in both problems.

Lo, L. L.↗

NASA NTRS1995

On the Nonlinear Evolution of a Stationary Cross-Flow Vortex in a Fully Three-Dimensional Boundary Layer Flow

We consider the nonlinear stability of a fully three-dimensional boundary layer flow in an incompressible fluid and derive an equation governing the nonlinear development of a stationary cross-flow vortex. The amplitude equation is a novel integro-differential equation which has spatial derivatives of the amplitude occurring in the kernal function. It is shown that the evolution of the cross-flow vortex is strongly coupled to the properties of an unsteady wall layer which is in fact driven by an unknown slip velocity, proportional to the amplitude of the cross-flow vortex. The work is extended to obtain the corresponding equation for rotating disk flow. A number of special cases are examined and the numerical solution for one of cases, and further analysis, demonstrates the existence of finite-distance as well as focussing type singularities. The numerical solutions also indicate the presence of a new type of nonlinear wave solution for a certain set of parameter values.

Gajjar, J. S. B.↗

NASA NTRS1982

Statics of wrinkling films

An analytical investigation of the equilibrium of wrinkling films is conducted. Zak (1979) has shown that wrinkling occurs in connection with the instability of a smooth film having no resistance to bending in the case of compression. The governing equation for the equilibrium of a film with possible regions of wrinkling is considered. The introduction of fictitious stress reduces the governing equation to a form which formally coincides with the governing equation for a string. Equilibrium conditions in the case of an absence of external forces are explored, taking into account the stretching of a semispherical film, the torsion of a convex film of revolution, and stress singularities. A study is conducted of the equilibrium under conditions in which external forces normal to the surface of a film are present. Attention is also given to the equilibrium in a potential field.

Zak, M.↗

NASA NTRS1987

Dusty cometary atmospheres

Governing equations describing the gas and dust production and the accelerating dusty gas flow in a cometary atmosphere are compiled and compared with the information available on the inner coma of the Halley comet. Cometary dusty gas flow calculations, based on the Probstein (1968) approach, yield a first order differential equation for the gas velocity with only one physical solution possible, which is the one of a transonic accelerating gas 'wind'. According to this model, the comet wind starts subsonically at the nucleus, goes through the singular sonic point, and then accelerates further.

Gombosi, Tamas I.↗

NASA NTRS1992

Confining potential in momentum space

A method is presented for the solution in momentum space of the bound state problem with a linear potential in r space. The potential is unbounded at large r leading to a singularity at small q. The singularity is integrable, when regulated by exponentially screening the r-space potential, and is removed by a subtraction technique. The limit of zero screening is taken analytically, and the numerical solution of the subtracted integral equation gives eigenvalues and wave functions in good agreement with position space calculations.

Norbury, John W.↗

NASA NTRS1993

On the 'delta-equations' for vortex sheet evolution

We use a set of equations, sometimes referred to as the 'delta-equations', to approximate the two-dimensional inviscid motion of an initially circular vortex sheet released from rest in a cross-flow. We present numerical solutions of these equations for the case with delta-square = 0 (for which the equations are exact) and for delta-square greater than 0. For small values of the smoothing parameter delta, a spectral filter must be used to eliminate spurious instabilities due to round-off error. Two singularities appear simultaneously in the vortex sheet when delta-square = 0 at a critical time t(c). After t(c), the solutions do not converge as the computational mesh is refined. With delta-square greater than 0, converged solutions were found for all values of delta-square when t is less than t(c), and for all but the two smallest values of delta-square used when t is greater than t(c). Our results show that, when delta-square is greater than 0, the vortex sheet deforms into two doubly branched spirals some time after t(c). The limiting solution as delta approaching 0 clearly exists and equals the delta = 0 solution when t is less than t(c).

Rottman, James W.↗

NASA NTRS1992

Confining potential in momentum space

Norbury, John W.↗

NASA NTRS1994

A Galerkin method for linear PDE systems in circular geometries with structural acoustic applications

A Galerkin method for systems of PDE's in circular geometries is presented with motivating problems being drawn from structural, acoustic, and structural acoustic applications. Depending upon the application under consideration, piecewise splines or Legendre polynomials are used when approximating the system dynamics with modifications included to incorporate the analytic solution decay near the coordinate singularity. This provides an efficient method which retains its accuracy throughout the circular domain without degradation at singularity. Because the problems under consideration are linear or weakly nonlinear with constant or piecewise constant coefficients, transform methods for the problems are not investigated. While the specific method is developed for the two dimensional wave equations on a circular domain and the equation of transverse motion for a thin circular plate, examples demonstrating the extension of the techniques to a fully coupled structural acoustic system are used to illustrate the flexibility of the method when approximating the dynamics of more complex systems.

Smith, Ralph C.↗

NASA NTRS1994

Influence of cross-flow on nonlinear Tollmien-Schlichting/vortex interaction

The transition of an incompressible three-dimensional boundary layer with strong cross-flow is considered theoretically and computationally in the context of vortex/wave interactions. Specifically the work centers on two lower-branch Tollmien-Schlichting waves which mutually interact nonlinearly to induce a longitudinal vortex flow. The vortex motion in turn gives rise to significant wave modulation via wall-shear forcing. The characteristic Reynolds number is large and, as a consequence, the waves' and the vortex motion are governed primarily by triple deck theory. The nonlinear interaction is captured by a viscous partial-differential system for the vortex coupled with a pair of amplitude equations for each wave pressure. Following analysis and computation over a wide range of parameters, three distinct responses are found to emerge in the nonlinear behavior of the flow solution downstream: an algebraic finite-distance singularity, far-downstream saturation or far-downstream wave decay leaving pure vortex flow. These depend on the input conditions, the wave angles and the size of the cross flow.

Davis, D. A. R.↗

NASA NTRS1998

The Real-Time Wall Interference Correction System of the NASA Ames 12-Foot Pressure Wind Tunnel

An improved version of the Wall Signature Method was developed to compute wall interference effects in three-dimensional subsonic wind tunnel testing of aircraft models in real-time. The method may be applied to a full-span or a semispan model. A simplified singularity representation of the aircraft model is used. Fuselage, support system, propulsion simulator, and separation wake volume blockage effects are represented by point sources and sinks. Lifting effects are represented by semi-infinite line doublets. The singularity representation of the test article is combined with the measurement of wind tunnel test reference conditions, wall pressure, lift force, thrust force, pitching moment, rolling moment, and pre-computed solutions of the subsonic potential equation to determine first order wall interference corrections. Second order wall interference corrections for pitching and rolling moment coefficient are also determined. A new procedure is presented that estimates a rolling moment coefficient correction for wings with non-symmetric lift distribution. Experimental data obtained during the calibration of the Ames Bipod model support system and during tests of two semispan models mounted on an image plane in the NASA Ames 12 ft. Pressure Wind Tunnel are used to demonstrate the application of the wall interference correction method.

Ulbrich, Norbert↗

NASA NTRS2001

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗

NASA NTRSDate not supplied

A New Method for the Accurate Calculation of Derivatives at Structured Multiblock Grid Singularities

The new Structured Mapping At Reduced dimension Topology elements (SMART) method for the accurate calculation of derivatives at multiblock structured grid block interfaces is presented. In the SMART method, a ’local’ curvilinear mapping is defined at each interface point, giving a consistent multidimensional definition of the derivatives. The resulting system of equations is overdetermined, and a least squares method is used to determine the value of the derivatives in the local coordinate system. In the special case of a nonsingular interface point, the SMART method automatically recovers the spatial differencing stencils used at the interior points in the neighboring grid blocks.

Computational Aeroacoustics↗

NASA NTRS1972

Lifting-surface theory for calculating the loading induced on a wing by a flap

A method is described for using lifting-surface theory to obtain the pressure distribution on a wing with a trailing-edge flap or control surface. The loading has a logarithmic singularity at the flap edges, which may be determined directly by the method of matched asymptotic expansions. Expressions are given for the singular flap loading for various flap hinge line and side edge geometries, both for steady and unsteady flap deflection. The regular part of the flap loading must be obtained by inverting the lifting-surface-theory integral equation relating the pressure and the downwash on the wing: procedures are described to accomplish this for a general wing and flap geometry. The method is applied to several example wings, and the results are compared with experimental data. Theory and test correlate well.

Johnson, W. A.↗

NASA NTRS1973

Optimal low thrust geocentric transfer

A computer code which will rapidly calculate time-optimal low thrust transfers is being developed as a mission analysis tool. The final program will apply to NEP or SEP missions and will include a variety of environmental effects. The current program assumes constant acceleration. The oblateness effect and shadowing may be included. Detailed state and costate equations are given for the thrust effect, oblateness effect, and shadowing. A simple but adequate model yields analytical formulas for power degradation due to the Van Allen radiation belts for SEP missions. The program avoids the classical singularities by the use of equinoctial orbital elements. Kryloff-Bogoliuboff averaging is used to facilitate rapid calculation. Results for selected cases using the current program are given.

Edelbaum, T. N.↗