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 487 records · Page 27

The F(N) method for the one-angle radiative transfer equation applied to plant canopies

The paper presents a semianalytical solution method, called the F(N) method, for the one-angle radiative transfer equation in slab geometry. The F(N) method is based on two integral equations specifying the intensities exiting the boundaries of the vegetation canopy; the solution is obtained through an expansion in a set of basis functions with expansion coefficients to be determined. The advantage of this method is that it avoids spatial truncation error entirely because it requires discretization only in the angular variable.

Ganapol, B. D.↗

Introduction to Generalized Functions with Applications in Aerodynamics and Aeroacoustics

Generalized functions have many applications in science and engineering. One useful aspect is that discontinuous functions can be handled as easily as continuous or differentiable functions and provide a powerful tool in formulating and solving many problems of aerodynamics and acoustics. Furthermore, generalized function theory elucidates and unifies many ad hoc mathematical approaches used by engineers and scientists. We define generalized functions as continuous linear functionals on the space of infinitely differentiable functions with compact support, then introduce the concept of generalized differentiation. Generalized differentiation is the most important concept in generalized function theory and the applications we present utilize mainly this concept. First, some results of classical analysis, are derived with the generalized function theory. Other applications of the generalized function theory in aerodynamics discussed here are the derivations of general transport theorems for deriving governing equations of fluid mechanics, the interpretation of the finite part of divergent integrals, the derivation of the Oswatitsch integral equation of transonic flow, and the analysis of velocity field discontinuities as sources of vorticity. Applications in aeroacoustics include the derivation of the Kirchhoff formula for moving surfaces, the noise from moving surfaces, and shock noise source strength based on the Ffowcs Williams-Hawkings equation.

Farassat, F.↗

Artificial Boundary Conditions Based on the Difference Potentials Method

While numerically solving a problem initially formulated on an unbounded domain, one typically truncates this domain, which necessitates setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The issue of setting the ABC's appears to be most significant in many areas of scientific computing, for example, in problems originating from acoustics, electrodynamics, solid mechanics, and fluid dynamics. In particular, in computational fluid dynamics (where external problems present a wide class of practically important formulations) the proper treatment of external boundaries may have a profound impact on the overall quality and performance of numerical algorithms. Most of the currently used techniques for setting the ABC's can basically be classified into two groups. The methods from the first group (global ABC's) usually provide high accuracy and robustness of the numerical procedure but often appear to be fairly cumbersome and (computationally) expensive. The methods from the second group (local ABC's) are, as a rule, algorithmically simple, numerically cheap, and geometrically universal; however, they usually lack accuracy of computations. In this paper we first present a survey and provide a comparative assessment of different existing methods for constructing the ABC's. Then, we describe a relatively new ABC's technique of ours and review the corresponding results. This new technique, in our opinion, is currently one of the most promising in the field. It enables one to construct such ABC's that combine the advantages relevant to the two aforementioned classes of existing methods. Our approach is based on application of the difference potentials method attributable to V. S. Ryaben'kii. This approach allows us to obtain highly accurate ABC's in the form of certain (nonlocal) boundary operator equations. The operators involved are analogous to the pseudodifferential boundary projections first introduced by A. P. Calderon and then also studied by R. T. Seeley. The apparatus of the boundary pseudodifferential equations, which has formerly been used mostly in the qualitative theory of integral equations and PDE'S, is now effectively employed for developing numerical methods in the different fields of scientific computing.

Tsynkov, Semyon V.↗

Higher Order Bases in a 2D Hybrid BEM/FEM Formulation

The advantages of using higher order, interpolatory basis functions are examined in the analysis of transverse electric (TE) plane wave scattering by homogeneous, dielectric cylinders. A boundary-element/finite-element (BEM/FEM) hybrid formulation is employed in which the interior dielectric region is modeled with the vector Helmholtz equation, and a radiation boundary condition is supplied by an Electric Field Integral Equation (EFIE). An efficient method of handling the singular self-term arising in the EFIE is presented. The iterative solution of the partially dense system of equations is obtained using the Quasi-Minimal Residual (QMR) algorithm with an Incomplete LU Threshold (ILUT) preconditioner. Numerical results are shown for the case of an incident wave impinging upon a square dielectric cylinder. The convergence of the solution is shown versus the number of unknowns as a function of the completeness order of the basis functions.

Fink, Patrick W.↗

An Overview of Gemma FY20 Verification Activities

Gemma verification activities for FY20 can be divided into three categories: the development of specialized quadrature rules, initial progress towards the development of manufactured solutions for code verification, and automated code-verification testing. In the method-of-moments implementation of the electric-field integral equation, the presence of a Green’s function in the four-dimensional integrals yields singularities in the integrand when two elements are nearby. To address these challenges, we have developed quadrature rules to integrate the functions through which the singularities can be characterized. Code verification is necessary to develop confidence in the implementation of the numerical methods in Gemma. Therefore, we have begun investigating the use of manufactured solutions to more thoroughly verify Gemma. Manufactured solutions provide greater flexibility for testing aspects of the code; however, the aforementioned singularities provide challenges, and existing work is limited in rigor and quantity. Finally, we have implemented automated code-verification testing using the VVTest framework to automate the mesh refinement and execution of a Gemma simulation to generate mesh convergence data. This infrastructure computes the observed order of accuracy from these data and compares it with the theoretical order of accuracy to either develop confidence in the implementation of the numerical methods or detect coding errors.

97 MATHEMATICS AND COMPUTING↗

Three-dimensional marginal separation

The three dimensional marginal separation of a boundary layer along a line of symmetry is considered. The key equation governing the displacement function is derived, and found to be a nonlinear integral equation in two space variables. This is solved iteratively using a pseudo-spectral approach, based partly in double Fourier space, and partly in physical space. Qualitatively, the results are similar to previously reported two dimensional results (which are also computed to test the accuracy of the numerical scheme); however quantitatively the three dimensional results are much different.

Duck, Peter W.↗

Hybrid Modeling for Complex Systems Health Management

The research work presents application of hybrid physics-informed machine learning to a representative electric powertrain for unmanned aerial vehicles. The model is composed of physics-derived and empirical equations, integrated with connected networks that are strategically placed within the model to substitute equations that are subject to large uncertainty. Polynomial fit driven by heuristics or empirical observations can be substituted by more flexible networks that can minimize the error between model predictions and observations without being restricted to a predefined functional form. This modeling strategy allows training of networks deep inside the model and unknown parameters in a single learning stage. The powertrain model consists of Li-ion batteries, electronic speed controller with pulse-width modulation, and brush-less DC motor with connected propeller. Results obtained from combination of laboratory and simulation tests are discussed in this work.

PINNS↗

Physics Informed Neural Nets for Systems Health Management

To facilitate and solve the prediction problem, awareness of the current health state of the system is key, since it is necessary to perform condition-based predictions. To accurately predict the future state of any system, it is required to possess knowledge of its current health state and future operational conditions. Development in data-driven algorithms in regression of complex nonlinear functions and classification tasks have generated a growing interest in artificial intelligence for industrial applications. Complex multi-physics models as well as digital twins, once purely built on physics and corresponding simplified lumped parameter iterations, can now benefit from machine learning algorithms to mitigate the lack of understanding of some complex behavior. The research work presents application of physics-informed neural nets application to a representative electric powertrain for unmanned aerial vehicles. The model is composed of physics-derived and empirical equations, integrated with connected networks that are strategically placed within the model to substitute equations that are subject to large uncertainty. Polynomial fit driven by heuristics or empirical observations can be substituted by more flexible networks that can minimize the error between model predictions and observations without being restricted to a predefined functional form. This modeling strategy allows training of networks deep inside the model and unknown parameters in a single learning stage.

Physics Informed↗

Model Based Approaches for Fault Detection, Prognostics, Decision Making in Complex Systems

The presentation discusses application of model based approaches to complex systems. The model is composed of physics-derived and empirical equations, integrated with connected networks that are strategically placed within the model to substitute equations that are subject to large uncertainty. Polynomial fit driven by heuristics or empirical observations can be substituted by more flexible networks that can minimize the error between model predictions and observations without being restricted to a predefined functional form. This modeling strategy allows training of networks deep inside the model and unknown parameters in a single learning stage.

Physics Informed↗

Stability in integrable nonlocal nonlinear equations

Recently a variety of nonlocal integrable systems has been introduced that besides fields located at particular space-time points simultaneously also contain fields that are located at different, but symmetrically related, points. Here we investigate different types of soliton solutions with regard to their stability against linear perturbations obtained for the nonlocal version of the Hirota/nonlinear Schrödinger equation and the so-called Alice and Bob versions of the Korteweg-de Vries and Bousinesq equations. In this work, we encounter different types of scenarios: Solition solutions that are linearly stable or unstable and also solutions that change their stability properties depending on the parameter regime they are in.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A new approach to the solution of large, full matrix equations: A two-dimensional potential flow feasibility study

An approach to the solution of matrix problems resulting from integral equations of mathematical physics is presented. Based on the inherent smoothness in such equations, the problem is reformulated using a set of orthogonal basis vectors, leading to an equivalent coefficient problem which can be of lower order without significantly impairing the accuracy of the solution. This approach was evaluated using a two-dimensional Neumann problem describing the inviscid, incompressible flow over an airfoil. Two different kinds of mode functions were investigated, namely eigenfunction series and Fourier series. The method using Fourier series was found preferable. It uses all of the coefficients from a Fast Fourier Transform algorithm in an approximate method which exploits the known structure of the transformed coefficient matrix and very promising results for the flow over a realistic airfoil are obtained. On the basis of the results presented here, an order of magnitude reduction in this computer time can be expected for such problems as compared with the time for a direct matrix solution.

James, R. M.↗

Application of Fourier transforms for microwave radiometric inversions

Existing microwave radiometer technology now provides a suitable method for remote determination of the ocean surface's absolute brightness temperature. To extract the brightness temperature of the water from the antenna temperature, an unstable Fredholm integral equation of the first kind is solved. Fourier transform techniques are used to invert the integral after it is placed into a cross correlation form. Application and verification of the methods to a two-dimensional modeling of a laboratory wave tank system are included. The instability of the ill-posed Fredholm equation is examined and a restoration procedure is included which smooths the resulting oscillations. With the recent availability and advances of fast Fourier transform (FFT) techniques, the method presented becomes very attractive in the evaluation of large quantities of data.

Holmes, J. J.↗

A similarity rule for compressibility and sidewall-boundary-layer effects in two-dimensional wind tunnels

The effect of the sidewall boundary layer on flow in two-dimensional wind tunnels is determined. The small-disturbance and isentropic approximations are made, and the sidewall-boundary-layer dynamics are modeled with the von Karman momentum-integral equation. The effects of the edge-velocity-gradient term in the sidewall momentum integral, which is usually dominant near the model, and the compressibility term are shown to be similar. It is shown that the effect of sidewall suction around the model is not similar to two-dimensional flow. Comparisons with experiment are made to verify the similarity rule.

Barnwell, R. W.↗

A lifting surface theory in rotational flow

The partial differential equation for small disturbance steady rotational flow in three dimensions is solved through an integral equation approach. The solution is obtained by using the method of weighted residuals. Specific applications are directed to wings in nonuniform subsonic parallel streams with velocity varying in vertical and spanwise directions and to airfoils in nonuniform freestream. Comparison with limited known results indicates that the present method is reasonably accurate. Numerical results for the lifting pressure of airfoil, lift, induced drag, and pitching moments of airfoil, lift, induced drag, and pitching moments of elliptic, rectangular, and delta wings in a jet, wake, or monotonic sheared stream are presented. It is shown that, in addition to the effect of local dynamic pressures, a positive velocity gradient tends to enhance the lift.

Shiau, M. J.↗

Shear buckling of specially orthotropic plates with centrally located cutouts

There is significant industry demand for a method of analyzing the shear buckling of composite plates with cutouts. This method should be able to easily accommodate frequent changes in model design; inflexibility being the major drawback to current finite element methods. The approach taken here is broken into two problems, prebuckling and buckling. To solve for the prebuckling stresses, complex variable equations are used in conjunction with boundary collocation. The least squares approach is utilized to improve the accuracy of the results. The buckling problem is solved using the Ritz method. A product of the Ritz method is a complicated integral equation which is solved using numerical integration. To date, the aforementioned method of determining the prebuckling stresses was verified against infinite plate theory and finite elements. Recent efforts have focused on reducing run time, improving accuracy and adding more boundary condition choices.

Klang, Eric C.↗

Shear buckling of specially orthotropic plates with centrally located cutouts

There is significant industry demand for a method of analyzing the shear buckling of composite plates with cutouts. This method should be able to easily accommodate frequent changes in model design; inflexibility is the major drawback to current finite element methods. The approach taken here is broken into two problems, prebuckling and buckling. To solve for the prebuckling stresses, complex variable equations are used in conjunction with boundary collocation. The least squares approach is utilized to improve the accuracy of the results. The buckling problem is solved using the Ritz method. A product of the Ritz method is a complicated integral equation which is solved using numerical integration. To date, the aforementioned method of determining the prebuckling stresses was verified against infinite plate theory and finite elements.

Owen, Vicki L.↗