Search NASASearch

SEARCH · Search NASA

Results for “consistent discretization”

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 181 records · Page 10

The contribution to the gamma-ray background from a superposition of 3C 273-like jet sources

Recent theories suggest that quasar jets may produce the gamma-ray emission observed in sources such as 3C 273. One of these theories is tested to see if it is consistent with observations of the gamma-ray background. Quasar jets are taken as the source of this radiation, each with the spectrum of the nearby quasar 3C 273, and with an evolution in time based on a model for the observed gamma-rays of 3C 273. In that model, a plasma jet interacts with the ambient intergalactic gas to produce (from nuclear collisions) secondary electrons which lead to isotropic bremsstrahlung gamma-rays, and ultimately to a double radio source. It is shown that the gamma-ray background can be modeled in general as a superposition of such discrete sources of various redshifts. The results are consistent with the observed background intensity and spectral form; of course, there are uncertain and adjustable parameters.

Sadun, Alberto C.

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian

Engine structures modeling software system: Computer code. User's manual

ESMOSS is a specialized software system for the construction of geometric descriptive and discrete analytical models of engine parts, components and substructures which can be transferred to finite element analysis programs such as NASTRAN. The software architecture of ESMOSS is designed in modular form with a central executive module through which the user controls and directs the development of the analytical model. Modules consist of a geometric shape generator, a library of discretization procedures, interfacing modules to join both geometric and discrete models, a deck generator to produce input for NASTRAN and a 'recipe' processor which generates geometric models from parametric definitions. ESMOSS can be executed both in interactive and batch modes. Interactive mode is considered to be the default mode and that mode will be assumed in the discussion in this document unless stated otherwise.

Source record

A solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method.

Rosenfeld, Moshe

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation

On High-Order/Low-Order and Micro-Macro Methods for Implicit Time-Stepping of the BGK Model

In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.

BGK model

A mesh gradient technique for numerical optimization

A class of successive-improvement optimization methods in which directions of descent are defined in the state space along each trial trajectory are considered. The given problem is first decomposed into two discrete levels by imposing mesh points. Level 1 consists of running optimal subarcs between each successive pair of mesh points. For normal systems, these optimal two-point boundary value problems can be solved by following a routine prescription if the mesh spacing is sufficiently close. A spacing criterion is given. Under appropriate conditions, the criterion value depends only on the coordinates of the mesh points, and its gradient with respect to those coordinates may be defined by interpreting the adjoint variables as partial derivatives of the criterion value function. In level 2, the gradient data is used to generate improvement steps or search directions in the state space which satisfy the boundary values and constraints of the given problem.

Willis, E. A., Jr.

Photoelectric spectrophotometry of OQ 172 and OH 471

Absolute spectral energy distributions for the large redshift quasars OQ 172 and OH 471 are discussed along with similar data for two other quasars 4C05.34 and PHL 957. Assuming cosmological redshifts, OQ 172 and OH 471 are not as luminous as PHL 957. If these quasars are basically similar and if radiative processes dominate, the strength of Ly alpha and the behavior of the continuum at the Lyman limit strongly suggest that these objects consist of a central ionizing source surrounded by discrete clouds, filaments or a gaseous structure such as a disk. This gaseous matter does not cover the whole solid angle surrounding the source.

Oke, J. B.

Photoelectric spectrophotometry of OQ 172 and OH 471

Absolute spectral energy distributions for the large-redshift quasars OQ 172 and OH 471 are discussed along with similar data for two other quasars: 4C 05.34 and PHL 957. If the redshift are assumed to be cosmological, OQ 172 and OH 471 are not as luminous as PHL 957. If these quasars are basically similar and if radiative processes dominate, the strength of alpha and the behavior of the continuum at the Lyman limit strongly suggest that these objects consist of a central ionizing source surrounded by discrete clouds, filaments, or a gaseous structure such as a disk. This gaseous matter does not cover the whole solid angle surrounding the source.

Oke, J. B.

Viking magnetic properties experiment - Extended mission results

The backhoe magnets on Viking Lander (VL) 2 were successfully cleaned, followed by a test involving successive insertions of the cleaned backhoe into the surface. Rapid saturation of the magnets confirmed evidence from primary mission results that the magnetic mineral in the Martian surface is widely distributed, most probably in the form of composite particles of magnetic and nonmagnetic minerals. An image of the VL 2 backhoe taken via the X4 magnifying mirror demonstrates the fine-grained nature of the attracted magnetic material. The presence of maghemite and its occurrence as a pigment in, or a thin coating on, all mineral particles or as discrete, finely divided and widely distributed crystallites, are consistent with data from the inorganic analysis experiments and with laboratory simulations of results of the biology experiments on Mars.

Hargraves, R. B.

X-ray spectra of active galactic nuclei

Characteristics of the X-ray spectra of seven Seyfert 1 galaxies, five BL Lac objects, and three quasars are discussed. It is found that all the Seyfert 1 spectra can be fitted to a single power law with an energy index of 0.6 to 0.7, that there is no obvious correlation between spectral index and luminosity, and that variability appears to be inversely correlated with luminosity. It is suggested that the spectra of all the BL Lac sources can be synthesized from two power-law forms, viz., a very hard component with an index of approximately 0.2 and a much softer component with an index of about 2. The extent to which the measured spectrum of the diffuse X-ray background is consistent with the spectra of candidate populations of discrete sources is examined.

Holt, S. S.

Acoustic pressures emanating from a turbomachine stage

A knowledge of the acoustic energy emission of each blade row of a turbomachine is useful for estimating the overall noise level of the machine and for determining its discrete frequency noise content. Because of the close spacing between the rotor and stator of a compressor stage, the strong aerodynamic interactions between them have to be included in obtaining the resultant flow field. A three dimensional theory for determining the discrete frequency noise content of an axial compressor consisting of a rotor and a stator each with a finite number of blades are outlined. The lifting surface theory and the linearized equation of an ideal, nonsteady compressible fluid motion are used for thin blades of arbitrary cross section. The combined pressure field at a point of the fluid is constructed by linear addition of the rotor and stator solutions together with an interference factor obtained by matching them for net zero vorticity behind the stage.

Ramachandra, S. M.

Acoustic pressures emanating from a turbomachine stage

A knowledge of the acoustic energy emission of each blade row of a turbomachine is useful for estimating the overall noise level of the machine and for determining its discrete frequency noise content. Because of the close spacing between the rotor and stator of a compressor stage, the strong aerodynamic interactions between them have to be included in obtaining the resultant flow field. This paper outlines a three-dimensional theory for determining the discrete frequency noise content of an axial compressor consisting of a rotor and a stator each with a finite number of blades. The lifting surface theory and the linearized equation of an ideal, nonsteady compressible fluid motion are used for thin blades of arbitrary cross section. The combined pressure field at a point of the fluid is constructed by linear addition of the rotor and stator solutions together with an interference factor obtained by matching them for net zero vorticity behind the stage.

Ramachandra, S. M.

Optical signal processing of spatially distributed sensor data in smart structures

Smart structures which contain dense two- or three-dimensional arrays of attached or embedded sensor elements inherently require signal multiplexing and processing capabilities to permit good spatial data resolution as well as the adequately short calculation times demanded by real time active feedback actuator drive circuitry. This paper reports the implementation of an in-line optical signal processor and its application in a structural sensing system which incorporates multiple discrete optical fiber sensor elements. The signal processor consists of an array of optical fiber couplers having tailored s-parameters and arranged to allow gray code amplitude scaling of sensor inputs. The use of this signal processor in systems designed to indicate the location of distributed strain and damage in composite materials, as well as to quantitatively characterize that damage, is described. Extension of similar signal processing methods to more complicated smart materials and structures applications are discussed.

Bennett, K. D.

Relative-Error-Covariance Algorithms

Two algorithms compute error covariance of difference between optimal estimates, based on data acquired during overlapping or disjoint intervals, of state of discrete linear system. Provides quantitative measure of mutual consistency or inconsistency of estimates of states. Relative-error-covariance concept applied, to determine degree of correlation between trajectories calculated from two overlapping sets of measurements and construct real-time test of consistency of state estimates based upon recently acquired data.

Bierman, Gerald J.

An aeroelastic analysis with a generalized dynamic wake

An aeroelastic model with generalized dynamic wake is developed for application in the integration of aerodynamic, dynamic, and structural optimization of a rotor blade. The investigation is carried out with special attention to efficiency and accuracy of aeroelastic modeling. Each blade is assumed to be an elastic beam undergoing flap bending, lead-lag bending, elastic twist and axial deflections. The nonuniform blade is discretized into finite beam elements, each of which consists of twelve degrees of freedom. Such important blade design variables as pretwist, and chordwise offsets of the blade center of gravity and of the aerodynamic center from the elastic axis have been included in the analysis. Aerodynamic loads are computed from unsteady blade element theory where the rotor three-dimensional unsteady wake is modeled using a generalized dynamic wake theory. The noncirculatory loads based on unsteady thin airfoil theory are also included.

He, Cheng J.

Improving Rotor-Stator Interaction Noise Code Through Analysis of Input Parameters

There are two major sources of aircraft noise. The first is from the airframe and the second is from the engines. The focus of the acoustics branch at NASA Glenn is on the engine noise sources. There are two major sources of engine noise; fan noise and jet noise. Fan noise, produced by rotating machinery of the engine, consists of both tonal noise, which occurs at discrete frequencies, and broadband noise, which occurs across a wide range of frequencies. The focus of my assignment is on the broadband noise generated by the interaction of fan flow turbulence and the stator blades. such as the sweep and stagger angles and blade count, as well as the flow parameters such as intensity of turbulence in the flow. The tool I employed in this work is a computer program that predicts broadband noise from fans. The program assumes that the complex shape of the curved blade can be represented as a single flat plate, allowing it to use fairly simple equations that can be solved in a reasonable amount of time. While the results from such representation provided reasonable estimates of the broadband noise levels, they did not usually represent the entire spectrum accurately. My investigation found that the discrepancy between data and theory can be improved if the leading edge and the trailing edge of the blade are treated separately. Using this approach, I reduced the maximum error in noise level from a high of 30% to less than 5% for the cases investigated. Detailed results of this investigation will be discussed at my presentation. The objective of this study is to investigate the influence of geometric parameters

Unton, Timothy J.

SEDS: The Spitzer Extended Deep Survey. Survey Design, Photometry, and Deep IRAC Source Counts

The Spitzer Extended Deep Survey (SEDS) is a very deep infrared survey within five well-known extragalactic science fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the Hubble Deep Field North, and the Extended Groth Strip. SEDS covers a total area of 1.46 deg(exp 2) to a depth of 26 AB mag (3sigma) in both of the warm Infrared Array Camera (IRAC) bands at 3.6 and 4.5 micron. Because of its uniform depth of coverage in so many widely-separated fields, SEDS is subject to roughly 25% smaller errors due to cosmic variance than a single-field survey of the same size. SEDS was designed to detect and characterize galaxies from intermediate to high redshifts (z = 2-7) with a built-in means of assessing the impact of cosmic variance on the individual fields. Because the full SEDS depth was accumulated in at least three separate visits to each field, typically with six-month intervals between visits, SEDS also furnishes an opportunity to assess the infrared variability of faint objects. This paper describes the SEDS survey design, processing, and publicly-available data products. Deep IRAC counts for the more than 300,000 galaxies detected by SEDS are consistent with models based on known galaxy populations. Discrete IRAC sources contribute 5.6 +/- 1.0 and 4.4 +/- 0.8 nW / square m/sr at 3.6 and 4.5 micron to the diffuse cosmic infrared background (CIB). IRAC sources cannot contribute more than half of the total CIB flux estimated from DIRBE data. Barring an unexpected error in the DIRBE flux estimates, half the CIB flux must therefore come from a diffuse component.

Spitzer