Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase interpolation”

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 127 records · Page 7

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗

Time-dependent wave-packet method for the complete determination of S-matrix elements for reactive molecular collisions in three dimensions

An alternative time-dependent wave-packet method for treating three-dimensional gas phase reactive atom-diatom collisions is presented. The method employs a nonreactive body-frame wave packet propagation procedure, made possible by judicious use of absorbing optical potentials, a novel scheme for interpolating the wave function from coordinates in one arrangement to those in another and the fact that the time-dependent Schroedinger equation is an initial-value problem. The last feature makes possible a computationally viable and accurate procedure for changing from one arrangement's coordinates to another. In addition, the method allows the determination of S-matrix elements over a wide range of energies from a single wave-packet propagation. The method is illustrated by carrying out detailed calculations of inelastic and reactive scattering in the H + H2 system using the Liu-Siegbahn-Truhlar-Horowitz potential surface.

Judson, Richard S.↗

Identification of Bare-Airframe Dynamics from Closed-Loop Data Using Multisine Inputs and Frequency Responses

Amethod is presented for computing multiple-input multiple-output frequency responses of bare-airframe dynamics for systems excited using orthogonal phase-optimized multisines and including correlated data arising from control mixing or feedback control. The estimation was posed as the solution to an underdetermined system of linear equations, for which additional information was supplied using interpolation of the frequency responses. A simulation model of the NASA T-2 aircraft having two inputs and two outputs was used to investigate the method in the open-loop configuration and under closed-loop control. The method was also applied to flight test data from the X-56A aeroelastic demonstrator having five inputs and ten outputs and flying under closed-loop control with additional control allocation mixing. Results demonstrated that the proposed method accurately estimates the bare airframe frequency responses in the presence of correlated data from control mixing and feedback control. Results also agreed with estimates obtained using different methods that are less sensitive to correlated inputs.

Grauer, Jared A.↗

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗

Dislocation dynamics in heterogeneous nanostructured materials

Crystalline materials can be strengthened by introducing dissimilar phases that impede dislocation glide. At the same time, the changes in microstructure and chemistry usually make the materials less ductile. One way to circumvent the strength–ductility dilemma is to take advantage of heterogeneous nanophases which simultaneously serve as dislocation barriers and sources. Owing to their superior mechanical properties, heterogeneous nanostructured materials (HNMs) have attracted a lot of attention worldwide. Nevertheless, it has been difficult to characterize dislocation dynamics in HNMs using classical continuum models, mainly due to the challenges in describing the elastic and plastic heterogeneity among the phases. Here, in this work, we advance a phase-field dislocation dynamics (PFDD) model to treat multi-phase materials, consisting of phases differing in composition, structural order, and size in the same system. We then apply the advanced PFDD model to exploring two important but divergent materials design problems in HNMs: dislocation/obstacle interactions and dislocation/interface interactions. Results show that the interactions between a dislocation and distribution of obstacles varying in structure and composition cannot be understood by simply interpolating from their individual interactions with a dislocation. It is also found that materials containing interfaces with nanoscale thicknesses and compositional gradients have a much higher dislocation bypass stress than those with sharp interfaces, providing an explanation for the simultaneous high strength and toughness of thick interface-containing nanolaminates as observed in recent experiments.

36 MATERIALS SCIENCE↗

Thermodynamic modeling of CsF with LiF-NaF-KF for molten fluoride-fueled reactors

Gibbs energy models were developed to describe the thermochemical behavior of CsF in molten FLiNaK (46.5LiF-11.5NaF-42KF mol%), a proposed molten salt reactor (MSR) fuel solvent and coolant, as cesium is of concern due to its high radiotoxicity and volatility. Initially, it was necessary to obtain a more accurate Gibbs energy function for CsF which required fitting parameters to reported vapor pressures over condensed phase CsF. The pseudo-binary systems CsF-LiF, CsF-NaF and CsF-KF were then evaluated utilizing phase equilibria and enthalpy of mixing (Δ mix H) values, together with original differential scanning calorimetry (DSC) measurements performed for the CsF-LiF and CsF-KF systems. The CsF-LiF-NaF, CsF-LiF-KF and CsF-NaF-KF pseudo-ternary system representations were obtained by interpolation of the constituent pseudo-binary systems, with DSC measurements performed for the CsF-LiF-NaF system to corroborate the calculated liquidus temperature. Ultimately, the pseudo-ternary systems were interpolated to obtain Gibbs energy models for the pseudo-quaternary CsF-LiF-NaF-KF system, supported by DSC measurements at low CsF compositions (1–10 mol%), yielding computed equilibria and cesium-containing vapor pressures that compare favorably with reported values. In conclusion, the Molten Salt Thermal Properties Database – Thermochemical (MSTDB-TC) was subsequently expanded to include these Gibbs energy models allowing description of the thermochemical behavior of the CsF-LiF-NaF-KF system.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

FLUID: A numerical interpolation procedure for obtaining thermodynamic and transport properties of fluids

A computer program subroutine, FLUID, was developed to calculate thermodynamic and transport properties of pure fluid substances. It provides for determining the thermodynamic state from assigned values for temperature-density, pressure-density, temperature-pressure, pressure-entropy, or pressure-enthalpy. Liquid or two-phase (liquid-gas) conditions are considered as well as the gas phase. A van der Waals model is used to obtain approximate state values; these values are then corrected for real gas effects by model-correction factors obtained from tables based on experimental data. Saturation conditions, specific heat, entropy, and enthalpy data are included in the tables for each gas. Since these tables are external to the FLUID subroutine itself, FLUID can implement any gas for which a set of tables has been generated. (A setup phase is used to establish pointers dynamically to the tables for a specific gas.) Data-table preparation is described. FLUID is available in both SFTRAN and FORTRAN

Fessler, T. E.↗

Tidal estimation in the Atlantic and Indian Oceans, 3 deg x 3 deg solution

An estimation technique was developed to extrapolate tidal amplitudes and phases over entire ocean basins using existing gauge data and the altimetric measurements provided by satellite oceanography. The technique was previously tested. Some results obtained by using a 3 deg by 3 deg grid are presented. The functions used in the interpolation are the eigenfunctions of the velocity (Proudman functions) which are computed numerically from a knowledge of the basin's bottom topography, the horizontal plan form and the necessary boundary conditions. These functions are characteristic of the particular basin. The gravitational normal modes of the basin are computed as part of the investigation; they are used to obtain the theoretical forced solutions for the tidal constituents. The latter can provide the simulated data for the testing of the method and serve as a guide in choosing the most energetic functions for the interpolation.

Sanchez, Braulio V.↗

Interfacial tension and vapor-liquid equilibria in the critical region of mixtures

In the critical region, the concept of two-scale-factor universality can be used to accurately predict the surface tension between near-critical vapor and liquid phases from the singularity in the thermodynamic properties of the bulk fluid. In the present work, this idea is generalized to binary mixtures and is illustrated using the data of Hsu et al. (1985) for CO2 + n-butane. The pressure-temperature-composition-density data for coexisting, near-critical phases of the mixtures are fitted with a thermodynamic potential comprised of a sum of a singular term and nonsingular terms. The nonuniversal amplitudes characterizing the singular term for the mixtures are obtained from the amplitudes for the pure components by interpolation in a space of thermodynamic 'field' variables. The interfacial tensions predicted for the mixtures from the singular term are within 10 percent of the data on three isotherms in the pressure range (Pc - P)/Pc of less than 0.5. This difference is comparable to the combined experimental and model errors.

Moldover, Michael R.↗

A Priori Subgrid Scale Modeling for a Droplet Laden Temporal Mixing Layer

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using a direct numerical simulation (DNS) database. The DNS is for a Reynolds number (based on initial vorticity thickness) of 600, with droplet mass loading of 0.2. The gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. Since Large Eddy Simulation (LES) of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be given by the filtered variables plus a correction based on the filtered standard deviation, which can be computed from the sub-grid scale (SGS) standard deviation. This model predicts unfiltered variables at droplet locations better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: Smagorinsky, gradient and scale-similarity. When properly calibrated, the gradient and scale-similarity methods give results in excellent agreement with the DNS.

Okongo, Nora↗

Apparently superluminal superfluids

We consider the superfluid phase of a specific renormalizable relativistic quantum field theory. We prove that, within the regime of validity of perturbation theory and of the superfluid effective theory, there are consistent and regular vortex solutions where the superfluid’s velocity field as traditionally defined smoothly interpolates between zero and arbitrarily large superluminal values. We show that this solution is free of instabilities and of superluminal excitations. We show that, in contrast, a generic vortex solution for an ordinary fluid does develop an instability if the velocity field becomes superluminal. All this questions the characterization of a superfluid velocity field as the actual velocity of “something”.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Holographic BCFT with a Defect on the End-of-the-World brane

In this paper, we propose a new gravity dual for a 2d BCFT with two conformal boundaries by introducing a defect that connects the two End-of-the-World branes. We demonstrate that the BCFT dual to this bulk model exhibits a richer lowest spectrum. The corresponding lowest energy eigenvalue can continuously interpolate between - πc 24 Δ x and 0 where Δ x is the distance between the boundaries. This range was inaccessible to the conventional AdS/BCFT model with distinct boundary conditions. We compute the holographic entanglement entropy and find that it exhibits three different phases, one of which breaks the time reflection symmetry. We also construct a wormhole saddle, analogous to a 3d replica wormhole, which connects different boundaries through the AdS bulk. This saddle is present only if the BCFT is non-unitary and is always subdominant compared to the disconnected saddle.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Sparsely Sampled Phase-Insensitive Ultrasonic Transducer Arrays

Three methods of interpretation of outputs from sparsely sampled two-dimensional array of receiving ultrasonic transducers used in transmission experiments investigated. Methods are: description of sampled beam in terms of first few spatial moments of sampled distribution of energy; use of signal-dependent cutoff to limit extent of effective receiver aperture; and use of spatial interpolation to increase apparent density of sampling during computation. Methods reduce errors in computations of shapes of ultrasonic beams.

Patrick H. Johnston↗

Performance evaluation of the JPL interim digital SAR processor

The performance of the Interim Digital SAR Processor (IDP) was evaluated. The IDP processor was originally developed for experimental processing of digital SEASAT SAR data. One phase of the system upgrade which features parallel processing in three peripheral array processors, automated estimation for Doppler parameters, and unsupervised image pixel location determination and registration was executed. The method to compensate for the target range curvature effect was improved. A four point interpolation scheme is implemented to replace the nearest neighbor scheme used in the original IDP. The processor still maintains its fast throughput speed. The current performance and capability of the processing modes now available on the IDP system are updated.

Wu, C.↗

Planetary Waves

The mean temperature field at various pressure levels for January and July for both hemispheres are given. In summer (July in the Northern Hemisphere, January in the Southern Hemisphere) the fields are nearly zonally symmetric, but in winter large longitudinal variations are evident. These are mainly of low wave number, i.e., they may be represented by Fourier analysis around the globe using just a few (one or two) waves. This is the basis for representing the climatology of longitudinal variation in terms of wave components, since the fields may be defined with fewer values than by using a grid in longitude that has a sufficiently small interval to adequately represent the smooth variations. The amplitude and phase of temperature and geopotential height for wave numbers one and two, with In(pressure) as the vertical coordinate are given. The fields were calculated and plotted at pressure intervals of 0.2 in In(pressure) and at latitude intervals of 4 deg, but were interpolated to intervals of 0.5 in In(pressure) (approximately 3.5 km) and 10 deg latitude for tabulation. Tables and figures give wave coefficients for the monthly mean temperature fields. Thus they represent the quasi-statinary planetary waves.

Barnett, J. J.↗

A Priori Subgrid Analysis of Temporal Mixing Layers with Evaporating Droplets

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using three sets of results from a Direct Numerical Simulation (DNS) database, with Reynolds numbers (based on initial vorticity thickness) as large as 600 and with droplet mass loadings as large as 0.5. In the DNS, the gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. The Large Eddy Simulation (LES) equations corresponding to the DNS are first derived, and key assumptions in deriving them are first confirmed by computing the terms using the DNS database. Since LES of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be the sum of the filtered variables and a correction based on the filtered standard deviation; this correction is then computed from the Subgrid Scale (SGS) standard deviation. This model predicts the unfiltered variables at droplet locations considerably better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: the Smagorinsky approach, the Gradient model and the Scale-Similarity formulation. When the proportionality constant inherent in the SGS models is properly calculated, the Gradient and Scale-Similarity methods give results in excellent agreement with the DNS.

Okongo, Nora↗

Real-Time Estimation of Bare-Airframe Frequency Responses from Closed-Loop Data and Multisine Inputs

A method is presented for computing frequency responses of multiple-input multiple-output bare-airframe dynamics from flight test data containing feedback control and/or mixing of control effectors. Orthogonal phase-optimized multisines are used to simultaneously excite each input with unique harmonic frequencies, at which frequency responses are computed as ratios of output-to-input Fourier transform data. The confounding effects of feedback and mixing for frequency response estimation are resolved by interpolating the frequency responses among all the harmonic frequencies. The method can be run in batch for post-flight analysis, or in real time as the aircraft is flying. The effectiveness of the method was verified using closed-loop simulations of the subscale NASA T-2 generic transport airplane. The method was also demonstrated using flight test data from the X-56A MUTT aeroelastic airplane, which was flown with feedback control and mixing.

Jared A Grauer↗