Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete phase model”

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 271 records · Page 15

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗

A ring-source model for jet noise

A model consisting of two ring sources was developed to study the direct radiation of jet noise in terms of correlation, coherence, and phase and also to aid in solving the inverse radiation problem of determining the noise source in terms of far-field measurements. The rings consist of discrete sources which are either monopoles or quadrupoles with Gaussian profiles. Only adjacent sources, both within the rings and between rings, are correlated. Results show that from the far-field information can be used to determine when the sources are compact or noncompact with respect to the acoustic wavelength and to distinguish between the types of sources. In addition, from the inverse radiation approach, the center of mass, the location and separation distance of the ring, and the diameters can be recovered.

Maestrello, L.↗

Internal Flow Simulation of Enhanced Performance Solid Rocket Booster for the Space Transportation System

An enhanced performance solid rocket booster concept for the space shuttle system has been proposed. The concept booster will have strong commonality with the existing, proven, reliable four-segment Space Shuttle Reusable Solid Rocket Motors (RSRM) with individual component design (nozzle, insulator, etc.) optimized for a five-segment configuration. Increased performance is desirable to further enhance safety/reliability and/or increase payload capability. Performance increase will be achieved by adding a fifth propellant segment to the current four-segment booster and opening the throat to accommodate the increased mass flow while maintaining current pressure levels. One development concept under consideration is the static test of a "standard" RSRM with a fifth propellant segment inserted and appropriate minimum motor modifications. Feasibility studies are being conducted to assess the potential for any significant departure in component performance/loading from the well-characterized RSRM. An area of concern is the aft motor (submerged nozzle inlet, aft dome, etc.) where the altered internal flow resulting from the performance enhancing features (25% increase in mass flow rate, higher Mach numbers, modified subsonic nozzle contour) may result in increased component erosion and char. To assess this issue and to define the minimum design changes required to successfully static test a fifth segment RSRM engineering test motor, internal flow studies have been initiated. Internal aero-thermal environments were quantified in terms of conventional convective heating and discrete phase alumina particle impact/concentration and accretion calculations via Computational Fluid Dynamics (CFD) simulation. Two sets of comparative CFD simulations of the RSRM and the five-segment (IBM) concept motor were conducted with CFD commercial code FLUENT. The first simulation involved a two-dimensional axi-symmetric model of the full motor, initial grain RSRM. The second set of analyses included three-dimensional models of the RSRM and FSM aft motors with four-degree vectored nozzles.

Ahmad, Rashid A.↗

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗

National facilities study. Volume 4: Space operations facilities task group

The principal objectives of the National Facilities Study (NFS) were to: (1) determine where U.S. facilities do not meet national aerospace needs; (2) define new facilities required to make U.S. capabilities 'world class' where such improvements are in the national interest; (3) define where consolidation and phase-out of existing facilities is appropriate; and (4) develop a long-term national plan for world-class facility acquisition and shared usage. The Space Operations Facilities Task Group defined discrete tasks to accomplish the above objectives within the scope of the study. An assessment of national space operations facilities was conducted to determine the nation's capability to meet the requirements of space operations during the next 30 years. The mission model used in the study to define facility requirements is described in Volume 3. Based on this model, the major focus of the Task Group was to identify any substantive overlap or underutilization of space operations facilities and to identify any facility shortfalls that would necessitate facility upgrades or new facilities. The focus of this initial study was directed toward facility recommendations related to consolidations, closures, enhancements, and upgrades considered necessary to efficiently and effectively support the baseline requirements model. Activities related to identifying facility needs or recommendations for enhancing U.S. international competitiveness and achieving world-class capability, where appropriate, were deferred to a subsequent study phase.

Source record↗

Ultramicrowave communications system, phase 3

The ultramicrowave communications system program investigated the feasibility of a solid state system that meets the projected space to space requirements, while using the advantages of the 100 to 200 GHz band. The program successfully demonstrated a laboratory model of a high frequency communications system operating between 100 to 200 GHz. In the process, vendor claims for performance specifications of discrete components were evaluated, and a window was provided into system design and integration problems.

Source record↗

Modeling Optical Properties of Mineral Aerosol Particles by Using Nonsymmetric Hexahedra

We explore the use of nonsymmetric geometries to simulate the single-scattering properties of airborne dust particles with complicated morphologies. Specifically, the shapes of irregular dust particles are assumed to be nonsymmetric hexahedra defined by using the Monte Carlo method. A combination of the discrete dipole approximation method and an improved geometric optics method is employed to compute the single-scattering properties of dust particles for size parameters ranging from 0.5 to 3000. The primary optical effect of eliminating the geometric symmetry of regular hexahedra is to smooth the scattering features in the phase function and to decrease the backscatter. The optical properties of the nonsymmetric hexahedra are used to mimic the laboratory measurements. It is demonstrated that a relatively close agreement can be achieved by using only one shape of nonsymmetric hexahedra. The agreement between the theoretical results and their measurement counterparts can be further improved by using a mixture of nonsymmetric hexahedra. It is also shown that the hexahedron model is much more appropriate than the "equivalent sphere" model for simulating the optical properties of dust particles, particularly, in the case of the elements of the phase matrix that associated with the polarization state of scattered light.

Bi, Lei↗

Lateral conduction effects on heat-transfer data obtained with the phase-change paint technique

A computerized tool, CAPE, (Conduction Analysis Program using Eigenvalues) has been developed to account for lateral heat conduction in wind tunnel models in the data reduction of the phase-change paint technique. The tool also accounts for the effects of finite thickness (thin wings) and surface curvature. A special reduction procedure using just one time of melt is also possible on leading edges. A novel iterative numerical scheme was used, with discretized spatial coordinates but analytic integration in time, to solve the inverse conduction problem involved in the data reduction. A yes-no chart is provided which tells the test engineer when various corrections are large enough so that CAPE should be used. The accuracy of the phase-change paint technique in the presence of finite thickness and lateral conduction is also investigated.

Maise, G.↗

Fluid-Kinetic Coupling: Advanced Discretizations for Simulations on Emerging Heterogeneous Architectures (LDRD FY20-0643)

Plasma physics simulations are vital for a host of Sandia mission concerns, for fundamental science, and for clean energy in the form of fusion power. Sandia's most mature plasma physics simulation capabilities come in the form of particle-in-cell (PIC) models and magnetohydrodynamics (MHD) models. MHD models for a plasma work well in denser plasma regimes when there is enough material that the plasma approximates a fluid. PIC models, on the other hand, work well in lower-density regimes, in which there is not too much to simulate; error in PIC scales as the square root of the number of particles, making high-accuracy simulations expensive. Real-world applications, however, almost always involve a transition region between the high-density regimes where MHD is appropriate, and the low-density regimes for PIC. In such a transition region, a direct discretization of Vlasov is appropriate. Such discretizations come with their own computational costs, however; the phase-space mesh for Vlasov can involve up to six dimensions (seven if time is included), and to apply appropriate homogeneous boundary conditions in velocity space requires meshing a substantial padding region to ensure that the distribution remains sufficiently close to zero at the velocity boundaries. Moreover, for collisional plasmas, the right-hand side of the Vlasov equation is a collision operator, which is non-local in velocity space, and which may dominate the cost of the Vlasov solver. The present LDRD project endeavors to develop modern, foundational tools for the development of continuum-kinetic Vlasov solvers, using the discontinuous Petrov-Galerkin (DPG) methodology, for discretization of Vlasov, and machine-learning (ML) models to enable efficient evaluation of collision operators. DPG affords several key advantages. First, it has a built-in, robust error indicator, allowing us to adapt the mesh in a very natural way, enabling a coarse velocity-space mesh near the homogeneous boundaries, and a fine mesh where the solution has fine features. Second, it is an inherently high-order, high-intensity method, requiring extra local computations to determine so-called optimal test functions, which makes it particularly suited to modern hardware in which floating-point throughput is increasing at a faster rate than memory bandwidth. Finally, DPG is a residual-minimizing method, which enables high-accuracy computation: in typical cases, the method delivers something very close to the $L^2$ projection of the exact solution. Meanwhile, the ML-based collision model we adopt affords a cost structure that scales as the square root of a standard direct evaluation. Moreover, we design our model to conserve mass, momentum, and energy by construction, and our approach to training is highly flexible, in that it can incorporate not only synthetic data from direct-simulation Monte Carlo (DSMC) codes, but also experimental data. We have developed two DPG formulations for Vlasov-Poisson: a time-marching, backward-Euler discretization and a space-time discretization. We have conducted a number of numerical experiments to verify the approach in a 1D1V setting. In this report, we detail these formulations and experiments. We also summarize some new theoretical results developed as part of this project (published as papers previously): some new analysis of DPG for the convection-reaction problem (of which the Vlasov equation is an instance), a new exponential integrator for DPG, and some numerical exploration of various DPG-based time-marching approaches to the heat equation. As part of this work, we have contributed extensively to the Camellia open-source library; we also describe the new capabilities and their usage. We have also developed a well-documented methodology for single-species collision operators, which we applied to argon and demonstrated with numerical experiments. We summarize those results here, as well as describing at a high level a design extending the methodology to multi-species operators. We have released a new open-source library, MLC, under a BSD license; we include a summary of its capabilities as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Principles of two-dimensional terahertz spectroscopy of collective excitations: The case of Josephson plasmons in layered superconductors

Two-dimensional terahertz spectroscopy (2DTS), a terahertz analog of nuclear magnetic resonance, is a new technique poised to address many open questions in complex condensed matter systems. The conventional theoretical framework used ubiquitously for interpreting multidimensional spectra of discrete quantum level systems is, however, insufficient for the continua of collective excitations in strongly correlated materials. Here, in this work, we develop a theory for 2DTS of a model collective excitation, the Josephson plasma resonance in layered superconductors. Starting from a mean-field approach at temperatures well below the superconducting phase transition, we obtain expressions for the multidimensional nonlinear responses that are amenable to intuition derived from the conventional single-mode scenario. We then consider temperatures near the superconducting critical temperature T c , where dynamics beyond mean-field become important and conventional intuition fails. As fluctuations proliferate near T c , the dominant contribution to nonlinear response comes from an optical parametric drive of counterpropagating Josephson plasmons, which gives rise to 2D spectra that are qualitatively different from the mean-field predictions. As such, and in contrast to one-dimensional spectroscopy techniques, such as third harmonic generation, 2DTS can be used to directly probe thermally excited finite-momentum plasmons and their interactions. Our theory can readily be tested in cuprates, and we discuss implications beyond the present context of Josephson plasmons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Preliminary results in using Deep Learning to emulate BLOB, a nuclear interaction model

Purpose: A reliable model to simulate nuclear interactions is fundamental for Ion-therapy. We already showed how BLOB (“Boltzmann-Langevin One Body”), a model developed to simulate heavy ion interactions up to few hundreds of MeV/u, could simulate also 12 C reactions in the same energy domain. However, its computation time is too long for any medical application. For this reason we present the possibility of emulating it with a Deep Learning algorithm. Methods: The BLOB final state is a Probability Density Function (PDF) of finding a nucleon in a position of the phase space. We discretised this PDF and trained a Variational Auto-Encoder (VAE) to reproduce such a discrete PDF. As a proof of concept, we developed and trained a VAE to emulate BLOB in simulating the interactions of 12 C with 12 C at 62 MeV/u. To have more control on the generation, we forced the VAE latent space to be organised with respect to the impact parameter (b) training a classifier of b jointly with the VAE. Results: In this work, the distributions obtained from the VAE are similar to the input ones and the computation time needed to use the VAE as a generator is negligible. Conclusions: We show that it is possible to use a Deep Learning approach to emulate a model developed to simulate nuclear reactions in the energy range of interest for Ion-therapy. We foresee the implementation of the generation part in C++ and to interface it with the most used Monte Carlo toolkit: Geant4.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Testing and Modeling of Functionally Graded Aluminum-Doped Zinc Oxide Using Spark Plasma Sintering and Discrete Powder Layers of Varying Composition

Functionally graded material (FGM) ZnO is made via spark plasma sintering/field-assisted sintering technique (SPS/FAST) by varying the Al dopant content along the z-axis or pressing direction. A wide range of Al content (0–5 wt%) is used by adding Al powder to nanosized ZnO powder. Thermoelectric (TE) measurements are done on the FGM and individual layers made separately. X-ray diffraction (XRD) shows two phases, ZnO and spinel phase (ZnAl 2 O 4 ). High-resolution Raman spectroscopy reveals doped ZnO and spinel phase (ZnAl 2 O 4 ) spatially and shows segregation in the layer with highest Al content. Electron backscatter diffraction (EBSD) reveals noticeable grain growth with decreasing Al content, and there is a common, random preferred orientation in all layers. The tested properties are used to simulate efficiency curves for a discretely graded, five-layer FGM as well as a homogenous material, where both graded structures provide an opportunity to widen the current density ranges and therefore the temperature range of useful energy conversion.

36 MATERIALS SCIENCE↗

Simulating Heisenberg interactions in the Ising model with strong drive fields

The time evolution of an Ising model with large driving fields over discrete time intervals is shown to be reproduced by an effective XXZ-Heisenberg model at leading order in the inverse field strength. For specific orientations of the drive field, the dynamics of the XXX-Heisenberg model is reproduced. Finally, these approximate equivalences, valid above a critical driving field strength set by dynamical phase transitions in the Ising model, are expected to enable quantum devices that natively evolve qubits according to the Ising model to simulate more complex systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Handbook of operation and maintenance instructions for Model 5208A bit synchronizer unit

Electrical/mechanical characteristics, operation, theory of operation, parts list, and electrical diagrams for the model 5208A bit synchronizer unit are provided. The unit is capable of performing the functions of bit synchronization and signal conditioning on demodulated wave trains containing amplitude, frequency, and phase distortions, as well as additive noise. Operating frequencies in the 1 to 27 MHz bit rate range are accommodated through the use of discrete frequency plug-in subunits. The unit provides 0 and 180 degree clocks and either a single signal hard decision, or a three bit quantized soft decision data output.

Source record↗

Comparison of Acoustic Models and Trajectory Generation Methods for an Acoustically-Aware Aircraft

Motivation - Noise management is one of the major barriers to Urban Air Mobility - Approaches to noise mitigation (non-exhaustive) - Vehicle configuration - Directivity control via propeller phase synchronization - Trajectory optimization Objective - Create framework for trajectory generation integrating location-based acoustic metrics and vehicle performance limitations - Multiple trajectory optimization methods and acoustic noise models - Mission-relevant constraints - Mission duration, airspace restrictions, ... - Vehicle dynamic constraints - Aircraft structural limitations, min/max airspeed, ... - Vehicle separation/obstacle avoidance - Acoustic constraints at a number of discrete observer locations

Kasey A. Ackerman↗

Comparison of Acoustic Models and Trajectory Generation Methods for an Acoustically-Aware Aircraft

Motivation - Noise management is one of the major barriers to Urban Air Mobility - Approaches to noise mitigation (non-exhaustive) - Vehicle configuration - Directivity control via propeller phase synchronization - Trajectory optimization Objective - Create framework for trajectory generation integrating location-based acoustic metrics and vehicle performance limitations - Multiple trajectory optimization methods and acoustic noise models - Mission-relevant constraints - Mission duration, airspace restrictions, ... - Vehicle dynamic constraints - Aircraft structural limitations, min/max airspeed, ... - Vehicle separation/obstacle avoidance - Acoustic constraints at a number of discrete observer locations

Kasey A Ackerman↗

Cold-Flow Testing of a Proposed Integrated Center-Body Diffuser/Steam Blocker Concept for Plum Brook Station's B-2 Test Facility

The center-body diffuser (CBD) steam blocker (SB) system is a concept that incorporates a set of secondary drive nozzles into the envelope of a CBD, such that both nozzle systems (i.e., the rocket engine and the steam blocking nozzles) utilize the same supersonic diffuser, and will operate either singularly or concurrently. In this manner, the SB performs as an exhaust system stage when the rocket engine is not operating, and virtually eliminates discharge flow on rocket engine shutdown. A 2.25-percent scale model of a proposed SB integrated into a diffuser for the Plum Brook B-2 facility was constructed and cold-flow tested for the purpose of evaluating performance characteristics of various design options. These specific design options addressed secondary drive nozzle design (method of steam injection), secondary drive nozzle location relative to CBD throat, and center-body throat length to diameter (L/D) ratios. The objective of the test program is to identify the desired configuration to carry forward should the next phase of design proceed. The tested scale model can provide data for various pressure ratios; however, its design is based on a proposed B-2 spray chamber (SC) operating pressure of 4.0 psia and a steam supply pressure of 165 psia. Evaluation of the test data acquired during these tests indicate that either the discrete axial or annular nozzle configuration integrated into a CBD, with an annular throat length of 1.5 L/D at the nominal injection position, would be suitable to carry forward from the SB's perspective. Selection between these two then becomes more a function of constructability and implementation than performance. L/D also has some flexibility, and final L/D selection can be a function of constructability issues within a limited range.

Edwards, Daryl A.↗

Spectrally-Invariant Approximation Within Atmospheric Radiative Transfer

Certain algebraic combinations of single scattering albedo and solar radiation reflected from, or transmitted through, vegetation canopies do not vary with wavelength. These "spectrally invariant relationships" are the consequence of wavelength independence of the extinction coefficient and scattering phase function in vegetation. In general, this wavelength independence does not hold in the atmosphere, but in cloud~dominated atmospheres the total extinction and total scattering phase function vary only weakly with wavelength. This paper identifies the atmospheric conditions under which the spectrally invariant approximation can accurately describe the extinction. and scattering properties of cloudy atmospheres. The validity of the assumptions and the accuracy of the approximation are tested with ID radiative transfer calculations using publicly available radiative transfer models: Discrete Ordinate Radiative Transfer (DISORT) and Santa Barbara DISORT Atmospheric Radiative Transfer (SBDART). It is shown for cloudy atmospheres with cloud optical depth above 3, and for spectral intervals that exclude strong water vapor absorption, that the spectrally invariant relationships found in vegetation canopy radiative transfer are valid to better than 5%. The physics behind this phenomenon, its mathematical basis, and possible applications to remote sensing and climate are discussed.

Marshak, A.↗