Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integral solution”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗
NASA NTRSDate not supplied

Molecular flow venting of a volume with an outgassing or desorbing source

A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution.

Michael Woronowicz↗
NASA NTRSDate not supplied

Molecular flow venting of a volume with an outgassing or desorbing source

"A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution."

Michael Woronowicz↗
NASA NTRS1980

A mixed-mode crack analysis of isotropic solids using conservation laws of elasticity

A simple and convenient method of analysis for studying two-dimensional mixed-mode crack problems is presented. The analysis is formulated on the basis of conservation laws of elasticity and of fundamental relationships in fracture mechanics. The problem is reduced to the determination of mixed-mode stress-intensity factor solutions in terms of conservation integrals involving known auxiliary solutions. One of the salient features of the present analysis is that the stress-intensity solutions can be determined directly by using information extracted in the far field. Several examples with solutions available in the literature are solved to examine the accuracy and other characteristics of the current approach. This method is demonstrated to be superior in its numerical simplicity and computational efficiency to other approaches. Solutions of more complicated and practical engineering fracture problems dealing with the crack emanating from a circular hole are presented also to illustrate the capacity of this method

Yau, J. F.↗
NASA NTRS1972

Multipath propagation study for L-band, over-ocean, satellite-aircraft communication link

The results are presented of a study of multipath propagation between an aeronautical relay satellite in synchronous earth orbit and an aircraft flying over the ocean. An analytical model for the reflected power is developed and the results of a computer solution are compared with an approximate solution obtained through integration by the method of steepest descent. The computer solution is shown to differ considerably from the approximate solution. Both methods are compared with the limited experimental results available.

Deal, J. H.↗
NASA NTRS1980

Generation of C-type cascade grids for viscous flow computation

A rapid procedure for generating C-type cascade grids suitable for viscous flow computations in turbomachinery blade rows is presented. The resulting mesh is periodic from one blade passage to the next, nearly orthogonal, and continuous across the wake downstream of a blade. The procedure employs a pair of conformal mappings that take the exterior of the cascade into the interior of an infinite strip with curved boundaries. The final transformation to a rectangular computational domain is accomplished numerically. The boundary values are obtained from a panel solution of an integral equation and the interior values by a rapid ADI solution of Laplace's equation. Examples of C-type grids are presented for both compressor and turbine blades and the extension of the procedure to three dimensions is briefly outlined.

Sockol, P. M.↗
NASA NTRS1971

FORMAC integration program: A special applications package used in developing techniques of orbital decay and long term ephemeris prediction for satellites in earth orbit

The symbolic manipulation capabilities of the FORMAC (Formula Manipulation Compiler) language are employed to expand and analytically evaluate integrals. The program integration is effected by expanding the integral(s) into a series of subintegrals and then substituting a pre-derived and pre-coded solution for that particular subintegral. Derivation of the integral solutions necessary for precoding is included, as is a discussion of the FORMAC system limitations encountered in the programming effort.

Rowe, C. K.↗
NASA NTRS1994

A fast numerical solution of scattering by a cylinder: Spectral method for the boundary integral equations

It is known that the exact analytic solutions of wave scattering by a circular cylinder, when they exist, are not in a closed form but in infinite series which converges slowly for high frequency waves. In this paper, we present a fast number solution for the scattering problem in which the boundary integral equations, reformulated from the Helmholtz equation, are solved using a Fourier spectral method. It is shown that the special geometry considered here allows the implementation of the spectral method to be simple and very efficient. The present method differs from previous approaches in that the singularities of the integral kernels are removed and dealt with accurately. The proposed method preserves the spectral accuracy and is shown to have an exponential rate of convergence. Aspects of efficient implementation using FFT are discussed. Moreover, the boundary integral equations of combined single and double-layer representation are used in the present paper. This ensures the uniqueness of the numerical solution for the scattering problem at all frequencies. Although a strongly singular kernel is encountered for the Neumann boundary conditions, we show that the hypersingularity can be handled easily in the spectral method. Numerical examples that demonstrate the validity of the method are also presented.

Hu, Fang Q.↗
NASA NTRS2011

Combined Weekly Coordinate Solutions from SLR and DORIS

In International Terrestrial Reference Frame (ITRF) 2005 and ITRF2008, the approach for the construction of solutions by the IERS has been for individual analysis centers of each technique to process geodetic tracking data, and for each technique to develop a solution (or contribution) that is integrated into the final ITRF solution by careful combination of the technique solutions. The connections between the geodetic networks are realized by the application of local ties. In an alternate approach, we may assure processing homogeneity by creating normal equations for different techniques with the same orbit determination software, using identically derived algorithms. Another derivative of this approach is to realize the ties between the techniques using satellites tracked with multiple techniques; in effect tieing the networks together using satellite dynamics. In this solution, we develop a time series and a set of cumulative solutions from Satellite Laser Ranging (SLR) & Doppler Orbitography and Radio-positioning Integrated by Satellite (DORIS) based on homogeneous processing with the NASA GEODYN precise orbit determination suite of programs, where we jointly combine weekly the SLR data to Lageos1, Lageos2, Starlette, and Stella with the DORIS data from SPOT2-SPOT5, as well as satellites that utilize both techniques (TOPEX/Poseidon, Envisat, Jason-2). We discuss the modeling that is applied including upgrades implemented since the submission of the GSC ITRF2008 contributions for IDS. Firstly, we compare the SLR-only solutions comprising four geodetic satellites with the standard approach of utilizing only Lageos1 & Lageos2. Secondly, we evaluate the impact on the DORIS coordinates of the joint analysis with the SLR data.

Lemoine, F.↗
NASA NTRS2008

High Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

Taylor series integration is implemented in a spacecraft trajectory analysis code-the Spacecraft N-body Analysis Program (SNAP) - and compared with the code s existing eighth-order Runge-Kutta Fehlberg time integration scheme. Nine trajectory problems, including near Earth, lunar, Mars and Europa missions, are analyzed. Head-to-head comparison at five different error tolerances shows that, on average, Taylor series is faster than Runge-Kutta Fehlberg by a factor of 15.8. Results further show that Taylor series has superior convergence properties. Taylor series integration proves that it can provide rapid, highly accurate solutions to spacecraft trajectory problems.

Scott, James R.↗
NASA NTRS2001

Coarse Grid Modeling of Turbine Film Cooling Flows Using Volumetric Source Terms

The recent trend in numerical modeling of turbine film cooling flows has been toward higher fidelity grids and more complex geometries. This trend has been enabled by the rapid increase in computing power available to researchers. However, the turbine design community requires fast turnaround time in its design computations, rendering these comprehensive simulations ineffective in the design cycle. The present study describes a methodology for implementing a volumetric source term distribution in a coarse grid calculation that can model the small-scale and three-dimensional effects present in turbine film cooling flows. This model could be implemented in turbine design codes or in multistage turbomachinery codes such as APNASA, where the computational grid size may be larger than the film hole size. Detailed computations of a single row of 35 deg round holes on a flat plate have been obtained for blowing ratios of 0.5, 0.8, and 1.0, and density ratios of 1.0 and 2.0 using a multiblock grid system to resolve the flows on both sides of the plate as well as inside the hole itself. These detailed flow fields were spatially averaged to generate a field of volumetric source terms for each conservative flow variable. Solutions were also obtained using three coarse grids having streamwise and spanwise grid spacings of 3d, 1d, and d/3. These coarse grid solutions used the integrated hole exit mass, momentum, energy, and turbulence quantities from the detailed solutions as volumetric source terms. It is shown that a uniform source term addition over a distance from the wall on the order of the hole diameter is able to predict adiabatic film effectiveness better than a near-wall source term model, while strictly enforcing correct values of integrated boundary layer quantities.

Heidmann, James D.↗
NASA NTRS1973

Application of matched asymptotic expansions to lunar and interplanetary trajectories. Volume 1: Technical discussion

Previously published asymptotic solutions for lunar and interplanetary trajectories have been modified and combined to formulate a general analytical solution to the problem on N-bodies. The earlier first-order solutions, derived by the method of matched asymptotic expansions, have been extended to second order for the purpose of obtaining increased accuracy. The derivation of the second-order solution is summarized by showing the essential steps, some in functional form. The general asymptotic solution has been used as a basis for formulating a number of analytical two-point boundary value solutions. These include earth-to-moon, one- and two-impulse moon-to-earth, and interplanetary solutions. The results show that the accuracies of the asymptotic solutions range from an order of magnitude better than conic approximations to that of numerical integration itself. Also, since no iterations are required, the asymptotic boundary value solutions are obtained in a fraction of the time required for comparable numerically integrated solutions. The subject of minimizing the second-order error is discussed, and recommendations made for further work directed toward achieving a uniform accuracy in all applications.

Lancaster, J. E.↗
NASA NTRS1971

Computational investigations of low-discrepancy point sets

The quasi-Monte Carlo method of integration offers an attractive solution to the problem of evaluating integrals in a large number of dimensions; however, the associated error bounds are difficult to obtain theoretically. Since these bounds are associated with the L2 discrepancy of the set of points used in the integration. Numerical calculations of the L2 discrepancy for several types of quasi-Monte Carlo formulae are presented.

Warnock, T. T.↗
NASA NTRS1971

Interactive program for analysis and design problems in advanced composites technology

During the past year an experimental program in the fracture of advanced fiber composites has been completed. The experimental program has given direction to additional experimental and theoretical work. A synthesis program for designing low weight multifastener joints in composites is proposed, based on extensive analytical background. A number of failed joints have been thoroughly analyzed to evaluate the failure hypothesis used in the synthesis procedure. Finally, a new solution is reported for isotropic and anisotropic laminates using the boundary-integral method. The solution method offers significant savings of computer core and time for important problems.

Cruse, T. A.↗
NASA NTRS1981

A solution for two-dimensional Fredholm integral equations of the second kind with periodic, semiperiodic, or nonperiodic kernels

A numerical scheme for solving two dimensional Fredholm integral equations of the second kind is developed. The proof of the convergence of the numerical scheme is shown for three cases: the case of periodic kernels, the case of semiperiodic kernels, and the case of nonperiodic kernels. Applications to the incompressible, stationary Navier-Stokes problem are of primary interest.

Gabrielsen, R. E.↗
NASA NTRS1984

Random equations in aerodynamics

Literature was reviewed to identify aerodynamic models which might be treated by probablistic methods. The numerical solution of some integral equations that arise in aerodynamical problems were investigated. On the basis of the numerical studies a qualitative theory of random integral equations was developed to provide information on the behavior of the solutions of these equations (in particular, boundary and asymptotic behavior, and stability) and their statistical properties without actually obtaining explicit solutions of the equations.

Bharucha-Reid, A. T.↗
NASA NTRS1985

A procedure for obtaining explicit solutions to a class of Fredholm integral equations

A Fredholm integral equation is solved in closed form by reducing it to a standard boundary value problem for a sectionally analytic function. The procedure maybe extended to more general classes of functions whose Fourier transforms only exist along certain lines in the complex k-plane. Examples are used to illustrate the results.

Goldstein, M. E.↗