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 163 records · Page 9

Generalized Fluid System Simulation Program, Version 5.0-Educational

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the point, drag and click method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids and 21 different resistance/source options are provided for modeling momentum sources or sinks in the branches. This Technical Memorandum illustrates the application and verification of the code through 12 demonstrated example problems.

Majumdar, A. K.↗

Generalized Fluid System Simulation Program, Version 5.0-Educational. Supplemental Information for NASA/TM-2011-216470

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the point, drag and click method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids and 21 different resistance/source options are provided for modeling momentum sources or sinks in the branches. This Technical Memorandum illustrates the application and verification of the code through 12 demonstrated example problems. This supplement gives the input and output data files for the examples.

Majumdar, A. K.↗

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗

Investigation for improving Global Positioning System (GPS) orbits using a discrete sequential estimator and stochastic models of selected physical processes

GEODYNII is a conventional batch least-squares differential corrector computer program with deterministic models of the physical environment. Conventional algorithms were used to process differenced phase and pseudorange data to determine eight-day Global Positioning system (GPS) orbits with several meter accuracy. However, random physical processes drive the errors whose magnitudes prevent improving the GPS orbit accuracy. To improve the orbit accuracy, these random processes should be modeled stochastically. The conventional batch least-squares algorithm cannot accommodate stochastic models, only a stochastic estimation algorithm is suitable, such as a sequential filter/smoother. Also, GEODYNII cannot currently model the correlation among data values. Differenced pseudorange, and especially differenced phase, are precise data types that can be used to improve the GPS orbit precision. To overcome these limitations and improve the accuracy of GPS orbits computed using GEODYNII, we proposed to develop a sequential stochastic filter/smoother processor by using GEODYNII as a type of trajectory preprocessor. Our proposed processor is now completed. It contains a correlated double difference range processing capability, first order Gauss Markov models for the solar radiation pressure scale coefficient and y-bias acceleration, and a random walk model for the tropospheric refraction correction. The development approach was to interface the standard GEODYNII output files (measurement partials and variationals) with software modules containing the stochastic estimator, the stochastic models, and a double differenced phase range processing routine. Thus, no modifications to the original GEODYNII software were required. A schematic of the development is shown. The observational data are edited in the preprocessor and the data are passed to GEODYNII as one of its standard data types. A reference orbit is determined using GEODYNII as a batch least-squares processor and the GEODYNII measurement partial (FTN90) and variational (FTN80, V-matrix) files are generated. These two files along with a control statement file and a satellite identification and mass file are passed to the filter/smoother to estimate time-varying parameter states at each epoch, improved satellite initial elements, and improved estimates of constant parameters.

Goad, Clyde C.↗

EUTERPE: A global gyrokinetic code for stellarator geometry

The current state of the EUTERPE code is described with emphasis on the implemented models and their numerical implementation. The code solves the multi-species electromagnetic gyrokinetic equations in the full volume of a three-dimensional domain. Noise reduction of the particle-in-cell method is achieved by using a δf-method and Fourier filters. The field equations are discretized with B-splines and the resulting system of equations is solved iteratively. For linear simulations a phase-factor transformation is applied in order to strongly reduce the necessary grid resolution. Apart from the full gyrokinetic model, other numerically less expensive hybrid models are also implemented. They are mainly tailored for comparison with fluid theory and for studying the interaction of the bulk plasma with fast particles. The code is parallelized for CPUs by particle and domain decomposition. Good scalability up to several thousand nodes is demonstrated.

97 MATHEMATICS AND COMPUTING↗

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗

Performance of a finite phase state bit-synchronization loop with and without sequential filters

A discrete phase state bit-synchronization loop is proposed for synchronization and detection of a binary nonreturn-to-zero (NRZ) process. The general loop model is developed for 11 and 17 phase states with optional sequential loop filter. Loop performance is considered in terms of the phase error density function, rms phase error, and mean time to bit slippage. It is shown that satisfactory loop performance can be achieved at -3 dB SNR through use of sequential loop filters.

Ransom, J. J.↗

A two-phase restricted equilibrium model for combustion of metalized solid propellants

An Eulerian-Lagrangian two-phase approach was adopted to model the multi-phase reacting internal flow in a solid rocket with a metalized propellant. An Eulerian description was used to analyze the motion of the continuous phase which includes the gas as well as the small (micron-sized) particulates, while a Lagrangian description is used for the analysis of the discrete phase which consists of the larger particulates in the motor chamber. The particulates consist of Al and Al2O3 such that the particulate composition is 100 percent Al at injection from the propellant surface with Al2O3 fraction increasing due to combustion along the particle trajectory. An empirical model is used to compute the combustion rate for agglomerates while the continuous phase chemistry is treated using chemical equilibrium. The computer code was used to simulate the reacting flow in a solid rocket motor with an AP/HTPB/Al propellant. The computed results show the existence of an extended combustion zone in the chamber rather than a thin reaction region. The presence of the extended combustion zone results in the chamber flow field and chemical being far from isothermal (as would be predicted by a surface combustion assumption). The temperature in the chamber increases from about 2600 K at the propellant surface to about 3350 K in the core. Similarly the chemical composition and the density of the propellant gas also show spatially non-uniform distribution in the chamber. The analysis developed under the present effort provides a more sophisticated tool for solid rocket internal flow predictions than is presently available, and can be useful in studying apparent anomalies and improving the simple correlations currently in use. The code can be used in the analysis of combustion efficiency, thermal load in the internal insulation, plume radiation, etc.

Sabnis, J. S.↗

The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation

We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.

42 ENGINEERING↗

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependant flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. This Technical Memorandum illustrates the application and verification of the code through 25 demonstrated example problems.

Majumdar, A. K.↗

Advanced ISDN satellite designs and experiments

The research performed by GTE Government Systems and the University of Colorado in support of the NASA Satellite Communications Applications Research (SCAR) Program is summarized. Two levels of research were undertaken. The first dealt with providing interim services Integrated Services Digital Network (ISDN) satellite (ISIS) capabilities that accented basic rate ISDN with a ground control similar to that of the Advanced Communications Technology Satellite (ACTS). The ISIS Network Model development represents satellite systems like the ACTS orbiting switch. The ultimate aim is to move these ACTS ground control functions on-board the next generation of ISDN communications satellite to provide full-service ISDN satellite (FSIS) capabilities. The technical and operational parameters for the advanced ISDN communications satellite design are obtainable from the simulation of ISIS and FSIS engineering software models of the major subsystems of the ISDN communications satellite architecture. Discrete event simulation experiments would generate data for analysis against NASA SCAR performance measure and the data obtained from the ISDN satellite terminal adapter hardware (ISTA) experiments, also developed in the program. The Basic and Option 1 phases of the program are also described and include the following: literature search, traffic mode, network model, scenario specifications, performance measures definitions, hardware experiment design, hardware experiment development, simulator design, and simulator development.

Pepin, Gerard R.↗

Pair distribution function analysis of discrete nanomaterials in PDFgui

Pair distribution functions (PDFs) are a leading tool for atomic structure analysis of nanomaterials. However, the most widely used programs for refining atomic structure against PDF data are based on extended crystallographic models, which cannot be applied to discrete, whole nanoparticles. Furthermore, this work describes a straightforward approach to simulate and refine atomistic models of discrete clusters and nanoparticles employing widely used PDF modelling programs such as PDFgui that utilize extended crystallographic models. In this approach, the whole particle to be modelled is contained within an expanded, and otherwise empty, unit cell that is sufficiently large to avoid correlations between atoms in neighbouring unit cells over the r range analyzed. The PDF of the particle is simulated as a composite using two conventional `phases': one that calculates the atom–atom correlations and one that approximates the local number density. This approach is first validated for large nanoparticles that are well modelled by a conventional shape factor model, and then applied to simulate the PDF of discrete particles and low-dimensional materials (graphene and MXene) and to model the experimental PDF data for single-layer FeS nanosheets. A comparison of this approach with the DiffPy-CMI program, which calculates the PDF of discrete species, shows that the composite modelling approach is equally or more accurate. Example input files for implementing this approach within PDFgui and TOPAS, and recommendations for selecting model parameters for reliable application of this refinement strategy, are provided.

36 MATERIALS SCIENCE↗

Signatures of fractionalization in the optical phonons of the hyperhoneycomb Kitaev magnet 𝛽−Li 2 ⁢IrO 3

Here, in this study, we propose that the signatures of spin fractionalization in quantum magnets can be identified through a detailed analysis of the temperature dependence of the asymmetric Fano lineshape of optical phonons overlapping with a continuum of spin excitations. We focus on the hyperhoneycomb magnet 𝛽−Li 2 ⁢IrO 3 , a promising candidate for being in proximity to a three-dimensional Kitaev quantum spin liquid. The Raman response in 𝛽−Li 2 ⁢IrO 3 notably displays a distinctive asymmetric Fano lineshape in the 24 meV Raman-active optical phonon. This asymmetry arises from the interaction between the discrete phonon mode and the spin excitation continuum, which could be fractionalized if the material is indeed near a quantum spin-liquid phase. Our theoretical model considers the coupling of this optical phonon to Majorana fermions in the Kitaev model on the hyperhoneycomb lattice. Our findings reveal that the temperature-dependent Fano lineshape is consistent with the fractionalization of spins into Majorana fermions and ℤ ⁢2 fluxes.

Kitaev model↗

Combustion chamber analysis code

A three-dimensional, time dependent, Favre averaged, finite volume Navier-Stokes code has been developed to model compressible and incompressible flows (with and without chemical reactions) in liquid rocket engines. The code has a non-staggered formulation with generalized body-fitted-coordinates (BFC) capability. Higher order differencing methodologies such as MUSCL and Osher-Chakravarthy schemes are available. Turbulent flows can be modeled using any of the five turbulent models present in the code. A two-phase, two-liquid, Lagrangian spray model has been incorporated into the code. Chemical equilibrium and finite rate reaction models are available to model chemically reacting flows. The discrete ordinate method is used to model effects of thermal radiation. The code has been validated extensively against benchmark experimental data and has been applied to model flows in several propulsion system components of the SSME and the STME.

Przekwas, A. J.↗

Beam Propagation Through Atmospheric Turbulence Using an Altitude-Dependent Structure Profile with Non-Uniformly Distributed Phase Screens

For free-space optical communication systems, numerical wave optics simulations provide a useful technique for modeling turbulence-induced beam degradation when the analytical theory is insufficient for characterizing the atmospheric channel. Motivated by such applications we use a split-step method modeling the turbulence as a series of random phase screens using the Hufnagel-Valley turbulence profile. We employ a space-to-ground case study to examine the irradiance and phase statistics for uniformly and non-uniformly located screens and find better agreement with theory using a non-uniform discretization minimizing the contribution of each screen to the total scintillation. In this poster, we summarize the method and the results of the case study including a comparison to layered models used in astronomical imaging applications.

Chahine, Yousef K.↗

Tools for unbinned unfolding

Machine learning has enabled differential cross section measurements that are not discretized. Going beyond the traditional histogram-based paradigm, these unbinned unfolding methods are rapidly being integrated into experimental workflows. Here, in order to enable widespread adaptation and standardization, we develop methods, benchmarks, and software for unbinned unfolding. For methodology, we demonstrate the utility of boosted decision trees for unfolding with a relatively small number of high-level features. This complements state-of-the-art deep learning models capable of unfolding the full phase space. To benchmark unbinned unfolding methods, we develop an extension of existing dataset to include acceptance effects, a necessary challenge for real measurements. Additionally, we directly compare binned and unbinned methods using discretized inputs for the latter in order to control for the binning itself. Lastly, we have assembled two software packages for the OmniFold unbinned unfolding method that should serve as the starting point for any future analyses using this technique. One package is based on the widely-used RooUnfold framework and the other is a standalone package available through the Python Package Index (PyPI).

47 OTHER INSTRUMENTATION↗

Continuous symmetry breaking in a two-dimensional Rydberg array

Spontaneous symmetry breaking underlies much of our classification of phases of matter and their associated transitions. The nature of the underlying symmetry being broken determines many of the qualitative properties of the phase; this is illustrated by the case of discrete versus continuous symmetry breaking. Indeed, in contrast to the discrete case, the breaking of a continuous symmetry leads to the emergence of gapless Goldstone modes controlling, for instance, the thermodynamic stability of the ordered phase. Here, we realize a two-dimensional dipolar XY model that shows a continuous spin-rotational symmetry using a programmable Rydberg quantum simulator. We demonstrate the adiabatic preparation of correlated low-temperature states of both the XY ferromagnet and the XY antiferromagnet. In the ferromagnetic case, we characterize the presence of a long-range XY order, a feature prohibited in the absence of long-range dipolar interaction. In conclusion, our exploration of the many-body physics of XY interactions complements recent works using the Rydberg-blockade mechanism to realize Ising-type interactions showing discrete spin rotation symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Strategies for simultaneous strengthening and toughening via nanoscopic intracrystalline defects in a biogenic ceramic

While many organisms synthesize robust skeletal composites consisting of spatially discrete organic and mineral (ceramic) phases, the intrinsic mechanical properties of the mineral phases are poorly understood. Using the shell of the marine bivalve Atrina rigida as a model system, and through a combination of multiscale structural and mechanical characterization in conjunction with theoretical and computational modeling, we uncover the underlying mechanical roles of a ubiquitous structural motif in biogenic calcite, their nanoscopic intracrystalline defects. These nanoscopic defects not only suppress the soft yielding of pure calcite through the classical precipitation strengthening mechanism, but also enhance energy dissipation through controlled nano- and micro-fracture, where the defects’ size, geometry, orientation, and distribution facilitate and guide crack initialization and propagation. These nano- and micro-scale cracks are further confined by larger scale intercrystalline organic interfaces, enabling further improved damage tolerance.

36 MATERIALS SCIENCE↗