Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “D CODES”

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 829 records · Page 46

Evaluation of DIII-D plasmas for Future Measurements of the Driven Current Density in the DIII-D High Field Side Lower Hybrid Experiment

High field side (HFS) lower hybrid current drive (LHCD) is a promising method for efficiently driving off-axis current in steady state tokamak power plants. The first test of this technology is underway at the DIII-D tokamak. The initial physics goal is to measure the LH-driven current density profile in order to validate the ray-tracing/Fokker-Plank codes GENRAY/CQL3D and gain confidence in their ability to predict LH deposition in future DIII-D experiments and fusion power plants. DIII-D’s motional Stark effect (MSE) diagnostic will be used to constrain the magnetic equilibrium reconstructions. From these reconstructions, the ohmic, bootstrap, and driven non-inductive components of the current density profile can be extracted. Several recent DIII-D plasmas have been identified in which GENRAY/CQL3D predict measurable amounts of current with only 100 kW of coupled power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toroidal phase of post-disruption runaway electron loss to wall in presence of applied 3-D fields in DIII-D

The toroidal distribution of runaway electrons (REs) striking the centre post in DIII-D post-disruption RE final loss events is simulated utilizing the MARS-F code combined with the REORBIT module, by tracing the guiding-centre drift orbits of test REs in the presence of external (applied) plus internal (intrinsic instability) 3-D fields. To better recover the experimental results, three different equilibria with different safety factor profiles are adopted. Dominant resistive magnetohydrodynamic instabilities are found to be m/n=2/1, 3/2 and 1/1 modes (where m and n are the poloidal and toroidal mode numbers). Externally applied n=1 resonant magnetic perturbations (RMPs), with both even and odd parity configuration, are considered in each equilibrium. Locking of MHD instability to the applied RMPs is simulated by minimizing the total perturbed magnetic energy in the plasma. Here, the simulated toroidal locations, where most REs strike in the presence of an edge m/n= 2/1 instability and the n=1 RMPs in even parity, are found to agree well with experiments. Modeling also captures the measured n = 2 RE loss pattern on the wall for experiments where an edge m/n=3/2 instability is dominant instead of the m/n=2/1 mode. The key overall findings include (a) the measured global RE toroidal impact pattern in experiments can be explained by a locked resistive kink instability, (b) the toroidal phase of the locked mode is well predicted by minimizing the total perturbed magnetic energy, and (c) this edge locked mode determines the toroidal phase of the peak RE wall impact.

Bai, Xue [University of California San Diego, La J↗

The 3-D inelastic analysis methods for hot section components (base program)

A series of computer codes are sought that would permit more accurate and efficient three-dimensional inelastic structural analysis of combustor liners, turbine blades, and turbine vanes. Each code embodies a progression of mathematical models for increasingly comprehensive representation of the geometrical features, loading conditions, and forms of nonlinear material response that distinguish these three groups of hot section components.

Todd, E. S.↗

Implementation of algebraic stress models in a general 3-D Navier-Stokes method (PAB3D)

A three-dimensional multiblock Navier-Stokes code, PAB3D, which was developed for propulsion integration and general aerodynamic analysis, has been used extensively by NASA Langley and other organizations to perform both internal (exhaust) and external flow analysis of complex aircraft configurations. This code was designed to solve the simplified Reynolds Averaged Navier-Stokes equations. A two-equation k-epsilon turbulence model has been used with considerable success, especially for attached flows. Accurate predicting of transonic shock wave location and pressure recovery in separated flow regions has been more difficult. Two algebraic Reynolds stress models (ASM) have been recently implemented in the code that greatly improved the code's ability to predict these difficult flow conditions. Good agreement with Direct Numerical Simulation (DNS) for a subsonic flat plate was achieved with ASM's developed by Shih, Zhu, and Lumley and Gatski and Speziale. Good predictions were also achieved at subsonic and transonic Mach numbers for shock location and trailing edge boattail pressure recovery on a single-engine afterbody/nozzle model.

Abdol-Hamid, Khaled S.↗

High sensitivity of simulated fog properties to parameterized aerosol activation in case studies from ParisFog

Aerosols influence fog properties such as visibility and lifetime by affecting fog droplet number concentrations (N d ). Numerical weather prediction (NWP) models often represent aerosol–fog interactions using highly simplified approaches. Incorporating prognostic size-resolved aerosol microphysics from climate models could allow them to simulate N d and aerosol–fog interactions without incurring excessive computational expense. However, microphysics code designed for coarse spatial resolution may struggle with sub-kilometer-scale grid spacings. Here, we test the ability of the UK Met Office Unified Model to simulate aerosol and fog properties during case studies from the ParisFog field campaign in 2011. We examine the sensitivity of fog properties to variations in N d caused by modifications to simulated aerosol activation. Our model, with a 500 m horizontal resolution and interactive aerosol and cloud microphysics, significantly underpredicts N d , although it only slightly underestimates the cloud condensation nuclei concentration. With an updated version of the Abdul-Razzak and Ghan (2000) activation scheme, we produce N d that are more consistent with those predicted by a cloud parcel model under fog-like conditions. We activate droplets only by adiabatic cooling. We incorporate more realistic hygroscopicities for sulfate and organic aerosols and explore the sensitivity of simulated N d to unresolved updrafts. We find that both N d and simulated fog liquid water content are very sensitive to the updated activation scheme but remain less affected by the update to hygroscopicities. Our improvements offer insights into the physical processes regulating N d in stable conditions, potentially laying foundations for improved operational fog forecasts that incorporate interactive aerosol simulations or aerosol climatologies.

Ghosh, Pratapaditya [Carnegie Mellon University, P↗

A 3-D Navier-Stokes CFD study of turbojet/ramjet nozzle plume interactions at Mach 3.0 and comparison with data

Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually consist of a single integrated turboramjet or as in this study, a turbojet housed in an upper bay with a separate ramjet housed in a lower bay. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. One nozzle concept under consideration has a common nozzle, where the plumes from the turbojet and ramjet interact with one another as they expand to ambient conditions. In this paper, the two plumes interact at the end of a common 2-D cowl, when they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. At this condition, the turbojet engine operated at a higher NPR than the ramjet, where the turbojet overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data used to compare with the PARC3D code with k-kl two equation turbulence model. The 2-D and 3-D centerline CFD solutions are in good agreement, but as the CFD solutions approach the outer sidewall, a slight variance occurs. The outer wall boundary layers are thin and do not present much of an interaction, however, where the confluence interaction shocks interact with the thin boundary layer on the outer wall, strong vortices run down each shock causing substantial disturbances in the boundary layer. These disturbances amplify somewhat as they propagate downstream axially from the confluence point. The nozzle coefficient (CFG) is reduced 1/2 percent as a result of this sidewall interaction, from 0.9850 to 0.9807. This three-dimensional reduction is in better agreement with the experimental value of 0.9790.

Chang, Ing↗

Research activities in nuclear astrophysics and related areas

NASA/GRO grant NAG 5-2081, at the University of Chicago, has provided support for a broad program of theoretical research in nuclear astrophysics and related areas, with regard to gamma-ray and hard X-ray emission from classical nova explosions. This research emphasized the possible detection of 22Na gamma-ray line emission from nearby novae involving ONeMg white dwarfs, the detailed examination of 26Al production in novae, and the possible detection of the predicted early gamma ray emission from novae that arises from the decay of the short lived, positron emitting isotopes of CNO elements. Studies of nova related problems have consumed an increasing fraction of the Principal Investigator's research efforts over the past decade. Current research addresses problems associated with the standard model for the outbursts of the classical novae: the occurrence of thermonuclear runaways (TNR) in the accreted hydrogen rich envelopes on white dwarfs in close binary systems (see, e.g., the reviews by Truran 1982; and Shara 1989). Research in progress and planned for the next three years has three main objectives: (1) to gain an improved understanding of the early evolution of the light curves of, particularly, the fastest novae; (2) to gain an improved understanding of the relative importance of the various possible mechanisms of envelope hydrogen depletion (e.g. winds, common envelope driven mass loss, and nuclear burning) to the long term evolution of novae in outburst; and (3) to seek to provide a somewhat more definitive statement of the role of classical novae in nucleosynthesis. Our proposed 2-D studies of convection during the early phases of the TNR and our systematic attempt to incorporate an improved treatment of radiation hydrodynamics into the hydrodynamic code utilized in our calculations, are particularly relevant to the first of these objectives. Further 2-D studies of the effects of common envelope evolution are intended to provide more realistic constraints on the mass depletion mechanisms. Finally, detailed calculations of the thermonuclear history of the matter ejected in novae will be carried out for representative nova configurations involving both carbon-oxygen (CO) and oxygen-neon-magnesium (ONeMg) white dwarfs.

Source record↗

Simulation of DIII-D disruption with argon pellet injection and runaway electron beam

Abstract The next generation of large tokamaks, including ITER, will be equipped with a disruption mitigation system (DMS) that can be activated if a disruption is deemed to be imminent. Introducing impurities by pellet (large or shattered) or massive gas injection has been shown to be an effective mitigation mechanism on many tokamaks. The goal of the mitigation is to lessen the thermal and electromagnetic loads from the disruption without generating enough high-energy (runaway) electrons to damage the device. Variations of this mitigation process with impurity injection are presently being tested on many experiments. We have modeled one such impurity injection experiment on DIII-D using the M3D-C1 nonlinear 3D extended MHD code (Jardin et al 2012 Comput. Sci. Discovery 6 014002), The model includes an argon large pellet injection and ablation model, impurity ionization, recombination, and radiation, and runaway electron formation and subsequent evolution, including both Dreicer and avalanche sources. We obtain reasonable agreement with the experimental results for the timescale of the thermal and current quench and for the magnitude of the runaway electron plateau formed during the mitigation. This is the first 3D full MHD simulation with pellets and REs to simulate the disruption process and it also provides a partial validation of the M3D-C1 DMS model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

HPB_strengthmodel

Python-implementation of the Hunter-Preston strength model, but with generalized drag coefficient B. For details, see D. N. Blaschke, A. Hunter, and D. L. Preston, Int. J. Plast. 131 (2020) 102750. This code was used to generate most of the figures in this paper.

Blaschke, Daniel N.↗

Data Repository for Multi-Objective Urban Observational Strategies: A risk-based framework for expanding flood sensor networks.

These data support the manuscript "Multi-Objective Urban Observational Strategies: A risk-based framework for expanding flood sensor networks." These data are generated to allow water managers to reason about optimal locations to expand a flood observation system from multiple perspectives, specifically focusing on flood hazards, and population exposure to flooding. The data included are a) a shapefile of individual sensor locations b) a shapefile of river reach catchments, c) raster of FEMA flood likelihood layers d) shapefile of population locations and population socioeconomic characteristics. The code is written in R and includes all files necessary to generate the figures for the associated manuscript. Interactive maps of the final calculated maps of hazard, vulnerability, exposure, and risk are also included as html files.

54 ENVIRONMENTAL SCIENCES↗

3-D reconstruction and microstructural characterization of neutron-irradiated U-10Zr fuel using FIB-SEM serial sectioning

Herein, focused ion beam-scanning electron microscopy serial sectioning was applied to characterize the three-dimensional (3-D) porosity and phase regions of a neutron-irradiated U-10 wt% Zr fuel. The specimen was removed from an intermediate radial region of a fuel pin irradiated to 5.7 at.% burn-up. Backscattered electron imaging and energy-dispersive spectroscopy were performed on each serial section, allowing for the characterization of microstructural morphology and composition. Porosity size followed a lognormal distribution, ranging from 1.46 × 10 –4 to 25.58 µm 3 with a total porosity volume fraction of 13.02%. Distinctive microstructural regions were identified by composition and porosity: (1) a Zr-rich region with an average composition of 28.4 wt% Zr and a local porosity fraction of 6.88%, and (2) a U-rich region with an average composition of 97.0 wt% U and a local porosity fraction of 14.11%; subdivided into U-rich—high porosity (16.68%) and U-rich—low porosity (8.04%) regions. The detailed 3-D compositional and porosity regions can improve nuclear fuel performance codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Seismic Source Modeling Software Enhancements (FY21)

Seismic source modeling allows researchers both to simulate how a source that induces seismic waves interacts with the Earth to produce observed seismograms and, inversely, to infer what the time histories, sizes, and force distributions were for a seismic source given observed seismograms. In this report, we discuss improvements made in FY21 to our software as applies to both the forward and inverse seismic source modeling problems. For the forward portion of the problem, we have added the ability to use full 3-D nonlinear simulations by implementing 3-D time varying boundary conditions within Sandia’s linear seismic code Parelasti. Secondly, on the inverse source modeling side, we have developed software that allows us to invert seismic gradiometer-derived observations in conjunction with standard translational motion seismic data to infer properties of the source that may improve characterization in certain circumstances. First, we describe the basic theory behind each software enhancement and then demonstrate the software in action with some simple examples.

58 GEOSCIENCES↗

Performance measures for transform data coding.

This paper develops performance criteria for evaluating transform data coding schemes under computational constraints. Computational constraints that conform with the proposed basis-restricted model give rise to suboptimal coding efficiency characterized by a rate-distortion relation R(D) similar in form to the theoretical rate-distortion function. Numerical examples of this performance measure are presented for Fourier, Walsh, Haar, and Karhunen-Loeve transforms.

Pearl, J.↗

Perspectives on the mixing of a row of jets with a confined crossflow

An interactive computer code, written for a microcomputer, is presented which displays 2-D and 3-D oblique plots of the temperature distribution downstream of jets mixing with a confined crossflow, for either single-side or opposed jet injection. Temperature profiles calculated with this routine are presented to show the effects of flow and geometric variables on the mixing. Examples are shown to illustrate the different perspectives on the mixing available by exercising various view options. In addition, the program is used to calculate profiles for opposed rows of jets with their centerlines in-line, by assuming that the confining effect of an opposite wall is equivalent to that of a plane of symmetry between opposed jets.

Holdeman, J. D.↗

Grid flexibility and patching techniques

The numerical determination of combustor flowfields is of great value to the combustor designer. An a priori knowledge of the flow behavior can speed the combustor design process and reduce the number of experimental test rigs required to arrive at an optimal design. Even 2-D steady incompressible isothermal flow predictions are of use; many codes of this kind are available, each employing different techniques to surmount the difficulties arising from the nonlinearity of the governing equations and from typically irregular combustor geometries. Mapping techniques (algebraic and elliptic PDE), and adaptive grid methods (both multi-grid and grid embedding) as applied to axisymmetric combustors are discussed.

Keith, T. G.↗

Low thrust chemical rocket technology

An on-going technology program to improve the performance of low thrust chemical rockets for spacecraft on-board propulsion applications is reviewed. Improved performance and lifetime is sought by the development of new predictive tools to understand the combustion and flow physics, introduction of high temperature materials and improved component designs to optimize performance, and use of higher performance propellants. Improved predictive technology is sought through the comparison of both local and global predictions with experimental data. Predictions are based on both the RPLUS Navier-Stokes code with finite rate kinetics and the JANNAF methodology. Data were obtained with laser-based diagnostics along with global performance measurements. Results indicate that the modeling of the injector and the combustion process needs improvement in these codes and flow visualization with a technique such as 2-D laser induced fluorescence (LIF) would aid in resolving issues of flow symmetry and shear layer combustion processes. High temperature material fabrication processes are under development and small rockets are being designed, fabricated, and tested using these new materials. Rhenium coated with iridium for oxidation protection was produced by the Chemical Vapor Deposition (CVD) process and enabled an 800 K increase in rocket operating temperature. Performance gains with this material in rockets using Earth storable propellants (nitrogen tetroxide and monomethylhydrazine or hydrazine) were obtained through component redesign to eliminate fuel film cooling and its associated combustion inefficiency while managing head end thermal soakback. Material interdiffusion and oxidation characteristics indicated that the requisite lifetimes of tens of hours were available for thruster applications. Rockets were designed, fabricated, and tested with thrusts of 22, 62, 440 and 550 N. Performance improvements of 10 to 20 seconds specific impulse were demonstrated. Higher performance propellants were evaluated: Space storable propellants, including liquid oxygen (LOX) as the oxidizer with nitrogen hydrides or hydrocarbon as fuels. Specifically, a LOX/hydrazine engine was designed, fabricated, and shown to have a 95 pct theoretical c-star which translates into a projected vacuum specific impulse of 345 seconds at an area ratio of 204:1. Further performance improvment can be obtained by the use of LOX/hydrogen propellants, especially for manned spacecraft applications, and specific designs must be developed and advanced through flight qualification.

Schneider, Steven J.↗