Search NASA⌕ Search

SEARCH · Search NASA

Results for “test case”

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 1,441 records · Page 80

Less is more: Sampling chemical space with active learning

We present the development of accurate and transferable machine learning (ML) potentials for predicting molecular energetics is a challenging task. The process of data generation to train such ML potentials is a task neither well understood nor researched in detail. In this work, we present a fully automated approach for the generation of datasets with the intent of training universal ML potentials. It is based on the concept of active learning (AL) via Query by Committee (QBC), which uses the disagreement between an ensemble of ML potentials to infer the reliability of the ensemble’s prediction. QBC allows the presented AL algorithm to automatically sample regions of chemical space where the ML potential fails to accurately predict the potential energy. AL improves the overall fitness of ANAKIN-ME (ANI) deep learning potentials in rigorous test cases by mitigating human biases in deciding what new training data to use. AL also reduces the training set size to a fraction of the data required when using naive random sampling techniques. To provide validation of our AL approach, we develop the COmprehensive Machine-learning Potential (COMP6) benchmark (publicly available on GitHub) which contains a diverse set of organic molecules. Active learning-based ANI potentials outperform the original random sampled ANI-1 potential with only 10% of the data, while the final active learning-based model vastly outperforms ANI-1 on the COMP6 benchmark after training to only 25% of the data. Finally, we show that our proposed AL technique develops a universal ANI potential (ANI-1x) that provides accurate energy and force predictions on the entire COMP6 benchmark. Finally, this universal ML potential achieves a level of accuracy on par with the best ML potentials for single molecules or materials, while remaining applicable to the general class of organic molecules composed of the elements CHNO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Highly parallelisable simulations of time-dependent viscoplastic fluid flow with structured adaptive mesh refinement

We present the extension of an efficient and highly parallelisable framework for incompressible fluid flow simulations to viscoplastic fluids. The system is governed by incompressible conservation of mass, the Cauchy momentum equation, and a generalised Newtonian constitutive law. In order to simulate a wide range of viscoplastic fluids, we employ the Herschel-Bulkley model for yield-stress fluids with nonlinear stress-strain dependency above the yield limit. We utilise Papanastasiou regularisation in our algorithm to deal with the singularity in apparent viscosity. The resulting system of partial differential equations is solved using the IAMR (Incompressible Adaptive Mesh Refinement) code, which uses second-order Godunov methodology for the advective terms and semi-implicit diffusion in the context of an approximate projection method to solve adaptively refined meshes. By augmenting the IAMR code with the ability to simulate regularised Herschel-Bulkley fluids, we obtain efficient numerical software for time-dependent viscoplastic flow in three dimensions, which can be used to investigate systems not considered previously due to computational expense. We validate results from simulations using this new capability against previously published data for Bingham plastics and power-law fluids in the two-dimensional lid-driven cavity. In doing so, we expand the range of Bingham and Reynolds numbers which have been considered in the benchmark tests. Moreover, extensions to time-dependent flow of Herschel-Bulkley fluids and three spatial dimensions offer new insights into the flow of viscoplastic fluids in this test case, and we provide missing benchmark results for these extensions.

42 ENGINEERING↗

Non-orthogonal determinants in multi-Slater-Jastrow trial wave functions for fixed-node diffusion Monte Carlo

The accuracy and efficiency of ab initio Quantum Monte Carlo (QMC) algorithms benefit greatly from compact variational trial wave functions that accurately reproduce ground state properties of a system. We investigate the possibility of using multi-Slater-Jastrow trial wave functions with non-orthogonal determinants by optimizing identical single particle orbitals independently in separate determinants. As a test case, we compute variational and fixed-node diffusion Monte Carlo (FN-DMC) energies of a C 2 molecule. For a given multi-determinant expansion, we find that this non-orthogonal orbital optimization results in a consistent improvement in the variational energy and the FN-DMC energy on the order of a few tenths of an eV. In some cases, fewer non-orthogonal determinants are required compared to orthogonal ones in order to achieve similar accuracy in FN-DMC. Our calculations indicate that trial wave functions with non-orthogonal determinants can improve computed energies in a QMC calculation when compared to their orthogonal counterparts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The use of strain to tailor electronic thermoelectric transport properties: A first principles study of 2H-phase CuAlO 2

Using first principles calculations, the use of strain to adjust electronic transport and the resultant thermoelectric (TE) properties is discussed using 2H phase CuAlO 2 as a test case. Transparent oxide materials, such as CuAlO 2 , a p-­type transparent conducting oxide (TCO), have recently been studied for high temperature thermoelectric power generators and coolers for waste heat. Given TCO materials with relative ease of fabrication, low cost of materials, and non-­toxicity, the ability to tailor them to specific temperature ranges, power needs, and size requirements, through the use of strain opens an interesting avenue. We find that strain can have a significant effect on these properties, in some cases detrimental and in others beneficial, including the potential for n-­type power factors larger than the highest p-type case. The physical reasons for this behavior are explained in the terms of the thermoelectric transport distribution and the Landauer distribution of modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comparison of Optical Modelling Tools for Sunshape and Surface Slope Error

The sunshape and mirror surface slope error are two main factors that influence the characteristics of incoming and reflected solar beams in a solar concentrator. Different optical modelling tools, however, address these two factors with different methods and simplifications. In this study, simulation methods are reviewed in five optical modelling tools, namely SolTrace, Tonatiuh, Tracer, Heliosim and SolarPILOT. A series of simple test cases are simulated in each tool to compare and verify the implementation of sunshape and slope error models separately. This strategy has successfully identified discrepancies among the tools and contributed to a better understanding of their reasons. It has also uncovered few glitches in some of those tools and contributed to improve them. It is our hope that this study will help to resolve common misconceptions, to build confidence amongst CSP researches on the accuracy of tools with which they have less experience, and to ensure better agreement in optical modelling work in future.

14 SOLAR ENERGY↗

Use of Cernox thermometers in AC specific heat measurements under pressure

We report on the resistance behavior of bare-chip Cernox thermometers under pressures up to 2 GPa, generated in a piston-cylinder pressure cell. Our results clearly show that Cernox thermometers, frequently used in low-temperature experiments due to their high sensitivity, remain highly sensitive even under applied pressure. We show that these thermometers are therefore ideally suited for measurements of heat capacity under pressure utilizing an ac oscillation technique up to at least 150 K. Our Cernox-based system is very accurate in determining changes in the specific heat as a function of pressure as demonstrated by measurements of the heat capacity on three different test cases: (i) the superconducting transition in elemental Pb ( T c = 7.2 K), (ii) the antiferromagnetic transition in the rare-earth compound GdNiGe 3 ( T N = 26 K), and (iii) the structural/magnetic transition in the iron-pnictide BaFe 2 As 2 ( T s,N = 130 K). In conclusion, the chosen examples demonstrate the versatility of our technique for measuring the specific heat under the pressure of various condensed-matter systems with very different transition temperatures as well as amounts of removed entropy.

47 OTHER INSTRUMENTATION↗

Saturation and nonlinear electromagnetic stabilization of ITG turbulence

Energy transfer in ion-temperature-gradient-driven (ITG) turbulence and its role in modeling transport are examined for finite normalized plasma pressure β for a number of test cases and experimental discharges. The analysis shows that like the zero-β case, finite-β ITG turbulence saturates by nonlinear energy transfer to stable modes mediated by a zonal flow. Electromagnetic effects reliably increase stable mode amplitudes but affect heat fluxes only at the ≈5% level. The most important change with increased β is an increase in the correlation time of the triplet interaction of the unstable mode, stable mode, and zonal flow, thus providing a heightened nonlinear energy transfer efficiency, which allows the instability to saturate at lower amplitude. The heat flux is examined in connection with nonlinear electromagnetic stabilization, the phenomenon where the flux falloff with β is more pronounced than the falloff predicted by quasilinear transport models. The inclusion of the triplet correlation time in the quasilinear model captures most of the nonlinearly enhanced stabilization for the configurations studied here.

Whelan, G. G.↗

ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field

We have implemented the Martini force field within Lawrence Livermore National Laboratory’s molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein–lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-color scattering for the measurement of neutrals at the edge of fusion devices

Laser two-color scattering (TCS) is proposed to detect the neutral species in the edge of fusion devices, namely, tokamaks. Here, TCS uses two wavelengths to probe both the laser Rayleigh scattering and Thomson scattering of the neutral-electron bath, with emphasis on neutral density measurements such as that of hydrogen and deuterium. Modeling of the Rayleigh scattering of tokamak neutral species under various plasma conditions (electron density and temperature) shows that, with an appropriate filtering of the Thomson signal and by going to ultraviolet-region wavelengths, identification of the Rayleigh signal can be achieved. Photon count and signal fractions are calculated in two test cases, one in the midplane region of the National Spherical Torus Experiment and one in the divertor region of DIII-D. An uncertainty analysis and discussion of the feasibility of the TCS diagnostic is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fully variational incremental CASSCF

We report the complete-active-space self-consistent field (CASSCF) method is a canonical electronic structure theory that holds a central place in conceptualizing and practicing first principles simulations. For application to realistic molecules, however, the CASSCF must be approximated to circumvent its exponentially scaling computational costs. Applying the many-body expansion - also known as the method of increments - to CASSCF (iCASSCF) has been shown to produce a polynomially scaling method that retains much of the accuracy of the parent theory and is capable of treating full valence active spaces. Due to an approximation made in the orbital gradient, the orbital parameters of the original iCASSCF formulation could not be variationally optimized, which limited the accuracy of its nuclear gradient. Herein, a variational iCASSCF is introduced and implemented, where all parameters are fully optimized during energy minimization. This method is able to recover electronic correlations from the full valence space in large systems, produce accurate gradients, and optimize stable geometries as well as transition states. Demonstrations on challenging test cases, such as the oxoMn(salen)Cl complex with 84 electrons in 84 orbitals and the automerization of cyclobutadiene, show that the fully variational iCASSCF is a powerful tool for describing challenging molecular chemistries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Instrumentation development and testing of a wind turbine blade for sub-scale wake studies

Wind turbine blade aerodynamics and the resulting wake flow are complex, and wind tunnel testing of these flows can provide critical insight. The data from such tests are also valuable for validation of numerical models. For experiments using sub-scale turbines to be useful, the blade aerodynamics and wakes must exhibit the important features observed in utility-scale turbine flows. In this work, ~1 m blades designed to produce physics relevant to larger scale turbines were manufactured and affixed to an existing turbine. One of the blades was instrumented for blade surface pressure measurements. Time-dependent surface pressure measurements coupled with instantaneous measurements of the inflow with Cobra probes and the wake with particle image velocimetry allowed for characterization of the inflow, blade flow, and near wake. The results demonstrate that the instrumentation was effective in characterizing the blade loading and the flow field. For the one test case discussed in this paper, the measurements of inflow, blade loading, and wake properties facilitate understanding of the wake's behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Verification of a fully implicit particle-in-cell method for the <!--${MathJax: TeX-AMS-MML_HTMLorMML}--> v &#x2225; -formalism of electromagnetic gyrokinetics in the XGC code

A fully implicit particle-in-cell method for handling the v ∥ -formalism of electromagnetic gyrokinetics has been implemented in XGC. By choosing the v ∥ -formalism, here we avoid introducing the nonphysical skin terms in Ampère's law, which are responsible for the well-known “cancellation problem” in the p ∥ -formalism. The v ∥ -formalism, however, is known to suffer from a numerical instability when explicit time integration schemes are used due to the appearance of a time derivative in the particle equations of motion from the inductive component of the electric field. Here, using the conventional δf scheme, we demonstrate that our implicitly discretized algorithm can provide numerically stable simulation results with accurate dispersive properties. We verify the algorithm using a test case for shear Alfvén wave propagation in addition to a case demonstrating the ion temperature gradient-kinetic ballooning mode (ITG-KBM) transition. The ITG-KBM transition case is compared to results obtained from other δf gyrokinetic codes/schemes, whose verification has already been archived in the literature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accurate frozen core approximation for all-electron density-functional theory

We implement and benchmark the frozen core approximation, a technique commonly adopted in electronic structure theory to reduce the computational cost by means of mathematically fixing the chemically inactive core electron states. The accuracy and efficiency of this approach are well controlled by a single parameter, the number of frozen orbitals. Explicit corrections for the frozen core orbitals and the unfrozen valence orbitals are introduced, safeguarding against seemingly minor numerical deviations from the assumed orthonormality conditions of the basis functions. A speedup of over twofold can be achieved for the diagonalization step in all-electron density-functional theory simulations containing heavy elements, without any accuracy degradation in terms of the electron density, total energy, and atomic forces. This is demonstrated in a benchmark study covering 103 materials across the Periodic Table and a large-scale simulation of CsPbBr3 with 2560 atoms. Our study provides a rigorous benchmark of the precision of the frozen core approximation (sub-meV per atom for frozen core orbitals below −200 eV) for a wide range of test cases and for chemical elements ranging from Li to Po. The algorithms discussed here are implemented in the open-source Electronic Structure Infrastructure software package.

Yu, Victor Wen-zhe↗

A high temperature cell for investigating interfacial structure on the molecular scale in molten salt/alloy systems

In this work, we describe the design and development of an in situ neutron reflectometry cell for high temperature investigations of structural changes occurring at the interface between inorganic salts, in their molten state up to 800 °C, and corrosion resistant alloys or other surfaces. In the cell, a molten salt is confined by an annular ring of single crystal sapphire constrained between the sample substrate and a sapphire plate using two gold O-rings, enclosing a liquid salt volume of 20 ml, along with a dynamic cell volume to accommodate expansion of the liquid with heating. Further, as a test case for the cell, we report on an in situ neutron reflectometry measurement of the interface between a eutectic salt mixture of MgCl 2 –KCl (32:68 molar ratio) and a single crystal sapphire substrate at 450 °C, resulting in the formation of a 60 Å layer having a scattering length density of 1.72 × 10 –6 Å –2 . While the origin of this layer is uncertain, it is likely to have resulted from the salt reacting with an existing impurity layer on the sapphire substrate.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Transient mixing behavior of a supercritical fluid injected into supercritical and subcritical environments

Transient, near-field mixing of a supercritical fluid injected into supercritical and subcritical environments is investigated experimentally. Fuel (propane) from a 100 μm diameter Bosch diesel injector is injected into a high-pressure test rig with a steady flow of heated, low-velocity nitrogen to control the ambient pressure and temperature. The fuel supply pressure is approximately 190 bar, allowing short injections lasting 3–4 ms to be performed in quick succession to acquire statistically significant datasets to represent the unsteady fuel–air mixing. Two test cases with the same ambient density and injector inlet pressure and temperature were chosen in an attempt to isolate effects of mixture phase behavior. Microscopic rainbow Schlieren deflectometry (RSD) is used for visualization and quantitative data acquisition. In the case of the supercritical ambient, the injector inlet and exit and ambient conditions are all supercritical. The RSD analysis reveals a mixing process similar to the single phase fluid mixing, though condensation is also observed locally, and can be explained by the adiabatic mixing theory presented in the paper. For the subcritical ambient case, the injector inlet is supercritical, but the injector exit and ambient conditions are both subcritical. The results show a rapid expansion in the near-field promoting mixing between propane and nitrogen. At downstream locations, significant condensation is also observed as predicted by the adiabatic mixing theory. As a result, the supercritical ambient case exhibited longer and narrower jet penetration compared to the subcritical ambient case.

02 PETROLEUM↗

Modeling ablator grain structure impacts in ICF implosions

High-density carbon is a leading ablator material for inertial confinement fusion (ICF). This and some other ablator materials have grain structure which is believed to introduce very small-scale (~nm) density inhomogeneity. In principle, such inhomogeneity can affect key ICF metrics like fuel compression and yield, by, for example, acting as a seed for instabilities and inducing mix between ablator and fuel. However, assessments of such effects are uncertain due to the difficulty of modeling this small-scale structure in ICF simulations, typically requiring reduced-resolution modeling that scales these features. We present a grain model and show both the impact of de-resolving grains and the complex mixing dynamics such structures can induce. We find that different methods for de-resolving grains can yield both different total deposition of kinetic energy perturbations and different fuel–ablator mixing. We then show a simple-to-implement approach for approximately conserving the deposition of perturbed kinetic energy and demonstrate that, for the present grain model and test cases, this approach yields a reasonably matched time history of mix width between less and more resolved grain models. The simulations here also demonstrate the complex interaction history between grain-induced mixing and instability around the fuel–ablator interface, showing, for example, that the grain-induced perturbations typically trigger instability of conduction-driven density gradients in the DT fuel, enhancing mix penetration early in the acceleration of the shell. Furthermore simulating both microscale and nanoscale grains, we find initial evidence for larger mixing in the microscale case of the present model, despite smaller deposited kinetic energy perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of ray-based cross-beam energy transfer models

Ray-based cross-beam energy transfer (CBET) models have become a common feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion experiments. Here, these models are necessary for achieving better agreement with experimental measurements, but their detailed implementation can vary widely between the codes and often rely on artificial multipliers. To address this, a series of 2D and 3D test cases has been developed with validated solutions from wave-based calculations. Comparisons of various ray-based CBET models to the wave-based calculations highlight the essential physics that is required for accurate ray-based CBET modeling. Quantitative comparison metrics and/or field data from the wave-based calculations have been made available for use in the validation of other ray-based CBET codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗