Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing continuum”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Galerkin approximations for dissipative magnetohydrodynamics

A Galerkin approximation scheme is proposed for voltage-driven, dissipative magnetohydrodynamics. The trial functions are exact eigenfunctions of the linearized continuum equations and represent helical deformations of the axisymmetric, zero-flow, driven steady state. The lowest nontrivial truncation is explored: one axisymmetric trial function and one helical trial function each for the magnetic and velocity fields. The system resembles the Lorenz approximation to Benard convection, but in the region of believed applicability, its dynamical behavior is rather different, including relaxation to a helically deformed state similar to those that have emerged in the much higher resolution computations of Dahlburg et al.

Chen, Hudong↗
NASA NTRS1991

Infrared cooling rate calculations in operational general circulation models - Comparisons with benchmark computations

The performance of several parameterized models is described with respect to numerical prediction and climate research at GFDL, NCAR, and GISS. The radiation codes of the models were compared to benchmark calculations and other codes for the intercomparison of radiation codes in climate models (ICRCCM). Cooling rates and fluxes calculated from the models are examined in terms of their application to established general circulation models (GCMs) from the three research institutions. The newest radiation parameterization techniques show the most significant agreement with the benchmark line-by-line (LBL) results. The LBL cooling rates correspond to cooling rate profiles from the models, but the parameterization of the water vapor continuum demonstrates uncertain results. These uncertainties affect the understanding of some lower tropospheric cooling, and therefore more accurate parameterization of the water vapor continuum, as well as the weaker absorption bands of CO2 and O3 is recommended.

Kiehl, J. T.↗
DOE OSTI2024

Fuzzy gauge theory for quantum computers

Continuous gauge theories, because of their bosonic degrees of freedom, have an infinite-dimensional local Hilbert space. Encoding these degrees of freedom on qubit-based hardware demands some sort of “qubitization” scheme, where one approximates the behavior of a theory while using only finitely many degrees of freedom. We propose a novel qubitization strategy for gauge theories, called “fuzzy gauge theory,” building on the success of the fuzzy σ -model in earlier work. We provide arguments that the fuzzy gauge theory lies in the same universality class as regular gauge theory, in which case its use would obviate the need of any further limit besides the usual spatial continuum limit. Furthermore, we demonstrate that these models are relatively resource-efficient for quantum simulations. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2023

Fierro Version 2.x

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three-dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2023

Fierro

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2021

GDSA PFLOTRAN Development (FY2021)

The Spent Fuel & Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is to develop a disposal system modeling and analysis capability for evaluating disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2021 advances of the PFLOTRAN Development group of the SFWST Campaign. The mission of this group is to develop a geologic disposal system modeling capability for nuclear waste that can be used to probabilistically assess the performance of generic disposal concepts. In FY 2021, development proceeded along three main thrusts: software infrastructure, code performance, and process model advancement. Software infrastructure improvements included implementing an Agile software development framework and making improvements to the QA Test Suite. Code performance improvements included development of advanced linear and nonlinear solvers as well as design of flexible smoothing algorithms for capillary pressure functions. Process modeling advancements included the addition of flexible thermal conductivity function definitions and refinement of multi-continuum reactive transport to support Sandia’s participation in DECOVALEX

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1972

Electrical conductivity and velocity of highly ionized plasma flows - Theory and experiment.

Use of an immersible, three-coil, magnetic-induction probe, previously tested in a low-density supersonic argon jet, to measure electrical conductivity and velocity profiles of a highly ionized high-density nitrogen jet in the continuum flow regime where effects due to probe bow shocks and boundary layers might not be negligible. Measured centerline values of electrical conductivity and velocity were compared with predictions based on a theoretical analysis previously developed to study the gas as it expanded adiabatically and inviscidly from an equilibrium sonic state to the nozzle exit. The resulting numerical exit plane values for electron density and electron temperature were then substituted into the Spitzer-Haerm conductivity formula to compute a theoretical conductivity which agreed within 40% of the measured conductivity, while the calculated and experimental velocity values differed by as much as 50%. The lack of agreement was attributed to the possible use of invalid assumptions and boundary conditions in the computer analysis or to the unknown effects of shocks on the probe data.

Vendell, E. W.↗
NASA NTRS1994

X-ray bumps, iron K-alpha lines, and X-ray suppression by obscuring tori in Seyfert galaxies

We investigate the X-ray spectral properties of unobscured type 1 and obscured type 2 Seyferts as predicted by the unified Seyfert scheme. We consider the reprocessing of X-ray photons by photoelectric absorption, iron fluorescence, and Compton downscattering in the obscuring tori surrounding these active nuclei, and compute by Monte Carlo methods the reprocessed spectra as a function of the viewing angle. Depending on the optical depth and shape of the torus, and on the viewing angle, the X-ray flux can be suppressed by substantial factors when our line of sight is obscured. We show that an immediate consequence of the existence of an obscuring thick torus is the production in the spectra of type 1 Seyfert galaxies of a bump in the continuum above 10-20 keV and an Fe K-alpha line with significant equivalent width. In those type 2 Seyferts for which the hard X-ray spectrum has been substantially suppressed, the equivalent width of the Fe K-alpha line in the transmitted spectrum can be very large.

Krolik, Julian H.↗
DOE OSTI2023

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗
DOE OSTI2022

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗
NASA NTRS1992

Aerodynamic characteristics of a hypersonic viscous optimized waverider at high altitudes

The present paper addresses the applicability of the basic concept of waveriding at high altitudes, and the extent to which the large viscous forces degrade the aerodynamic performance of waveriders. The waverider under consideration was designed using a continuum flow methodology. It is shown that the lift-to-drag ratio of high-altitude/high-Knudsen-number waveriders can be expected to be significantly lower than their low altitude/low Knudsen number counterparts. The aerodynamic performance of a representative waverider which was optimized for a 90-km, Mach-25 application is studied for altitudes ranging from 97 km to 145 km and incidence angles of 0 to 30 deg. Typical values of the lift-to-drag ratio were computed to be in the range of 0 to 0.3. Friction forces are mostly responsible for this poor performance. Friction forces account for more than 93 percent of the drag and significantly reduce lift.

Rault, Didier F. G.↗
DOE OSTI2023

Geometry of conformal manifolds and the inversion formula

Families of conformal field theories are naturally endowed with a Riemannian geometry which is locally encoded by correlation functions of exactly marginal operators. We show that the curvature of such conformal manifolds can be computed using Euclidean and Lorentzian inversion formulae, which combine the operator content of the conformal field theory into an analytic function. Analogously, operators of fixed dimension define bundles over the conformal manifold whose curvatures can also be computed using inversion formulae. These results relate curvatures to integrated four-point correlation functions which are sensitive only to the behavior of the theory at separated points. We apply these inversion formulae to derive convergent sum rules expressing the curvature in terms of the spectrum of local operators and their three-point function coefficients. We further show that the curvature can smoothly diverge only if a conserved current appears in the spectrum, or if the theory develops a continuum. We verify our results explicitly in 2d examples. In particular, for 2d (2,2) superconformal field theories we derive a lower bound on the scalar curvature, which is saturated by free theories when the central charge is a multiple of three.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The holar wind

String theory in AdS 3 with purely NS-NS fluxes and vanishing RR moduli has a continuum of winding string excitations in radial plane wave states. BTZ black holes can emit such strings, which then flow out toward the AdS 3 boundary as a stream of massive quanta, and form a black hole analogue of the solar wind. The winding string sector thus provides a decay channel for the black hole to evaporate without having either to couple the system to an external reservoir or to match the AdS 3 throat onto an asymptotically flat region. We compute the emission amplitude of this “holar wind” in the semi-classical approximation, and consider the associated version of the black hole information paradox.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

On the accuracy of modelling the dynamics of large space structures

Proposed space missions will require large scale, light weight, space based structural systems. Large space structure technology (LSST) systems will have to accommodate (among others): ocean data systems; electronic mail systems; large multibeam antenna systems; and, space based solar power systems. The structures are to be delivered into orbit by the space shuttle. Because of their inherent size, modelling techniques and scaling algorithms must be developed so that system performance can be predicted accurately prior to launch and assembly. When the size and weight-to-area ratio of proposed LSST systems dictate that the entire system be considered flexible, there are two basic modeling methods which can be used. The first is a continuum approach, a mathematical formulation for predicting the motion of a general orbiting flexible body, in which elastic deformations are considered small compared with characteristic body dimensions. This approach is based on an a priori knowledge of the frequencies and shape functions of all modes included within the system model. Alternatively, finite element techniques can be used to model the entire structure as a system of lumped masses connected by a series of (restoring) springs and possibly dampers. In addition, a computational algorithm was developed to evaluate the coefficients of the various coupling terms in the equations of motion as applied to the finite element model of the Hoop/Column.

Diarra, C. M.↗
NASA NTRS1985

On the accuracy of modelling the dynamics of large space structures

Proposed space missions will require large scale, light weight, space based structural systems. Large space structure technology (LSST) systems will have to accommodate (among others): ocean data systems; electronic mail systems; large multibeam antenna systems; and, space based solar power systems. The structures are to be delivered into orbit by the Space Shuttle. Because of their inherent size, modelling techniques and scaling algorithms must be developed so that system performance can be predicted accurately prior to launch and assembly. When the size and weight-to-area ratio of proposed LSST systems dictate that the entire system be considered flexible, there are two basic modelling methods which can be used. The first is a continuum approach, a mathematical formulation for predicting the motion of a general orbiting flexible body, in which elastic deformations are considered small compared with characteristic body dimensions. This approach is based on an a priori knowledge of the frequencies and shape functions of all modes included within the system model. Alternatively, finite element techniques can be used to model the entire structure as a system of lumped masses connected by a series of (restoring) springs and possibly dampers. In addition, a computational algorithm was developed to evaluate the coefficients of the various coupling terms in the equations of motion as applied to the finite element model of the Hoop/Column.

Diarra, C. M.↗
NASA NTRSDate not supplied

DRF: A Software Architecture for a Data Marketplace to Support Advanced Air Mobility

Advanced Air Mobility is a new aviation vision, where unmanned aerial systems will trans- port passengers and cargo across urban and rural areas. Critical to the realization of this vision is the development of a digital marketplace, which allows service providers and consumers operating in the airspace ecosystem to securely exchange data and reasoning insights. In this paper, we present the architecture of a decentralized data marketplace that connects data and reasoning service providers to vehicles and other service consumers along the cloud-to-edge continuum. We also present two example use cases to demonstrate the value of our approach.

Autonomy↗
DOE OSTI2022

Mesoflow: An Open-Source Reacting Flow Solver for Catalysis at Mesoscale

We present the capabilities and software performance metrics of our open-source continuum solver for catalysis, Mesoflow, developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex catalyst surface morphologies directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous catalyst interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPUs) compared to a single processor for representative problem sizes (2 million cell mesh). We will also present a brief introduction on how to build and use this software for application problems pertaining to catalytic upgrading and gas transport within porous catalyst particles.

adaptive meshing↗