Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical methods and 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 37 records · Page 2

Large eddy simulation - A predictive approach to turbulent flow computation

The filtering approach (Leonard, 1974) to large eddy simulation is described. Large eddy simulation is then surveyed with reference to subgrid scale models, numerical methods, channel flow, and free shear flows. It is concluded that the premise of LES (i.e., that the small eddies are much easier to model than the large ones) seems to be borne out by computational experience to date. The greatest obstacle in LES at present is the lack of computer power which effectively denies the ability to simulate all the significant eddies in flows of interest.

Ferziger, J. H.↗

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Combustion hot section technology

The Turbine Engine Hot Section Technology Combustion Program is briefly described. The overall objective of the project is to develop and verify improved and more accurate analysis methods for increasing the ability to design with confidence the combustion system for advanced aircraft turbine engines. The approach is to first assess and evaluate existing combustor aerothermal analysis models by means of a contracted effort initiated during FY-82. The program also includes both analytical and experimental research efforts in the areas of aerothermal modeling and liner cyclic life. It is expected that the combustor model development effort will generate improved understanding in the areas of high pressure flame radiation characteristics, model numerical methods and solution schemes, complex geometrical boundary conditions, fuel spray - flow field interactions, combustion kinetics, flow and mixing of dilution jets, turbulence and heat transfer, and soot and carbon formation.

Ercegovic, D. B.↗

Performance on HPC Platforms Is Possible Without C++

Computing at large scales has become extremely challenging due to increasing heterogeneity in both hardware and software. More and more scientific workflows must tackle a range of scales and use machine learning and AI intertwined with more traditional numerical modeling methods, placing more demands on computational platforms. These constraints indicate a need to fundamentally rethink the way computational science is done and the tools that are needed to enable these complex workflows. The current set of C++-based solutions may not suffice, and relying exclusively upon C++ may not be the best option, especially because several newer languages and boutique solutions offer more robust design features to tackle the challenges of heterogeneity. In June 2023, we held a mini symposium that explored the use of newer languages and heterogeneity solutions that are not tied to C++ and that offer options beyond template metaprogramming and Parallel. For for performance and portability. In conclusion, we describe some of the presentations and discussion from the mini symposium in this article.

97 MATHEMATICS AND COMPUTING↗

SAM User’s Guide

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory for advanced non-LWR safety analysis. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for Sodium-cooled Fast Reactors (SFR), Lead-cooled Fast Reactors (LFR), and Molten Salt Reactors (MSR) or Fluoride-cooled High-temperature Reactors (FHR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage the modern advanced software environments and numerical methods. This document provides a user’s guide, which will help users understand the input description and core capabilities of the SAM code. A brief overview of the code is presented, as well as how to obtain and run it. The input syntax for various parts of the code is provided. Additionally, a number of example problems, starting with simple unit component problems to problems with increasing complexity, are provided. Because the code is still under active development, this SAM User’s Guide will evolve with periodic updates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SAM User's Guide

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory for advanced non-LWR safety analysis. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for Sodium-cooled Fast Reactors (SFR), Lead-cooled Fast Reactors (LFR), and Molten Salt Reactors (MSR) or Fluoride-cooled High-temperature Reactors (FHR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage the modern advanced software environments and numerical methods. This document provides a user’s guide, which will help users understand the input description and core capabilities of the SAM code. A brief overview of the code is presented, as well as how to obtain and run it. The input syntax for various parts of the code is provided. Additionally, a number of example problems, starting with simple unit component problems to problems with increasing complexity, are provided. Because the code is still under active development, this SAM User’s Guide will evolve with periodic updates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Heat Transfer and Fluid Mechanics Institute, Meeting, 25th, University of California, Davis, Calif., June 21-23, 1976, Proceedings

Major research advances in heat transfer and fluid dynamics are outlined, with particular reference to relevant energy problems. Of significant importance are such topics as synthetic fuels in combustion, turbulence models, combustion modeling, numerical methods for interacting boundary layers, and light-scattering diagnostics for gases. The discussion covers thermal convection, two-phase flow and boiling heat transfer, turbulent flows, combustion, and aerospace heat transfer problems. Other areas discussed include compressible flows, fluid mechanics and drag, and heat exchangers. Featured topics comprise heat and salt transfer in double-diffusive systems, limits of boiling heat transfer in a liquid-filled enclosure, investigation of buoyancy-induced flow stratification in a cylindrical plenum, and digital algorithms for dynamic analysis of a heat exchanger. Individual items are announced in this issue.

Mckillop, A. A.↗

NASA Agricultural Aircraft Research Program in the Langley Vortex Research Facility and the Langley Full Scale Wind Tunnel

The current status of aerial applications technology research at the Langley's Vortex Research Facility and Full-Scale Wind Tunnel is reviewed. Efforts have been directed mainly toward developing and validating the required experimental and theoretical research tools. A capability to simulate aerial dispersal of materials from agricultural airplanes with small-scale airplane models, numerical methods, and dynamically scaled test particles was demonstrated. Tests on wake modification concepts have proved the feasibility of tailoring wake properties aerodynamically to produce favorable changes in deposition and to provide drift control. An aerodynamic evaluation of the Thrush Commander 800 agricultural airplane with various dispersal systems installed is described. A number of modifications intended to provide system improvement to both airplane and dispersal system are examined, and a technique for documenting near-field spray characteristics is evaluated.

Jordan, F. L., Jr.↗

Higher-level simulations of turbulent flows

The fundamentals of large eddy simulation are considered and the approaches to it are compared. Subgrid scale models and the development of models for the Reynolds-averaged equations are discussed as well as the use of full simulation in testing these models. Numerical methods used in simulating large eddies, the simulation of homogeneous flows, and results from full and large scale eddy simulations of such flows are examined. Free shear flows are considered with emphasis on the mixing layer and wake simulation. Wall-bounded flow (channel flow) and recent work on the boundary layer are also discussed. Applications of large eddy simulation and full simulation in meteorological and environmental contexts are included along with a look at the direction in which work is proceeding and what can be expected from higher-level simulation in the future.

Ferziger, J. H.↗

Higher level simulations of turbulent flows

The fundamentals of large eddy simulation are considered and various approaches to this simulation are compared. The subgrid model required by large eddy simulation is discussed as well as the use of this type of simulation in the development of models for the Reynolds-averaged equations and the application of direct simulation to the testing of both subgrid scale and Reynolds-averaged models. Numerical methods used in large eddy and direct simulation are described with emphasis on special purpose methods. Topics covered include the simulation of homogeneous flow, free shear flows, the mixing layer, wakes, and wall-bounded flows including channel flow and the boundary layer. Applications of large eddy simulation in the laboratory as in meteorological and other environmental flows are examined. Directions in which the work is proceeding and what can be expected from higher levels simulated on are examined.

Ferziger, J. H.↗

Nonlinear Constitutive Relations for High Temperature Application, 1984

Nonlinear constitutive relations for high temperature applications were discussed. The state of the art in nonlinear constitutive modeling of high temperature materials was reviewed and the need for future research and development efforts in this area was identified. Considerable research efforts are urgently needed in the development of nonlinear constitutive relations for high temperature applications prompted by recent advances in high temperature materials technology and new demands on material and component performance. Topics discussed include: constitutive modeling, numerical methods, material testing, and structural applications.

Source record↗

Opportunities and Challenges of Computational Electromagnetics Methods for Superconducting Circuit Quantum Device Modeling: A Practical Review

High-fidelity numerical methods that model the physical layout of a device are essential for the design of many technologies. For methods that characterize electromagnetic effects, these numerical methods are referred to as computational electromagnetics (CEM) methods. Although the CEM research field is mature, emerging applications can still stress the capabilities of the techniques in use today. The design of superconducting circuit quantum devices falls in this category due to the unconventional material properties and important features of the devices covering nanometer to centimeter scales. Such multiscale devices can stress the fundamental properties of CEM tools which can lead to an increase in simulation times, a loss in accuracy, or even cause no solution to be reliably found. While these challenges are being investigated by CEM researchers, knowledge about them is limited in the broader community of users of these CEM tools. This review is meant to serve as a practical introduction to the fundamental aspects of the major CEM techniques that a researcher may need to choose between to model a device, as well as provide insight into what steps they may take to alleviate some of their challenges. Our focus is on highlighting the main concepts without rigorously deriving all the details, which can be found in many textbooks and articles. After covering the fundamentals, we discuss more advanced topics related to the challenges of modeling multiscale devices with specific examples from superconducting circuit quantum devices. We conclude with a discussion on future research directions that will be valuable for improving the ability to successfully design increasingly more sophisticated superconducting circuit quantum devices. Although our focus and examples are taken from this area, researchers from other fields will still benefit from the details discussed here.

Elkin, Samuel T. [Virginia Tech., Blacksburg]↗

Magnetohydrodynamic (MHD) modelling of solar active phenomena via numerical methods

Numerical ideal MHD models for the study of solar active phenomena are summarized. Particular attention is given to the following physical phenomena: (1) local heating of a coronal loop in an isothermal and stratified atmosphere, and (2) the coronal dynamic responses due to magnetic field movement. The results suggest that local heating of a magnetic loop will lead to the enhancement of the density of the neighboring loops through MHD wave compression. It is noted that field lines can be pinched off and may form a self-contained magnetized plasma blob that may move outward into interplanetary space.

Wu, S. T.↗

Analysis of the distribution of pitch angles in model galactic disks - Numerical methods and algorithms

An automated mathematical method capable of successfully isolating the many different features in prototype and observed spiral galaxies and of accurately measuring the pitch angles and lengths of these individual features is developed. The method is applied to analyze the evolution of specific features in a prototype galaxy exhibiting flocculent spiral structure. The mathematical-computational method was separated into two components. Initially, the galaxy was partitioned into dense regions constituting features using two different methods. The results obtained using these two partitioning algorithms were very similar, from which it is inferred that no numerical biasing was evident and that capturing of the features was consistent. Standard least-squares methods underestimated the true slope of the cloud distribution and were incapable of approximating an orientation of 45 deg. The problems were overcome by introducing a superior fit least-squares method, developed with the intention of calculating true orientation rather than a regression line.

Russell, William S.↗

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗

A numerical study of the effects of wind tunnel wall proximity on an airfoil model

A procedure was developed for modeling wind tunnel flows using computational fluid dynamics. Using this method, a numerical study was undertaken to explore the effects of solid wind tunnel wall proximity and Reynolds number on a two-dimensional airfoil model at low speed. Wind tunnel walls are located at varying wind tunnel height to airfoil chord ratios and the results are compared with freestream flow in the absence of wind tunnel walls. Discrepancies between the constrained and unconstrained flows can be attributed to the presence of the walls. Results are for a Mach Number of 0.25 at angles of attack through stall. A typical wind tunnel Reynolds number of 1,200,000 and full-scale flight Reynolds number of 6,000,000 were investigated. At this low Mach number, wind tunnel wall corrections to Mach number and angle of attack are supported. Reynolds number effects are seen to be a consideration in wind tunnel testing and wall interference correction methods. An unstructured grid Navier-Stokes code is used with a Baldwin-Lomax turbulence model. The numerical method is described since unstructured flow solvers present several difficulties and fundamental differences from structured grid codes, especially in the area of turbulence modeling and grid generation.

Potsdam, Mark↗

Coupled Neutron Transport for HZETRN

Exposure estimates inside space vehicles, surface habitats, and high altitude aircrafts exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETC-HEDS, FLUKA, and MCNPX, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light particle transport procedure adversely affects low energy light ion fluence spectra and exposure quantities. A first order correction is presented to resolve the problem, and it is shown to be both accurate and efficient.

Slaba, Tony C.↗

NETL Plastic Pipes Project - June 2022 Progress Report

Plastic or composite pipelines have been the bane of the utility locating industry because conduits are neither conductive nor magnetic, which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipeline where the resistive cover allows for adequate penetration of radar energy. However, GPR cannot be used in areas where the soil cover is conductive due to significant clay and/or salt content. This study takes a comprehensive look at near-surface geophysical methods that potentially are useful for locating buried plastic/composite pipelines, either singly or in combination with other geophysical methods. Specifically, this modeling study uses computational numerical methods to forward model the response of GPR, resistivity, seismic, and gravity gradiometry methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft to 4 ft; (3) various degrees in contrast in physical properties (i.e., permittivity, elasticity, electrical resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar frequencies, seismic source frequency, electrode spacing).

47 OTHER INSTRUMENTATION↗