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 325 records · Page 18

Separated flow over a body of revolution

A method is developed to determine the flowfield of a body of revolution in separated flow. The technique employed is the use of the computer to integrate various solutions and solution properties of the sub-flowfields which make up the entire flowfield without resorting to a finite difference solution to the complete Navier-Stokes equations. The technique entails the use of the unsteady cross flow analogy and a new solution to the required two-dimensional unsteady separated flow problem based upon an unsteady discrete-vorticity wake. Data for the forces and moments on bodies of revolution at high angle of attack (outside the range of linear inviscid theories) such that the flow is substantially separated are produced which compare well with experimental data at low speeds. In addition, three-dimensional steady separation regions and wake vortex patterns are determined.

Marshall, F. J.↗

Liquid immersion apparatus for minute articles

Apparatus is disclosed for immersing minute integrated circuit chips in an etching solution in manufacturing integrated circuits during research and development. The apparatus includes a holder, having a handle and basket support for carrying a removable unitary basket and lid structure where fluid flow-through passages are formed, and wherein graduated openings in the handle provide for adjustably supporting the basket in a breaker at a desired level.

Belcher, J. G., Jr.↗

PANEL Code for Planar Cascades

Plane cascade flow solved using improved surface-singularity methods. Current external aerodynamic integral-equation techniques adapted for use in internal flow calculation. Inherent computational speed and flexibility of integral equation solutions make them useful for design calculations.

Mcfarland, E. R.↗

Aerothermal loads analysis for high speed flow over a quilted surface configuration

Attention is given to hypersonic laminar flow over a quilted surface configuration that simulates an array of Space Shuttle Thermal Protection System panels bowed in a spherical shape as a result of thermal gradient through the panel thickness. Pressure and heating loads to the surface are determined. The flow field over the configuration was mathematically modeled by means of time-dependent, three-dimensional conservation of mass, momentum, and energy equations. A boundary mapping technique was then used to obtain a rectangular, parallel piped computational domain, and an explicit MacCormack (1972) explicit time-split predictor corrector finite difference algorithm was used to obtain steady state solutions. Total integrated heating loads vary linearly with bowed height when this value does not exceed the local boundary layer thickness.

Olsen, G. C.↗

Aerothermal loads analysis for high speed flow over a quilted surface configuration

Attention is given to hypersonic laminar flow over a quilted surface configuration that simulates an array of Space Shuttle Thermal Protection System panels bowed in a spherical shape as a result of thermal gradients through the panel thickness. Pressure and heating loads to the surface are determined. The flow field over the configuration was mathematically modeled by means of time-dependent, three-dimensional conservation of mass, momentum, and energy equations. A boundary mapping technique was then used to obtain a rectangular, parallelepiped computational domain, and an explicit MacCormack (1972) explicit time-split predictor-corrector finite difference algorithm was used to obtain steady state solutions. Total integrated heating loads vary linearly with bowed height when this value does not exceed the local boundary layer thickness.

Olsen, G. C.↗

Improving the Accuracy of the Boundary Integral Method Based on the Helmholtz Integral

Several recent papers in the literature have been based on various forms of the Helmholtz integral to compute the radiation fields of vibrating bodies. The surface integral form is given. The symbols of P,R micron, rho,G,R,V, and S micron are acoustic pressure, source coordinate, angular frequency, fluid density, Green function, field coordinate, surface velocity and body surface respectively. A discretized form of the surface integral is also given. Solutions to the surface integral are complicated with the singularity of the Green function at R=R micron and with the uniqueness problem at interior eigen frequencies of the enclosed space. The use of the interior integral circumvents the singularity problem since the field points are chosen in the interior space of the vibrating body where a zero pressure condition exists. The interior integral form is given. The method to improve the accuracy is detailed. Examples of the method is presented for a variety of radiators.

Koopmann, G. H.↗

Aeroelastic tailoring of composite wings with external stores

The use of an integrating matrix solution procedure allows the aeroelastic stability of a large number of wing/store combinations to be studied with relative speed and accuracy. Results are presented which show the effects of tailoring the orientation of the plies of laminated composite wings when external stores are present. The stores utilized take the form of concentrated masses that are varied in terms of size and location, and combinations of stores are also studied. It is seen that there are cases for which small changes in the laminate can produce large changes in the aeroelastic stability boundary, and for such situations care should be used in tailoring a wing. For a given wing, the choice fiber angle to optimize the aeroelastic performance is quite sensitive to both the size and location of the store.

Green, J. A.↗

Electromagnetic fields backscattered from an s-shaped inlet cavity with an absorber coating on its inner walls

The EM backscatter from a two-dimensional S-shaped inlet cavity is analyzed using three different techniques, namely a hybrid combination of asymptotic high frequency and modal methods, an integral equation method, and the geometrical optics ray method, respectively. This inlet has a thin absorber coating on its perfectly conducting inner walls and the planar interior termination is made perfectly conducting. The effect of the absorber on the inner wall is treated via a perturbation scheme in the hybrid approach where it is assumed that the loss is sufficiently small for the method to be valid. The results are compared with the backscatter from a straight inlet cavity to evaluate the effect of offsetting the termination in the S-bend configuration such that it is not visible from the open end of the inlet. The envelope of the backscatter pattern for the straight inlet is always seen to peak around the forward axis due to the large return from the directly visible termination, and the pattern envelope tapers off away from the forward axis. Offsetting the termination causes the envelope of the backscatter pattern to flatten out, thereby reducing the return near the forward axis by several dB. The absorber coating reduces the pattern level of the straight inlet in directions away from the forward axis but has little effect on the peak near the axis; furthermore, the absorber coating is seen to consistently reduce the backscatter from the S-bend inlet for almost all incidence angles. The hybrid method gives excellent agreement with experimental data and with the integral equation solution, whereas, the geometrical optics ray tracing method is able to generally predict the average of the bachscatter pattern but not the pattern details.

Burkholder, R. J.↗

I-BIEM, an iterative boundary integral equation method for computer solutions of current distribution problems with complex boundaries: A new algorithm. I - Theoretical

A new algorithm for an iterative computation of solutions of Laplace's or Poisson's equations in two dimensions, using Green's second identity, is presented. This algorithm converges strongly and geometrically and can be applied to curved, irregular, or moving boundaries with nonlinear and/or discontinuous boundary conditions. It has been implemented in Pascal on a number of micro- and minicomputers and applied to several geometries. Cases with known analytic solutions have been tested. Convergence to within 0.1 percent to 0.01 percent of the theoretical values are obtained in a few minutes on a microcomputer.

Cahan, B. D.↗

Dust formation around M-type stars

IRAS LRS spectra of M Mira variable have shown variation in the appearance of the 9.7 micron silicate feature, which is correlated with the shape of light curve. The LRS spectra of 100 Mira variables have been studied using simple dust shell models containing mixtures of silicate and aluminum oxide dust grains. It has been shown that the aluminum oxide grains account for the observed broad feature around 12 microns and that the variation of the spectra can be interpreted in terms of the variation of the temperature at the inner boundary of silicate dust shell. It has been proposed that silicate mantle growth on aluminum oxide grains is a possible explanation for the results. In this report, the model spectra are calculated taking account of silicate mantle growth, and the physical parameters which may determine the appearance of the 9.7 micron feature in M Mira variables are investigated. In the model calculation it is assumed that aluminum oxide grains are already formed at the bottom of the circumstellar envelope because of their high condensation temperature. The growth of silicate mantle and the motion of gas and grains from r=r(sub 0), where the mantle growth starts, are investigated. Sticking and sputtering processes due to the relative motion of grain to the ambient gas are taken into account. The thermal velocity is assumed to be negligible to the drift velocity. Acceleration by radiation pressure is considered in the gas motion equation. The formal solution is integrated to obtain the emergent spectra. Physical conditions inside r(sub 0) are regarded as boundary conditions. Observed spectra are compared to model spectra to investigate the conditions at the bottom of circumstellar envelope. In modeling the envelope, a parameter C(sub l) is introduced to take account of the density fluctuation of the envelope phenomenologically.

Onaka, Takashi↗

Applications of the conjugate gradient FFT method in scattering and radiation including simulations with impedance boundary conditions

The theoretical and computational aspects related to the application of the Conjugate Gradient FFT (CGFFT) method in computational electromagnetics are examined. The advantages of applying the CGFFT method to a class of large scale scattering and radiation problems are outlined. The main advantages of the method stem from its iterative nature which eliminates a need to form the system matrix (thus reducing the computer memory allocation requirements) and guarantees convergence to the true solution in a finite number of steps. Results are presented for various radiators and scatterers including thin cylindrical dipole antennas, thin conductive and resistive strips and plates, as well as dielectric cylinders. Solutions of integral equations derived on the basis of generalized impedance boundary conditions (GIBC) are also examined. The boundary conditions can be used to replace the profile of a material coating by an impedance sheet or insert, thus, eliminating the need to introduce unknown polarization currents within the volume of the layer. A general full wave analysis of 2-D and 3-D rectangular grooves and cavities is presented which will also serve as a reference for future work.

Barkeshli, Kasra↗

Stretched three-dimensional plasma equilibria with field-aligned flow

A three-dimensional equilibrium theory for stretched plasma configurations, such as the earth's magnetotail, is extended to include the effects of field-aligned flow. The MHD equations for this case can be solved in a general way by reduction to a set of ordinary differential equations and an ordinary integral. The solutions represent lowest-order solutions of an asymptotic expansion of the MHD equations for small electric field and weak time dependence. Simplified equations are presented for two-dimensional equilibria and for incompressible flow. Possible magnetospheric applications include the configuration of the geotail near and beyond the termination of the closed field line region, the steady motion of a plasmoid (a plasma bubble severed from the earth) through the distant geotail, and configurations at the magnetopause, the interface between the magnetosphere and the shocked solar wind plasma. For illustration, solutions for the steady motion of a two-dimensional plasmoid through the distant magnetotail are presented.

Birn, J.↗

Wavelets and electromagnetics

Wavelets are an exciting new topic in applied mathematics and signal processing. This paper will provide a brief review of wavelets which are also known as families of functions with an emphasis on interpretation rather than rigor. We will derive an indirect use of wavelets for the solution of integral equations based techniques adapted from image processing. Examples for resistive strips will be given illustrating the effect of these techniques as well as their promise in reducing dramatically the requirement in order to solve an integral equation for large bodies. We also will present a direct implementation of wavelets to solve an integral equation. Both methods suggest future research topics and may hold promise for a variety of uses in computational electromagnetics.

Kempel, Leo C.↗

Antenna analysis using neural networks

Conventional computing schemes have long been used to analyze problems in electromagnetics (EM). The vast majority of EM applications require computationally intensive algorithms involving numerical integration and solutions to large systems of equations. The feasibility of using neural network computing algorithms for antenna analysis is investigated. The ultimate goal is to use a trained neural network algorithm to reduce the computational demands of existing reflector surface error compensation techniques. Neural networks are computational algorithms based on neurobiological systems. Neural nets consist of massively parallel interconnected nonlinear computational elements. They are often employed in pattern recognition and image processing problems. Recently, neural network analysis has been applied in the electromagnetics area for the design of frequency selective surfaces and beam forming networks. The backpropagation training algorithm was employed to simulate classical antenna array synthesis techniques. The Woodward-Lawson (W-L) and Dolph-Chebyshev (D-C) array pattern synthesis techniques were used to train the neural network. The inputs to the network were samples of the desired synthesis pattern. The outputs are the array element excitations required to synthesize the desired pattern. Once trained, the network is used to simulate the W-L or D-C techniques. Various sector patterns and cosecant-type patterns (27 total) generated using W-L synthesis were used to train the network. Desired pattern samples were then fed to the neural network. The outputs of the network were the simulated W-L excitations. A 20 element linear array was used. There were 41 input pattern samples with 40 output excitations (20 real parts, 20 imaginary). A comparison between the simulated and actual W-L techniques is shown for a triangular-shaped pattern. Dolph-Chebyshev is a different class of synthesis technique in that D-C is used for side lobe control as opposed to pattern shaping. The interesting thing about D-C synthesis is that the side lobes have the same amplitude. Five-element arrays were used. Again, 41 pattern samples were used for the input. Nine actual D-C patterns ranging from -10 dB to -30 dB side lobe levels were used to train the network. A comparison between simulated and actual D-C techniques for a pattern with -22 dB side lobe level is shown. The goal for this research was to evaluate the performance of neural network computing with antennas. Future applications will employ the backpropagation training algorithm to drastically reduce the computational complexity involved in performing EM compensation for surface errors in large space reflector antennas.

Smith, William T.↗

Quasi-optical grids with thin rectangular patch/aperture elements

Theoretical analysis is presented for an efficient and accurate performance evaluation of quasi-optical grids comprised of thin rectangular patch/aperture elements with/without a dielectric substrate/superstrate. The convergence rate of this efficient technique is improved by an order of magnitude with the approximate edge conditions incorporated in the basis functions of the integral equation solution. Also presented are the interesting applications of this efficient analytical technique to the design and performance evaluation of the coupling grids and beam splitters in the optical systems as well as thermal protection sunshields used in the communication systems of satellites and spacecrafts.

Wu, Te-Kao↗

Development of Fast Algorithms Using Recursion, Nesting and Iterations for Computational Electromagnetics

In the first phase of our work, we have concentrated on laying the foundation to develop fast algorithms, including the use of recursive structure like the recursive aggregate interaction matrix algorithm (RAIMA), the nested equivalence principle algorithm (NEPAL), the ray-propagation fast multipole algorithm (RPFMA), and the multi-level fast multipole algorithm (MLFMA). We have also investigated the use of curvilinear patches to build a basic method of moments code where these acceleration techniques can be used later. In the second phase, which is mainly reported on here, we have concentrated on implementing three-dimensional NEPAL on a massively parallel machine, the Connection Machine CM-5, and have been able to obtain some 3D scattering results. In order to understand the parallelization of codes on the Connection Machine, we have also studied the parallelization of 3D finite-difference time-domain (FDTD) code with PML material absorbing boundary condition (ABC). We found that simple algorithms like the FDTD with material ABC can be parallelized very well allowing us to solve within a minute a problem of over a million nodes. In addition, we have studied the use of the fast multipole method and the ray-propagation fast multipole algorithm to expedite matrix-vector multiplication in a conjugate-gradient solution to integral equations of scattering. We find that these methods are faster than LU decomposition for one incident angle, but are slower than LU decomposition when many incident angles are needed as in the monostatic RCS calculations.

Chew, W. C.↗

SCARLET I: Mechanization solutions for deployable concentrator optics integrated with rigid array technology

The SCARLET I (Solar Concentrator Army with Refractive Linear Element Technology) solar array wing was designed and built to demonstrate, in flight, the feasibility of integrating deployable concentrator optics within the design envelope of typical rigid array technology. Innovative mechanism designs were used throughout the array, and a full series of qualification tests were successfully performed in anticipation of a flight on the Multiple Experiment Transporter to Earth Orbit and Return (METEOR) spacecraft. Even though the Conestoga launch vehicle was unable to place the spacecraft in orbit, the program effort was successful in achieving the milestones of analytical and design development functional validation, and flight qualification, thus leading to a future flight evaluation for the SCARLET technology.

Wachholz, James J.↗

Numerical Simulation of Delamination Growth in Composite Materials

The use of decohesion elements for the simulation of delamination in composite materials is reviewed. The test methods available to measure the interfacial fracture toughness used in the formulation of decohesion elements are described initially. After a brief presentation of the virtual crack closure technique, the technique most widely used to simulate delamination growth, the formulation of interfacial decohesion elements is described. Problems related with decohesion element constitutive equations, mixed-mode crack growth, element numerical integration and solution procedures are discussed. Based on these investigations, it is concluded that the use of interfacial decohesion elements is a promising technique that avoids the need for a pre-existing crack and pre-defined crack paths, and that these elements can be used to simulate both delamination onset and growth.

Camanho, P. P.↗