Search NASA⌕ Search

SEARCH · Search NASA

Results for “CHARACTERISTIC METHOD”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dancoff-based Wigner-Seitz approximation for the subgroup resonance self-shielding in the VERA neutronic simulator MPACT

The MPACT neutronics module of the Virtual Environment for Reactor Analysis (VERA) has used the subgroup method for resonance self-shielding calculation, for which two-dimensional (2D) fixed-source transport calculations are performed using the method of characteristics for resonance energy groups. When considering thermal feedbacks, the subgroup calculation must be performed at each outer iteration. Therefore, the computing time for cross section processing is a significant burden for computational efficiency. The Dancoff-based Wigner-Seitz approximation (DWA) capability has been implemented into MPACT in conjunction with the subgroup method, which has been used in SCALE/XSProc since SCALE version 6.0 and has recently been called an equivalent Dancoff-factor cell (EDC) method. The issue of relatively large reactivity bias in DWA for the gadolinia rods was resolved by introducing multiple Dancoff factors. Benchmark results for the VERA pressurized and boiling water reactor benchmark suites show that the DWA capability would significantly enhance computational efficiency with comparable accuracy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1974

Nickel-cadium batteries for Apollo telescope mount

The operational testing and evaluation program is presented which was conducted on 20-ampere-hour nickel-cadmium (Ni-Cd) batteries for use on the Apollo telescope mount (ATM). The test program was initiated in 1967 to determine if the batteries could meet ATM mission requirements and to determine operating characteristics and methods. The ATM system power and charging power for the Ni-Cd secondary batteries is provided by a solar array during the 58-minute daylight portion of the orbit; during the 36-minute night portion of the orbit, the Ni-Cd secondary batteries will supply ATM system power. The test results reflect battery operating characteristics and parameters relative to simulated ATM orbital test conditions. Maximum voltage, charge requirements, capacity, temperature, and cyclic characteristics are presented.

Kirsch, W. W.↗
DOE OSTI2021

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1990

A review of instability and noise propagation in supersonic flows

Originally analytical and numerical models were to be developed for noise production in supersonic jets, wakes and free shear layers. While the effort was concentrated initially on these aspects, other topics were also pursued, most were of interest to the Jet Noise Group of the Aeroacoustics Branch. An overview is given of subjects reviewed and the investigations that were carried out. A significant effort was devoted to numerically predicting the flow field of a turbulent supersonic wall jet. This information is necessary for computing the pressure in the far field. The wall jet was selected because it represents a generic flow that can be associated with plug nozzle in supersonic engines. It combines the characteristic of a boundary layer with that of a free shear flow. The spatially evolving flow obtained using Dash's code would form the input for the stability analysis program. This analysis would determine the large scale instability wave within the flow. The far field pressure can be computed from the shape of the evolving large scale structure by asymptotic methods. Flow characteristics obtained from a program that analyses the turbulent downstream supersonic flow in a nozzle are described and compared with experimental results.

Daudpota, Q. Isa↗
NASA NTRS2005

Method for Estimating the Sonic-Boom Characteristics of Lifting Canard-Wing Aircraft Concepts

A method for estimating the sonic-boom overpressures from a conceptual aircraft where the lift is carried by both a canard and a wing during supersonic cruise is presented and discussed. Computer codes used for the prediction of the aerodynamic performance of the wing, the canard-wing interference, the nacelle-wing interference, and the sonic-boom overpressures are identified and discussed as the procedures in the method are discussed. A canard-wing supersonic-cruise concept was used as an example to demonstrate the application of the method.

Mack, Robert J.↗
NASA NTRS1973

Analysis of SRM model nozzle calibration test data in support of IA12B, IA12C and IA36 space shuttle launch vehicle aerodynamics tests

Variations of nozzle performance characteristics of the model nozzles used in the Space Shuttle IA12B, IA12C, IA36 power-on launch vehicle test series are shown by comparison between experimental and analytical data. The experimental data are nozzle wall pressure distributions and schlieren photographs of the exhaust plume shapes. The exhaust plume shapes were simulated experimentally with cold flow while the analytical data were generated using a method-of-characteristics solution. Exhaust plume boundaries, boundary shockwave locations and nozzle wall pressure measurements calculated analytically agree favorably with the experimental data from the IA12C and IA36 test series. For the IA12B test series condensation was suspected in the exhaust plumes at the higher pressure ratios required to simulate the prototype plume shapes. Nozzle calibration tests for the series were conducted at pressure ratios where condensation either did not occur or if present did not produce a noticeable effect on the plume shapes. However, at the pressure ratios required in the power-on launch vehicle tests condensation probably occurs and could significantly affect the exhaust plume shapes.

Baker, L. R., Jr.↗
NASA NTRS2012

Method for Predicting the Energy Characteristics of Li-Ion Cells Designed for High Specific Energy

Novel electrode materials with increased specific capacity and voltage performance are critical to the NASA goals for developing Li-ion batteries with increased specific energy and energy density. Although performance metrics of the individual electrodes are critically important, a fundamental understanding of the interactions of electrodes in a full cell is essential to achieving the desired performance, and for establishing meaningful goals for electrode performance in the first place. This paper presents design considerations for matching positive and negative electrodes in a viable design. Methods for predicting cell-level performance, based on laboratory data for individual electrodes, are presented and discussed.

Bennett, William, R.↗
NASA NTRSDate not supplied

Validated Reduced Computational Methods for Realistic RDC Injection Modeling

This paper outlines a simplified computational method without chemistry for single and multi-injector performance for rotating detonation combustion (RDC) environment as well as relevant performance metrics. For non-reactive cases, an unsteady pressure profile from a Method of Characteristics (MOC) simulation was used to simulate an RDC combustion wave travelling around the annulus. In reactive cases, a pre-mixed stoichiometric mixture of hydrogen and air is injected into a 2D unwrapped combustion chamber with a specified injector area. The performance is evaluated based on several performance metrics including unsteady discharge coefficient, average mass flux coming into the combustor, and percentage time obstructed. Finally, the computational method was validated through an experimental RDC test setup with a single element injector in a hydrogen-air detonation driver environment.

Propulsion↗
DOE OSTI2022

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

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

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Laminar Burning Speed Measurements of Hydrogen/Natural Gas Mixtures

Abstract To address the shortage of fossil fuels and the environmental impact of exhaust emissions, research efforts in combustion seek to identify candidates for sustainable fuels. Hydrogen is an example, from which water is produced as the main combustion product. However, when used as neat hydrogen, it is very difficult to control hydrogen combustion due to its high flame speed, pre-ignition propensity, and flashback characteristics. One method to control hydrogen combustion is to blend it with a well-established fuel. Natural gas is a common fuel for use in engines or gas turbines. Natural gas is characterized by slow flame speeds and poor lean-burn ability; thus, engine power is decreased in the lean-burn region. Consequently, natural gas engines are typically operated at stoichiometric conditions to maximize output. Blending hydrogen to natural gas has several advantages: (i) Control of hydrogen’s ignition characteristics, (ii) extension of Natural gas’s lower flammability limit — making it easier to ignite, and (iii) reduced greenhouse gas emissions. In this study, the laminar burning speed (LBS) for neat natural gas and natural gas/hydrogen mixtures at 5 atm and 296 K were measured. Schlieren optical imaging was utilized to validate laminar flame conditions and further investigate the structure of natural gas/hydrogen flames. Results indicate that replacing 50% of natural gas with hydrogen increases the peak laminar burning velocity by a factor of ∼1.6 times. LBS simulations were conducted with NUI 1.1 and UCF NG/H2 mechanisms; Predictions by both mechanisms were found to be satisfactory for neat NG mixtures. However, for 50% NG/ 50% H2 mixture, the true LBS for fuel-lean conditions was under-predicted by both models, while predictions for fuel-rich conditions were satisfactory. Additionally, a literature review and similar experimental analysis will be performed for methane, ammonia, and hydrogen fuel blends to better understand the combustion phenomena of natural gas fuels.

Yovino, Louis↗
DOE OSTI2022

MPACT Theory Manual (V.4.3)

This manual presents the theory underlying the three-dimensional (3D) whole-core, pin-resolved neutron transport calculation methodologies employed in the MPACT code. MPACT’s primary goal is to provide accurate sub-pin power distributions in a computationally efficient manner. To accomplish this, MPACT offers several different transport method options. The 2D/1D method, in which 3D problems are decomposed into an axial stack of radial “slices,” is currently the most commonly used. In this option, two-dimensional (2D) planar solutions are provided by the method of characteristics (MOC), and axial solutions are provided via one-dimensional (1D) approximate diffusion, or P3 solutions. The radial and axial solutions are coupled by (i) axial and radial transverse leakages, and (ii) a global 3D coarse mesh finite difference (CMFD) solve, which provides both acceleration and stability to the solution iteration scheme. The subsequent chapters of this manual present a range of topics, including MOC, CMFD, axial nodal transport solvers, 2D/1D, self-shielding, depletion, thermal-hydraulics, and transient methods. Development of the underlying theory of the 2D/1D method is presented, diagrams are included to highlight important algorithmic flow, and important concepts are discussed as appropriate. This manual is intended to be self-sufficient, but references to published articles and other materials are included for further reading. MPACT is a relatively new code, with new capabilities and many computational methods that did not exist until recently. This manual is an even newer document, with chapters written by several different code contributors working under time constraints. For these reasons, the manual is not yet complete and finalized. Future versions of this manual will address current deficiencies.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1973

Propagation of plane waves in flow through a variable area duct

Duct geometry is an important factor that influences the transmission of sound in a duct. Nonconstant area produces variations in steady flow quantities through the duct, which can cause reflections of acoustic disturbances and the creation of standing wave patterns, and which can attenuate and disperse propagating waves. To study the problem, a second-order-accurate numerical method has been developed. The method of characteristics is used to solve the acoustic equations with the assumption of quasi-one-dimensional flow in the duct. Numerical results are presented over a broad range of frequencies. Low-frequency results are compared with approximate solutions for compact nozzles and high-frequency calculations are compared with results from the short-wave theory. The numerical results are also compared with an analytical solution of the acoustic equations for exponential horn sections with zero mean flow to further establish the validity of the numerical method.

King, L. S.↗
NASA NTRS1976

Numerical solution of periodic transonic flow through a fan stage

A numerical method of solution of the inviscid, compressible, two-dimensional unsteady flow on a blade-to-blade stream surface through a stage (rotor and stator), or a single blade row, of an axial flow compressor or fan is described. A cyclic procedure has been developed for representation of adjacent blade-to-blade passages, which asymptotically achieves the correct phase between all passages of a stage. A shock-capturing finite-difference method is employed in the interior of the passage, and a method-of-characteristics technique is used at the boundaries. The blade slipstreams form two of the passage boundaries, and are treated as moving contact surfaces capable of supporting jumps in entropy and tangential velocity. The Kutta condition is imposed by requiring the slipstreams to originate at the trailing edges, which are assumed to be sharp. Results are presented for several transonic fan rotors, and compared with available experimental data consisting of holographic observations of shock structure and pressure contour maps. A subcritical stator solution is also compared with results from a relaxation method. Finally, a periodic solution for a stage consisting of 44 rotor blades and 46 stator blades is discussed.

Erdos, J. I.↗
NASA NTRS1977

Computation of unsteady transonic flows through rotating and stationary cascades. 1: Method of analysis

A numerical method of solution of the inviscid, compressible, two-dimensional unsteady flow on a blade-to-blade stream surface through a stage (rotor and stator) or a single blade row of an axial flow compressor or fan is described. A cyclic procedure has been developed for representation of adjacent blade-to-blade passages which asymptotically achieves the correct phase between all passages of a stage. A shock-capturing finite difference method is employed in the interior of the passage, and a method of characteristics technique is used at the boundaries. The blade slipstreams form two of the passage boundaries and are treated as moving contact surfaces capable of supporting jumps in entropy and tangential velocity. The Kutta condition is imposed by requiring the slipstreams to originate at the trailing edges, which are assumed to be sharp. Results are presented for several transonic fan rotors and compared with available experimental data, consisting of holographic observations of shock structure and pressure contour maps. A subcritical stator solution is also compared with results from a relaxation method. Finally, a periodic solution for a stage consisting of 44 rotor blades and 46 stator blades is discussed.

Erdos, J. I.↗
NASA NTRS1952

Transonic flow past a wedge profile with detached bow wave

A theoretical study has been made of the aerodynamic characteristics at zero angle of attack of a thin, doubly symmetrical double-wedge profile in the range of supersonic flight speed in which the bow wave is detached. The analysis utilizes the equations of the transonic small-disturbance theory and involves no assumptions beyond those implicit in this theory. The mixed flow about the front half of the profile is calculated by relaxation solution of boundary conditions along the shock polar and sonic line. The purely subsonic flow about the rear of the profile is found by means of the method of characteristics specialized to the transonic small-disturbance theory. Complete calculations were made for four values of the transonic similarity parameter. These were found sufficient to bridge the gap between the previous results of Guderley and Yoshihara at a Mach number of 1 and the results which are readily obtained when the bow wave is attached and the flow is completely supersonic.

Vincenti, Walter G↗
NASA NTRS1974

Attenuation characteristics of nonlinear pressure waves propagating in pipes

A series of experiments was conducted to investigate temporal and spatial velocity distributions of fluid flow in 3-in. open-end pipes of various lengths up to 210 ft, produced by the propagation of nonlinear pressure waves of various intensities. Velocity profiles across each of five sections along the pipes were measured as a function of time with the use of hot-film and hot-wire anemometers for two pressure waves produced by a piston. Peculiar configurations of the velocity profiles across the pipe section were noted, which are uncommon for steady pipe flow. Theoretical consideration was given to this phenomenon of higher velocity near the pipe wall for qualitative confirmation. Experimentally time-dependent velocity distributions along the pipe axis were compared with one-dimensional theoretical results obtained by the method of characteristics with or without diffusion term for the purpose of determining the attenuation characteristics of the nonlinear wave propagation in the pipes.

Shih, C. C.↗
NASA NTRS1972

Initial development of a hypersonic free mixing layer

A preliminary experimental investigation to establish some of the characteristics and further the understanding of the initial development of a turbulent free mixing layer for hypersonic speeds has been conducted. Mean profile data at about 6 inches downstream of the exit of a hypersonic nozzle have been obtained in nitrogen for a nominal Mach number of 19.5, total temperature of about 1670 K and Reynolds number range from about 50,000 to 110,000 per foot and have been compared with profiles upstream of the nozzle exit. Static pressure varied across the shear layer for the present tests. The outer 80 percent of the high-velocity portion of the free shear layer can be calculated by a rotational method of characteristics. However, turbulent mixing is evidently important in the low-velocity region, and effects of eddy viscosity and eddy conductivity should be included in a theoretical analysis.

Harvey, W. D.↗