Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 937 records · Page 52

Selective decay hypothesis at high mechanical and magnetic Reynolds numbers

Implications of certain applications of turbulence theory to two-dimensional turbulence and magnetohydrodynamic flow are discussed. It is shown that the use of the Navier-Stokes equation (NSE) for measurements of turbulent fluctuations has been effective only for three-dimensional flows. For two-dimensional flows, used for the study of large-scale motions in the atmosphere or ocean, enstrophy is cascaded to high wave numbers and dissipated at a finite rate even at infinite Re. MHD flows are numerically calculated for the two-dimensional case and analytically for the three-dimensional case, for which discrepancies in the relative rates of energy and cross helicity decay lead to a recommendation that numerical calculations for the three-dimensional case be carried out to determine the precise decayed states.

Matthaeus, W. H.↗

Current-Phase Relations of Superconducting PtSi Constriction Josephson Junctions

Here, we study the current-phase relations of superconductor–constriction–superconductor Josephson junctions made from platinum silicide thin films by fabricating dc-superconducting quantum interference devices incorporating pairs of junctions and measuring their magnetic-field-dependent electrical transport. By comparing the supercurrent interference patterns with numerical Ginzburg–Landau simulations, we extract the current-phase relations of individual junctions and quantify the degree of nonlinearity. These measurements show that while the constrictions themselves sustain substantial nonlinearity, the effect is moderated by the kinetic inductance of the device leads, which is an important consideration for practical superconducting circuit applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Research highlights: June 1990 - May 1991

Linear instability calculations at MSFC have suggested that the Geophysical Fluid Flow Cell (GFFC) should exhibit classic baroclinic instability at accessible parameter settings. Interest was in the mechanisms of transition to temporal chaos and the evolution of spatio-temporal chaos. In order to understand more about such transitions, high resolution numerical experiments for the physically simplest model of two layer baroclinic instability were conducted. This model has the advantage that the numerical code is exponentially convergent and can be efficiently run for very long times, enabling the study of chaotic attractors without the often devastating effects of low-order trunction found in many previous studies. Numerical algorithms for implementing an empirical orthogonal function (EOF) analysis of the high resolution numerical results were completed. Under conditions of rapid rotation and relatively low differential heating, convection in a spherical shell takes place as columnar banana cells wrapped around the annular gap, but with axes oriented along the axis of rotation; these were clearly evident in the GFFC experiments. The results of recent numerical simulations of columnar convection and future research plans are presented.

Source record↗

A numerical study of transient heat and mass transfer in crystal growth

A numerical analysis of transient heat and solute transport across a rectangular cavity is performed. Five nonlinear partial differential equations which govern the conservation of mass, momentum, energy and solute concentration related to crystal growth in solution, are simultaneously integrated by a numerical method based on the SIMPLE algorithm. Numerical results showed that the flow, temperature and solute fields are dependent on thermal and solutal Grashoff number, Prandtl number, Schmidt number and aspect ratio. The average Nusselt and Sherwood numbers evaluated at the center of the cavity decrease markedly when the solutal buoyancy force acts in the opposite direction to the thermal buoyancy force. When the solutal and thermal buoyancy forces act in the same direction, however, Sherwood number increases significantly and yet Nusselt number decreases. Overall effects of convection on the crystal growth are seen to be an enhancement of growth rate as expected but with highly nonuniform spatial growth variations.

Han, Samuel Bang-Moo↗

Multilateration - A nondegenerate method of obtaining station coordinates and satellite ephemerides

A technique for the determination of three-dimensional station coordinates and satellite ephemerides is developed which is based on the principle of multilateration. The method makes use of a system of six ground stations which simultaneously measure the slant range between each station and one or two satellites. It is demonstrated that a minimum of six stations is required in order to yield a system which will be free of mathematical degeneracies. It will be seen that the method of multilateration is not dependent upon the position of the satellite or any other dynamical considerations in the equations used to determine the relative station coordinates. In fact, the satellite coordinates are obtained as a direct by-product of the method. Numerical results are presented which indicate that the method of multilateration can determine the relative three-dimensional station coordinates with an accuracy that is limited only by the hardware measurement system. If a highly accurate laser ranging system is used, then accuracies in the 1-cm range can be expected.

Ong, K. M.↗

A Lagrangian dual method for two-stage robust optimization with binary uncertainties

This report presents a new exact method to calculate worst-case parameter realizations in two-stage robust optimization problems with categorical or binary-valued uncertain data. Traditional exact algorithms for these problems, notably Benders decomposition and column-and-constraint generation, compute worst-case parameter realizations by solving mixed-integer bilinear optimization subproblems. However, their numerical solution can be computationally expensive not only due to their resulting large size after reformulating the bilinear terms, but also because decision-independent bounds on their variables are typically unknown. We propose an alternative Lagrangian dual method that circumvents these difficulties and is readily integrated in either algorithm. We specialize the method to problems where the binary parameters switch on or off constraints as these are commonly encountered in applications, and discuss extensions to problems that lack relatively complete recourse and to those with integer recourse. Numerical experiments provide evidence of significant computational improvements over existing methods.

42 ENGINEERING↗

Spatial Intermittency of Particle Distribution in Relativistic Plasma Turbulence

Abstract Relativistic magnetically dominated turbulence is an efficient engine for particle acceleration in a collisionless plasma. Ultrarelativistic particles accelerated by interactions with turbulent fluctuations form nonthermal power-law distribution functions in the momentum (or energy) space, f ( γ ) d γ ∝ γ − α d γ , where γ is the Lorenz factor. We argue that in addition to exhibiting non-Gaussian distributions over energies, particles energized by relativistic turbulence also become highly intermittent in space. Based on particle-in-cell numerical simulations and phenomenological modeling, we propose that the bulk plasma density has lognormal statistics, while the density of the accelerated particles, n , has a power-law distribution function, P ( n ) dn ∝ n − β dn . We argue that the scaling exponents are related as β ≈ α + 1, which is broadly consistent with numerical simulations. Non-space-filling, intermittent distributions of plasma density and energy fluctuations may have implications for plasma heating and for radiation produced by relativistic turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

A Numerical Investigation of Turbine Noise Source Hierarchy and Its Acoustic Transmission Characteristics

Understanding the relative importance of the various turbine noise generation mechanisms and the characteristics of the turbine acoustic transmission loss are essential ingredients in developing robust reduced-order models for predicting the turbine noise signature. A computationally based investigation has been undertaken to help guide the development of a turbine noise prediction capability that does not rely on empiricism. The investigation relies on highly detailed numerical simulations of the unsteady flowfield inside a modern high-pressure turbine (HPT). The simulations are developed using TURBO, which is an unsteady Reynolds-averaged Navier-Stokes (URANS) code capable of multi-stage simulations. The purpose of this study is twofold. First, to determine an estimate of the relative importance of the contributions to the coherent part of the acoustic signature of a turbine from the three potential sources of turbine noise generation, namely, blade-row viscous interaction, potential field interaction, and entropic source associated with the interaction of the blade rows with the temperature nonuniformities caused by the incomplete mixing of the hot fluid and the cooling flow. Second, to develop an understanding of the turbine acoustic transmission characteristics and to assess the applicability of existing empirical and analytical transmission loss models to realistic geometries and flow conditions for modern turbine designs. The investigation so far has concentrated on two simulations: (1) a single-stage HPT and (2) a two-stage HPT and the associated inter-turbine duct/strut segment. The simulations are designed to resolve up to the second harmonic of the blade passing frequency tone in accordance with accepted rules for second order solvers like TURBO. The calculations include blade and vane cooling flows and a radial profile of pressure and temperature at the turbine inlet. The calculation can be modified later to include the combustor pattern factor at the turbine inlet to include that contribution to turbine noise. We shall present preliminary analysis of the results obtained so far in order to assess the validity of such an approach and to seek feedback on improving the approach. This work addresses both Area 1 (Turbine Tone Noise) and Area 5 (Influence of the Turbine on Combustor Noise) topics.

VanZante, Dale↗

Atmospheric Chemical Transport Based on High Resolution Model- Derived Winds: A Case Study

Flight 10 of NASA's Subsonic Assessment (SASS) Ozone and Nitrogen Oxide Experiment (SONEX) extended southwest of Lajes, Azores. A variety of chemical signatures were encountered. These signatures are examined in detail, relating them to meteorological data from a high resolution numerical model having horizontal grid spacing of 30 and 90 km and 26 vertical levels. The meteorological output at hourly intervals is used to create backward trajectories from the locations of the chemical signatures. Four major categories of chemical signatures are discussed-stratospheric, lightning, continental pollution, and a transition layer. The strong stratospheric signal is encountered just south of the Azores in a region of depressed tropopause height. Three chemical signatures at different altitudes in the upper troposphere are attributed to lightning. Backward trajectories arriving at locations of these signatures are related to locations of cloud-to-ground lightning. Results show that the trajectories pass through regions of lightning 1-2 days earlier over the eastern Gulf of Mexico and off the southeast coast of the United States. The lowest leg of the flight exhibits a chemical signature consistent with continental pollution. Trajectories arriving at this signature are found to pass over the highly populated Northeast Corridor of the United States. Surface based pollution apparently is lofted to the altitudes of the trajectories by convective clouds along the East Coast that did not contain lightning. Finally, a chemical transition layer is described. Its chemical signature is intermediate to those of lightning and continental pollution. Trajectories arriving in this layer pass between the trajectories of the lightning and pollution signatures. Thus, they probably are impacted by both sources.

Hannan, John R.↗

Application of linear prolongation to coarse mesh finite difference acceleration in CASMO5

The Coarse-Mesh Finite Difference (CMFD) method has been used for over a decade to accelerate the convergence of the Method of Characteristics (MOC) solution to the two- dimensional particle transport equation in CASMO5. Numerical testing, along with widespread use in production-level calculations, have shown that the current CMFD implementation provides stability and robustness for a wide range of realistic reactor physics problems. However, the recent development of linear prolongation has attracted attention from the community as a way to further improve the performance and stability of CMFD. Two interpolation methods for linear prolongation are presented in this work and implemented into a test version of CASMO5. The performance of the proposed interpolations, referred to as the bilinear and linear directional schemes, is evaluated in terms of runtime relative to the default constant or uniform scaling update. Numerical results indicate that the use of linear prolongation can reduce the transport solver runtime on average by approximately 10% when tested with two hundred randomly selected cases. The new directional linear interpolation, combined with default constant boundary updates, is found to provide the highest reduction in runtime for the cases analyzed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Relations between Shot Noise, Gain Bandwidth, and Saturation of Instabilities

There are numerous instabilities present in charged particle beams that undergo exponential growth and reach saturation. In various applications, such as free-electron lasers or micro-bunching light sources, achieving saturation is desirable. Conversely, there are applications where these instabilities are utilized as linear broad-band amplifiers for signals embedded in the charged beam. In the latter scenario, the saturation of an instability induces non-linear distortions in the imprinted signal, thereby limiting the useful range of such amplifiers. Accurate evaluation of these instabilities necessitates a complete and comprehensive modeling approach that includes shot noise within the beam. Unfortunately, such modeling is not always feasible or practical. In this paper, we introduce a methodology utilizing the frequency and bandwidth of the instability as key parameters. Through this, we derive an estimation for the range of linear instability growth. Our derivation is conducted in a model-independent manner, making it applicable to a broad spectrum of instabilities. To validate our approach, we employ established and thoroughly benchmarked simulations with a free electron laser (FEL) code as well as self-consistent 3-dimensional simulation of plasma-cascade instability using code SPACE.

43 PARTICLE ACCELERATORS↗

Numerical solution of compressible viscous flows

The presented review is concerned with the problem of calculating compressible viscous flows. Basic numerical considerations and problems associated with calculating viscous flows are examined and current numerical approaches toward the solution of the Navier-Stokes equations are discussed. It is pointed out that the numerical solution of the full time-dependent equations for turbulent flow is not practical with present computers. Therefore, turbulence effects must be accounted for by modeling. Developments related to turbulence modeling are described. In connection with a discussion of numerical methods for solving viscous flow equations, attention is given to numerical domains of dependence of typical explicit and implicit methods, the diffusion problem, the convection-diffusion problem, and the split-hybrid method.

Maccormack, R. W.↗

Modeling creep behavior of fiber composites

A micromechanical model for the creep behavior of fiber composites is developed based on a typical cell consisting of a fiber and the surrounding matrix. The fiber is assumed to be linearly elastic and the matrix nonlinearly viscous. The creep strain rate in the matrix is assumed to be a function of stress. The nominal stress-strain relations are derived in the form of differential equations which are solved numerically for off-axis specimens under uniaxial loading. A potential function and the associated effective stress and effective creep strain rates are introduced to simplify the orthotropic relations.

Chen, J. L.↗

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗

Scaling and Systems Considerations in Pulsed Inductive Thrusters

Performance scaling in pulsed inductive thrusters is discussed in the context of previous experimental studies and modeling results. Two processes, propellant ionization and acceleration, are interconnected where overall thruster performance and operation are concerned, but they are separated here to gain physical insight into each process and arrive at quantitative criteria that should be met to address or mitigate inherent inductive thruster difficulties. The effects of preionization in lowering the discharge energy requirements relative to a case where no preionization is employed, and in influencing the location of the initial current sheet, are described. The relevant performance scaling parameters for the acceleration stage are reviewed, emphasizing their physical importance and the numerical values required for efficient acceleration. The scaling parameters are then related to the design of the pulsed power train providing current to the acceleration stage. The impact of various choices in pulsed power train and circuit topology selection are reviewed, paying special attention to how these choices mitigate or exacerbate switching, lifetime, and power consumption issues.

Polzin, Kurt A.↗

Shear Instabilities as a Probe of Jupiter's Atmosphere

Linear wave patterns in Jupiter clouds with wavelengths strongly clustered around 300 km are commonly observed in the planet's equatorial atmosphere. We propose that the preferred wavelength is related to the thickness of an unstable shear layer within the clouds. We numerically analyze the linear stability of wavelike disturbances that have nonzero horizontal phase speeds in Jupiter's atmosphere and find that. if the static stability in the strongly clustered around 300 km are commonly observed in the planet's equatorial atmosphere. We propose that the preferred wavelength is related to the thickness of an unstable shear layer within the clouds. We numerically analyze the linear stability of wavelike disturbances that have nonzero horizontal phase speeds in Jupiter's atmosphere and find that. if the static stability in the shear layer is very low (but still nonnegative), a deep vertical shear layer like the one measured by the Galileo probe can generate the instabilities. The fastest growing waves grow exponentially within an hour, and their wavelengths match the observations. Close to zero values of static stability that permit the growth of instabilities are within the range of values measured by the Galileo probe in a hot spot. Our model probes Jupiter's equatorial atmosphere below the cloud deck and suggests that thick regions of wind shear and low static stability exist outside hot spots.

Bosak, Tanja↗

On an origin of numerical diffusion: Violation of invariance under space-time inversion

The invariant properties of the convection equation du/dt + adu/dx = 0 (where d is the partial differential operator) with respect to spatial reflection, time reversal, and space-time inversion are studied. Generally, a finite-difference analog of this equation may possess some or none of these properties. It is shown that, under certain conditions, the von Neumann amplification factor of an analog satisfies a special relation for each invariant property this analog possesses. Particularly, an analog is neutrally stable and thus free of numerical diffusion if it possesses the invariant property related to space-time inversion. It is also explained why generally (1) an upwind scheme possesses neither the invariant property related to spatial reflection nor that related to space-time inversion, and (2) an explicit scheme possesses neither the invariant property related to time reversal nor that related to space-time inversion. Extension to the viscous case and a remarkable connection between the current work and a new numerical framework for solving conservation laws are also discussed.

Chang, Sin-Chung↗

A Systems Biology Approach to Energy Flow in H 2 Producing Microbial Communities: Microbial Mats

This report describes the results of work conducted at NASA Ames Research Center, in collaboration with investigators at Lawrence Livermore National Laboratory, during the time period 12/15/2014 to 12/15/2017. The work was an extension of a collaboration with LLNL that began in 2008. Work this reporting period focused primarily on tasks related to microbial mats from both the Elkhorn Slough and Baja California field sites, including the analysis and publication of data produced prior to 12/15/2014. During this reporting period, in addition to the microbial mat work, we also conducted a series of experiments with a particular algal-bacterial association (C. sorokiniana and A. brasilense). A number of aspects of this well-characterized algal-bacterial association enabled us not only to further our understanding of interactions between these particular organisms, but to apply that knowledge and techniques to our collection of microorganisms isolated from microbial mats. Microbial mat experiments focused both on intact mat communities and on cultures of organisms isolated from the mats (the Microbial Mat Living Library). In a series of co-culture experiments, we documented both positive and negative interactions amongst heterotrophic and autotrophic isolates of microbial mat microorganisms. One particularly important and novel mat cyanobacterium (ESFC-1) was studied in detail. We developed techniques to pursue an investigation of functional redundancy (the idea that closely related organisms perform identical roles in a community) in cyanobacterial autotrophs in the Elkhorn Slough microbial mat. In those experiments, we documented numerous instances of increases in growth rate between co-cultures of two organisms (relative to their growth rate in isolation) performing what would seem to be identical functional roles. Simplified microbial mat communities consisting of organisms from the Microbial Mat Living Library were constructed and studied in laboratory experiments. The “constructed mats” were grown in standard laboratory glassware and plasticware, but cultures were also combined with artificial materials in an effort to re-create the three-dimensional structure of naturally occurring microbial mats. These constructed mats re-created a number of functions attributable to natural microbial mats (N-fixation, H2-production) and simplified their study using molecular ecological techniques due to the smaller number of microbial “species” present. Cyanobacteria which migrated in response to vertical variations in light intensity (irradiance) obtained characteristic positions in constructed mats that were similar to their positions in natural mats. A detailed investigation of nitrogen cycling in both the Baja California and Elkhorn Slough microbial mats was initiated during this reporting period. The nitrogen cycle in the microbial mats has not been studied in depth, despite the fact that nitrogen usually limits growth in marine environments. The organisms conducting nitrogen cycling in both microbial mats were studied using the genomes available from earlier work on these mats, as well as with a targeted approach using molecular ecological methods developed here. Quantitative measurements of rates of nitrogen transformations were made using stable isotopically labeled nitrogen. The processes quantified and examined using molecular ecological methods included: nitrogen fixation, ammonification, ammonium and nitrate assimilation, denitrification, anaerobic oxidation of ammonium, and dissimilatory nitrate reduction to ammonium. In a series of experiments designed to examine a potential application for microbial mats, we demonstrated that both Elkhorn Slough and Baja California microbial mats were very effective in removing large quantities of both nitrate and ammonium from waste streams containing these compounds at levels where they would be considered to be pollutants. Considerable effort was expended during this reporting period to analyze, and prepare for publication, metagenomes and metatranscriptomes collected from two experimental manipulations of microbial mats. The first, a diel study, was conducted in November 2011, while the second, an investigation of the effects of sulfate on community composition, was conducted during this reporting period. For the earlier study, reference-based and reference-free methods were used to assess microbial ecology and genetic partitioning in these complex microbial systems, and that work is available as a pre-print. Four metagenomes from the same experiment were published. In the second, longer (nearly 200 day) manipulation, rRNA iTag libraries for cDNA and DNA were analyzed to determine the effects of sulfate concentration on microbial mat community composition over the course of the experiment. Work on the algal-bacterial association (Chlorella sorokiniana and Azospirillum brasilense) focused initially on demonstrating the transfer of carbon and nitrogen between the partners, using stable isotopic labeling and the LLNL NanoSIMS. Subsequently, the potential role of diazotrophy (nitrogen fixation) by the bacterial partner in supporting the growth of the association was assessed in a series of experiments. It appears to be the case that, although the bacterial partner demonstrably increases the fitness of the algae under nitrogen limitation, and cells of C. sorokiniana receive nitrogen originally fixed from atmospheric nitrogen by A. brasilense, the amount of nitrogen obtained by diazotrophy is too small to explain the increase in fitness.

08 HYDROGEN↗