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 271 records · Page 15

Vortex methods for separated flows

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented in an elementary fashion and includes the relationship with traditional point-vortex studies, the convergence to smooth solutions of the Euler equations, and the essential differences between two- and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. The overlap with the excellent review articles available is kept to a minimum and more emphasis is placed on the area of expertise, namely two-dimensional flows around bluff bodies. When solid walls are present, complete mathematical models are not available and a more heuristic attitude must be adopted. The imposition of inviscid and viscous boundary conditions without conformal mappings or image vortices and the creation of vorticity along solid walls are examined in detail. Methods for boundary-layer treatment and the question of the Kutta condition are discussed. Practical aspects and tips helpful in creating a method that really works are explained. The topics include the robustness of the method and the assessment of accuracy, vortex-core profiles, timemarching schemes, numerical dissipation, and efficient programming. Calculations of flows past streamlined or bluff bodies are used as examples when appropriate.

Spalart, Philippe R.↗

The structure of synchronously rotating close binaries built on polytropic model nu = 3

Mathematical procedures are developed which can provide an exact solution for the structure of synchronously rotating close binaries built on the polytropic model with nu equal to 3. Results are discussed for separations of centers of 10.0, 4.0, 3.0, and 2.8 solar radii. For the present discussion, both components are assumed to have the solar mass and, in the zeroth approximation, the solar radius. The present results are compared with those obtained by Chandrasekhar (1933) using first-order perturbation theory.

Green, L. C.↗

On estimating gravity anomalies: A comparison of least squares collocation with least squares techniques

The least squares collocation algorithm for estimating gravity anomalies from geodetic data is shown to be an application of the well known regression equations which provide the mean and covariance of a random vector (gravity anomalies) given a realization of a correlated random vector (geodetic data). It is also shown that the collocation solution for gravity anomalies is equivalent to the conventional least-squares-Stokes' function solution when the conventional solution utilizes properly weighted zero a priori estimates. The mathematical and physical assumptions underlying the least squares collocation estimator are described, and its numerical properties are compared with the numerical properties of the conventional least squares estimator.

Argentiero, P.↗

On estimating gravity anomalies - A comparison of least squares collocation with conventional least squares techniques

The least squares collocation algorithm for estimating gravity anomalies from geodetic data is shown to be an application of the well known regression equations which provide the mean and covariance of a random vector (gravity anomalies) given a realization of a correlated random vector (geodetic data). It is also shown that the collocation solution for gravity anomalies is equivalent to the conventional least-squares-Stokes' function solution when the conventional solution utilizes properly weighted zero a priori estimates. The mathematical and physical assumptions underlying the least squares collocation estimator are described.

Argentiero, P.↗

Impedance characteristics of coaxial and planar magnetoplasma capacitors

A theory has been developed for the impedance of a homogeneous magnetoplasma enclosed between two specular reflecting coaxial electrodes, with a static magnetic field parallel to the electrode axes. The parallel-plate magnetoplasma capacitor is treated as a sub-case. Starting with the Vlasov equation, an integral equation is derived for the electric field. Solving this equation, and integrating to obtain the voltage, gives the capacitor impedance. This includes a capacitive component, and a resistive component expressing the Landau damping associated with the open orbits of electrons reflected at the electrodes. A direct numerical solution of the field integral equation has been carried out for a range of values of magnetic field, plasma density, and signal frequency. The values of impedance so obtained are compared with the predictions of macroscopic theory, and of an approximate microscopic theory in which open orbits are ignored and solutions are obtained using finite Fourier transform methods. The mathematical relations between these theories are demonstrated.

Harker, K. J.↗

End effects in scarf joints

Scarf joints with small scarf angles are especially sensitive to stiffness mismatch between adherends and to adherend tip bluntness. Pre-assembly breakage of an adherend tip where it is only a few microns thick can cause significant reduction in joint strength. Mathematically, the reason for such sensitivity is that the solutions to the governing differential equation develop boundary layer character when the scarf angles are small. The boundary layers are regions with large adhesive stresses. Experimental strength data for laminated composite adherends agree with the results of this analysis.

Adkins, D. W.↗

Inertial oscillation of a vertical rotating draft with application to a supercell storm

An analytic model (vertical rotating draft) which includes the gross features of a supercell storm on an f-plane, undergoes an inertial oscillation that appears to have been overlooked in previous analytic and numerical models. The oscillation is nonlinear and consists of a long quiescent phase and a short intense phase. During the intense phase, the rotating draft has the following features of a supercell: the diameter of the core contracts as it spins up and expands as it spins down; if vertical wind shear is included, the track of the rotating draft turns to the right (an anticyclonic rotating draft turns to the left); this turning point is followed by a predominantly upward flow; and the horizontal pressure gradient is very small (a property of most tornadoless supercells). The rapid spin-up during the intense phase and the high Rossby numbers obtainable establish the ability of the Coriolis force to spin up single cyclonic or anticyclonic supercells by means of this inertial oscillation. This surprising result has implications for numerical supercell simulations, which generally do not rely on the Coriolis force as a source of rotation. The physics and mathematics of the inertial oscillation are given, and the solution is applied to a documented supercell.

Supercell thunderstorm↗

Failure of the method of slowly varying amplitude and phase for non-linear, singular oscillators

It is shown that the method of slowly varying amplitude and phase yields erroneous results in the study of the mathematical properties of nonlinear singular oscillator systems. The analytical solution is described in which the phase function is constant and for which a special limiting behavior exists when the wavelength is zero. The previous method based on the condition of boundedness cannot be satisfied for nonlinear singular characteristics, and the erroneous designation of the expansion parameter is identified.

Mickens, R. E.↗

Magnetohydrodynamic wave propagation in one-dimensional nonhomogeneous, self-gravitating clouds

We study the propagation of magnetohydrodynamic (MHD) waves through nonhomogeneous, selfgravitating, magnetic media representative of molecular cloud environments and focus on the issues of cloud support and line profiles. Since the general treatment of this topic is burdened by a complex mathematical formalism, we consider simplifying geometries which yield analytical solutions in the linear wave limit. In particular, we study both magnetoacoustic and Alfven wave propagation along the density gradient in a one-imensional slab. Of specific relevance to molecular clouds, we find that the back reaction of the Alfven waves can provide a pressure along the direction of the magnetic field lines; this pressure can help support a density enhancement against gravitational collapse. Furthermore, we find that the velocity amplitudes of these waves increase as the density decreases, in rough agreement with observational estimates of line width versus density relations.

Fatuzzo, Marco↗

Development of a nondestructive vibration technique for bond assessment of Space Shuttle tiles

This final report describes the achievements of the above titled project. The project is funded by NASA-KSC (Grant No. NAG 10-0117) for the period of 1 Jan. to 31 Dec. 1993. The purpose of this project was to develop a nondestructive, noncontact technique based on 'vibration signature' of tile systems to quantify the bond conditions of the thermal protection system) tiles of Space Shuttle orbiters. The technique uses a laser rapid scan system, modal measurements, and finite element modeling. Finite element models were developed for tiles bonded to both clamped and deformable integrated skin-stringer orbiter mid-fuselage. Results showed that the size and location of a disbonded tile can be determined from frequency and mode shape information. Moreover, a frequency response survey was used to quickly identify the disbonded tiles. The finite element results were compared with experimentally determined frequency responses of a 17-tile test panel, where a rapidscan laser system was employed. An excellent degree of correlation between the mathematical simulation and experimental results was realized. An inverse solution for single-tile assemblies was also derived and is being implemented into a computer program that can interact with the modal testing software. The output of the program displays the size and location of disbond. This program has been tested with simulated input (i.e., finite element data), and excellent agreement between predicted and simulated disbonds was shown. Finally, laser vibration imaging and acoustic emission techniques were shown to be well suited for detecting and monitoring the progressive damage in Graphite/Epoxy composite materials.

Moslehy, Faissal A.↗

Improved Tracking of Targets by Cameras on a Mars Rover

A paper describes a method devised to increase the robustness and accuracy of tracking of targets by means of three stereoscopic pairs of video cameras on a Mars-rover-type exploratory robotic vehicle. Two of the camera pairs are mounted on a mast that can be adjusted in pan and tilt; the third camera pair is mounted on the main vehicle body. Elements of the method include a mast calibration, a camera-pointing algorithm, and a purely geometric technique for handing off tracking between different camera pairs at critical distances as the rover approaches a target of interest. The mast calibration is an extension of camera calibration in which the camera images of calibration targets at known positions are collected at various pan and tilt angles. In the camerapointing algorithm, pan and tilt angles are computed by a closed-form, non-iterative solution of inverse kinematics of the mast combined with mathematical models of the cameras. The purely geometric camera-handoff technique involves the use of stereoscopic views of a target of interest in conjunction with the mast calibration.

Kim, Won↗

Drag and Lift Forces Between a Rotating Conductive Sphere and a Cylindrical Magnet

Modeling the interaction between a non-uniform magnetic field and a rotating conductive object allows study of the drag force which is used in applications such as eddy current braking and linear induction motors as well as the transition to a repulsive force that is the basis for magnetic levitation systems. Here, we study the interaction between a non-uniform field generated by a cylindrical magnet and a rotating conductive sphere. Each eddy current in the sphere generates a magnetic field which in turn generates another eddy current, eventually feeding back on itself. A two step mathematics process is developed to find a closed form solution in terms of only two eddy currents. However, the complete solution requires decomposition of the magnetic field into a summation of spherical harmonics, making it more suitable for a graduate level electromagnetism lecture or lab. Finally, the forces associated with these currents are calculated and then verified experimentally.

Nurge, Mark A.↗

Radiative transfer calculated from a Markov chain formalism

The theory of Markov chains is used to formulate the radiative transport problem in a general way by modeling the successive interactions of a photon as a stochastic process. Under the minimal requirement that the stochastic process is a Markov chain, the determination of the diffuse reflection or transmission from a scattering atmosphere is equivalent to the solution of a system of linear equations. This treatment is mathematically equivalent to, and thus has many of the advantages of, Monte Carlo methods, but can be considerably more rapid than Monte Carlo algorithms for numerical calculations in particular applications. We have verified the speed and accuracy of this formalism for the standard problem of finding the intensity of scattered light from a homogeneous plane-parallel atmosphere with an arbitrary phase function for scattering. Accurate results over a wide range of parameters were obtained with computation times comparable to those of a standard 'doubling' routine. The generality of this formalism thus allows fast, direct solutions to problems that were previously soluble only by Monte Carlo methods. Some comparisons are made with respect to integral equation methods.

Esposito, L. W.↗

On designing for quality

The problem of ensuring the required quality of products and/or technological processes often becomes more difficult due to the fact that there is not general theory of determining the optimal sets of value of the primary factors, i.e., of the output parameters of the parts and units comprising an object and ensuring the correspondence of the object's parameters to the quality requirements. This is the main reason for the amount of time taken to finish complex vital article. To create this theory, one has to overcome a number of difficulties and to solve the following tasks: the creation of reliable and stable mathematical models showing the influence of the primary factors on the output parameters; finding a new technique of assigning tolerances for primary factors with regard to economical, technological, and other criteria, the technique being based on the solution of the main problem; well reasoned assignment of nominal values for primary factors which serve as the basis for creating tolerances. Each of the above listed tasks is of independent importance. An attempt is made to give solutions for this problem. The above problem dealing with quality ensuring an mathematically formalized aspect is called the multiple inverse problem.

Vajingortin, L. D.↗

Improved Timing Scheme for Spaceborne Precipitation Radar

An improved timing scheme has been conceived for operation of a scanning satellite-borne rain-measuring radar system. The scheme allows a real-time-generated solution, which is required for auto targeting. The current timing scheme used in radar satellites involves pre-computing a solution that allows the instrument to catch all transmitted pulses without transmitting and receiving at the same time. Satellite altitude requires many pulses in flight at any time, and the timing solution to prevent transmit and receive operations from colliding is usually found iteratively. The proposed satellite has a large number of scanning beams each with a different range to target and few pulses per beam. Furthermore, the satellite will be self-targeting, so the selection of which beams are used will change from sweep to sweep. The proposed timing solution guarantees no echo collisions, can be generated using simple FPGA-based hardware in real time, and can be mathematically shown to deliver the maximum number of pulses per second, given the timing constraints. The timing solution is computed every sweep, and consists of three phases: (1) a build-up phase, (2) a feedback phase, and (3) a build-down phase. Before the build-up phase can begin, the beams to be transmitted are sorted in numerical order. The numerical order of the beams is also the order from shortest range to longest range. Sorting the list guarantees no pulse collisions. The build-up phase begins by transmitting the first pulse from the first beam on the list. Transmission of this pulse starts a delay counter, which stores the beam number and the time delay to the beginning of the receive window for that beam. The timing generator waits just long enough to complete the transmit pulse plus one receive window, then sends out the second pulse. The second pulse starts a second delay counter, which stores its beam number and time delay. This process continues until an output from the first timer indicates there is less than one transmit pulse width until the start of the next receive event. This blocks future transmit pulses in the build-up phase. The feedback phase begins with the first timer paying off and starting the first receive window. When the first receive window is complete, the timing generator transmits the next beam from the list. When the second timer pays off, the second receive event is started. Following the second receive event, the timing generator will transmit the next beam on the list and start an additional timer. The timers work in a circular buffer fashion so there only need to be enough to cover the maximum number of echoes in flight.

Berkun, Andrew↗

Polynomial approximation of functions of matrices and its application to the solution of a general system of linear equations

During the process of solving a mathematical model numerically, there is often a need to operate on a vector v by an operator which can be expressed as f(A) while A is NxN matrix (ex: exp(A), sin(A), A sup -1). Except for very simple matrices, it is impractical to construct the matrix f(A) explicitly. Usually an approximation to it is used. In the present research, an algorithm is developed which uses a polynomial approximation to f(A). It is reduced to a problem of approximating f(z) by a polynomial in z while z belongs to the domain D in the complex plane which includes all the eigenvalues of A. This problem of approximation is approached by interpolating the function f(z) in a certain set of points which is known to have some maximal properties. The approximation thus achieved is almost best. Implementing the algorithm to some practical problem is described. Since a solution to a linear system Ax = b is x= A sup -1 b, an iterative solution to it can be regarded as a polynomial approximation to f(A) = A sup -1. Implementing the algorithm in this case is also described.

Tal-Ezer, Hillel↗