SEARCH · Engineering Papers

Results for “computational physics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Embedded Model Error Representation for Bayesian Model Calibration

Model error estimation remains one of the key challenges in uncertainty quantification and predictive science. For computational models of complex physical systems, model error, also known as structural error or model inadequacy, is often the largest contributor to the overall predictive uncertainty. This work builds on a recently developed framework of embedded, internal model correction, in order to represent and quantify structural errors, together with model parameters,within a Bayesian inference context. We focus specifically on a Polynomial Chaos representation with additive modification of existing model parameters, enabling a non-intrusive procedure for efficient approximate likelihood construction, model error estimation, and disambiguation of model and data errors’ contributions to predictive uncertainty. The framework is demonstrated on several synthetic examples, as well as on a chemical ignition problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

An Adaptive Mesh Algorithm: Mesh Structure and Generation

The purpose of Adaptive Mesh Refinement is to minimize spatial errors over the computational space not to minimize the number of computational elements. The additional result of the technique is that it may reduce the number of computational elements needed to retain a given level of spatial accuracy. Adaptive mesh refinement is a computational technique used to dynamically select, over a region of space, a set of computational elements designed to minimize spatial error in the computational model of a physical process. The fundamental idea is to increase the mesh resolution in regions where the physical variables are represented by a broad spectrum of modes in k-space, hence increasing the effective global spectral coverage of those physical variables. In addition, the selection of the spatially distributed elements is done dynamically by cyclically adjusting the mesh to follow the spectral evolution of the system. Over the years three types of AMR schemes have evolved; block, patch and locally refined AMR. In block and patch AMR logical blocks of various grid sizes are overlaid to span the physical space of interest, whereas in locally refined AMR no logical blocks are employed but locally nested mesh levels are used to span the physical space. The distinction between block and patch AMR is that in block AMR the original blocks refine and coarsen entirely in time, whereas in patch AMR the patches change location and zone size with time. The type of AMR described herein is a locally refi ned AMR. In the algorithm described, at any point in physical space only one zone exists at whatever level of mesh that is appropriate for that physical location. The dynamic creation of a locally refi ned computational mesh is made practical by a judicious selection of mesh rules. With these rules the mesh is evolved via a mesh potential designed to concentrate the nest mesh in regions where the physics is modally dense, and coarsen zones in regions where the physics is modally sparse.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

Impurity precipitation in atomized particles evidenced by nano x-ray diffraction computed tomography

Performances and physical properties of high technology materials are influenced or even determined by their initial microstructure and by the behavior of impurity phases. Characterizing these impurities and their relations with the surrounding matrix is therefore of primary importance but it unfortunately often requires a destructive approach, with the risk of misinterpreting the observations. The improvement we have done in high resolution X-ray diffraction computed tomography combined with the use of an X-ray nanoprobe allows non-destructive crystallographic description of materials with microscopic heterogeneous microstructure (with a grain size between 10 nm and 10 μm). In this study, the grain localization in a 2D slice of a 20 μm solidified atomized γU-Mo particle is shown and a minority U(C,O) phase (1 wt. %) with sub-micrometer sized grains was characterized inside. Evidence is presented showing that the onset of U(C,O) grain crystallization can be described by a precipitation mechanism since one single U-Mo grain has direct orientation relationship with more than one surrounding U(C,O) grains.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2015

VERAIn

CASL's modeling and simulation technology, the Virtual Environment for Reactor Applications (VERA), incorporates coupled physics and science-based models, state-of-the-art numerical methods, modern computational science, integrated uncertainty quantification (UQ) and validation against data from operating pressurized water reactors (PWRs), single-effect experiments, and integral tests. The computational simulation component of VERA is the VERA Core Simulator (VERA-CS). The core simulator is the specific collection of multi-physics computer codes used to model and deplete a LWR core over multiple cycles. The core simulator has a single common input file that drives all of the different physics codes. The parser code, VERAIn, converts VERA Input into an XML file that is used as input to different VERA codes.

Simunovic, Srdjan↗

DOE OSTI2015

VERAIn

Simunovic, Srdjan↗

NASA NTRS1991

Implementation of a three degree of freedom, motor/brake hybrid force output device for virtual environment control tasks

The advent of high resolution, physical model based computer graphics has left a gap in the design of input/output technology appropriate for interacting with such complex virtual world models. Since virtual worlds consist of physical models, it is appropriate to output the inherent force information necessary for the simulation to the user. The detailed design, construction, and control of a three degree freedom force output joystick will be presented. A novel kinematic design allows all three axes to be uncoupled, so that the system inertia matrix is diagonal. The two planar axes are actuated through an offset gimbal, and the third through a sleeved cable. To compensate for friction and inertia effects, this transmission is controlled by a force feedforward and a closed force feedback proportional loop. Workspace volume is a cone of 512 cubic inches, and the device bandwidth is maximized at 60 Hz for the two planar and 30 Hz for the third axis. Each axis is controlled by a motor/proportional magnetic particle brake combination fixed to the base. The innovative use of motors and brakes allows objects with high resistive torque requirements to be simulated without the stability and related safety issues involved with high torque, energy storing motors alone. Position, velocity, and applied endpoint force are sensed directly. Different control strategies are discussed and implemented, with an emphasis on how virtual environment force information, generated by the MIT Media Lab Computer Graphics and Animation Group BOLIO system, is transmitted to the device controller. The design constraints for a kinesthetic force feedback device can be summarized as: How can the symbiosis between the sense of presence in the virtual environment be maximized without compromising the interaction task under the constraints of the mechanical device limitations? Research in this field will yield insights to the optimal human sensory feedback mix for a wide spectrum of control and interaction problems. A flexible research tool that is designed as an easily reproducible product prototype has been constructed to explore the variety of possible force interaction.

Russo, Massimo↗

DOE OSTI2014

Charmed and light pseudoscalar meson decay constants from four-flavor lattice QCD with physical light quarks

We compute the leptonic decay constants f D+ , f Ds , and f K+ , and the quark-mass ratios m c =m s and m s =m l in unquenched lattice QCD using the experimentally determined value of f π+ for normalization. We use the MILC Highly Improved Staggered Quark (HISQ) ensembles with four dynamical quark flavors -- up, down, strange, and charm -- and with both physical and unphysical values of the light sea-quark masses. The use of physical pions removes the need for a chiral extrapolation, thereby eliminating a significant source of uncertainty in previous calculations. Four different lattice spacing ranging from a ≈ 0:06 fm to 0:15 fm are included in the analysis to control the extrapolation to the continuum limit. Our primary results are f D+ = 212:6(0:4)$(^{+1.0}_{-1.2})$ MeV, f Ds = 249:0(0:3)$(^{+1.1}_{-1.5})$ MeV, and f Ds /f D+ = 1:1712(10)$(^{+29}_{-32})$, where the errors are statistical and total systematic, respectively. The errors on our results for the charm decay constants and their ratio are approximately two to four times smaller than those of the most precise previous lattice calculations. We also obtain f K+ / f π+ = 1:1956(10)$(^{+26}_{-18})$, updating our previous result, and determine the quark-mass ratios m s /m l = 27:35(5)$(^{+10}_{-7})$ and m c /m s = 11:747(19)$(^{+59}_{-43})$. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the CKM matrix elements !Vus! = 0:22487(51)(29)(20)(5), !Vcd! = 0:217(1)(5)(1) and !Vcs! = 1:010(5)(18)(6), where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of !Vus!, the uncertainty in |Vud|, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

The Parareal Method for Solving Space-Dependent Neutron Kinetics Equations in Hexagonal-z Geometry

Neutron kinetics plays an important role in reactor safety and analysis. The backward Euler method is the most widely used time integration method in the calculation of space-dependent nuclear reactor kinetics. Xia (2006) used this method to discretize the derivative of flux, and adopted the analytic integration to discretize the delayed neutron precursor concentrations. Diagonally Implicit Runge-Kutta (DIRK) method owns high accuracy and excellent stability, and it was applied in the rectangle core for calculating core the nuclear reactor kinetics by Wenbo Zhao (2013). Recently, DIRK was applied to the neutron kinetics for hexagonal-z geometry application (CAI, 2016). As solving the neutron kinetics equations is very time-consuming and the number of available cores continues to increase with parallel architectures evolving, parallel algorithms need to be designed in order to effectively utilize the available resources. However, it is difficult to parallel in time axis because the later moment is strongly dependent on the previous moment. Some methods have been proposed to parallel in time algorithm over the last decades. One of the famous methods to realize time parallelism is the Parareal algorithm (Lions, 2001), which offers a promising approach to increase the degree of concurrency. More than one decade from this method first proposed, this method has been studied and applied in the many fields of applied mathematics, physics, chemistry, parallel computing, etc. Fischer (2014 applied the Parareal method in the field of the Navier-Stokes Equations. Baudron (2014) combined the Parareal method with finite element method to solve the neutron kinetics equations, where the Parareal method was implemented based on the implicit Euler method or method. In this paper, the Parareal method is studied in the hexagonal-z geometry with the help of DIRK developed for neutron kinetics equations (CAI, 2016)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Predicting metal–metal interactions. II. Accelerating generalized schemes through physical insights

Operando-computational frameworks that integrate descriptors for catalyst stability within catalyst screening paradigms enable predictions of rates and selectivity on chemically faithful representations of nanoparticles under reaction conditions. These catalyst stability descriptors can be efficiently predicted by density functional theory (DFT)-based models. The alloy stability model, for example, predicts the stability of metal atoms in nanoparticles with site-by-site resolution. In this work, we use physical insights to present accelerated approaches of parameterizing this recently introduced alloy-stability model. These accelerated approaches meld quadratic functions for the energy of metal atoms in terms of the coordination number with linear correlations between model parameters and the cohesive energies of bulk metals. By interpolating across both the coordination number and chemical space, these accelerated approaches shrink the training set size for 12 fcc p- and d-block metals from 204 to as few as 24 DFT calculated total energies without sacrificing the accuracy of our model. We validate the accelerated approaches by predicting adsorption energies of metal atoms on extended surfaces and 147 atom cuboctahedral nanoparticles with mean absolute errors of 0.10 eV and 0.24 eV, respectively. This efficiency boost will enable a rapid and exhaustive exploration of the vast material space of transition metal alloys for catalytic applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1980

An analytical transformation technique for generating uniformly spaced computational mesh

An analytical transformation method which can map arbitrary physical coordinate grid distribution into desired computational coordinate with uniform grid distribution is derived. The transformation function and its higher derivatives are differentiable. Salient features include; (1) precise control of grid sizes; (2) more than one location of clustered grids; (3) exact positioning of particular computational nodes in the physical plane; and (4) ensuring several patches of uniformly spaced grids in the physical plane for the higher accuracies (such as at the boundaries), etc., while keeping the variation of grid spacing continuous to avoid numerical instability.

Oh, Y. H.↗

NASA NTRS1995

NASA Tech Briefs, March 1995

This issue contains articles with a special focus on Computer-Aided design and engineering amd a research report on the Ames Research Center. Other subjects in this issue are: Electronic Components and Circuits, Electronic Systems, Physical Sciences, Materials, Computer Programs, Mechanics, Machinery, Manufacturing/Fabrication, Mathematics and Information Sciences and Life Sciences

Source record↗

NASA NTRS1973

Generation and physical characteristics of the ERTS MSS system corrected computer compatible tapes

The generation and format are discussed of the ERTS system corrected multispectral scanner computer compatible tapes. The discussion includes spacecraft sensors, scene characteristics, data transmission, and conversion of data to computer compatible tapes at the NASA Data Processing Facility. Geometeric and radiometric corrections, tape formats, and the physical characteristics of the tapes are also included.

Thomas, V. L.↗

DOE OSTI2025

Defect Generation and Evolution in Irradiated Epitaxial Films and Heterostructures of Fe 3 O 4 and Cr 2 O 3

Abstract The functionality of nuclear structural materials, sensors, and microelectronics in harsh environments such as radiation relies on understanding defect generation and evolution processes in oxide layers. The initial radiation response of epitaxial thin films of Fe 3 O 4 (111), Cr 2 O 3 (0001), and Fe 3 O 4 (111)/Cr 2 O 3 (0001) heterostructures deposited on Al 2 O 3 (0001) by oxygen‐assisted molecular beam epitaxy and irradiated with 200 keV He + is characterized. X‐ray diffraction and X‐ray absorption near edge spectroscopy showed that the Cr 2 O 3 layers underwent significant lattice expansion and disordering under irradiation, whereas the Fe 3 O 4 layers do not exhibit noticeable changes. In contrast, positron annihilation spectroscopy revealed an evolution of cation vacancy point defects in the Fe 3 O 4 layers into larger vacancy clusters with increasing irradiation, while the cation vacancies in Cr 2 O 3 remained primarily as single vacancies and small clusters. The results suggest that the Fe 3 O 4 lattice can utilize the free volume of the larger vacancy clusters to relax but the small vacancies in the Cr 2 O 3 lattice do not facilitate relaxation. Comparing defect concentrations in the single layer films versus the heterostructure suggests that point defects may cross the interface from Fe 3 O 4 into Cr 2 O 3 . Together, these results enhance the understanding of the initial defect evolution mechanisms in oxide layers in harsh irradiation environments.

36 MATERIALS SCIENCE↗

NASA NTRS1997

Dynamics of Numerics & Spurious Behaviors in CFD Computations

The global nonlinear behavior of finite discretizations for constant time steps and fixed or adaptive grid spacings is studied using tools from dynamical systems theory. Detailed analysis of commonly used temporal and spatial discretizations for simple model problems is presented. The role of dynamics in the understanding of long time behavior of numerical integration and the nonlinear stability, convergence, and reliability of using time-marching approaches for obtaining steady-state numerical solutions in computational fluid dynamics (CFD) is explored. The study is complemented with examples of spurious behavior observed in steady and unsteady CFD computations. The CFD examples were chosen to illustrate non-apparent spurious behavior that was difficult to detect without extensive grid and temporal refinement studies and some knowledge from dynamical systems theory. Studies revealed the various possible dangers of misinterpreting numerical simulation of realistic complex flows that are constrained by available computing power. In large scale computations where the physics of the problem under study is not well understood and numerical simulations are the only viable means of solution, extreme care must be taken in both computation and interpretation of the numerical data. The goal of this paper is to explore the important role that dynamical systems theory can play in the understanding of the global nonlinear behavior of numerical algorithms and to aid the identification of the sources of numerical uncertainties in CFD.

Yee, Helen C.↗

DOE OSTI2021

Towards Predictive Plasma Science and Engineering through Revolutionary Multi-Scale Algorithms and Models (Final Report)

This report describes the high-level accomplishments from the Plasma Science and Engineering Grand Challenge LDRD at Sandia National Laboratories. The Laboratory has a need to demonstrate predictive capabilities to model plasma phenomena in order to rapidly accelerate engineering development in several mission areas. The purpose of this Grand Challenge LDRD was to advance the fundamental models, methods, and algorithms along with supporting electrode science foundation to enable a revolutionary shift towards predictive plasma engineering design principles. This project integrated the SNL knowledge base in computer science, plasma physics, materials science, applied mathematics, and relevant application engineering to establish new cross-laboratory collaborations on these topics. As an initial exemplar, this project focused efforts on improving multi-scale modeling capabilities that are utilized to predict the electrical power delivery on large-scale pulsed power accelerators. Specifically, this LDRD was structured into three primary research thrusts that, when integrated, enable complex simulations of these devices: (1) the exploration of multi-scale models describing the desorption of contaminants from pulsed power electrodes, (2) the development of improved algorithms and code technologies to treat the multi-physics phenomena required to predict device performance, and (3) the creation of a rigorous verification and validation infrastructure to evaluate the codes and models across a range of challenge problems. These components were integrated into initial demonstrations of the largest simulations of multi-level vacuum power flow completed to-date, executed on the leading HPC computing machines available in the NNSA complex today. These preliminary studies indicate relevant pulsed power engineering design simulations can now be completed in (of order) several days, a significant improvement over pre-LDRD levels of performance.

42 ENGINEERING↗

DOE OSTI2019

Accomplishments and Year-End Performance Report: Wind Energy Program, Fiscal Year 2019

Offering broad-based technical expertise and world-class capabilities and facilities, NREL leverages these assets to provide the wind industry with a better understanding of fundamental physics, high-performance computing-enabled simulation tools, and the physical validation necessary to significantly lower the cost of wind energy. This report provides an overview of the many achievements NREL delivered on behalf of DOE's Wind Energy Technologies Office (WETO) and other partners during Fiscal Year (FY) 2019.

17 WIND ENERGY↗

DOE OSTI2023

Integration of a physics-based direct normal irradiance ($\mathrm{DNI}$) model to enhance the National Solar Radiation Database ($\mathrm{NSRDB}$)

The National Solar Radiation Database (NSRDB) is an extensively used dataset that furnishes satellite-retrieved solar resource data across the United States and an expanding list of other countries. Although the NSRDB uses a physical model to compute global horizontal irradiance (GHI), it currently employs an empirical approach based on surface observations to estimate cloudy-sky direct normal irradiance (DNI). Recently, a new physics-based approach, known as the Fast All-sky Radiation Model for Solar applications with DNI (FARMS-DNI), was developed to improve DNI forecasting. FARMS-DNI integrates direct and scattered solar radiances within the circumsolar region, resulting in improved day-ahead forecasting of DNI by incorporating it into the Weather Research and Forecasting model with Solar extensions (WRF-Solar). This study incorporates FARMS-DNI into the NSRDB algorithm to produce high-spatiotemporal-resolution DNI data from satellite data. The accuracy of the NSRDB based on FARMS-DNI is analyzed using surface observations from 19 sites situated within the National Oceanic and Atmospheric Administration (NOAA) Surface Radiation Budget (SURFRAD) and Solar Radiation (SOLRAD) networks, the University of Oregon (UO) network, the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) network, and at the National Renewable Energy Laboratory (NREL). The results demonstrate that FARMS-DNI reduces the significant overestimation of DNI in the conventional NSRDB at all surface sites, particularly in cloud overcast conditions classified using both satellite retrievals and surface observations. Consequently, this new model can effectively improve the overall accuracy of the NSRDB. The results also suggest that further improvement of DNI estimates at individual time steps, however, requires advanced satellite techniques and precise identification of clouds and retrieval of cloud properties.

14 SOLAR ENERGY↗

DOE OSTI2017

Sub genome anchored physical frameworks of the allotetraploid Upland cotton (Gossypium hirsutum L.) genome, and an approach toward reference-grade assemblies of polyploids

Like those of many agricultural crops, the cultivated cotton is an allotetraploid and has a large genome (~2.5 gigabase pairs). The two sub genomes, A and D, are highly similar but unequally sized and repeat-rich, which pose significant challenges for accurate genome reconstruction using standard approaches. Here we report the development of BAC libraries, sub genome specific physical maps, and a new-generation sequencing approach that will lead to a reference-grade genome assembly for Upland cotton. Three BAC libraries were constructed, fingerprinted, and integrated with BAC-end sequences (BES) to produce a de novo whole-genome physical map. The BAC map was partitioned by sub genomes through alignment to the diploid progenitor D-genome reference sequence with densely spaced BES anchor points and computational filtering. The physical maps were validated with FISH and genetic mapping of SNP markers derived from BES. Two pairs of homeologous chromosomes, A11/D11 and A12/D12, were used to assess multiplex sequencing approaches for completeness and scalability. The results represent the first sub genome anchored physical maps of Upland cotton, and a new-generation approach to the whole-genome sequencing, which will lead to the reference-grade assembly of allopolyploid cotton and serve as a general strategy for sequencing other polyploid species.

59 BASIC BIOLOGICAL SCIENCES↗