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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The piecewise-linear predictor-corrector code - A Lagrangian-remap method for astrophysical flows

We describe a time-explicit finite-difference algorithm for solving the nonlinear fluid equations. The method is similar to existing Eulerian schemes in its use of operator-splitting and artificial viscosity, except that we solve the Lagrangian equations of motion with a predictor-corrector and then remap onto a fixed Eulerian grid. The remap is formulated to eliminate errors associated with coordinate singularities, with a general prescription for remaps of arbitrary order. We perform a comprehensive series of tests on standard problems. Self-convergence tests show that the code has a second-order rate of convergence in smooth, two-dimensional flow, with pressure forces, gravity, and curvilinear geometry included. While not as accurate on idealized problems as high-order Riemann-solving schemes, the predictor-corrector Lagrangian-remap code has great flexibility for application to a variety of astrophysical problems.

Lufkin, Eric A.↗
NASA NTRS1995

The history of re-connection and the concept of the solar wind plasma with relatively small electrical conductivity

Petschek's 're-connection' model, aspiring to be universal, treated as a boundary problem meets unresolvable difficulties connected with impossibility to specify correctly boundary and initial conditions. This problem was incorrectly formulated. Hence, ineradicable logarithmic singularities occurred on the boundary surfaces. Attempts to eliminate them by incorporating the finite electrical conductivity are incorrect. This should lead to the change in the equation type, boundary condition type and in consequence to the change in solutions. Besides, the slow mode shocks cannot be driven by small internal source. As an alternative a new plasma concept is suggested. The state of fully ionized plasma in space depends completely on the entropy of the plasma heating source and on the process in which plasma is involved. The presumptive source of the solar wind creation - the induction electric field of the solar origin - has very low entropy. The state of plasma should be very far from the thermodynamic equilibrium. Debye's screening is not complete. The excitation of the powerful resonant self-consistent electric fields in plasma provides low electric conductivity. The MHD problems should be treated in frameworks of dissipative theories.

Chertkov, A. D.↗
NASA NTRS1955

On the Kernel function of the integral equation relating the lift and downwash distributions of oscillating finite wings in subsonic flow

This report treats the Kernel function of an integral equation that relates a known prescribed downwash distribution to an unknown lift distribution for a harmonically oscillating finite wing in compressible subsonic flow. The Kernel function is reduced to a form that can be accurately evaluated by separating the Kernel function into two parts: a part in which the singularities are isolated and analytically expressed and a nonsingular part which may be tabulated. The form of the Kernel function for the sonic case (Mach number 1) is treated separately. In addition, results for the special cases of Mach number of 0 (incompressible case) and frequency of 0 (steady case) are given. The derivation of the integral equation which involves this Kernel function is reproduced as an appendix. Another appendix gives the reduction of the form of the Kernel function obtained herein for the three-dimensional case to a known result of Possio for two-dimensional flow. A third appendix contains some remarks on the evaluation of the Kernel function, and a fourth appendix presents an alternate form of expression for the Kernel function.

Watkins, Charles E↗
NASA NTRS1988

Dynamic structural correlation via nonlinear programming techniques

A solution to the correlation between structural dynamic test results and finite element analyses of the same components is presented in this paper. Basically, the method can be categorized as a Levenberg-Marquardt type Gauss-Newton method which requires only the differences between FE modal analyses and test results and their first derivatives with respect to preassigned design variables. With proper variable normalization and equation scaling, the method has been made numerically better-conditioned and the inclusion of the Levenberg-Marquardt technique overcomes any remaining difficulty encountered in inverting singular or near-singular matrices. An important feature is that each iteration requires only one function evaluation along with the associated design sensitivity analysis and so the procedure is computationally efficient.

Ting, T.↗
NASA NTRS1994

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

Vehicle control for fuel optimization

The problem of predicting the minimum fuel trajectory for a six degree of freedom vehicle which has a motion characterized by the first order differential equations of translational and rotational dynamics is considered. The thrust direction and center of gravity of the vehicle are assumed to be fixed with respect to the vehicle. Thrust magnitude and the control moment are used as control variables and appear linearly in the equation of motion. Pontryagin's Maximum Principle is used to solve the variational problem. With this formulation, the extremal controls are bang-bang with the exception of the singular case.

Dannenberg, K. D.↗
NASA NTRS1977

Unsteady supersonic aerodynamic theory by the method of potential gradient

A generalized solution of the hyperbolic wave equation has been derived. The method used has been further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals have been evaluated in closed form, while numerical integration methods are suggested for treating more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient because line integrations along the other two sides need only be considered. Furthermore, all the integrands vanish on the Mach cone and the need for determining the hyperbolic curves of intersection of the cone with the lifting surface is avoided. Generalized aerodynamic coefficients for three AGARD planforms have been calculated and compared with the available results.

Jones, W. P.↗
NASA NTRS1983

Solution of Poisson equations for 3-dimensional grid generations

A method for generating three dimensional, finite difference grids about complicated geometries by using Poisson equations is developed. The inhomogenous terms are automatically chosen such that orthogonality and spacing restrictions at the body surface are satisfied. Spherical variables are used to avoid the axis singularity, and an alternating-direction-implicit (ADI) solution scheme is used to accelerate the computations. Computed results are presented that show the capability of the method. Since most of the results presented have been used as grids for flow-field computations, this is indicative that the method is a useful tool for generating three-dimensional grids about complicated geometries.

Fujii, K.↗
NASA NTRS1992

Use of the maximum entropy method to retrieve the vertical atmospheric ozone profile and predict atmospheric ozone content

A method is developed for prediction of ozone levels in planetary atmospheres. This method is formulated in terms of error covariance matrices, and is associated with both direct measurements, a priori first guess profiles, and a weighting function matrix. This is described by the following linearized equation: y = A(matrix) x X + eta, where A is the weighting matrix and eta is noise. The problems to this approach are: (1) the A matrix is near singularity; (2) the number of unknowns in the profile exceeds the number of data points, therefore, the solution may not be unique; and (3) even if a unique solution exists, eta may cause the solution to be ill conditioned.

Turner, B. Curtis↗
NASA NTRS2001

Finite Element Methods for real-time Haptic Feedback of Soft-Tissue Models in Virtual Reality Simulators

We have applied the linear elastic finite element method to compute haptic force feedback and domain deformations of soft tissue models for use in virtual reality simulators. Our results show that, for virtual object models of high-resolution 3D data (>10,000 nodes), haptic real time computations (>500 Hz) are not currently possible using traditional methods. Current research efforts are focused in the following areas: 1) efficient implementation of fully adaptive multi-resolution methods and 2) multi-resolution methods with specialized basis functions to capture the singularity at the haptic interface (point loading). To achieve real time computations, we propose parallel processing of a Jacobi preconditioned conjugate gradient method applied to a reduced system of equations resulting from surface domain decomposition. This can effectively be achieved using reconfigurable computing systems such as field programmable gate arrays (FPGA), thereby providing a flexible solution that allows for new FPGA implementations as improved algorithms become available. The resulting soft tissue simulation system would meet NASA Virtual Glovebox requirements and, at the same time, provide a generalized simulation engine for any immersive environment application, such as biomedical/surgical procedures or interactive scientific applications.

Frank, Andreas O.↗
NASA NTRS2004

Hamilton's Equations with Euler Parameters for Rigid Body Dynamics Modeling

A combination of Euler parameter kinematics and Hamiltonian mechanics provides a rigid body dynamics model well suited for use in strongly nonlinear problems involving arbitrarily large rotations. The model is unconstrained, free of singularities, includes a general potential energy function and a minimum set of momentum variables, and takes an explicit state space form convenient for numerical implementation. The general formulation may be specialized to address particular applications, as illustrated in several three dimensional example problems.

Shivarama, Ravishankar↗
NASA NTRS1972

Stresses in bonded materials with a crack perpendicular to the interface.

The problem of two elastic bonded half planes containing a crack perpendicular to the interface is considered. First the solution for the semi-infinite crack under concentrated wedge loading is given. Then the problem of a finite crack fully imbedded in one of the half planes of terminating at the interface is considered. The Mellin transform in conjunction with the dislocations is used to formulate the problem and to derive the integral equation. The integral equation is solved; the stress intensity factors, the crack surface displacement, and the stresses around the crack tip terminating at the interface are obtained. Noting that the power of the singularity for the interface tip of the crack is not -1/2, a tentative fracture criterion dealing with the initiation of fracture propagation is proposed.

Cook, T. S.↗
NASA NTRS1986

Theoretical three-and four-axis gimbal robot wrists

In high-performance flight simulations, a four-axis gimbal system allows all possible rotations with acceptable gimbal angle rates while it avoids the so-callled 'gimbal lock' that occurs when gimbal rotational axes are colinear. In this paper, pertinent equations (including quaternions) are assembled for a hypothetical robot wrist, functionally equivalent to this four-axis gimbal system, and also for a true three-axis gimbal robot wrist. These equations are used to simulate the rotation of a robot hand by the robot wrist in response to operator rotational velocity commands to the robot hand. Near gimbal lock (wrist singularity), excessive rotational rates occur. Scaling the rates, which is necessary for the three-gimbal robot wrist to prevent rate limiting, introduces an undesirable time delay in the robot hand rotation with respect to the commanded rotation. However, the merit of the four-gimbal robot wrist is that the fourth gimbal angle keeps the robot wrist away from the singularity so that the robot hand moves exactly as commanded. It appears that in a 'worst-type' maneuver of the robot hand, the fourth gimbal angle can be defined so that none of the gimbal angle rates exceed about twice the commanded rates.

Barker, L. K.↗
NASA NTRS1989

Blow-up of unsteady two-dimensional Euler and Navier-Stokes solutions having stagnation-point form

The time-dependent form of the classic, two-dimensional stagnation-point solution of the Navier-Stokes equations is considered. If the viscosity is zero, a class of solutions of the initial-value problem can be found in closed form using Lagrangian coordinates. These solutions exhibit singular behavior in finite time, because of the infinite domain and unbounded initial vorticity. Thus, the blow-up found by Stuart in three dimensions using the stagnation-point form, also occurs in two. The singularity vanishes under a discrete, finite-dimensional 'point vortex' approximation, but is recovered as the number of vortices tends to infinity. We find that a small positive viscosity does not arrest the breakdown, but does strongly alter its form. Similar results are summarized for certain Boussinesq stratified flows.

Childress, S.↗
NASA NTRS1998

Nonlinear Interaction of Detuned Instability Waves in Boundary-Layer Transition: Amplitude Equations

The non-equilibrium critical-layer analysis of a system of frequency-detuned resonant-triads is presented. In this part of the analysis, the system of partial differential critical-layer equations derived in Part I is solved analytically to yield the amplitude equations which are analyzed using a combination of asymptotic and numerical methods. Numerical solutions of the inviscid non-equilibrium oblique-mode amplitude equations show that the frequency-detuned self-interaction enhances the growth of the lower-frequency oblique modes more than the higher-frequency ones. All amplitudes become singular at the same finite downstream position. The frequency detuning delays the occurrence of the singularity. The spanwise-periodic mean-flow distortion and low-frequency nonlinear modes are generated by the critical-layer interaction between frequency-detuned oblique modes. The nonlinear mean flow and higher harmonics as well as the primary instabilities become as large as the base mean flow in the inviscid wall layer in the downstream region where the distance from the singularity is of the order of the wavelength scale.

Lee, Sang Soo↗
NASA NTRS1979

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

Out-of-Core Solutions of Complex Sparse Linear Equations

ETCLIB is library of subroutines for obtaining out-of-core solutions of complex sparse linear equations. Routines apply to dense and sparse matrices too large to be stored in core. Useful for solving any set of linear equations, but particularly useful in cases where coefficient matrix has no special properties that guarantee convergence with any of interative processes. The only assumption made is that coefficient matrix is not singular.

Yip, E. L.↗