Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integral 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1971

The spectrum of steady state turbulent convection.

Based on Heisenberg's statistical theory of turbulence, a model for steady state turbulent convection is herein proposed, and on the basis of this model, equations for the energy spectrum for steady state turbulent convection are derived. The spectrum is obtained from the solution of a nonlinear integral equation. After the integral equation is brought into a universally valid nondimensional form, it is transformed into a nonlinear first order differential equation to be solved numerically, with the Rayleigh number appearing as the only parameter. The energy spectrum has a substantial deviation from the Kolmogoroff law, as a result of the buoyancy force acting on the rising and falling eddies. The presented theory may be applicable to convection in planetary and stellar atmospheres wherein the radiative heat transport is small.

Winterberg, F.↗
DOE OSTI2021

Meshless discretization of the discrete-ordinates transport equation with integration based on Voronoi cells

The time-dependent, gray, linear radiation transport equation is discretized using the meshless local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The discretization includes streamline-upwind Petrov-Galerkin stabilization to prevent oscillations and improve numerical conditioning. The angular quadrature is selectively refineable to increase angular resolution in chosen directions. The time discretization is done using backward Euler. The transport solve for each direction and the solve for the scattering source are both done using Krylov iterative methods. The results indicate first-order convergence in time and second-order convergence in space for linear reproducing kernels.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Profiles of emission lines in Be stars.

Study of the broadening functions resulting from a gaseous ring in circular motion around a star according to Kepler's law. When the distribution in the gaseous ring has a circular symmetry in the equatorial plane, the broadening profile is related to the surface density distribution along the radius by an integral equation which can be transformed into the Abel integral equation and solved analytically. Profiles corresponding to gaseous rings with a uniform density distribution but different widths are used to illustrate the general properties of the profile broadened by the gaseous ring in circular motion. The emission profile has also been studied for cases in which the circular motion does not follow Kepler's law.

Huang, S. S.↗
NASA NTRS1989

Resonance in cylindrical-rectangular and wraparound microstrip structures

A rigorous analysis of the resonance frequency problem of both the cylindrical-rectangular and the wraparound microstrip structure is presented. The problem is formulated in terms of a set of vector integral equations. Using Galerkin's method to solve the integral equations, the complex resonance frequencies are studied with sinusoidal basis functions which incorporate the edge singularity. The complex resonance frequencies are computed using a perturbation approach. Modes suitable for resonator or antenna applications are investigated. The edge singularity of the patch current is shown to have no significant effect on the accuracy of the results. It is shown that the HE10 modes of the cylindrical-rectangular and wraparound patches are more appropriate for resonator applications. The HE01 and TE01 modes of the cylindrical-rectangular and wraparound patches, respectively, are efficient radiating modes.

Ali, Sami M.↗
NASA NTRS1996

A Hybrid Boundary Element-Finite Volume Method for Unsteady Transonic Airfoil Flows

A hybrid boundary element finite volume method for unsteady transonic flow computation has been developed. In this method, the unsteady Euler equations in a moving frame of reference are solved in a small embedded domain (inner domain) around the airfoil using an implicit finite volume scheme. The unsteady full-potential equation, written in the same frame of reference and in the form of the Poisson equation. is solved in the outer domain using the integral equation boundary element method to provide the boundary conditions for the inner Euler domain. The solution procedure is a time-accurate stepping procedure, where the outer boundary conditions for the inner domain are updated using the integral equation -- boundary element solution over the outer domain. The method is applied to unsteady transonic flows around the NACA0012 airfoil undergoing pitching oscillation and ramp motion. The results are compared with those of an implicit Euler equation solver, which is used throughout a large computational domain, and experimental data.

Hu, Hong↗
DOE OSTI2020

Method of Green’s Potentials for Elliptic PDEs in Domains with Random Apertures

Problems with topological uncertainties appear in many fields ranging from nano-device engineering to the design of bridges. In many of such problems, a part of the domains boundaries is subjected to random perturbations making inefficient conventional schemes that rely on discretization of the whole domain. Here, we study elliptic PDEs in domains with boundaries comprised of a deterministic part and random apertures, and apply the method of modified potentials with Green’s kernels defined on the deterministic part of the domain. This approach allows to reduce the dimension of the original differential problem by reformulating it as a boundary integral equation posed on the random apertures only. The multilevel Monte Carlo method is then applied to this modified integral equation and its optimal ϵ -2 asymptotical complexity is shown. Finally, we provide the qualitative analysis of the proposed technique and support it with numerical results.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2011

Excitation of a Parallel Plate Waveguide by an Array of Rectangular Waveguides

This work addresses the problem of excitation of a parallel plate waveguide by an array of rectangular waveguides that arises in applications such as the continuous transverse stub (CTS) antenna and dual-polarized parabolic cylindrical reflector antennas excited by a scanning line source. In order to design the junction region between the parallel plate waveguide and the linear array of rectangular waveguides, waveguide sizes have to be chosen so that the input match is adequate for the range of scan angles for both polarizations. Electromagnetic wave scattered by the junction of a parallel plate waveguide by an array of rectangular waveguides is analyzed by formulating coupled integral equations for the aperture electric field at the junction. The integral equations are solved by the method of moments. In order to make the computational process efficient and accurate, the method of weighted averaging was used to evaluate rapidly oscillating integrals encountered in the moment matrix. In addition, the real axis spectral integral is evaluated in a deformed contour for speed and accuracy. The MoM results for a large finite array have been validated by comparing its reflection coefficients with corresponding results for an infinite array generated by the commercial finite element code, HFSS. Once the aperture electric field is determined by MoM, the input reflection coefficients at each waveguide port, and coupling for each polarization over the range of useful scan angles, are easily obtained. Results for the input impedance and coupling characteristics for both the vertical and horizontal polarizations are presented over a range of scan angles. It is shown that the scan range is limited to about 35 for both polarizations and therefore the optimum waveguide is a square of size equal to about 0.62 free space wavelength.

Rengarajan, Sembiam↗
DOE OSTI2025

Evaluation of Maximum Allowable Working Pressure and Svensson Burst Pressure Recommended in API 579-1 2021 Edition

Abstract API 579-1/ASME FFS-1 2021 Edition provides the minimum wall thickness, the maximum allowable working pressure (MAWP), and the membrane stress equations for thin and thick-walled cylindrical shells subject to internal pressure in Section 2C.3.3.1 of Appendix 2C – Thickness, MAWP, and Stress Equations for an FFS Assessment. The minimum wall thickness and MAWP are determined using the hoop stress and the Tresca yield criterion. Section 2C.7 – Estimation of Burst Pressure newly added the Svensson method for calculating burst pressure of cylindrical shells under internal pressure, where the plastic yielding is characterized by the von Mises yield criterion. For thin-walled cylinders, the von Mises flow solution of burst pressure in Equation (2C.179) was recommended. For thick-walled cylinders, an implicit burst pressure solution in an integral equation (2C.176) was recommended. But this integral equation is inconvenient to use in practice. It is well known that the classic plasticity theory includes the Tresca and von Mises yield criteria, with the Tresca criterion predicting a lower bound solution and the von Mises criteria predicting an upper bound solution. In addition, the present author developed an average shear stress yield criterion that can determine more accurate limit and burst pressures for thin and thick-walled cylinders. This work uses these three yield criteria to evaluate the minimum required wall thickness, MAWP and Svensson burst pressure recommended in the API 579 code.

burst pressure↗
DOE OSTI2025

Evaluation of Maximum Allowable Working Pressure and Svensson Burst Pressure Recommended in API 579-1 2021 Edition

ABSTRACT API 579-1/ASME FFS-1 2021 Edition provides the minimum wall thickness, the maximum allowable working pressure (MAWP), and the membrane stress equations for thin and thick-walled cylindrical shells subject to internal pressure in Section 2C.3.3.1 of Appendix 2C – Thickness, MAWP, and Stress Equations for an FFS Assessment. The minimum wall thickness and MAWP are determined using the hoop stress and the Tresca yield criterion. Section 2C.7 – Estimation of Burst Pressure newly added the Svensson method for calculating burst pressure of cylindrical shells under internal pressure, where the plastic yielding is characterized by the von Mises yield criterion. For thin-walled cylinders, the von Mises flow solution of burst pressure in Equation (2C.179) was recommended. For thick-walled cylinders, an implicit burst pressure solution in an integral equation (2C.176) was recommended. But this integral equation is inconvenient to use in practice. It is well known that the classic plasticity theory includes the Tresca and von Mises yield criteria, with the Tresca criterion predicting a lower bound solution and the von Mises criteria predicting an upper bound solution. In addition, the present author developed an average shear stress yield criterion that can determine more accurate limit and burst pressures for thin and thick-walled cylinders. This work uses these three yield criteria to evaluate the minimum required wall thickness, MAWP and Svensson burst pressure recommended in the API 579 code.

burst pressure↗
NASA NTRS1984

The influence of a wind tunnel on helicopter rotational noise: Formulation of analysis

An analytical model is discussed that can be used to examine the effects of wind tunnel walls on helicopter rotational noise. A complete physical model of an acoustic source in a wind tunnel is described and a simplified version is then developed. This simplified model retains the important physical processes involved, yet it is more amenable to analysis. The simplified physical model is then modeled as a mathematical problem. An inhomogeneous partial differential equation with mixed boundary conditions is set up and then transformed into an integral equation. Details of generating a suitable Green's function and integral equation are included and the equation is discussed and also given for a two-dimensional case.

Mosher, M.↗
NASA NTRS1997

Frictionless Contact of Multilayered Composite Half Planes Containing Layers With Complex Eigenvalues

A previously developed local-global stiffness matrix methodology for the response of a composite half plane, arbitrarily layered with isotropic, orthotropic or monoclinic plies, to indentation by a rigid parabolic punch is further extended to accommodate the presence of layers with complex eigenvalues (e.g., honeycomb or piezoelectric layers). First, a generalized plane deformation solution for the displacement field in an orthotropic layer or half plane characterized by complex eigenvalues is obtained using Fourier transforms. A local stiffness matrix in the transform domain is subsequently constructed for this class of layers and half planes, which is then assembled into a global stiffness matrix for the entire multilayered half plane by enforcing continuity conditions along the interfaces. Application of the mixed boundary condition on the top surface of the half plane indented by a rigid punch results in an integral equation for the unknown pressure in the contact region. The integral possesses a divergent kernel which is decomposed into Cauchy-type and regular parts using the asymptotic properties of the local stiffness matrix and a relationship between Fourier and finite Hilbert transform of the contact pressure. The solution of the resulting singular integral equation is obtained using a collocation technique based on the properties of orthogonal polynomials developed by Erdogan and Gupta. Examples are presented that illustrate the important influence of low transverse properties of layers with complex eigenvalues, such as those exhibited by honeycomb, on the load versus contact length response and contact pressure distributions for half planes containing typical composite materials.

Zhang, Wang↗
NASA NTRS1981

General relaxation schemes in multigrid algorithms for higher order singularity methods

Relaxation schemes based on approximate and incomplete factorization technique (AF) are described. The AF schemes allow construction of a fast multigrid method for solving integral equations of the second and first kind. The smoothing factors for integral equations of the first kind, and comparison with similar results from the second kind of equations are a novel item. Application of the MD algorithm shows convergence to the level of truncation error of a second order accurate panel method.

Oskam, B.↗
NASA NTRS1987

Scattering from a microstrip patch

A solution to the problem of plane wave scattering by a rectangular microstrip patch on a grounded dielectric substrate is presented. The model does not include the microstrip feed, and thus doe not include the so-called 'antenna mode' component of the scattering. The solution begins by formulating an electric field integral equation for the surface current density on the microstrip patch. The integral equation is solved using the method of moments. Computed data for the patch radar cross section (RCS) is found to be in close agreement with measurements over a broad frequency range. The microstrip RCS versus frequency consists of a number of large peaks which are identified as impedance or pattern factor resonance peaks.

Newman, Edward H.↗
NASA NTRS1993

Computer code for scattering from impedance bodies of revolution. Part 3: Surface impedance with s and phi variation. Analytical and numerical results

The third phase of the development of the computer codes for scattering by coated bodies that has been part of an ongoing effort in the Electromagnetics Laboratory of the Electrical Engineering and Computer Science Department at the University of Illinois at Chicago is described. The work reported discusses the analytical and numerical results for the scattering of an obliquely incident plane wave by impedance bodies of revolution with phi variation of the surface impedance. Integral equation formulation of the problem is considered. All three types of integral equations, electric field, magnetic field, and combined field, are considered. These equations are solved numerically via the method of moments with parametric elements. Both TE and TM polarization of the incident plane wave are considered. The surface impedance is allowed to vary along both the profile of the scatterer and in the phi direction. Computer code developed for this purpose determines the electric surface current as well as the bistatic radar cross section. The results obtained with this code were validated by comparing the results with available results for specific scatterers such as the perfectly conducting sphere. Results for the cone-sphere and cone-cylinder-sphere for the case of an axially incident plane were validated by comparing the results with the results with those obtained in the first phase of this project. Results for body of revolution scatterers with an abrupt change in the surface impedance along both the profile of the scatterer and the phi direction are presented.

Uslenghi, Piergiorgio L. E.↗
NASA NTRS2005

Diffraction of Harmonic Flexural Waves in a Cracked Elastic Plate Carrying Electrical Current

The scattering effect of harmonic flexural waves at a through crack in an elastic plate carrying electrical current is investigated. In this context, the Kirchhoffean bending plate theory is extended as to include magnetoelastic interactions. An incident wave giving rise to bending moments symmetric about the longitudinal z-axis of the crack is applied. Fourier transform technique reduces the problem to dual integral equations, which are then cast to a system of two singular integral equations. Efficient numerical computation is implemented to get the bending moment intensity factor for arbitrary frequency of the incident wave and of arbitrary electrical current intensity. The asymptotic behaviour of the bending moment intensity factor is analysed and parametric studies are conducted.

Ambur, Damodar R.↗
DOE OSTI2024

Axisymmetric hybrid Vlasov equilibria with applications to tokamak plasmas

We derive axisymmetric equilibrium equations in the context of the hybrid Vlasov model with kinetic ions and massless fluid electrons, assuming isothermal electrons and deformed Maxwellian distribution functions for the kinetic ions. The equilibrium system comprises a Grad–Shafranov partial differential equation and an integral equation. These equations can be utilized to calculate the equilibrium magnetic field and ion distribution function, respectively, for given particle density or given ion and electron toroidal current density profiles. The resulting solutions describe states characterized by toroidal plasma rotation and toroidal electric current density. Additionally, due to the presence of fluid electrons, these equilibria also exhibit a poloidal current density component. This is in contrast to the fully kinetic Vlasov model, where axisymmetric Jeans equilibria can only accommodate toroidal currents and flows, given the absence of a third integral of the microscopic motion.

Physics↗