Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Relativity”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Two dimensional sound propagation over a ridge

The sound propagation over an arbitrarily-shaped ridge situated on an infinite plane is treated using the Boundary Element Method, for which the Helmholtz integral equation in two-dimensional half space is formulated. A previously developed computational scheme for the acoustic problem at wave length shorther than the characteristic dimension of the radiator or scatterer is applied to the propagation of relatively low frequency sound over large terrain features. The numerical results using the Boundary Element Method are compared to those using algorithms for the Matched Asymptotic Expansion technique.

Eversman, W.↗
NASA NTRS1992

CFD lends the government a hand

The present survey of important and novel CFD applications being developed and implemented by U.S. Government contractors gives attention to naval vessel flow-modeling, Army ballistic and rotary wing aerodynamics, and NASA hypersonic vehicle related applications of CFD. CFD-generated knowledge of numerical algorithms, fluid motion, and supercomputer use is being incorporated into such additional areas as computational electromagnetics and acoustics. Attention is presently given to CFD methods' development status in such fields as submarine boundary layers, hypersonic kinetic energy projectile shock structures, helicopter main rotor tip flows, and National Aerospace Plane aerothermodynamics.

Lekoudis, Spiro↗
NASA NTRS1993

A hybrid structured-unstructured grid method for unsteady turbomachinery flow computations

A hybrid grid technique for the solution of 2D, unsteady flows is developed. This technique is capable of handling complex, multiple component geometries in relative motion, such as those encountered in turbomachinery. The numerical approach utilizes a mixed structured-unstructured zonal grid topology along with modeling equations and solution methods that are most appropriate in the individual domains, therefore combining the advantages of both structured and unstructured grid techniques.

Mathur, Sanjay R.↗
NASA NTRS1993

Turbulent structure of a wingtip vortex in the near field

The turbulent rollup of a vortex generated by a rectangular wing has been investigated. Extensive mean and turbulence measurements of the flowfield on a wingtip and in the near field have been completed. Velocity fluctuation measurements show that the near-field core is not laminar. A large axial velocity excess was found to exist in the core of the vortex. A momentum balance in the near-field of the wingtip showed that the magnitude of the core Reynolds-stress gradient terms are the same order as the largest terms in the governing equations. Navier-Stokes computations of the identical configuration, including wind tunnel walls and using measured inflow and outflow boundary conditions, reproduced many of the features of the experiment. Inherent limitations of the Baldwin-Barth turbulence model combined with limited grid resolution caused the computed vortex core to be more diffuse than desired. The momentum balance also demonstrated that the level of numerically generated false diffusion in the vortex core is relatively high.

Zilliac, Gregory G.↗
NASA NTRS1995

Viscous-flow analysis of a subsonic transport aircraft high-lift system and correlation with flight data

High-lift system aerodynamics has been gaining attention in recent years. In an effort to improve aircraft performance, comprehensive studies of multi-element airfoil systems are being undertaken in wind-tunnel and flight experiments. Recent developments in Computational Fluid Dynamics (CFD) offer a relatively inexpensive alternative for studying complex viscous flows by numerically solving the Navier-Stokes (N-S) equations. Current limitations in computer resources restrict practical high-lift N-S computations to two dimensions, but CFD predictions can yield tremendous insight into flow structure, interactions between airfoil elements, and effects of changes in airfoil geometry or free-stream conditions. These codes are very accurate when compared to strictly 2D data provided by wind-tunnel testing, as will be shown here. Yet, additional challenges must be faced in the analysis of a production aircraft wing section, such as that of the NASA Langley Transport Systems Research Vehicle (TSRV). A primary issue is the sweep theory used to correlate 2D predictions with 3D flight results, accounting for sweep, taper, and finite wing effects. Other computational issues addressed here include the effects of surface roughness of the geometry, cove shape modeling, grid topology, and transition specification. The sensitivity of the flow to changing free-stream conditions is investigated. In addition, the effects of Gurney flaps on the aerodynamic characteristics of the airfoil system are predicted.

Potter, R. C.↗
NASA NTRS1998

Bi-Level Integrated System Synthesis (BLISS)

BLISS is a method for optimization of engineering systems by decomposition. It separates the system level optimization, having a relatively small number of design variables, from the potentially numerous subsystem optimizations that may each have a large number of local design variables. The subsystem optimizations are autonomous and may be conducted concurrently. Subsystem and system optimizations alternate, linked by sensitivity data, producing a design improvement in each iteration. Starting from a best guess initial design, the method improves that design in iterative cycles, each cycle comprised of two steps. In step one, the system level variables are frozen and the improvement is achieved by separate, concurrent, and autonomous optimizations in the local variable subdomains. In step two, further improvement is sought in the space of the system level variables. Optimum sensitivity data link the second step to the first. The method prototype was implemented using MATLAB and iSIGHT programming software and tested on a simplified, conceptual level supersonic business jet design, and a detailed design of an electronic device. Satisfactory convergence and favorable agreement with the benchmark results were observed. Modularity of the method is intended to fit the human organization and map well on the computing technology of concurrent processing.

Sobieszczanski-Sobieski, Jaroslaw↗
NASA NTRS1999

Potential Flow Interactions With Directional Solidification

The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.

Buddhavarapu, Sudhir S.↗
NASA NTRS2001

Effect of Under-Resolved Grids on High Order Methods

There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.

Yee, H. C.↗
NASA NTRS2003

Dynamic Turbulence Modelling in Large-eddy Simulations of the Cloud-topped Atmospheric Boundary Layer

The use of large eddy simulation, or LES, to study the atmospheric boundary layer dates back to the early 1970s when Deardor (1972) used a three-dimensional simulation to determine velocity and temperature scales in the convective boundary layer. In 1974 he applied LES to the problem of mixing layer entrainment (Deardor 1974) and in 1980 to the cloud-topped boundary layer (Deardor 1980b). Since that time the LES approach has been applied to atmospheric boundary layer problems by numerous authors. While LES has been shown to be relatively robust for simple cases such as a clear, convective boundary layer (Mason 1989), simulation of the cloud-topped boundary layer has proved more of a challenge. The combination of small length scales and anisotropic turbulence coupled with cloud microphysics and radiation effects places a heavy burden on the turbulence model, especially in the cloud-top region. Consequently, over the past few decades considerable effort has been devoted to developing turbulence models that are better able to parameterize these processes. Much of this work has involved taking parameterizations developed for neutral boundary layers and deriving corrections to account for buoyancy effects associated with the background stratification and local buoyancy sources due to radiative and latent heat transfer within the cloud (see Lilly 1962; Deardor 1980a; Mason 1989; MacVean & Mason 1990, for example). In this paper we hope to contribute to this effort by presenting a number of turbulence models in which the model coefficients are calculated dynamically during the simulation rather than being prescribed a priori.

Kirkpatrick, M. P.↗
NASA NTRS1998

Bi-Level Integrated System Synthesis (BLISS)

BLISS is a method for optimization of engineering systems by decomposition. It separates the system level optimization, having a relatively small number of design variables, from the potentially numerous sub-system optimizations that may each have a large number of local design variables. The subsystem optimizations are autonomous and may be conducted concurrently. Subsystem and system optimizations alternate, linked by sensitivity data, producing a design improvement in each iteration. Starting from a best guess initial design, the method improves that design in iterative cycles, each cycle comprised of two steps. In step one, the system level variables are frozen and the improvement is achieved by separate, concurrent, and autonomous optimizations in the local variable subdomains. In step two, further improvement is sought in the space of the system level variables. Optimum sensitivity data link the second step to the first. The method prototype was implemented using MATLAB and iSIGHT programming software and tested on a simplified, conceptual level supersonic business jet design, and a detailed design of an electronic device.

Sobieszczanski-Sobieski, Jaroslaw↗
NASA NTRS2013

Practical UAV Optical Sensor Bench with Minimal Adjustability

A multiple-pass optical platform eliminates essentially all optical alignment degrees of freedom, save one. A four-pass absorption spectrometer architecture is made rigid by firmly mounting dielectric-coated mirror prisms with no alignment capability to the platform. The laser diode beam is collimated by a small, custom-developed lens, which has only a rotational degree of freedom along the standard optical "z" axis. This degree is itself eliminated by adhesive after laser collimation. Only one degree of freedom is preserved by allowing the laser diode chip and mount subassembly to move relative to the collimating lens by using over-sized mounting holes. This allows full 360 deg motion of a few millimeters relative to the lens, which, due to the high numerical aperture of the lens, provides wide directional steering of the collimated laser beam.

Pilgrim, Jeffrey↗
NASA NTRS2013

Discrete Kinetic Eigenmode Spectra of Electron Plasma Oscillations in Weakly Collisional Plasma: A Numerical Study

It has been demonstrated that in the presence of weak collisions, described by the Lenard-Bernstein collision operator, the Landau-damped solutions become true eigenmodes of the system and constitute a complete set. We present numerical results from an Eulerian Vlasov code that incorporates the Lenard-Bernstein collision operator. The effect of the collisions on the numerical recursion phenomenon seen in Vlasov codes is discussed. The code is benchmarked against exact linear eigenmode solutions in the presence of weak collisions, and a spectrum of Landau-damped solutions is determined within the limits of numerical resolution. Tests of the orthogonality and the completeness relation are presented.

Kinetic Eigenmode Spectra↗
NASA NTRS2016

Retrieval of Daily Maps of PM2.5 Aerosol in the Problematic California Valleys: Bright, Speckled Reflectances, Thin AOT, but High Pollution

The San Joaquin Valley suffers from severe episodes of respirable aerosol (PM2.5) in wintertime.We provide maps of aerosol episodes using daily snapshots of PM2.5 and its changing features despite numerous difficulties inherent to sampling the region. Linear relationships relating aerosol optical thickness (AOT) to PM2.5 give an explained variance of approximately 3.The GEO-CAPE mission has as a goal the provision of relevant measures of respirable aerosol to the community,but has not formulated a science goal beyond the limited goal of retrieval of AOT, bringing the usefulness of GEO-CAPE into doubt.Our special focus was on the DISCOVER-AQ period, Jan-Feb 2013, which had many supporting measurements.Both high pollution and retrieval difficulties tend to occur in many Mediterranean agricultural regions like the San Joauin. One difficulty is the relatively bright surfaces with considerable exposed soil. NASAs MAIAC and MODIS Deep Blue retrieval techniques are shown to have considerable skill even at low aerosol optical thickness (AOT) values, as evaluated by concurrent AERONET sunphotometer measurements.More significantly, these AOT values can correspond to high daytime PM2.5 since aerosol mixed layer depth is thin and variable, 200m 600 m. The thin layers derive from typical subsidence of dry air between more stormy periods. This situation provides an advantage: water vapor column is also almost completely limited to a similar mixed layer depth, and can thus serve as a measure of aerosol dilution.Using the MAIAC Water Vapor Column:In order to make the maps below, we used the MAIAC data but subtracted partial water-vapor columns estimated from MERRA Reanalysis Data availabe from the GSFC GMAO using kriging. We did not use the mixed-layer estimates from MERRA, since such analyses were found problematic during our forecasting exercises for DISCOVER-AQ. Observations from the aircraft soundings suggested that this overlying moisture was mostly due to larger scale flows, not ML venting.However, the specific humidity at the surface and a nearly well-mixed ML was analyzed by kriging from the surface network (MesoWest, University of Utah). These were thought to be truer, uninfluenced by physical process modeling that combines with data observations. (TBD: How different are they?) Procedure: Subtract overlying partial water columns from MAIAC column water and divide this by a surface value of water vapor. (MAIAC column is expressed in cm of water, i.e., water vapor at surface conditions.This method appears to bring out useful details in the distribution of submicron particles in the very problematic Wintertime San Joaquin Valley, and allow analysis of pollution episodes throughout the valley, rather than long-term averages.

san joaquin valley↗
NASA NTRS2019

Assessment of Mixer-Ejector Nozzle with Thermal Acoustic Shield for Jet Noise Reduction

A tendency for excessive exhaust jet mixing noise from low bypass ratio turbofan engines is recognized as a key challenge in the design of commercial supersonic aircraft. In this work we investigate a unique combination of two noise mitigation methods as a novel strategy to reduce jet mixing noise. First, a thermal acoustic shield (TAS) is used to reflect high frequency acoustic waves at small angles to the jet axis; second, a mixer-ejector (ME) nozzle is used to mechanically shield noise propagating at large angles to the axis. The ME shroud also provides a convenient location for a TAS nozzle and improves TAS effectiveness by limiting the downstream extent of high frequency noise generation. In an additional benefit for a velocity-matched TAS stream, the ME allows a reduction in strength of the TAS outer shear layer which could serve as a secondary noise source. The present work provides a quantitative assessment of the ME-TAS concept, using a combination of RANS CFD simulations, acoustic analogy calculations for the farfield Green's function, and surrogate-based modeling and parameter space exploration. We first evaluate a subscale configuration, then use scaling arguments to apply subscale results to the systems-level analysis of a flight configuration; the latter configuration includes a generic low bypass ratio turbofan engine with an engine-driven electric generator for supplementary heating of the TAS stream. Additional RANS CFD calculations are performed for a notional ME-TAS geometry based on the full scale configuration, and various modeling assumptions and operational characteristics are evaluated. The ME-TAS concept is shown to provide effective shielding for high frequency jet noise, and should enable comparable noise suppression to a stand-alone ME of considerably greater length, weight and drag. In addition to investigating the integrated ME-TAS system, the present work differs from previous research into TAS and related fluidic shield concepts through the inclusion of modern numerical analysis tools and the systematic numerical examination of various design parameters.

thermal acoustic shield↗
DOE OSTI2014

Investigations of Bubbly Flows using Direct Numerical Simulations

Direct numerical simulations are used to examine the behavior of relatively high void fraction bubbly flows. The full Navier-Stokes equations are solved for both the bubbles and the ambient liquid and the flow near the moving and deforming phase boundaries is resolved fully. This yields detailed data for the evolution of the system that can be used to extract information about the basic bubble-bubble and bubble-flow interactions as well as about the averaged properties of the system. Such information are critical for the generation of averaged engineering models of industrial size systems as well as for novel applications of bubbly flows where experiments are difficult. The project has yielded detailed information about many aspects of bubbly flows. The results for flows in vertical channels are particularly significant and have yielded an essentially complete understanding of many aspects of the flows. The important insight is that the flow has a particularly simple structure, consisting of a homogeneous core in hydrostatic balance and a wall layer whose composition depends of whether the liquid is flowing upward or downward. These results hold for both laminar and turbulent flows and have allowed us to derive an analytical model to predict the flow. We have also shown that bubble deformability results in a completely different flow structure, but also giving rise to relatively simple predictive models. Preliminary studies of transient migration of a non-homogeneous bubble distribution and coalescence induced flow regime transitions in vertical channels were also done.

42 ENGINEERING↗
DOE OSTI2018

A Workshop on the History of Plutonium Processing at Hanford Washington as the Foundation for Modern-day Plutonium Process Science

Numerous national laboratories in the United States have roots related to the Manhattan Project that contributed to the end of World War II. The pioneering science and engineering that took place over the decades that followed continues to influence modern-day actinide-related research at those institutions, leading to unique strengths and cutting-edge ideas and perspectives. These same national laboratories also face a challenge: figuring out how to capture the expertise of a retiring workforce who had first-hand experience in these pioneering nuclear efforts. Recently, at PNNL a “Train-the-Trainer” approach has been piloted where mentors and mentees are paired together to develop and deliver interactive workshops on historical topics that relate to modern-day research, such as Plutonium Processing or Environmental Management of Legacy Nuclear Waste. This approach was used to develop a workshop that beings with the Historic Hanford process plutonium chemistry and engineering, describes its influence on current plutonium separation and process science, and then details the chemistry of plutonium as it progresses from a minor constituent in irradiated fuel to purified solid product. This workshop is targeted to the practitioner of applied science and technology wishing to gain a better understanding of the history and process-related chemistry of element 94 so that they may better direct the tools of their specific craft toward the present and future needs of plutonium science.

Smith, Frances N.↗
DOE OSTI2020

Diversified Therapeutic Phage Cocktails from Close Relatives of the Target Bacterium

This project tackles the antibiotic resistance crisis, developing a new method for discovering numerous efficacious bacteriophages for therapeutic cocktails against bacterial pathogens. The phage therapy approach to infectious disease, recently rekindled in U.S. medicine, requires numerous phages for each bacterial pathogen. Our approach 1) uses Sandia-unique software to identify dormant phages (prophages) integrated into bacterial chromosomes, 2) identifies prophage-laden bacteria that are close relatives of the target pathogenic strain to be killed, and 3) engineers away properties of these phages that are undesirable for therapy. We have perfected our phage-finding software, implemented our phage therapy strategy by targeting the pathogen Pseudomonas aeruginosa, and prepared new software to assist the phage engineering. We then turned toward Burkholderia pathogens, aiming to overcome the difficulty to transform these bacteria with a novel phage conjugation approach. Our work demonstrates the validity of a new approach to phage therapy for killing antibiotic resistant pathogens.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

On the equivalence of dynamically orthogonal and bi-orthogonal methods: Theory and numerical simulations

The Karhunen–Lòeve (KL) decomposition provides a low-dimensional representation for random fields as it is optimal in the mean square sense. Although for many stochastic systems of practical interest, described by stochastic partial differential equations (SPDEs), solutions possess this low-dimensional character, they also have a strongly time-dependent form and to this end a fixed-in-time basis may not describe the solution in an efficient way. Motivated by this limitation of standard KL expansion, Sapsis and Lermusiaux (2009) [26] developed the dynamically orthogonal (DO) field equations which allow for the simultaneous evolution of both the spatial basis where uncertainty ‘lives’ but also the stochastic characteristics of uncertainty. Recently, Cheng et al. (2013) [28] introduced an alternative approach, the bi-orthogonal (BO) method, which performs the exact same tasks, i.e. it evolves the spatial basis and the stochastic characteristics of uncertainty. In the current work we examine the relation of the two approaches and we prove theoretically and illustrate numerically their equivalence, in the sense that one method is an exact reformulation of the other. We show this by deriving a linear and invertible transformation matrix described by a matrix differential equation that connects the BO and the DO solutions. We also examine a pathology of the BO equations that occurs when two eigenvalues of the solution cross, resulting in an instantaneous, infinite-speed, internal rotation of the computed spatial basis. We demonstrate that despite the instantaneous duration of the singularity this has important implications on the numerical performance of the BO approach. On the other hand, it is observed that the BO is more stable in nonlinear problems involving a relatively large number of modes. Several examples, linear and nonlinear, are presented to illustrate the DO and BO methods as well as their equivalence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗