SEARCH · Search NASA
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 review of various methods of calculating turbulent chemically reacting flow such as the Green Function, Navier-Stokes equation, and others is presented. Nonequilibrium degrees of freedom were employed to study the mixing behavior of a multiscale turbulence field. Classical and modern theories are discussed.
In the area of computational acoustics, procedures which accurately predict the far-field sound radiation are much sought after. A systematic development of such procedures are found in a sequence of papers by Atassi. The method presented here is an alternate approach to predicting far field sound based on simple layer potential theoretic methods. The main advantages of this method are: it requires only a simple free space Green's function, it can accommodate arbitrary shapes of Kirchoff surfaces, and is readily extendable to three-dimensional problems. Moreover, the procedure presented here, though tested for unsteady lifting airfoil problems, can easily be adapted to other areas of interest, such as jet noise radiation problems. Results are presented for lifting airfoil problems and comparisons are made with the results reported by Atassi. Direct comparisons are also made for the flat plate case.
This paper presents the results and methods of computing the high-frequency radiation emitted by freely falling particles moving in circular geodesic orbits in a spherically symmetric gravitational field. The high-frequency radiation, to which the methods of this paper apply, is the principal part of radiated energy only in the case of a particle moving in a highly relativistic, and therefore unstable, circular geodesic. The geodesic synchrotron radiation emitted in this case shows excitation of high-frequency harmonics and a narrow angular distribution. A Green's-function solution of the scalar wave equation is obtained using WKB methods. For application to relativistic circular orbits, a parabolic WKB approximation is required and yields solutions in terms of parabolic cylinder functions.
We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.
We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.
The linear acoustic propagation in rigid planar sharp bends and bifurcation ducts is analyzed using a Green's-function integral technique often known as the surface element method. The acoustic characteristics of the sharp bend differ substantially from those of a circular bend with identical turning angle, width, and centerline length. The acoustic loading resulting from a duct bifurcation is highly two dimensional beyond the first cutoff frequency.
A full-wave space-domain integral equation analysis of aperture coupled shielded microstrip lines is presented based on the Equivalence Principle. The formulation includes the capability to model multilayered substrates through the derivation of the associated dyadic Green's functions which represent the layers through impedance boundary conditions. The method of moments is used to solve for the line currents and slot voltage with even and odd mode excitations which are then interpreted through transmission line analysis to determine the two-port scattering parameters. A parametric study together with experimental data is presented which demonstrates the behavior of the coupler and the accuracy of the technique.
Ffowcs Williams and Hawkings (1969) based their derivation of the fundamental equation of the sound generated by arbitrarily moving aerodynamic surfaces on the study of mass and momentum balance of a control volume imbedding a mathematical surface(s) exactly corresponding to real surface(s). These investigators also sketched an alternative method, employing generalized functions, for its derivation. This latter method, which was later developed by Farassat (1975), is purely mathematical and formal. Goldstein (1976) used the free-space Green function to produce an implicit derivation of the Ffowcs Williams and Hawkings equation. In the study presented here, Lowson's (1965) concept of moving point singularities is generalized to moving surface singularities, and a new derivation is given of the fundamental equation. The derivation is based on topological considerations of the underlying space, the fluid medium, and the integral properties of the Dirac delta function.
Stochastic radiation track structure codes are of great interest for space radiation studies and hadron therapy in medicine. These codes are used for a many purposes, notably for microdosimetry and DNA damage studies. In the last two decades, they were also used with the Independent Reaction Times (IRT) method in the simulation of chemical reactions, to calculate the yield of various radiolytic species produced during the radiolysis of water and in chemical dosimeters. Recently, we have developed a Green's function based code to simulate reversible chemical reactions with an intermediate state, which yielded results in excellent agreement with those obtained by using the IRT method. This code was also used to simulate and the interaction of particles with membrane receptors. We are in the process of including this program for use with the Monte‐Carlo track structure code Relativistic Ion Tracks (RITRACKS). This recent addition should greatly expand the capabilities of RITRACKS, notably to simulate DNA damage by both the direct and indirect effect.
Jet flows interacting with nearby surfaces exhibit a complex behavior in which acoustic and aerodynamic characteristics are altered. The physical understanding and prediction of these characteristics are essential to designing future low noise aircraft. A new approach is created for predicting scattered jet mixing noise that utilizes an acoustic analogy and steady Reynolds-averaged Navier-Stokes solutions. A tailored Green's function accounts for the propagation of mixing noise about the airframe and is calculated numerically using a newly developed ray tracing method. The steady aerodynamic statistics, associated unsteady sound source, and acoustic intensity are examined as jet conditions are varied about a large flat plate. A non-dimensional number is proposed to estimate the effect of the aerodynamic noise source relative to jet operating condition and airframe position.The steady Reynolds-averaged Navier-Stokes solutions, acoustic analogy, tailored Green's function, non-dimensional number, and predicted noise are validated with a wide variety of measurements. The combination of the developed theory, ray tracing method, and careful implementation in a stand-alone computer program result in an approach that is more first principles oriented than alternatives, computationally efficient, and captures the relevant physics of fluid-structure interaction.
In this paper, an efficient numerical solution for the scattering problem of inhomogeneous dielectric rough surfaces is presented. The inhomogeneous dielectric random surface represents a bare soil surface and is considered to be comprised of a large number of randomly positioned dielectric humps of different sizes, shapes, and dielectric constants above an impedance surface. Clods with nonuniform moisture content and rocks are modeled by inhomogeneous dielectric humps and the underlying smooth wet soil surface is modeled by an impedance surface. In this technique, an efficient numerical solution for the constituent dielectric humps over an impedance surface is obtained using Green's function derived by the exact image theory in conjunction with the method of moments. The scattered field from a sample of the rough surface is obtained by summing the scattered fields from all the individual humps of the surface coherently ignoring the effect of multiple scattering between the humps. The statistical behavior of the scattering coefficient sigma(sup 0) is obtained from the calculation of scattered fields of many different realizations of the surface. Numerical results are presented for several different roughnesses and dielectric constants of the random surfaces. The numerical technique is verified by comparing the numerical solution with the solution based on the small perturbation method and the physical optics model for homogeneous rough surfaces. This technique can be used to study the behavior of scattering coefficient and phase difference statistics of rough soil surfaces for which no analytical solution exists.
Jet flows interacting with nearby surfaces exhibit a complex behavior in which acoustic and aerodynamic characteristics are altered. The physical understanding and prediction of these characteristics are essential to designing future low noise aircraft. A new approach is created for predicting scattered jet mixing noise that utilizes an acoustic analogy and steady Reynolds-averaged Navier-Stokes solutions. A tailored Green's function accounts for the propagation of mixing noise about the air-frame and is calculated numerically using a newly developed ray tracing method. The steady aerodynamic statistics, associated unsteady sound source, and acoustic intensity are examined as jet conditions are varied about a large at plate. A non-dimensional number is proposed to estimate the effect of the aerodynamic noise source relative to jet operating condition and airframe position. The steady Reynolds-averaged Navier-Stokes solutions, acoustic analogy, tailored Green's function, non- dimensional number, and predicted noise are validated with a wide variety of measurements. The combination of the developed theory, ray tracing method, and careful implementation in a stand-alone computer program result in an approach that is more first principles oriented than alternatives, computationally efficient, and captures the relevant physics of fluid-structure interaction.
Higher order divergence- and curl-conforming bases have been shown to provide significant benefits, in both convergence rate and accuracy, in the 2D hybrid finite element/boundary element formulation (P. Fink and D. Wilton, National Radio Science Meeting, Boulder, CO, Jan. 2000). A critical issue in achieving the potential for accuracy of the approach is the accurate evaluation of all matrix elements. These involve products of high order polynomials and, in some instances, singular Green's functions. In the 2D formulation, the use of a generalized Gaussian quadrature method was found to greatly facilitate the computation and to improve the accuracy of the boundary integral equation self-terms. In this paper, a 3D, hybrid electric field formulation employing higher order bases and higher order elements is presented. The improvements in convergence rate and accuracy, compared to those resulting from lower order modeling, are established. Techniques developed to facilitate the computation of the boundary integral self-terms are also shown to improve the accuracy of these terms. Finally, simple preconditioning techniques are used in conjunction with iterative solution procedures to solve the resulting linear system efficiently. In order to handle the boundary integral singularities in the 3D formulation, the parent element- either a triangle or rectangle-is subdivided into a set of sub-triangles with a common vertex at the singularity. The contribution to the integral from each of the sub-triangles is computed using the Duffy transformation to remove the singularity. This method is shown to greatly facilitate t'pe self-term computation when the bases are of higher order. In addition, the sub-triangles can be further divided to achieve near arbitrary accuracy in the self-term computation. An efficient method for subdividing the parent element is presented. The accuracy obtained using higher order bases is compared to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained using higher order bases are compared to those obtained with lower order bases for selected sample
A method is presented to analyze a circular microstrip disk antenna loaded with a narrow, inclined slot at its center and fed by a coaxial probe. Making use of the exact dyadic Green's functions for a grounded dielectric slab due to electric and magnetic currents, Richmond's reaction integral equation in conjunction with the method of moments (Galerkin), a generalized h-parameter equivalent network for the loaded patch is obtained. An expression for the input impedance of a slot loaded microstrip patch is obtained from the equivalent parameters. The theoretical results are compared with the experiment.
Compositions analogous to lunar green, organge, and brown glasses were synthesized under consistent conditions, then quenched into a variety of different media when the samples were removed from the furnace. Iron valence and coordination are a direct function of quench media used, spanning the range from brine/ice (most effective quench), water, butyl phthalate, silicone oil, liquid nitrogen, highly reducing CO-CO2 gas, to air (least efficient quench). In the green and brown glasses, Fe(3+) in four-fold and six-fold coordination is observed in the slowest-quenched samples; Fe(2+) coordination varies directly with quench efficiency. Less pronounced changes were observed in the Ti-rich orange glass. Therefore the remote-sensed spectrum of a glass-bearing regolith on the Moon may be influenced by the process by which the glass cooled, and extreme caution must be used when comparing spectra of synthetic glass analogs with real lunar glasses.
Compositions analogous to lunar green, orange, and brown glasses were synthesized under consistent conditions, then quenched into a variety of different media when the samples were removed from the furnace. Iron valence and coordination are a direct function of quench media used, spanning the range from brine/ice (most effective quench), water, butyl phthalate, silicone oil, liquid nitrogen, highly reducing CO-CO2 gas, to air (least efficient quench). In the green and brown glasses, Fe(3+) in four-fold and six-fold coordination is observed in the slowest-quenched samples; Fe(2+) coordination varies directly with quench efficiency. Less pronounced changes were observed in the Ti-rich orange glass. Therefore the remote-sensed spectrum of a glass-bearing regolith on the moon may be influenced by the process by which the glass cooled, and extreme caution must be used when comparing spectra of synthetic glass analogs with real lunar glasses.
A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport. The code is established to operate on the Langley Research Center nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code is highly efficient and compares well with the perturbation approximations.
The properties of finite-amplitude thermal convection for a Boussinesq fluid contained in a spherical shell are investigated. All nonlinear terms are retained in the equations, and both axisymmetric and nonaxisymmetric solutions are studied. The velocity is expanded in terms of poloidal and toroidal vectors. Spherical surface harmonics resolve the horizontal structure of the flow, but finite differences are used in the vertical. With a few modifications, the transform method developed by Orszag (1970) is used to calculate the nonlinear terms, while Green's function techniques are applied to the poloidal equation and diffusion terms.