Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATHEMATICAL SOLUTIONS”

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 289 records · Page 16

Constraints on the physical properties of optical bullets in SS 433

The present study discusses possible mechanisms for continuously heating the H-alpha emitting bullets of SS 433 out to distances of 5 x 10 to the 14th cm and for turning off this emission at 10 to the 15th. Various observational contraints are used to establish bounds on permissible solutions in terms of the two key bullet parameters, mass, and angular radius seen from the central source. The analysis is carried out for the mathematically simplest case of uniform spherical bullets. For radiative heat of such bullets by starlight, solutions are found to exist only for the very massive bullets with about 0.03 radians, which are highly implausible on the grounds of the large implied mean kinetic luminosity of about 10 to the 41st ergs/s. It is concluded that collisional interaction is the most likely mechanism for heating the optical bullets of SS 433. The effects of these constraints being fragmented rather than uniform, and being elongated rather than spherical are discussed.

Brown, John C.↗

Segregation effects during solidification in weightless melts

The generalized problem of determining the temperature and solute concentration profiles during directional solidification of binary alloys with surface evaporation was mathematically formulated. Realistic initial and boundary conditions were defined, and a computer program was developed and checked out. The programs computes the positions of two moving boundaries, evaporation and solidification, and their velocities. Temperature and solute concentration profiles in the semiinfinite material body at selected instances of time are also computed.

Li, C.↗

Application of the Parker flash method to high temperature thermoelectrics

Mathematical aspects of thermal diffusivity measurements are reviewed and applied to develop an experimental setup at JPL. Certain approximations commonly appearing in the literature are examined in some detail and rejected based upon sensitive slope variations having significant impact on the results. The final mathematical formulation used was one developed by Parker et al. (1961), which evaluated the diffusivity in a way that could be simply evaluated by examining an oscilloscope trace. A computer algorithm was written to examine solution sensitivity to the number of terms taken in the infinite series. The mathematical technique is being used at JPL to evaluate diffusivities of graphite, Boron Carbon and Lanthanum Sulphur compounds in the temperature range of 300 K to 1200 K.

Noon, E. L.↗

Program Computes Flows Of Fluids And Heat

SINDA'85/FLUINT incorporates lumped-parameter-network and one-dimensional-flow mathematical models. System enables analysis of mutual influences of thermal and flow phenomena. Offers two finite-difference numerical solution techniques: forward-difference explicit approximation and Crank-Nicholson approximation. Enables simulation of nonuniform heating and facilitates mathematical modeling of thin-walled heat exchangers. Ability to model nonequilibrium behavior within two-phase volumes included. Recent changes in program improve modeling of real evaporator pumps and other capillary-assist evaporators. Written in FORTRAN 77.

Cullimore, Brent↗

The problem of exact interior solutions for rotating rigid bodies in general relativity

The (3 + 1) dyadic formalism for timelike congruences is applied to derive interior solutions for stationary, axisymmetric, rigidly rotating bodies. In this approach the mathematics is formulated in terms of three-space-covariant, first-order, vector-dyadic, differential equations for a and Omega, the acceleration and angular velocity three-vectors of the rigid body; for T, the stress dyadic of the matter; and for A and B, the 'electric' and 'magnetic' Weyl curvature dyadics which describe the gravitational field. It is shown how an appropriate ansatz for the forms of these dyadics can be used to discover exact rotating interior solutions such as the perfect fluid solution first published in 1968. By incorporating anisotropic stresses, a generalization is found of that previous solution and, in addition, a very simple new solution that can only exist in toroidal configurations.

Wahlquist, H. D.↗

A finite element formulation of Euler equations for the solution of steady transonic flows

The main objective of the considered investigation is related to the development of a relaxation scheme for the analysis of inviscid, rotational, transonic flow problems. To formulate the equations of motion for inviscid flows in a fixed coordinate system, an Eulerian type variational principle is required. The derivation of an Eulerian variational principle which is employed in the finite element formulation is discussed. The presented numerical method describes the mathematical formulation and the application of a numerical process for the direct solution of steady Euler equations. The development of the procedure as an extension of existing potential flow formulations provides the applicability of previous procedures, e.g., proper application of the artificial viscosity for supersonic elements, and the accurate modeling of the shock.

Ecer, A.↗

Modeling and estimation of a low degree geopotential model from terrestrial gravity data

The development of appropriate modeling and adjustment procedures for the estimation of harmonic coefficients of the geopotential, from surface gravity data was studied, in order to provide an optimum way of utilizing the terrestrial gravity information in combination solutions currently developed at NASA/Goddard Space Flight Center, for use in the TOPEX/POSEIDON mission. The mathematical modeling was based on the fundamental boundary condition of the linearized Molodensky boundary value problem. Atmospheric and ellipsoidal corrections were applied to the surface anomalies. Terrestrial gravity solutions were found to be in good agreement with the satellite ones over areas which are well surveyed (gravimetrically), such as North America or Australia. However, systematic differences between the terrestrial only models and GEMT1, over extended regions in Africa, the Soviet Union, and China were found. In Africa, gravity anomaly differences on the order of 20 mgals and undulation differences on the order of 15 meters, over regions extending 2000 km in diameter, occur. Comparisons of the GEMT1 implied undulations with 32 well distributed Doppler derived undulations gave an RMS difference of 2.6 m, while corresponding comparison with undulations implied by the terrestrial solution gave RMS difference on the order of 15 m, which implies that the terrestrial data in that region are substantially in error.

Pavlis, Nikolaos K.↗

Finite element (MARC) solution technologies for viscoplastic analyses

A need for development of realistic constitutive models for structural components operating at high temperatures, accompanied by appropriate solution technologies for stress/life analyses of these components is studied. Viscoplastic models provide a better description of inelastic behavior of materials, but their mathematical structure is very complex. The highly nonlinear and stiff nature of the constitutive equations makes analytical solutions difficult. Therefore, suitable solution, finite element or other numerical, technologies must be developed to make these models adaptable for better and rational designs of components. NASA-Lewis has developed several solution technologies and successfully applied them to the solution of a number of uniaxial and multiaxial problems. Some of these solution technologies are described along with the models and representative results. The solution technologies developed and presented encompass a wide range of models, such as, isotropic, anisotropic, metal matrix composites, and single crystal models.

Arya, V. K.↗

Automating FEA programming

In this paper we briefly describe a combined symbolic and numeric approach for solving mathematical models on parallel computers. An experimental software system, PIER, is being developed in Common Lisp to synthesize computationally intensive and domain formulation dependent phases of finite element analysis (FEA) solution methods. Quantities for domain formulation like shape functions, element stiffness matrices, etc., are automatically derived using symbolic mathematical computations. The problem specific information and derived formulae are then used to generate (parallel) numerical code for FEA solution steps. A constructive approach to specify a numerical program design is taken. The code generator compiles application oriented input specifications into (parallel) FORTRAN77 routines with the help of built-in knowledge of the particular problem, numerical solution methods and the target computer.

Sharma, Naveen↗

A theory of supercritical wing sections, with computer programs and examples.

Mathematical methods for the design of supercritical wings, which depend on the numerical solution of the partial differential equations of two-dimensional gas dynamics, are developed. The main contribution is a computer program for the design of shockless transonic airfoils using the hodograph transformation and analytic continuation into the complex domain. The mathematical theory is described, and a manual for users of the programs is provided. Numerical examples are given and computational results are discussed, and the computer programs themselves are listed. The analysis routine can be used to ascertain whether the profiles behave well at off-design conditions, or to smooth coordinates and obtain a desirable shape more quickly when perfectly shockless flow is not essential.

Bauer, F.↗

The effect of elastic boundary conditions on the dynamic response of rectangular plates

Natural frequencies and forced steady-state harmonic response for the vibration of uniform rectangular plates with edges elastically restrained against rotation and transverse translation are addressed. A single mode Rayleigh-Ritz solution is derived using functions that describe the normal modes of vibration of a beam whose ends are elastically restrained. The finite element solution is obtained for comparison. MACSYMA symbolic manipulation system is implemented as an aid to the mathematical rigor of the Ritz approach, and NASTRAN finite element code is used to model the mechanical system. Comparisons are made to published results and the solutions of this study are found to give lower frequencies for some values of boundary restraint. Steady-state harmonic amplitudes of displacement and acceleration are found to agree favorably for the two solutions. Low predictions of steady-state strain from NASTRAN result in some cases when compared to the Ritz values. Finally, a subjective assessment is made about the merit of using MACSYMA and NASTRAN.

Brewer, Terry K.↗

Finite elements of nonlinear continua.

The finite element method is extended to a broad class of practical nonlinear problems, treating both theory and applications from a general and unifying point of view. The thermomechanical principles of continuous media and the properties of the finite element method are outlined, and are brought together to produce discrete physical models of nonlinear continua. The mathematical properties of the models are analyzed, and the numerical solution of the equations governing the discrete models is examined. The application of the models to nonlinear problems in finite elasticity, viscoelasticity, heat conduction, and thermoviscoelasticity is discussed. Other specific topics include the topological properties of finite element models, applications to linear and nonlinear boundary value problems, convergence, continuum thermodynamics, finite elasticity, solutions to nonlinear partial differential equations, and discrete models of the nonlinear thermomechanical behavior of dissipative media.

Oden, J. T.↗

Spherical means of solutions of partial differential equations in a conical region

The spherical means of the solutions of a linear partial differential equation Lu = f in a conical region are studied. The conical region is bounded by a surface generated by curvilinear ti surfaces. The spherical mean is the average of u over a constant ti surface. The conditions on the linear differential operator, L, and on the orthogonal coordinates (ti, eta, zeta) are established so that the spherical mean of the solution subjected to the appropriate boundary and initial conditions can be determined directly as a problem with only space variable. Conditions are then established so that the spherical mean of the solution in one concial region will be proportional to that of a known solution in another conical region. Applications to various problems of mathematical physics and their physical interpretations are presented.

Ting, L.↗

Spherical means of solutions of partial differential equations in a conical region

The spherical means of the solutions of a linear partial differential equation Lu = f in a conical region are studied. The conical region is bounded by a surface generated by curvilinear xi lines and by two truncating xi surfaces. The spherical mean is the average of u over a constant xi surface. Conditions on the linear differential operator, L, and on the orthogonal coordinates xi, eta, and zeta are established so that the problem for the determination of the spherical mean of the solution subjected to the appropriate boundary and initial conditions can be reduced to a problem with only one space variable. Conditions are then established so that the spherical mean of the solution in one conical region will be proportional to that of a known solution in another conical region. Applications to various problems of mathematical physics and their physical interpretations are presented.

Ting, L.↗