Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D Modeling”

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,639 records · Page 91

A new multi-energy neutrino radiation-hydrodynamics code in full general relativity and its application to the gravitational collapse of massive stars

We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and momentum based on a truncated moment formalism. Regarding neutrino opacities, we take into account a baseline set in state-of-the-art simulations, in which inelastic neutrino–electron scattering, thermal neutrino production via pair annihilation, and nucleon–nucleon bremsstrahlung are included. While the Einstein field equations and the spatial advection terms in the radiation-hydrodynamics equations are evolved explicitly, the source terms due to neutrino–matter interactions and energy shift in the radiation moment equations are integrated implicitly by an iteration method. To verify our code, we first perform a series of standard radiation tests with analytical solutions that include the check of gravitational redshift and Doppler shift. A good agreement in these tests supports the reliability of the GR multi-energy neutrino transport scheme. We then conduct several test simulations of core-collapse, bounce, and shock stall of a 15M{sub ⊙} star in the Cartesian coordinates and make a detailed comparison with published results. Our code performs quite well to reproduce the results of full Boltzmann neutrino transport especially before bounce. In the postbounce phase, our code basically performs well, however, there are several differences that are most likely to come from the insufficient spatial resolution in our current 3D-GR models. For clarifying the resolution dependence and extending the code comparison in the late postbounce phase, we discuss that next-generation Exaflops class supercomputers are needed at least.

79 ASTRONOMY AND ASTROPHYSICS↗

A versatile technique to enable sub-milli-Kelvin instrument stability for precise radial velocity measurements: tests with the habitable-zone planet finder

Insufficient instrument thermomechanical stability is one of the many roadblocks for achieving 10 cm s{sup −1} Doppler radial velocity precision, the precision needed to detect Earth-twins orbiting solar-type stars. Highly temperature and pressure stabilized spectrographs allow us to better calibrate out instrumental drifts, thereby helping in distinguishing instrumental noise from astrophysical stellar signals. We present the design and performance of the Environmental Control System (ECS) for the Habitable-zone Planet Finder (HPF), a high-resolution (R = 50,000) fiber-fed near-infrared (NIR) spectrograph for the 10 m Hobby–Eberly Telescope at McDonald Observatory. HPF will operate at 180 K, driven by the choice of an H2RG NIR detector array with a 1.7 μm cutoff. This ECS has demonstrated 0.6 mK rms stability over 15 days at both 180 and 300 K, and maintained high-quality vacuum (<10{sup −7} Torr) over months, during long-term stability tests conducted without a planned passive thermal enclosure surrounding the vacuum chamber. This control scheme is versatile and can be applied as a blueprint to stabilize future NIR and optical high-precision Doppler instruments over a wide temperature range from ∼77 K to elevated room temperatures. A similar ECS is being implemented to stabilize NEID, the NASA/NSF NN-EXPLORE spectrograph for the 3.5 m WIYN telescope at Kitt Peak, operating at 300 K. A [full SolidWorks 3D-CAD model] and a comprehensive parts list of the HPF ECS are included with this manuscript to facilitate the adaptation of this versatile environmental control scheme in the broader astronomical community.

79 ASTRONOMY AND ASTROPHYSICS↗

A STUDY OF SOLAR PHOTOSPHERIC TEMPERATURE GRADIENT VARIATION USING LIMB DARKENING MEASUREMENTS

The variation in area of quiet magnetic network measured over the sunspot cycle should modulate the spatially averaged photospheric temperature gradient, since temperature declines with optical depth more gradually in magnetic flux tube atmospheres. Yet, limb darkening measurements show no dependence upon activity level, even at an rms precision of 0.04%. We study the sensitivity of limb darkening to changes in area filling factor using a 3D MHD model of the magnetized photosphere. The limb darkening change expected from the measured 11-year area variation lies below the level of measured limb darkening variations, for a reasonable range of magnetic flux in quiet network and internetwork regions. So the remarkably constant limb darkening observed over the solar activity cycle is not inconsistent with the measured 11-year change in area of quiet magnetic network. Our findings offer an independent constraint on photospheric temperature gradient changes reported from measurements of the solar spectral irradiance from the Spectral Irradiance Monitor, and recently, from wavelength-differential spectrophotometry using the Solar Optical Telescope aboard the HINODE spacecraft.

79 ASTRONOMY AND ASTROPHYSICS↗

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗

THE SPACE WEATHER OF PROXIMA CENTAURI b

A planet orbiting in the “habitable zone” of our closest neighboring star, Proxima Centauri, has recently been discovered, and the next natural question is whether or not Proxima b is “habitable.” Stellar winds are likely a source of atmospheric erosion that could be particularly severe in the case of M dwarf habitable zone planets that reside close to their parent star. Here, we study the stellar wind conditions that Proxima b experiences over its orbit. We construct 3D MHD models of the wind and magnetic field around Proxima Centauri using a surface magnetic field map for a star of the same spectral type and scaled to match the observed ∼600 G surface magnetic field strength of Proxima. We examine the wind conditions and dynamic pressure over different plausible orbits that sample the constrained parameters of the orbit of Proxima b. For all the parameter space explored, the planet is subject to stellar wind pressures of more than 2000 times those experienced by Earth from the solar wind. During an orbit, Proxima b is also subject to pressure changes of 1–3 orders of magnitude on timescales of a day. Its magnetopause standoff distance consequently undergoes sudden and periodic changes by a factor of 2–5. Proxima b will traverse the interplanetary current sheet twice each orbit, and likely crosses into regions of subsonic wind quite frequently. These effects should be taken into account in any physically realistic assessment or prediction of its atmospheric reservoir, characteristics, and loss.

79 ASTRONOMY AND ASTROPHYSICS↗

The Threatening Magnetic and Plasma Environment of the TRAPPIST-1 Planets

Recently, four additional Earth-mass planets were discovered orbiting the nearby ultracool M8 dwarf, TRAPPIST-1, making a remarkable total of seven planets with equilibrium temperatures compatible with the presence of liquid water on their surface. Temperate terrestrial planets around an M-dwarf orbit close to their parent star, rendering their atmospheres vulnerable to erosion by the stellar wind and energetic electromagnetic and particle radiation. Here, we use state-of-the-art 3D magnetohydrodynamic models to simulate the wind around TRAPPIST-1 and study the conditions at each planetary orbit. All planets experience a stellar wind pressure between 10{sup 3} and 10{sup 5} times the solar wind pressure on Earth. All orbits pass through wind pressure changes of an order of magnitude and most planets spend a large fraction of their orbital period in the sub-Alfvénic regime. For plausible planetary magnetic field strengths, all magnetospheres are greatly compressed and undergo much more dynamic change than that of the Earth. The planetary magnetic fields connect with the stellar radial field over much of the planetary surface, allowing the direct flow of stellar wind particles onto the planetary atmosphere. These conditions could result in strong atmospheric stripping and evaporation and should be taken into account for any realistic assessment of the evolution and habitability of the TRAPPIST-1 planets.

79 ASTRONOMY AND ASTROPHYSICS↗

Physical Correlations and Predictions Emerging from Modern Core-collapse Supernova Theory

Abstract In this paper, we derive correlations between core-collapse supernova observables and progenitor core structures that emerge from our suite of 20 state-of-the-art 3D core-collapse supernova simulations carried to late times. This is the largest such collection of 3D supernova models ever generated and allows one to witness and derive testable patterns that might otherwise be obscured when studying one or a few models in isolation. From this panoramic perspective, we have discovered correlations between explosion energy, neutron star gravitational birth masses, 56 Ni and α -rich freezeout yields, and pulsar kicks and theoretically important correlations with the compactness parameter of progenitor structure. We find a correlation between explosion energy and progenitor mantle binding energy, suggesting that such explosions are self-regulating. We also find a testable correlation between explosion energy and measures of explosion asymmetry, such as the ejecta energy and mass dipoles. While the correlations between two observables are roughly independent of the progenitor zero-age main-sequence (ZAMS) mass, the many correlations we derive with compactness cannot unambiguously be tied to a particular progenitor ZAMS mass. This relationship depends on the compactness/ZAMS mass mapping associated with the massive star progenitor models employed. Therefore, our derived correlations between compactness and observables may be more robust than with ZAMS mass but can nevertheless be used in the future once massive star modeling has converged.

Astronomy & Astrophysics↗

Assessing the Impact of a Novel TBC Material on Heat Transfer in a Spark Ignition Engine through 3D CFD-FEA Co-Simulation Routine

Thermal barrier coatings (TBCs) have been of interest since the 1970s for application in internal combustion (IC) engines. Thin TBCs exhibit a temperature swing phenomenon wherein wall temperatures dynamically respond to the transient working-gas temperature throughout the engine cycle, thus reducing the temperature difference driving the heat transfer. Determining these varying wall temperatures is necessary to evaluate and study the effect of coatings on wall heat transfer. This study focuses on developing a 3D computational fluid dynamics (CFD)-finite element analysis (FEA) coupled simulation, or co-simulation, routine to determine the wall temperatures of a piston coated with a thin TBC layer subject to spark ignition combustion heat flux. A CONVERGE 3D-CFD model was used to simulate the combustion process in a single-cylinder, light-duty experimental spark ignition (SI) engine. Transient piston heat transfer analysis was conducted using ABAQUS, a FEA package, under the simulated combustion heat flux load. The effect of the temperature swing phenomenon due to this TBC layer was observed in a CFD simulation by implementing the FEA results as the piston thermal boundary conditions. The boundary conditions were passed between the CFD and FEA tools until a quasi-steady state solution was achieved. Furthermore, a reduction in wall heat transfer was observed due to a reduced temperature difference between the wall and the working gas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Embedding a Critical Point in a Hadron to Quark–Gluon Crossover Equation of State

It is shown how to embed a critical point in a smooth background equation of state so as to yield the critical exponents and critical amplitude ratios expected of a transition in the same universality class as the liquid–gas phase transition and the 3D Ising model. There are only two independent critical exponents; the relations α + 2β + $\gamma$ = 2 and β(δ - 1) = $\gamma$ arise automatically, as does a new relation between the two critical amplitudes. The resulting equation of state has parameters which may be inferred by hydrodynamic modeling of heavy-ion collisions in the Beam Energy Scan II at the Relativistic Heavy Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Tracing the Emergence of Collectivity Phenomena in Small Systems

We study initial state momentum correlations and event-by-event geometry in p+Pb collisions at $\sqrt{2}$ = 5.02 TeV by following the approach of extending the IP-Glasma model to 3D using JIMWLK rapidity evolution. On examining the non-trivial rapidity dependence of the observables, we find that the geometry is correlated over large rapidity intervals, while the initial state momentum correlations have a relatively short range in rapidity. Based on our results, we discuss implications for the relevance of both effects in explaining the origin of collective phenomena in small systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Shattered pellet penetration in low and high energy plasmas on DIII-D

Shattered pellet injection (SPI) has been adopted as the baseline disruption mitigation system for ITER, as the radiative payload penetration into DIII-D plasmas from SPI is superior to those using the massive gas injection (MGI) method. Because of the substantial differences in the energy content of ITER plasma and those in present experiments, reliable 3D MHD modeling, benchmarked against present experiments is needed to project to ITER plasmas. In support of these needs, the depth of SPI fragment penetration in DIII-D plasmas was investigated by injecting SPI into two discharges with vastly different energy content and pedestal height. 400 Torr-L pure Ne fragmented pellets at a velocity of about 200 m s−1 were injected into a 0.2 MJ L-mode discharge and a 2 MJ super H-mode discharge. Results show deep penetration of SPI fragments into low-energy plasmas in DIII-D. SPI fragment penetration is reduced as the plasma energy content increases, with some discharges exhibiting penetration that is confined to the outer regions of the plasma. The injected SPI fragments are also spread out over a distance of about 20 cm, which results in some fragments arriving near the end of or after the thermal quench is over.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advanced manufacturing development of a composite empennage component for l-1011 aircraft

Tooling concepts were developed which would permit co-couring of the hat stiffeners to the skin to form the cover assembly in a single autoclave cycle. These tooling concepts include the use of solid rubber mandrels, foam mandrels, and formed elastometric bladders. A simplification of the root end design of the cover hat stiffeners was accomplished in order to facilitate fabrication. The conversion of the 3D NASTRAN model from level 15 to level 16 was completed and a successful check run accomplished. A detailed analysis of the thermal load requirement for the environmental chambers was carried out. Based on the thermal analysis, best function requirements, load inputs and ease of access, a system involving four chambers, two for the covers containing 6 and 4 specimens, respectively, and two for the spares containing 6 and 4 specimens, respectively, evolved.

Source record↗

Effects of stratospheric lapse rate on numerically simulated thunderstorm cloud top structure

The study of overshooting midlatitude severe thunderstorm tops has led to the detection of some intriguing patterns. This paper provides a summary of an investigation of the effects of stratospheric stability on thunderstorm cloud top structure, taking into account the utilization of a fully three-dimensional (3D) anelastic model discussed by Schlesinger (1975, 1978, 1984). The layout of model experiments is discussed, giving attention to initial environmental soundings, and initial forcing. The time variation of the vertical velocity extreme is examined, and a description of the mature storm structure is presented. Constant-altitude plane views throughout storm depth are considered along with cloud top height fields. With the stratospheric lapse, the warm region is small and breaks down into separate spots with secondary cold spots downshear of them.

Schlesinger, R. E.↗

Analysis of the low gravity tolerance of Bridgman-Stockbarger crystal growth. I - Steady and impulse accelerations

The effects of steady and impulse-type residual accelerations on dopant distributions during directional solidification in 2D and 3D 'generic' models of the Bridgman-Stockbarger technique are investigated using numerical methods. The calculations are based on the thermophysical properties of molten germanium doped with a low concentration of gallium. A Chebyshev collocation pseudospectral method is used for the solution of the governing momentum-, mass-, species-, and heat-transfer equations. Only convection caused by temperature gradients is considered. It is found that lateral nonuniformity in composition is very sensitive to the orientation of the steady component of the residual gravity vector and to the particular operating conditions under consideration. It is also found that laterally or radially averaged composition profiles are alone insufficient to describe the extent of residual convection in a spacecraft environment. The effects of impulse-type disturbances can be severe and can extend for times on the order of 1000 sec after the termination of the impulse.

Alexander, J. Iwan D.↗

Space shuttle main engine computed tomography applications

For the past two years the potential applications of computed tomography to the fabrication and overhaul of the Space Shuttle Main Engine were evaluated. Application tests were performed at various government and manufacturer facilities with equipment produced by four different manufacturers. The hardware scanned varied in size and complexity from a small temperature sensor and turbine blades to an assembled heat exchanger and main injector oxidizer inlet manifold. The evaluation of capabilities included the ability to identify and locate internal flaws, measure the depth of surface cracks, measure wall thickness, compare manifold design contours to actual part contours, perform automatic dimensional inspections, generate 3D computer models of actual parts, and image the relationship of the details in a complex assembly. The capabilities evaluated, with the exception of measuring the depth of surface flaws, demonstrated the existing and potential ability to perform many beneficial Space Shuttle Main Engine applications.

Sporny, Richard F.↗

On dipolarization and its relation to the substorm current wedge

The near-earth properties of magnetospheric substorms are investigated using the Hesse and Birn (1991) self-consistent 3D MHD model of magnetotail activity. The simulation shows that an increase of Bz, a decrease of the cross-tail current, and an increase of region-1 field-aligned currents occurs roughly simultaneously, consistent with the empirical picture described by McPherron et al. (1973). It was also found that the region where the cross-tail current is reduced and diverted is located tailward from the dipolarization region. The results indicate that the apparent tailward propagation of the dipolarization effects does not support an earlier suggestion that the underlying process is initiated further earthward.

Hesse, Michael↗

Rate calculations for the simultaneous vibrational relaxation and dissociation of nitrogen

Experimental data on vibrational relaxation and dissociation rates in highly nonequilibrium gases are scarce. Recently there have been several attempts at generating these rates by direct calculation of inelastic collision cross-sections. Some issues that need to be considered if such calculations are to be used as a basis for analysis and design are examined. In particular, the following are examined: (1) collinear vs 3D scattering models, (2) accurate intermolecular potentials, (3) closed channels and acceleration schemes, (4) thermal averaging and thermal rates, and (5) code validation. It has been necessary to confront these issues during the course of N2-N2 and N2-N inelastic cross-section calculations.

Gonzales, David A.↗

The VRT gas turbine combustor - Phase II

An innovative annular combustor configuration is being developed for aircraft and other gas turbine engines. This design has the potential of permitting higher turbine inlet temperatures by reducing the pattern factor and providing a major reduction in NO(x) emission. The design concept is based on a Variable Residence Time (VRT) technique which allows large fuel particles adequate time to completely burn in the circumferentially mixed primary zone. High durability of the combustor is achieved by dual-function use of the incoming air. In Phase I, the feasibility of the concept was demonstrated by water analogue tests and 3D computer modeling. The flow pattern within the combustor was as predicted. The VRT combustor uses only half the number of fuel nozzles of the conventional configuration. In Phase II, hardware was designed, procured, and tested under conditions simulating typical supersonic civil aircraft cruise conditions to the limits of the rig. The test results confirmed many of the superior performance predictions of the VRT concept. The Hastelloy X liner showed no signs of distress after nearly six hours of tests using JP5 fuel.

Melconian, Jerry O.↗