On the numerical solution of Poisson's equation
Numerical integration of Poisson equation for computer simulated plasma
COMPUTER SIMULATION↗
SEARCH · Engineering Papers
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.
A theoretical and experimental investigation of the hypersonic turbulent boundary layer subject to normal and longitudinal pressure gradients and cross flow along a windward plane of symmetry.
A theoretical and experimental investigation for the hypersonic turbulent boundary layer undergoing both normal and longitudinal pressure gradients and cross flow along the plane of symmetry on the windward side of a body of revolution is presented. The purpose of this investigation is to provide a method for calculating boundary layers that occur on the centerline of symmetry of an inlet. The three-dimensional compressible integral equations are written along the symmetry plane and integrated numerically, and the effect of the cross flow on the behavior of the boundary layer is also studied. The experiments were performed at Mach number 6 and Reynolds number of 44,000,000/ft for the axisymmetric case and for the case with cross flow. The velocity profiles and the boundary layer parameters agree relatively well with theoretical values with and without cross flow effects. The presence of the cross flow alters significantly the value of the profiles. Total temperature-velocity correlations deviate from the linear Crocco distribution. Heat transfer rates agree with theoretical values and values obtained by the flat plate reference enthalpy method. Cross flow influence on the heat transfer rates is of minor relevance.
Zakkay, V.↗
On the use of co-ordinate stretching in the numeral computation of high frequency scattering
The scattering of the sound of a jet engine by an airplane fuselage is modeled by solving the axially symmetric Helmholtz equation exterior to a long thin ellipsoid. The integral equation method based on the single layer potential formulation is used. A family of coordinate systems on the body is introduced and an algorithm is presented to determine the optimal coordinate system. Numerical results verify that the optimal choice enables the solution to be computed with a grid that is coarse relative to the wavelength.
Bayliss, A.↗
Electromagnetic Interactions GenERalized (EIGER): Applications at Johnson Space Center
This slide presentation reviews the software framework for the analysis and design of complex electromagnetic systems. The system called Electromagnetic Interactions GenERalized (EIGER) is a hybrid finite element solution for wave equations and deru=ives boundary element solutions for integral equation formulations. Written in Fortran 90 Eiger uses object-oriented design methods to abstract key analysis components. EIGER has been applied to higher order modeling for analysis of antennas to assist in the design of the antennas of a mini-AERCam. EIGER might also be used to model the Ion Cyclotron Resonant heating stage for the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR).
Dobbins, Justin A.↗
CALCULATION OF COLLISION INTEGRALS IN THE MOMENT EQUATIONS
Gas dynamics - collision integral in moment equations
RAREFIED GAS DYNAMICS↗
Numerical analysis of Weyl's method for integrating boundary layer equations
A fast method for accurate numerical integration of Blasius equation is proposed. It is based on the limit interchange in Weyl's fixed point method formulated as an iterated limit process. Each inner limit represents convergence to a discrete solution. It is shown that the error in a discrete solution admits asymptotic expansion in even powers of step size. An extrapolation process is set up to operate on a sequence of discrete solutions to reach the outer limit. Finally, this method is extended to related boundary layer equations.
Najfeld, I.↗
Changing stepsize in the integration of differential equations using modified divided differences
Multistep methods for solving differential equations based on numerical integration formulas or numerical differentiation formulas (for stiff equations) require special provision for changing the stepsize. New algorithms are given which make the use of modified divided differences an attractive way to carry out the change in stepsize for such methods. Error estimation and some of the important factors in stepsize selection and the selection of integration order are also considered.
Krogh, F. T.↗
The solution of the relaxation problem for the Boltzmann equation by the integral iteration method
The Boltzmann equation is considered in terms of the problem of relaxation of some initial distribution function which depends only on velocities, to Maxwell's distribution function. The Boltzmann equation is given for the relaxation problem in which the distribution function f(t, u, v) is time dependent and is also dependent on two other variables u and v (the velocities of rigid spherical molecules). An iteration process is discussed in which the velocity space u, v is subdivided into squares, the distribution function in each square being approximated by the second-order surface from the values of the distribution function at nine points. The set of all of these points forms a network of u, v values at the nodes of which the distribution function can be found.
Limar, Y. F.↗
Integration of control equations and the problem of small time constants
Numerical integration of control equations, and problem of small time constants
TIME CONSTANT↗
Treatment of body forces in boundary element design sensitivity analysis
The inclusion of body forces has received a good deal of attention in boundary element research. The consideration of such forces is essential in the desgin of high performance components such as fan and turbine disks in a gas turbine engine. Due to their critical performance requirements, optimal shapes are often desired for these components. The boundary element method (BEM) offers the possibility of being an efficient method for such iterative analysis as shape optimization. The implicit-differentiation of the boundary integral equations is performed to obtain the sensitivity equations. The body forces are accounted for by either the particular integrals for uniform body forces or by a surface integration for non-uniform body forces. The corresponding sensitivity equations for both these cases are presented. The validity of present formulations is established through a close agreement with exact analytical results.
Saigal, Sunil↗
Integration of the three dimensional heat conduction equation
Numerical integration of steady state heat conduction equations representing diurnal variations in upper atmosphere
Harris, I.↗
A finite element formulation for supersonic flows around complex configurations
The problem of small perturbation potential supersonic flow around complex configurations is considered. This problem requires the solution of an integral equation relating the values of the potential on the surface of the body to the values of the normal derivative, which is known from the small perturbation boundary conditions. The surface of the body is divided into small (hyperboloidal quadrilateral) surface elements which are described in terms of the Cartesian components of the four corner points. The values of the potential (and its normal derivative) within each element are assumed to be constant and equal to its value at the centroid of the element. This yields a set of linear algebraic equations whose coefficients are given by source and doublet integrals over the surface elements. Closed form evaluations of the integrals are presented.
Morino, L.↗
Integral representation of field variables for the finite element solution of viscous flow problems
The kinematics aspects of compressible and incompressible viscous flow problems are recasted into an integral representation for the velocity vector. For the incompressible flow problem, the kinetics aspect is similarly recasted into an integral representation. The result is a system of integral equations ideally suited for the finite element method. The integral representations are shown to permit the confinement of the computation field to the region of nonnegligible vorticity and dilatation which, for the incompressible flow, is identified with the viscous region. This drastic reduction in computation field results in a superior solution speed which is further improved by the use of a flowfield segmentation technique. The new approach is further shown to remove the difficulties associated with previous methods in specifying the far field and extraneous boundary conditions.
Wu, J. C.↗
Thermal microwave emission from half-space random media
Brightness temperatures resulting from microwave thermal emission from a half-space random medium are calculated. The random medium has a nonuniform temperature profile and is characterized by correlation functions that possess both vertical and lateral variations. Radiative transfer equations are derived. They are solved with an iterative integral equation approach for small scattering albedo and with a numerical approach for general cases. New results are illustrated, discussed, and compared with various special cases. Phenomena caused by resonant scattering, which is absent under Rayleigh approximations, are revealed. Scattering in the lateral direction contributes to the decrease of the equivalent reflectivity and the increase of the brightness temperature.
Tsang, L.↗
Vibrations of incompressible liquid-plastic tank combinations using a series representation of the liquid
A series representation is presented of the dynamic behavior of incompressible non-viscous liquids contained in partially filled elastic shells of revolution. Using a series in which each term satisfies exactly Laplace's equation, the liquid portion of the problem becomes an integral equation on the wetted wall of the shell. Physical insight required in selecting the terms of the series is such that selection may be readily programed on a computer. Liquid-shell vibration interaction is investigated by modifying an existing shell of revolution vibration program with the series expansion of the liquid pressure loading on the wetted shell wall. Application of the modified program is demonstrated on a hemispherical tank and a realistic scaled model of a launch vehicle tank having a more complex geometry. Comparison of the present analysis with a NASTRAN finite element analysis for a hemispherical tank, and comparison with test results for a launch vehicle tank shows excellent agreement and rapid convergence of the series representation of the liquid to the true solution.
Housner, J. M.↗
Analysis of a unidirectional composite containing broken fibers and matrix damage
An analytical solution is developed for the determination of the stresses and displacements in a unidirectional fiber-reinforced composite containing an arbitrary number of broken fibers as well as longitudinal yielding and splitting of the matrix. The solution is developed using a materials-modeling approach which is based on a shear-lag stress transfer mechanism. The equilibrium equation in the axial direction gives a pair of integral equations which are solved numerically. Excellent agreement is shown to exist between the solution and experimental results for notched unidirectional boron/aluminum laminates without splitting. For brittle matrix composites (i.e. epoxy) equally good results are indicated for both matrix yielding and splitting. For yielding without splitting the fracture strength depends on crack length while for large splitting it is crack length independent.
George, J. G.↗
Analysis of a unidirectional composite containing broken fibers and matrix damage
An analytical solution is developed for the determination of the stresses and displacements in a unidirectional fiber-reinforced composite containing an arbitrary number of broken fibers as well as longitudinal yielding and splitting of the matrix. The solution is developed using a 'materials-modeling' approach which is based on a shear-lag stress transfer mechanism. The equilibrium equation in the axial direction gives a pair of integral equations which are solved numerically. Excellent agreement is shown to exist between the solution and experimental results for notched unidirectional boron/aluminum laminates without splitting. For brittle matrix composites (i.e., epoxy) equally good results are indicated for both matrix yielding and splitting. For yielding without splitting the fracture strength is found to depend on crack length while for large splitting it is crack length independent.
Goree, J. G.↗
Estimating location parameters in a mixture
The problem of estimating the parameters in a finite mixture is considered. The approach is based on an integral equation formulation of the form h sub t (x) = integral (limits b and a) f(x-y) g sut t (y)dy where h sub t is a smoothed version of h and g sub t is a prior function that tends to be concentrated on the translation values. A solution for g sub t that uses the method of regularization and one based on a posterior operator approach is considered. Numerical simulations are presented to bring out some of the estimation and numerical problems of these approaches.
Heydorn, R. P.↗