Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Flow Physics at a Large Rotor Tip Gap in a Multi-Stage Axial Compressor

The flow physics in a large rotor tip gap in a 1.5-stage axial compressor is investigated in the current study. The flow structure in the rotor tip region is complex with several dominant vortical structures of opposite rotation, resulting in inhomogeneous and highly anisotropic turbulence. Earlier measurements show that eddy viscosity is negative over large parts of the tip region and eddy viscosity varies among stress/strain components. The present study aims to understand how the complex nature of rotor tip leakage flow affects compressor performance when the tip gap size is greater than 4–5% of the rotor span, which is typical of advanced small core engines. Unsteady Reynolds-averaged Navier-Stokes (URANS) and Large Eddy Simulation (LES) techniques are applied to study flow physics in a large rotor tip gap (5.5% of rotor span) in a 1.5-stage axial compressor. Calculated flow fields from the two different approaches are compared with available measurements and examined in detail. LES calculates the pressure rise in the present compressor fairly well, while URANS with a standard two-equation turbulence closure underpredicts the pressure rise by 15-20% of the measured values. The current study shows that URANS with the current turbulence closure produces much higher all-positive eddy viscosity in the tip-gap region compared to measurements and LES. The distribution of eddy viscosity in the URANS simulation is also wrong. Consequently, the flow in the tip region is highly damped with significantly larger blockage generation, which results in the tip leakage vortex (TLV) staying closer to the blade suction side compared to the measurement. When the TLV stays closer to the blade, both flow turning and the pressure rise across the compressor are reduced compared to the measurements. It appears that this effect is amplified by a large rotor tip gap.

tip gap↗
DOE OSTI2023

Molecular Symmetry in VQE: A Dual Approach for Trapped-Ion Simulations of Benzene

Understanding complex chemical systems—such as biomolecules, catalysts, and novel materials—is a central goal of quantum simulations. Near-term strategies hinge on the use of variational quantum eigensolver (VQE) algorithms combined with a suitable ansatz. However, straightforward application of many chemically-inspired ansatze yields prohibitively deep circuits. In this work, we employ several circuit optimization methods tailored for trapped-ion quantum devices to enhance the feasibility of intricate chemical simulations. The techniques aim to lessen the depth of the unitary coupled cluster with singles and doubles (uCCSD) ansatz's circuit compilation, a considerable challenge on current noisy quantum devices. Furthermore, we use symmetry-inspired classical post-selection methods to further refine the outcomes and minimize errors in energy measurements, without adding quantum overhead. Our strategies encompass optimal mapping from orbital to qubit, term reordering to minimize entangling gates, and the exploitation of molecular spin and point group symmetry to eliminate redundant parameters. The inclusion of error mitigation via post-selection based on known molecular symmetries improves the results to near milli-Hartree accuracy. These methods, when applied to a benzene molecule simulation, enabled the construction of an 8-qubit circuit with 69 two-qubit entangling operations, pushing the limits for variational quantum eigensolver (VQE) circuits executed on quantum hardware to date. 1 1 This manuscript has been authored in part by UT-Battelle, LLC, under contract DE-AC05-000R22725 with the US Department of Energy (DOE). The publisher acknowledges the US government license to provide public access under the DOE Public Access Plan(https://energy.gov/doe-public-access-plan).

Goings, Joshua↗
DOE OSTI2023

Does Carrier Envelope Phase Affect the Ionization Site in a Neutral Diatomic Molecule?

A recent work shows how to extract the ionization site of a neutral diatomic molecule by comparing Quantum Trajectory Monte Carlo (QTMC) simulations with experimental measurements of the final electron momenta distribution. This method was applied to an experiment using a 40-femtosecond infrared pulse, finding that a downfield atom is roughly twice as likely to be ionized as an upfield atom in a neutral nitrogen molecule. However, an open question remains as to whether an assumption of the zero carrier envelope phase (CEP) used in the above work is still valid for short, few-cycle pulses where the CEP can play a large role. Given experimentalists’ limited control over the CEP and its dramatic effect on electron momenta after ionization, it is desirable to see what influence the CEP may have in determining the ionization site. In this paper, we employ QTMC techniques to simulate strong-field ionization and electron propagation from neutral N2 using an intense 6-cycle laser pulse with various CEP values. Comparing simulated electron momenta to experimental data indicates that the ratio of down-to-upfield ions remains roughly 2:1 regardless of the CEP. This confirms that the ionization site of a neutral molecule is determined predominantly by the laser frequency and intensity, as well as the ground-state molecular wavefunction, and is largely independent of the CEP.

Schimmoller, Alex↗
DOE OSTI2021

Distribution of the Spectral Response of Cells in Silicon Modules - Mechanisms and Implications

Measurements of the spectral response of cells in a PV module can provide important insight into the physics of the device and are also needed for the accurate calibration of the performance of these devices using solar simulators. While techniques for measuring the spectral response of PV cells are well established, measurement of the spectral response of all cells in a module is not straightforward since in this case one has electrical access to a string of cells and not to cells individually. Here we present a new approach for testing the spectral response of individual cells in silicon modules that allows for the fast mapping of the response at different locations of the module. We show that in some cases the spectral response of different cells in a module can vary considerably and propose an underlying mechanism for this variation. We also discuss the implications of this observation for the uncertainty of the spectral correction and of the I SC calibration using a solar simulator.

modules↗
NASA NTRS2001

Technique for Very High Order Nonlinear Simulation and Validation

Finding the sources of sound in large nonlinear fields via direct simulation currently requires excessive computational cost. This paper describes a simple technique for efficiently solving the multidimensional nonlinear Euler equations that significantly reduces this cost and demonstrates a useful approach for validating high order nonlinear methods. Up to 15th order accuracy in space and time methods were compared and it is shown that an algorithm with a fixed design accuracy approaches its maximal utility and then its usefulness exponentially decays unless higher accuracy is used. It is concluded that at least a 7th order method is required to efficiently propagate a harmonic wave using the nonlinear Euler equations to a distance of 5 wavelengths while maintaining an overall error tolerance that is low enough to capture both the mean flow and the acoustics.

Dyson, Rodger W.↗
NASA NTRS2014

A Numerical and Experimental Study of Damage Growth in a Composite Laminate

The present study has three goals: (1) perform an experiment where a simple laminate damage process can be characterized in high detail; (2) evaluate the performance of existing commercially available laminate damage simulation tools by modeling the experiment; (3) observe and understand the underlying physics of damage in a composite honeycomb sandwich structure subjected to low-velocity impact. A quasi-static indentation experiment has been devised to provide detailed information about a simple mixed-mode damage growth process. The test specimens consist of an aluminum honeycomb core with a cross-ply laminate facesheet supported on a stiff uniform surface. When the sample is subjected to an indentation load, the honeycomb core provides support to the facesheet resulting in a gradual and stable damage growth process in the skin. This enables real time observation as a matrix crack forms, propagates through a ply, and then causes a delamination. Finite element analyses were conducted in ABAQUS/Explicit(TradeMark) 6.13 that used continuum and cohesive modeling techniques to simulate facesheet damage and a geometric and material nonlinear model to simulate core crushing. The high fidelity of the experimental data allows a detailed investigation and discussion of the accuracy of each numerical modeling approach.

McElroy, Mark↗
NASA NTRS2014

TPSAS-NF1676L-19616-DND

The present study has three goals: (1) perform an experiment where a simple laminate damage process can be characterized in high detail; (2) evaluate the performance of existing commercially available laminate damage simulation tools by modeling the experiment; (3) observe and understand the underlying physics of damage in a composite honeycomb sandwich structure subjected to low-velocity impact. A quasi-static indentation experiment has been devised to provide detailed information about a simple mixed-mode damage growth process. The test specimens consist of an aluminum honeycomb core with a cross-ply laminate facesheet supported on a stiff uniform surface. When the sample is subjected to an indentation load, the honeycomb core provides support to the facesheet resulting in a gradual and stable damage growth process in the skin. This enables real time observation as a matrix crack forms, propagates through a ply, and then causes a delamination. Finite element analyses were conducted in ABAQUS/Explicit(TradeMark) 6.13 that used continuum and cohesive modeling techniques to simulate facesheet damage and a geometric and material nonlinear model to simulate core crushing. The high fidelity of the experimental data allows a detailed investigation and discussion of the accuracy of each numerical modeling approach.

Mark McElroy↗
DOE OSTI2023

Review of Carbon Support Coordination Environments for Single Metal Atom Electrocatalysts (SACS)

This topical review focuses on the distinct role of carbon support coordination environment of single-atom catalysts (SACs) for electrocatalysis. The article begins with an overview of atomic coordination configurations in SACs, including a discussion of the advanced characterization techniques and simulation used for understanding the active sites. A summary of key electrocatalysis applications is then provided. These processes are oxygen reduction reaction (ORR), oxygen evolution reaction (OER), hydrogen evolution reaction (HER), nitrogen reduction reaction (NRR), and carbon dioxide reduction reaction (CO 2 RR). The review then shifts to modulation of the metal atom-carbon coordination environments, focusing on nitrogen and other non-metal coordination through modulation at the first coordination shell and modulation in the second and higher coordination shells. Representative case studies are provided, starting with the classic four-nitrogen-coordinated single metal atom (M$-$N 4 ) based SACs. Bimetallic coordination models including homo-paired and hetero-paired active sites are also discussed, being categorized as emerging approaches. The theme of the discussions is the correlation between synthesis methods for selective doping, the carbon structure–electron configuration changes associated with the doping, the analytical techniques used to ascertain these changes, and the resultant electrocatalysis performance. In conclusion, critical unanswered questions as well as promising underexplored research directions are identified.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Small-scale properties from exascale computations of turbulence on a $\mathbf{32\,768^3}$ periodic cube

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number R λ = 2500 on a 32 768 3 three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the 32 768 3 data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

intermittency↗
DOE OSTI2022

Differentiable Preisach Modeling for Characterization and Optimization of Particle Accelerator Systems with Hysteresis

Future improvements in particle accelerator performance are predicated on increasingly accurate online modeling of accelerators. Hysteresis effects in magnetic, mechanical, and material components of accelerators are often neglected in online accelerator models used to inform control algorithms, even though reproducibility errors from systems exhibiting hysteresis are not negligible in high precision accelerators. Here, we combine the classical Preisach model of hysteresis with machine learning techniques to efficiently create nonparametric, high-fidelity models of arbitrary systems exhibiting hysteresis. We experimentally demonstrate how these methods can be used in situ, where a hysteresis model of an accelerator magnet is combined with a Bayesian statistical model of the beam response, allowing characterization of magnetic hysteresis solely from beam-based measurements. Finally, we explore how using these joint hysteresis-Bayesian statistical models allows us to overcome optimization performance limitations that arise when hysteresis effects are ignored.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Path Integral Monte Carlo Simulations of Warm Dense Matter and Plasmas

New path integral Monte Carlo simulation (PIMC) techniques will be developed and applied to derive the equation of state (EOS) for the regime of warm dense matter and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 10,000 K, this method is not applicable at much higher temperature where electronic excitations become important because the number of partially occupied electronic orbitals reaches intractably large numbers and, more importantly, the use of zero-temperature exchange-correlation functionals introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it has commonly been assumed such methods can only be applied to elements without core electrons like hydrogen and helium, we recently showed how to extend PIMC to heavier elements by performing the first PIMC simulations of carbon and water plasmas [Driver, Militzer, Phys. Rev. Lett. 108 (2012) 115502]. Here we propose to continue this important development to extend the reach of PIMC simulations to yet heavier elements and also lower temperatures. The goal is to provide a robust first-principles simulation method that can accurately and efficiently study materials with excited electrons at solid-state densities in order to access parts of the phase diagram such the regime of warm dense matter and plasmas where so far only more approximate, semi-analytical methods could be applied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Scalable qubit representations of neutrino mixing matrices

Oscillating neutrino beams exhibit quantum coherence over distances of thousands of kilometers. Their unambiguously quantum nature suggests an appealing test system for direct quantum simulation. Such techniques may enable presently analytically intractable calculations involving multi-neutrino entanglements, such as collective neutrino oscillations in supernovae, but only once oscillation phenomenology is properly re-expressed in the language of quantum circuits. Here we resolve outstanding conceptual issues regarding encoding of arbitrarily mixed neutrino flavor states in the Hilbert space of an n-qubit quantum computer. We introduce algorithms to encode mixing and oscillation of any number of flavor-mixed neutrinos, both with and without CP-violation, with an efficient number of prescriptive input parameters in terms of sub-rotations of the PMNS matrix in standard form. As a result, examples encoded for an IBM-Q quantum computer are shown to converge to analytic predictions both with and without CP-violation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Procedure Parsing: A Method for Parsing Handwritten Documents into Computer-Based Procedures

The nuclear industry is heavily procedure driven, where almost everything has a step-by-step instruction that is expected to be followed in detail. Historically, these procedures were printed on paper copies. Recently, the industry transitioned towards electronic copies (i.e., PDFs on tablets). One major drive for this transition is the introduction of human error and loss of situation awareness when using paper copies. However, electronic copies of documents inherently have the same error traps as their paper cousins. Therefore, there is an increased interest in a way to utilize the information in the step-by-step guidance, but to present it in a dynamic manner that guides the user and adapts to any encountered conditions. Researchers at Idaho National Laboratory propose a flexible, automated method based on document parsing and augmented by natural language processing (NLP) techniques, to address these shortcomings and capitalize on these recent advancements in machine learning. The proposed method provides a cost-effective solution for computer-assisted procedure parsing of hand-written control room procedures, originally authored in Word or PDF formats, into instructions that can be displayed as computer-based procedures (CBP) in a modern graphical user interface. The researchers devised, implemented and demonstrated the Operating Procedure Extender for Novel Systems (OPENS) method in 2020. The key to OPENS is to map the original procedure text into a context-free grammar, tying content to equipment, locations, and other steps, actions, etc. This formal grammar is then used to isolate and define keywords and actions verbs, such as “measure” or “evaluate” and tie them to specific equipment referenced within that step or located in other steps, substeps, actions, subactions and tables throughout the procedure. OPENS generates an abstract syntax tree from the document which it uses to store a copy of this information in the open-standard, machine-readable and human-readable file formats XML and JSON. The XML is useful to preserve the relational aspects of the procedure for referencing tables and branching information so the user can be directed to the next appropriate active step based on the values entered for that step and previous steps. The JSON is useful for storing and exchanging data objects used to track responses to previous steps and state changes in simulated environments. In future iterations, these formats can also be used for storing more detailed information about input during plant operation or simulation. The techniques the researcher developed could further be improved by integration of recent advancements in machine learning. NLP methods could standardize documents, correct for grammatical error, and provide automated semantic validation. The researcher expects that self-supervised techniques applied to collections of natural language instructions could strengthen the model with broader context. All these methods together give us a practical way to automatically extract protocols from documents and user interactions, empowering researchers, procedure writers and nuclear operators while moving the industry forward.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS1978

Data Management Symposium, Huntsville, Ala., October 18, 19, 1977, Proceedings

The transfer, processing, and use of satellite data are discussed. Topics relating to simulation and processing techniques include computer design assessments through simulation, data system dynamic simulation, and future timing accuracy requirements and procedures for data processing. Subjects relevant to data base management systems and users and their needs include an image-based information system architecture for correlation satellite and topological data base, efficient searching and sorting applications using an associative array processor, analysis of user's needs for a large technical data base, and variable length data formats. The technology development outlook was considered with attention to data systems technology outlook for supporting NASA programs, a user-oriented interactive information extraction system, and flow control and sorting techniques for telemetry packets.

Castelli, A.↗
DOE OSTI2019

Influence of Electropulsing Treatment on the Initial Solidification of Molten Steel During Continuous Casting

An electropulsing-assisted mold simulator (EPMS) technique was developed to investigate the effect of pulsed electric current on the initial solidification behavior of molten steel during the process of continuous casting. The results indicated that the variation of mold hot surface temperatures in the meniscus area decreased from 375.8 K to 394 K for the case without electropulsing treatment (EPT) to 363.5 K to 377.3 K for the one with EPT. The mold surface heat fluxes fluctuated around the baseline of 0.83 and 2.42 MW/m{sup 2} in the cases with and without EPT, respectively, which indicated that a higher thermal resistance between the initial shell and copper mold was introduced by employing the EPT. The crystal fraction of slag films increased from 69.8 to 77.9 pct, whereas the grain size of mold flux film decreased, when compared the case of EPT treatment with the benchmark one, which suggested that the pulsed electric current promotes the movement of molten clusters and particles, leading to the enhancement of mold flux crystallization and the improvement of thermal resistance. Moreover, the solidification structure of steel along the direction of solidification (from the mold to melt) includes a fine chilled grain layer close to the mold, a columnar dendrite layer in the middle, and equiaxed grains next to the melt core, in which a compacted microstructure with smaller secondary dendrite arm spacing was observed in the case of the EPT. The results demonstrate that the EPT could refine the as-cast solidification microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Planning of distributed data production for High Energy and Nuclear Physics

Modern experiments in High Energy and Nuclear Physics heavily rely on distributed computations using multiple computational facilities across the world. One of the essential types of the computations is a distributed data production where petabytes of raw files from a single source has to be processed once (per production campaign) using thousands of CPUs at distant locations and the output has to be transferred back to that source. The data distribution over a large system does not necessary match the distribution of storage, network and CPU capacity. Therefore, bottlenecks may appear and lead to increased latency and degraded performance. In this paper we propose a new scheduling approach for distributed data production which is based on the network flow maximization model. In our approach a central planner defines how much input and output data should be transferred over each network link in order to maximize the computational throughput. Such plans are created periodically for a fixed planning time interval using up-to-date information on network, storage and CPU resources. The centrally created plans are executed in a distributed manner by dedicated services running at participating sites. In conclusion, our simulations based on the log records from the data production framework of the experiment STAR (Solenoid Tracker at RHIC) have shown that the proposed model systematically provides a better performance compared to the simulated traditional techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Surface Effects on the Dynamic Behavior of Vortices in Type II Superconducting Strips with Periodic and Conformal Pinning Arrays

Using molecular dynamics techniques, we simulate the vortex behavior in a type II superconducting strip in the presence of triangular and two types of conformal pinning arrays, one with a pinning gradient perpendicular to the driving force (C1) and the other parallel (C2), at zero temperature. A transport force is applied in the infinite direction of the strip, and the magnetic field is increased until the rate between the density of vortices (n{sub v}) and pinning (n{sub p}) reaches n{sub v}/n{sub p} = 1.5. Our results show a monotonic increase, by steps, of the vortex density with the applied magnetic field. Besides, each pinning lattice presents a different vortex penetration delay. For the triangular pinning array, different than the case of infinite films, here the system exhibits only one pronounced depinning force peak at n{sub v}/n{sub p} = 1. However, the depinning force peak is present for only one value of field, in the range of fields where n{sub v}/n{sub p} = 1 is stable. For the case of conformal pinning arrays, in analogy to what is observed in infinite films, no commensurability depinning force peaks were found. In the present case, the C1 array is more efficient at low fields, all arrays are equivalent in the intermediate fields, and the C2 array is more efficient for high fields. We also show that for the C1 array at high fields, vortices depin following the conformal arches, from the edge to the center. For the C2 array, the depinning force is higher when compared to that of C1, because this particular conformal structure prevents the formation of easy vortex flow channels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Reaction of amorphous/crystalline SiOC/Fe interfaces by thermal annealing

The development of revolutionary new alloys and composites is crucial to meeting materials requirements for next generation nuclear reactors. The newly developed amorphous silicon oxycarbide (SiOC) and crystalline Fe composite system has shown radiation tolerance over a wide range of temperatures. To advance understanding of this new composite, we investigate the structure and thermal stability of the interface between amorphous SiOC and crystalline Fe by combining various experimental techniques and simulation methods. We show that the SiOC/Fe interface is thermally stable up to at least 400 °C. When the annealing temperature reaches 600 °C, an intermixed region forms at this interface. This region appears to be a crystalline phase that forms an incoherent interface with the Fe layer. Density functional theory (DFT) Molecular dynamics (MD) is performed on the homogeneous SiFeOC phase to study the early stages of 2 formation of the intermixed layer. Both experimental and simulation results suggest this phase has the fayalite crystal structure. As a result, the physical processes involved in the formation of the intermixed region are discussed.

36 MATERIALS SCIENCE↗