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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

SYNTHESIS OF NOVEL ALL-DIELECTRIC GRATING FILTERS USING GENETIC ALGORITHMS

We are concerned with the design of inhomogeneous, all dielectric (lossless) periodic structures which act as filters. Dielectric filters made as stacks of inhomogeneous gratings and layers of materials are being used in optical technology, but are not common at microwave frequencies. The problem is then finding the periodic cell's geometric configuration and permittivity values which correspond to a specified reflectivity/transmittivity response as a function of frequency/illumination angle. This type of design can be thought of as an inverse-source problem, since it entails finding a distribution of sources which produce fields (or quantities derived from them) of given characteristics. Electromagnetic sources (electric and magnetic current densities) in a volume are related to the outside fields by a well known linear integral equation. Additionally, the sources are related to the fields inside the volume by a constitutive equation, involving the material properties. Then, the relationship linking the fields outside the source region to those inside is non-linear, in terms of material properties such as permittivity, permeability and conductivity. The solution of the non-linear inverse problem is cast here as a combination of two linear steps, by explicitly introducing the electromagnetic sources in the computational volume as a set of unknowns in addition to the material unknowns. This allows to solve for material parameters and related electric fields in the source volume which are consistent with Maxwell's equations. Solutions are obtained iteratively by decoupling the two steps. First, we invert for the permittivity only in the minimization of a cost function and second, given the materials, we find the corresponding electric fields through direct solution of the integral equation in the source volume. The sources thus computed are used to generate the far fields and the synthesized triter response. The cost function is obtained by calculating the deviation between the synthesized value of reflectivity/transmittivity and the desired one. Solution geometries for the periodic cell are sought as gratings (ensembles of columns of different heights and widths), or combinations of homogeneous layers of different dielectric materials and gratings. Hence the explicit unknowns of the inversion step are the material permittivities and the relative boundaries separating homogeneous parcels of the periodic cell.

Zuffada, Cinzia↗
NASA NTRS1974

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.↗
NASA NTRS1949

Response of a Rotating Propeller to Aerodynamic Excitation

The flexural vibration of a rotating propeller blade with clamped shank is analyzed with the object of presenting, in matrix form, equations for the elastic bending moments in forced vibration resulting from aerodynamic forces applied at a fixed multiple of rotational speed. Matrix equations are also derived which define the critical speeds end mode shapes for any excitation order and the relation between critical speed and blade angle. Reference is given to standard works on the numerical solution of matrix equations of the forms derived. The use of a segmented blade as an approximation to a continuous blade provides a simple means for obtaining the matrix solution from the integral equation of equilibrium, so that, in the numerical application of the method presented, the several matrix arrays of the basic physical characteristics of the propeller blade are of simple form, end their simplicity is preserved until, with the solution in sight, numerical manipulations well-known in matrix algebra yield the desired critical speeds and mode shapes frame which the vibration at any operating condition may be synthesized. A close correspondence between the familiar Stodola method and the matrix method is pointed out, indicating that any features of novelty are characteristic not of the analytical procedure but only of the abbreviation, condensation, and efficient organization of the numerical procedure made possible by the use of classical matrix theory.

Arnoldi, Walter E.↗
NASA NTRS1987

Damage growth in composite laminates with interleaves

The influence of placing interleaves between fiber reinforced plies in multilayered composite laminates is investigated. The geometry of the composite is idealized as two dimensional, isotropic, linearly elastic media made of a damaged layer bonded between two half planes and separated by thin interleaves of low extensional and shear moduli. The damage in the layer is taken in the form of a symmetric crack perpendicular to the interface and may extend up to the interface. The case of an H-shaped crack in the form of a broken layer with delamination along the interface is also analyzed. The interleaves are modeled as distributed shear and tension springs. Fourier integral transform techniques are used to develop solutions in terms of singular integral equations. An asymptotic analysis of the integral equations based on Muskhelishvili's techniques reveals logarithmically singular axial stresses in the half plane at the crack tips for the broken layer. For the H shaped crack, similar singularities are found to exist in the axial stresses at the interface crack tips in the layer and the half plane. The solution of the equations is found numerically for the stresses and displacements by using the Hadamard's concept of direct differentiation of Cauchy integrals as well as Gaussian integration techniques.

Goree, James G.↗
NASA NTRS1992

Calculation of stress intensity factors in an isotropic multicracked plate. Part 1: Theoretical development

An essential part of describing the damage state and predicting the damage growth in a multicracked plate is the accurate calculation of stress intensity factors (SIF's). Here, a methodology and rigorous solution formulation for SIF's of a multicracked plate, with fully interacting cracks, subjected to a far-field arbitrary stress state is presented. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. This analytical derivation and numerical solution are obtained by using intelligent application of symbolic computations and automatic FORTRAN generation capabilities (described in the second part of this paper). As a result, a symbolic/FORTRAN package, named SYMFRAC, that is capable of providing accurate SIF's at each crack tip was developed and validated.

Binienda, W. K.↗
NASA NTRS1975

Solution of elastoplastic torsion problem by boundary integral method

The boundary integral method was applied to the elastoplastic analysis of the torsion of prismatic bars, and the results are compared with those obtained by the finite difference method. Although fewer unknowns were used, very good accuracy was obtained with the boundary integral method. Both simply and multiply connected bodies can be handled with equal ease.

Mendelson, A.↗
NASA NTRS2003

Significance of Strain in Formulation in Theory of Solid Mechanics

The basic theory of solid mechanics was deemed complete circa 1860 when St. Venant provided the strain formulation or the field compatibility condition. The strain formulation was incomplete. The missing portion has been formulated and identified as the boundary compatibility condition (BCC). The BCC, derived through a variational formulation, has been verified through integral theorem and solution of problems. The BCC, unlike the field counterpart, do not trivialize when expressed in displacements. Navier s method and the stiffness formulation have to account for the extra conditions especially at the inter-element boundaries in a finite element model. Completion of the strain formulation has led to the revival of the direct force calculation methods: the Integrated Force Method (IFM) and its dual (IFMD) for finite element analysis, and the completed Beltrami-Michell formulation (CBMF) in elasticity. The benefits from the new methods in elasticity, in finite element analysis, and in design optimization are discussed. Existing solutions and computer codes may have to be adjusted for the compliance of the new conditions. Complacency because the discipline is over a century old and computer codes have been developed for half a century can lead to stagnation of the discipline.

Patnaik, Surya N.↗
NASA NTRS2010

The Emergence of Agent-Based Technology as an Architectural Component of Serious Games

The evolution of games as an alternative to traditional simulations in the military context has been gathering momentum over the past five years, even though the exploration of their use in the serious sense has been ongoing since the mid-nineties. Much of the focus has been on the aesthetics of the visuals provided by the core game engine as well as the artistry provided by talented development teams to produce not only breathtaking artwork, but highly immersive game play. Consideration of game technology is now so much a part of the modeling and simulation landscape that it is becoming difficult to distinguish traditional simulation solutions from game-based approaches. But games have yet to provide the much needed interactive free play that has been the domain of semi-autonomous forces (SAF). The component-based middleware architecture that game engines provide promises a great deal in terms of options for the integration of agent solutions to support the development of non-player characters that engage the human player without the deterministic nature of scripted behaviors. However, there are a number of hard-learned lessons on the modeling and simulation side of the equation that game developers have yet to learn, such as: correlation of heterogeneous systems, scalability of both terrain and numbers of non-player entities, and the bi-directional nature of simulation to game interaction provided by Distributed Interactive Simulation (DIS) and High Level Architecture (HLA).

Phillips, Mark↗
NASA NTRS2014

The Importance of Data Visualization: Incorporating Storytelling into the Scientific Presentation

From its inception in 2000, one of the primary tasks of the Biomedical Data Reduction Analysis (BDRA) group has been translation of large amounts of data into information that is relevant to the audience receiving it. BDRA helps translate data into an integrated model that supports both operational and research activities. This data integrated model and subsequent visual data presentations have contributed to BDRA's success in delivering the message (i.e., the story) that its customers have needed to communicate. This success has led to additional collaborations among groups that had previously not felt they had much in common until they worked together to develop solutions in an integrated fashion. As more emphasis is placed on working with "big data" and on showing how NASA's efforts contribute to the greater good of the American people and of the world, it becomes imperative to visualize the story of our data to communicate the greater message we need to share. METHODS To create and expand its data integrated model, BDRA has incorporated data from many different collaborating partner labs and other sources. Data are compiled from the repositories of the Lifetime Surveillance of Astronaut Health and the Life Sciences Data Archive, and from the individual laboratories at Johnson Space Center that support collection of data from medical testing, environmental monitoring, and countermeasures, as designated in the Medical Requirements Integration Documents. Ongoing communication with the participating collaborators is maintained to ensure that the message and story of the data are retained as data are translated into information and visual data presentations are delivered in different venues and to different audiences. RESULTS We will describe the importance of storytelling through an integrated model and of subsequent data visualizations in today's scientific presentations and discuss the collaborative methods used. We will illustrate the discussion with examples of graphs from BDRA's past work supporting operations and/or research efforts.

Babiak-Vazquez, A.↗
NASA NTRS1995

Driving Force Analysis in an Infinite Anisotropic Plate with Multiple Crack Interactions

The methodology and a rigorous solution formulation are presented for stress intensity factors (SIF's, k) and total strain energy release rates (SERR, G(sub T)) of a multicracked plate, that has fully interacting cracks and is subjected to a far-field arbitrary stress state. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. Parametric studies are conducted for two, three, and four crack problems. The sensitivity and characteristics of the model is demonstrated.

Binienda, Wieslaw K.↗
NASA NTRS1986

Viscous wing theory development. Volume 1: Analysis, method and results

Viscous transonic flows at large Reynolds numbers over 3-D wings were analyzed using a zonal viscid-inviscid interaction approach. A new numerical AFZ scheme was developed in conjunction with the finite volume formulation for the solution of the inviscid full-potential equation. A special far-field asymptotic boundary condition was developed and a second-order artificial viscosity included for an improved inviscid solution methodology. The integral method was used for the laminar/turbulent boundary layer and 3-D viscous wake calculation. The interaction calculation included the coupling conditions of the source flux due to the wing surface boundary layer, the flux jump due to the viscous wake, and the wake curvature effect. A method was also devised incorporating the 2-D trailing edge strong interaction solution for the normal pressure correction near the trailing edge region. A fully automated computer program was developed to perform the proposed method with one scalar version to be used on an IBM-3081 and two vectorized versions on Cray-1 and Cyber-205 computers.

Chow, R. R.↗
NASA NTRS1986

An exterior Poisson solver using fast direct methods and boundary integral equations with applications to nonlinear potential flow

A general method is developed combining fast direct methods and boundary integral equation methods to solve Poisson's equation on irregular exterior regions. The method requires O(N log N) operations where N is the number of grid points. Error estimates are given that hold for regions with corners and other boundary irregularities. Computational results are given in the context of computational aerodynamics for a two-dimensional lifting airfoil. Solutions of boundary integral equations for lifting and nonlifting aerodynamic configurations using preconditioned conjugate gradient are examined for varying degrees of thinness.

Young, D. P.↗
NASA NTRS1987

On the theory of the synthesis of single and dual offset shaped reflector antennas

Two reflectors shaped so as to transforn a geometrical optics (GO) feed power pattern into a uniform amplitude and phase aperture distribution are indicated schematically. Single- and dual-shaped geometrical optics problems are solved via exact solutions obtained by progressive integration. It is concluded that an infinite set of such solutions exists.

Galindo-Israel, Victor↗
NASA NTRS1995

Large-scale structure of the solar corona and inner heliosphere

The large-scale structure of the solar corona influences solar activity particularly coronal mass ejections (CMEs). The helmet streamers that are observed to dominate the structure of the inner corona are formed by the interaction of the solar wind with coronal magnetic fields. We have simulated this interaction in three dimensions using the magnetohydrodynamic (MHD) equations. In order to create a realistic model, we use the magnetic field that is observed at the Sun's surface (deduced from daily Wilcox Solar Observatory magnetograms) as input, in combination with specified density and temperature profiles at the surface. A self-consistent 3D solar-wind solution is developed by integrating the MHD equations in time to steady state. Such solutions can reproduce the observed structures that are seen in coronagraph images and eclipse photographs of the corona. This model allows us to accurately determine the position of the heliospheric current sheet. We will compare the results obtained from our model with Ulysses observations during the period May-June 1993, and with an eclipse photograph of the corona on November 3, 1994.

Mikic, Z.↗
NASA NTRS2000

Radiation Heat Transfer Between Diffuse-Gray Surfaces Using Higher Order Finite Elements

This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.

Gould, Dana C.↗
NASA NTRS2012

Quasi 1D Modeling of Mixed Compression Supersonic Inlets

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of the 2-dimensional bifurcated mixed compression supersonic inlet is being developed. The model utilizes computational fluid dynamics for both the supersonic and subsonic diffusers. The oblique shocks are modeled utilizing compressible flow equations. This model also implements variable geometry required to control the normal shock position. The model is flexible and can also be utilized to simulate other mixed compression supersonic inlet designs. The model was validated both in time and in the frequency domain against the legacy LArge Perturbation INlet code, which has been previously verified using test data. This legacy code written in FORTRAN is quite extensive and complex in terms of the amount of software and number of subroutines. Further, the legacy code is not suitable for closed loop feedback controls design, and the simulation environment is not amenable to systems integration. Therefore, a solution is to develop an innovative, more simplified, mixed compression inlet model with the same steady state and dynamic performance as the legacy code that also can be used for controls design. The new nonlinear dynamic model is implemented in MATLAB Simulink. This environment allows easier development of linear models for controls design for shock positioning. The new model is also well suited for integration with a propulsion system model to study inlet/propulsion system performance, and integration with an aero-servo-elastic system model to study integrated vehicle ride quality, vehicle stability, and efficiency.

Kopasakis, George↗
NASA NTRS1975

Second-order small-disturbance solutions for hypersonic flow over power-law bodies

Similarity solutions were found which give the adiabatic flow of an ideal gas about two-dimensional and axisymmetric power-law bodies at infinite Mach number to second order in the body slenderness parameter. The flow variables were expressed as a sum of zero-order and perturbation similarity functions for which the axial variations in the flow equations separated out. The resulting similarity equations were integrated numerically. The solutions, which are universal functions, are presented in graphic and tabular form. To avoid a singularity in the calculations, the results are limited to body power-law exponents greater than about 0.85 for the two-dimensional case and 0.75 for the axisymmetric case. Because of the entropy layer induced by the nose bluntness (for power-law bodies other than cones and wedges), only the pressure function is valid at the body surface. The similarity results give excellent agreement with the exact solutions for inviscid flow over wedges and cones having half-angles up to about 20 deg. They give good agreement with experimental shock-wave shapes and surface-pressure distributions for 3/4-power axisymmetric bodies, considering that Mach number and boundary-layer displacement effects are not included in the theory.

Townsend, J. C.↗