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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

INVESTIGATING THE FEASIBILITY AND STABILITY FOR MODELING ACOUSTIC WAVE SCATTERING USING A TIME-DOMAIN BOUNDARY INTEGRAL EQUATION WITH IMPEDANCE BOUNDARY CONDITION

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet and environmentally friendly aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction and are characterized by a frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. This work considers using either an impedance or an admittance (inverse of impedance) boundary condition to allow for acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces. Three acoustic liner models are discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and the Broadband Impedance Model. In both the Helmholtz and Three-Parameter models, liner impedance is specified at a given frequency, whereas the Broadband model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for acoustic liners are derived for each model and coupled with a time-domain boundary integral equation. The scattering solution is obtained iteratively using a boundary element method with constant spatial and third-order temporal basis functions. Time-domain boundary integral equations are unfortunately prone to numerical instabilities due to resonant frequencies resulting from non-trivial solutions in the interior domain. When reformulated with the Burton-Miller method, the instabilities are eliminated. Using a Burton-Miller reformulation, the stability of the boundary element method assuming a liner boundary condition is assessed using eigenvalue analysis. The stability of each liner model is discussed, and it is shown that the Three-Parameter and Broadband models are sufficient for modeling an acoustic liner on the surface of scattering bodies. The Helmholtz model demonstrates strict limitations for stability, whereas the Three-Parameter and Broadband models are stable for most cases. Also included in this work is an assessment of the spatial accuracy of the time-domain boundary element method with respect to the surface element basis functions, as well as a performance study of the numerical algorithm.

computational mathematics↗
NASA NTRS1976

Multistep methods of numerical integration using back-corrections

A class of linear multistep methods is proposed for the solution of the equations of motion of certain dynamical systems encountered in celestial mechanics and astrodynamics. These methods are distinguished from the classical predictor-corrector methods in that they permit 'back-corrections' of the solution to be made. As the integration advances in time, the numerical solution is corrected or improved at certain points in the past. The enhanced numerical stability of these methods allows the meaningful application of high-order algorithms. Consequently, step sizes larger than those attainable with the classical methods may be adopted, and greater overall efficiency may be realized. These methods are applied to the problem of determining the orbit of an artificial satellite, and the results are compared with those obtained using classical methods.

Feagin, T.↗
NASA NTRS1985

Computation of transonic vortex flows past delta wings Integral equation approach

The steady full-potential equation is written in the form of Poisson's equation, and the solution of the velocity field is expressed in terms of an integral equation. The solution consists of a surface integral of vorticity distribution on the wing and its free-vortex sheets and a volume integral of source distribution within a volume around the wing and its free-vortex sheets. The solution is obtained through successive iteration cycles. The source distribution is computed by using a mixed finite-difference scheme of the Murman-Cole type. The method is applied to delta wings. Numerical examples show that a conical shock is captured on the suction side of the wing. It is attached to the lower surface of the leading-edge vortex but does not necessarily reach to the wing surface.

Kandil, O. A.↗
NASA NTRS1984

Obtaining Runge-Kutta Solutions Between Time Steps

New interpolation method used with existing Runge-Kutta algorithms. Algorithm evaluates solution at intermediate point within integration step. Only few additional computations required to produce intermediate solution data. Runge-Kutta method provides accurate solution with larger time steps than allowable in other methods.

Horn, M. K.↗
NASA NTRS1992

Analytic solution of the thickness problem of a rectangular wing in steady subsonic flow

An analytic solution of the thickness problem of a rectangular wing with parabolic airfoil section in three-dimensional flow is presented. The free-air solution is obtained by integrating the equation of the axial perturbation velocity. The Prandtl-Glauert rule can be used to derive the subsonic solution. Parts of the free-air solution are verified by taking the limit of the axial peturbation velocity on the model surface as the wing span goes to infinity. Pressure coefficients are studied for a selected wing geometry in the flow field. The solution of the thickness problem of a rectangular wing in a rectangular wind tunnel is derived using the free-air solution and the method of images.

Ulbrich, N.↗
NASA NTRS1990

Boundary element analysis of sound scattered by a moving surface

A solution for the acoustic field scattered from a uniformly moving rigid body in the presence of a harmonic incident source has been obtained using a boundary integral method. A derivation of the Kirchhoff formula given by Farassat and Myers (1988) for moving surfaces forms the basis for the analysis, and the development of a boundary integral method for the solution of scattering problems from moving rigid bodies is described. Finite elements are used in conjunction with the Galerkin method in order to solve the integral equation that results from the Kirchhoff formula when the observer point is placed on the moving body surface. Once appropriate surface field values are known they are inserted back into the formula in order to predict the field scattered off the body. Tests, including the so called superposition method, are carried out in order to validate the technique and to establish some confidence in its accuracy. Application of the superposition method to moving bodies is presented, and results of the two approaches are discussed. Sample calculations of scattering from a simple body are presented to illustrate the effects of variations in relevant parameters.

Myers, M. K.↗
NASA NTRS1983

Calculations of atmospheric refraction for spacecraft remote-sensing applications

Analytical solutions to the refraction integrals appropriate for ray trajectories along slant paths through the atmosphere are derived in this paper. This type of geometry is commonly encountered in remote-sensing applications utilizing an occultation technique. The solutions are obtained by evaluating higher-order terms from expansion of the refraction integral and are dependent on the vertical temperature distributions. Refraction parameters such as total refraction angles, air masses, and path lengths can be accurately computed. It is also shown that the method can be used for computing refraction parameters in astronomical refraction geometry for large zenith angles.

Chu, W. P.↗
NASA NTRS1987

Investigation of the feasibility of an analytical method of accounting for the effects of atmospheric drag on satellite motion

An analytic technique for accounting for the joint effects of Earth oblateness and atmospheric drag on close-Earth satellites is investigated. The technique is analytic in the sense that explicit solutions to the Lagrange planetary equations are given; consequently, no numerical integrations are required in the solution process. The atmospheric density in the technique described is represented by a rotating spherical exponential model with superposed effects of the oblate atmosphere and the diurnal variations. A computer program implementing the process is discussed and sample output is compared with output from program NSEP (Numerical Satellite Ephemeris Program). NSEP uses a numerical integration technique to account for atmospheric drag effects.

Bozeman, Robert E.↗
NASA NTRS1973

Program for solution of ordinary differential equations

A program for the solution of linear and nonlinear first order ordinary differential equations is described and user instructions are included. The program contains a new integration algorithm for the solution of initial value problems which is particularly efficient for the solution of differential equations with a wide range of eigenvalues. The program in its present form handles up to ten state variables, but expansion to handle up to fifty state variables is being investigated.

Sloate, H.↗
NASA NTRS1996

Analytic corrections to CFD heating predictions accounting for changes in surface catalysis

Integral boundary-layer solution techniques applicable to the problem of determining aerodynamic heating rates of hypersonic vehicles in the vicinity of stagnation points and windward centerlines are briefly summarized. A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (but time intensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations is described. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. It allows a designer at a workstation to address two questions, given a single CFD solution. (1) How much does the heating change for a thermal protection system with different catalytic properties than was used in the original CFD solution? (2) How does the heating change at the interface of two different TPS materials with an abrupt change in catalytic efficiency? The answer to the second question is particularly important, because abrupt changes from low to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption.

Gnoffo, Peter A.↗
NASA NTRSDate not supplied

Wall-Modeled Large-Eddy Simulations for High-Lift Configurations using FUN3D

Wall-modeled large-eddy simulation (WMLES) capability has recently been implemented into FUN3D, an unstructured, node-centered, finite-volume solver developed at the NASA Langley Research Center. In this paper, WMLES is assessed for two configurations that are representative for high-lift applications. The first configuration is a nominal two-dimensional multielement airfoil that has been extensively studied in the literature. WMLES solutions are computed for four angles of attack and compared with previously reported solutions. Good agreement of integrated forces, surface pressures, and boundary-layer velocity profiles is shown with available experimental data especially at lower angles of attack. WMLES solutions are also computed for the NASA High-Lift Common Research Model over a large range of angles of attack. Forces, pitching moments, and pressure distributions are favorably compared with the experimental data up to the maximum lift, including the angle of attack where the maximum lift is obtained experimentally. Eddy visualization techniques of q-criterion and density-gradient magnitude illustrate the resolved content.

High-Lift flows↗
DOE OSTI2026

Field Validation of a Grid-Interactive Efficient Building Software Solution

The U.S. General Services Administration's (GSA's) Green Proving Ground (GPG) program, in partnership with the National Laboratory of the Rockies (NLR), completed a field study of a Grid-Interactive Efficient Buildings (GEB) software solution. The study focused on a single testbed facility to test the GEB functionality of the software solution, along with other features. The testbed facility - a courthouse - is a common building type in GSA's vast building portfolio, offering potentially impactful findings on a scalable level. The study evaluated Prescriptive Data's technology, Nantum OS, a connected building operating system ("GEB Solution") which aggregates multiple sources of previously siloed building data and combines that data with external sources, such as weather information or utility signals, into a single integrated platform. A GEB Solution is a type of Energy Management Information System (EMIS). EMIS is defined as a system of devices, data services, and software applications that communicates with any building system or third-party data source to aggregate and transform data into new capabilities to aid in the optimization of energy use at the building, campus, or agency level. This specific GEB Solution is an EMIS with ASO, automated system optimization, offering supervisory control of certain aspects of the Building Automation System (BAS). Multiple features were evaluated including, but not limited to, Continuous Demand Management to avoid setting new monthly kilowatt (kW) peaks, energy efficiency for reduction of kilowatt hours (kWh) and natural gas consumption, and automated demand response (ADR) for purposes of lowering demand during a utility called Demand Response (DR) event. The testbed facility was the Foley Federal Building and US Courthouse ("Foley Federal Building") located in Las Vegas, NV. This is a 209,496 sq. ft. building constructed in the 1960s with major renovations in 2004. The facility was a good candidate due to the large prevalence of office and courthouse spaces in the GSA portfolio of buildings. It also has many features which allow integration into and control of the building and a strong facilities team to assist with the study. Quantitative and qualitative performance objectives were developed using GSA's GPG GEB project template along with input from the vendor and building facility staff; these are outlined in Table 1. The quantitative performance objectives focused on continuous demand management, energy efficiency, and automated demand response. The qualitative performance objectives focused on the ease of installation and commissioning as well as the operability of the GEB solution. Other performance metrics that are reported on include carbon reduction, cost effectiveness, and occupant acceptance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1975

A general algorithm for relating ground trajectory distance, elapsed flight time, and aircraft airspeed and its application to 4-D guidance

A general solution using an elliptic integral approximation which relates flight time, aircraft airspeed, and ground distance on straight-line and circular-arc trajectory segments is developed. The solution procedure is applicable to both constant and accelerating aircraft flight. In addition, wind shear including both magnitude and heading change is incorporated in the solution. The solution equations are used in a four-dimensional (4-D) control algorithm where both flight time and final airspeed are specified. The results show that the algorithm converges rapidly and accurately.

Foudriat, E. C.↗
NASA NTRS1981

New approach to the solution of large, full matrix equations

A new approach to the solution of matrix equations resulting from integral equations is presented and applied to the solution of two-dimensional Neumann problems describing the inviscid, incompressible flow past an airfoil. The problem is reformulated in terms of a preselected set of mode functions giving an equivalent matrix equation to be solved for the mode-function expansion coefficients. Because of the inherent smoothness of the original problem, the coefficient problem can be solved approximately without significantly affecting the accuracy of the final solution. Very promising two-dimensional results are obtained and the extension of the method to three-dimensional problems is investigated. On the basis of these results it is shown that the computing time for the matrix solution for a large three-dimensional panel method calculation could be reduced by an order of magnitude compared with that required for a direct solution.

Clark, R. W.↗
NASA NTRS1984

The role of optimization in structural model refinement

To evaluate the role that optimization can play in structural model refinement, it is necessary to examine the existing environment for the structural design/structural modification process. The traditional approach to design, analysis, and modification is illustrated. Typically, a cyclical path is followed in evaluating and refining a structural system, with parallel paths existing between the real system and the analytical model of the system. The major failing of the existing approach is the rather weak link of communication between the cycle for the real system and the cycle for the analytical model. Only at the expense of much human effort can data sharing and comparative evaluation be enhanced for the two parallel cycles. Much of the difficulty can be traced to the lack of a user-friendly, rapidly reconfigurable engineering software environment for facilitating data and information exchange. Until this type of software environment becomes readily available to the majority of the engineering community, the role of optimization will not be able to reach its full potential and engineering productivity will continue to suffer. A key issue in current engineering design, analysis, and test is the definition and development of an integrated engineering software support capability. The data and solution flow for this type of integrated engineering analysis/refinement system is shown.

Lehman, L. L.↗
NASA NTRS1973

A note on relative motion in the general three-body problem.

It is shown that the equations of the general three-body problem take on a very symmetric form when one considers only their relative positions, rather than position vectors relative to some given coordinate system. From these equations one quickly surmises some well known classical properties of the three-body problem, such as the first integrals and the equilateral triangle solutions. Some new Lagrangians with relative coordinates are also obtained. Numerical integration of the new equations of motion is about 10% faster than with barycentric or heliocentric coordinates.

Broucke, R.↗
NASA NTRS1990

CFIE MM solution for TE and TM incidence on a 2-D conducting body with dielectric filled cavity

The problem of determining the scattering cross section of an arbitrarily shaped two-dimensional conducting body with an arbitrarily shaped dielectric-filled cavity is considered. The problem is solved using a method-of-moments (MM) solution for the combined field integral equations (CFIE). The particular form of the method of moments solution used has here uses a minimum number of expansion coefficients. Results are given for TE and TM incident waves.

Goggans, Paul M.↗
NASA NTRS1985

Ergonomic problems regarding the interactive touch input via screens in onboard and ground-based flight control

A significant problem concerning the integration of display and switching functions is related to the fact that numerous informative data which have to be processed by man must be read from only a few display devices. A satisfactory ergonomic design of integrated display devices and keyboards is in many cases difficult, because not all functions which can be displayed and selected are simultaneously available. A technical solution which provides an integration of display and functional elements on the basis of the highest flexibility is obtained by using a cathode ray tube with a touch-sensitive screen. The employment of an integrated data input/output system is demonstrated for the cases of onboard and ground-based flight control. Ergonomic studies conducted to investigate the suitability of an employment of touch-sensitive screens are also discussed.

Holzhausen, K. P.↗