Kinetic theory of dense-fluid mixtures. IV.
Kinetic prediction theory for viscosity, thermal conductivity and diffusivity of binary liquid mixtures composed of molecules interacting with square-well potential
VISCOUS FLUID↗
SEARCH · Engineering Papers
Results for “potential theory”
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.
Turbulence in electrostatic ion acoustic shocks
Three types of collisionless electrostatic ion acoustic shocks are investigated using a double plasma (DP) device: (1) laminar shocks; (2) small amplitude turbulent shocks in which the turbulence is confined to be upstream of the shock potential jump; and (3) large amplitude turbulent shocks in which the wave turbulence occurs throughout the shock transition. The wave turbulence is generated by ions which are reflected from the shock potential; linear theory spatial growth increments agree with experimental values. The experimental relationship between the shock Mach number and the shock potential is shown to be inconsistent with theoretical shock models which assume that the electrons are isothermal. Theoretical calculations which assume a trapped electron equation of a state and a turbulently flattened velocity distrubution function for the reflected ions yields a Mach number vs potential relationship in agreement with experiment.
Means, R. W.↗
Turbulence in electrostatic ion-acoustic shocks
Three types of collisionless electrostatic ion-acoustic shocks are investigated using the University of California, Los Angeles, double plasma device: (1) laminar shocks; (2) small amplitude turbulent shocks in which the turbulence is confined to be upstream of the shock potential jump; and (3) large amplitude turbulent shocks in which the wave turbulence occurs throughout the shock transition. The wave turbulence is generated by ions which are reflected from the shock potential; linear theory spatial growth increments agree with experimental values. The experimental relationship between the shock Mach number and the shock potential is shown to be inconsistent with theoretical shock models which assume that the electrons are isothermal. Theoretical calculations which assume a trapped electron equation of a state and a turbulently flattened velocity distribution function for the reflected ions yields a Mach number vs potential relationship in agreement with experiment.
Means, R. W.↗
Overview of advanced wing design
Examples of experiment theory correlation are presented to give an indication of the capabilities and limitations of wing design and analysis for transonic applications by potential flow theory. The examples include correlations of experimental pressure distributions with theoretical results from isolated wing codes and wing-body codes. Both conservative and non conservative differencing as well as body and boundary layer corrections are considered. A full potential isolated wing code correlates well with data from an isolated wing test but may give poor prediction of the aerodynamic characteristics of some wing-body configurations. Potential flow wing body codes were found to improve the correlation for the wing-body configurations considered.
Hicks, R. M.↗
Theoretical prediction of thick wing and pylon-fuselage-fanpod-nacelle aerodynamic characteristics at subcritical speeds. Part 1: Theory and results
The theoretical development and the comparison of results with data of a thick wing and pylon-fuselage-fanpod-nacelle analysis are presented. The analysis utilizes potential flow theory to compute the surface velocities and pressures, section lift and center of pressure, and the total configuration lift, moment, and vortex drag. The skin friction drag is also estimated in the analysis. The perturbation velocities induced by the wing and pylon, fuselage and fanpod, and nacelle are represented by source and vortex lattices, quadrilateral vortices, and source frustums, respectively. The strengths of these singularities are solved for simultaneously including all interference effects. The wing and pylon planforms, twists, cambers, and thickness distributions, and the fuselage and fanpod geometries can be arbitrary in shape, provided the surface gradients are smooth. The flow through nacelle is assumed to be axisymmetric. An axisymmetric center engine hub can also be included. The pylon and nacelle can be attached to the wing, fuselage, or fanpod.
Tulinius, J. R.↗
Computational methods for the problems of the tip vortex
The tip vortex problem for helicopter rotors in transonic flow was studied. The first model was that of potential flow past a single blade in parallel transonic flow. The incidence on the wing of a vortex sheet trailing from a blade in front was studied by modeling the incident sheet by a single vortex filament. A 3-D potential flow code was written to incorporate the usual boundary conditions at infinity together with the model of such a vortex filament and to analyze the resulting flow over the wing. The code was used to study the effect of the vortex filament on the production of shocks at the tip of the wing. Another problem studied within the framework of potential flow theory was the rollup of the vortex sheet behind the wing. A comparison was made of flow models described by the potential equation, by the Euler equations for steady flow, and by the Navier-Stokes equations. The question was also studied as to whether there is an advantage in solving the Navier-Stokes equations rather than the Euler equations for steady flow.
Garabedian, Paul R.↗
Multi-site reaction dynamics through multi-fragment density matrix embedding
The practical description of disordered chemical reactions, where the reactions involve multiple species at multiple sites, is presently a challenge using correlated electronic structure methods due to their high computational cost and steep scaling. Here, we describe the gradient theory of multi-fragment density matrix embedding theory, which potentially provides a minimal computational framework to model such processes at the correlated electron level. We present the derivation and implementation of the gradient theory, its validation on model systems and chemical reactions using density matrix embedding, and its application to a molecular dynamics simulation of proton transport in a small water cluster, a simple example of multi-site reaction dynamics.
Chemistry↗
Optimizing Power Generation of a Bottom-Raised Oscillating Surge Wave Energy Converter Using a Theoretical Model: Preprint
Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.
ACE↗
Optimizing Power Generation of a Bottomraised Oscillating Surge Wave Energy Converter Using a Theoretical Model
Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.
ACE↗
Nucleon-nucleon potentials from Δ-full chiral effective-field-theory and implications
We closely investigate NN potentials based upon the Δ-full version of chiral effective-field-theory. We find that recently constructed NN potentials of this kind, which (when applied together with three-nucleon forces) were presented as predicting accurate binding energies and radii for a range of nuclei from A = 16 to A = 132 and providing accurate equations of state for nuclear matter, yield a χ 2 /datum of 60 for the reproduction of the pp data below 100 MeV laboratory energy. This χ 2 is more than three times what the Hamada-Johnston potential of the year of 1962 achieved already some 60 years ago. We perceive this historical fact as concerning in view of the current emphasis on precision. We are able to trace the very large χ 2 as well as the apparent success of the potentials in nuclear structure to unrealistic predictions for P-wave states, in which the Δ-full next-to-next-to-leading order (NNLO) potentials are off by up to 40 times the NNLO truncation errors. In fact, we show that the worse the description of the P-wave states, the better the predictions in nuclear structure. Thus, these potentials cannot be seen as the solution to the outstanding problems in current microscopic nuclear structure physics.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Steady, subsonic, lifting surface theory for wings with swept, partial span, trailing edge control surfaces
A method for computing the lifting pressure distribution on a wing with partial span, swept control surfaces is presented. This method is valid within the framework of linearized, steady, potential flow theory and consists of using conventional lifting surface theory in conjuction with a flap pressure mode. The cause of a numerical instability that can occur during the quadrature of the flap pressure mode is discussed, and an efficient technique to eliminate the instability is derived. This technique is valid for both the flap pressure mode and regular pressure modes and could be used to improve existing lifting surface methods. Examples of the use of the flap pressure mode and comparisons among this method, other theoretical methods, and experiments are given. Discrepancies with experiment are indicated and candidate causes are presented. It is concluded that the method can lead to an efficient and accurate solution of the mathematical problem when a partial span, trailing edge flap is involved.
Medan, R. T.↗
ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials
Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.
boron↗
Potential vorticity transport in weakly and strongly magnetized plasmas
Tangled magnetic fields, often coexisting with an ordered mean field, have a major impact on turbulence and momentum transport in many plasmas, including those found in the solar tachocline and magnetic confinement devices. We present a novel mean field theory of potential vorticity mixing in β-plane magnetohydrodynamic (MHD) and drift wave turbulence. Our results show that mean square stochastic fields strongly reduce Reynolds stress coherence. This decoherence of potential vorticity flux due to stochastic field scattering leads to suppression of momentum transport and zonal flow formation. A simple calculation suggests that the breaking of the shear-eddy tilting feedback loop by stochastic fields is the key underlying physics mechanism. Furthermore, a dimensionless parameter that quantifies the increment in power threshold is identified and used to assess the impact of stochastic field on the L-H transition. We discuss a model of stochastic fields as a resisto-elastic network.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Potential flow in two-dimensional cascades
Potential flow theory for two-dimensional cascades, and application to fluid flow through axial flow compressors
TWO-DIMENSIONAL FLOW↗
Photodissociation Dynamics of CH 2 OO on Multiple Potential Energy Surfaces: Experiment and Theory
UV excitation of the CH 2 OO Criegee intermediate across most of the broad span of the (B 1A') – (X 1A') spectrum results in prompt dissociation to two energetically accessible asymptotes: O ( 1 D) + H 2 CO (X 1 A 1 ) and O ( 3 P) + H 2 CO (a 3A"). Dissociation proceeds on multiple singlet potential energy surfaces that are coupled by two regions of conical intersection (CoIn). Velocity map imaging (VMI) studies reveal a bimodal total kinetic energy (TKER) distribution for the O ( 1 D) + H 2 CO (X 1 A 1 ) products with the major and minor components accounting for ca. 40% and ca. 20% on average of the available energy (E avl ), respectively. The unexpected low TKER component corresponds to highly internally excited H 2 CO (X 1 A 1 ) products accommodating ca. 80% of E avl . Full dimensional trajectory calculations suggest that the bimodal TKER distribution of the O ( 1 D) + H 2 CO (X 1 A 1 ) products originates from two different dynamical pathways: a primary pathway (69%) evolving through one CoIn region to products and a smaller component (20%) sampling both CoIn regions enroute to products. Those that access both CoIn regions likely give rise to the more highly internally excited H 2 CO (X 1 A 1 ) products. The remaining trajectories (11%) dissociate to O ( 3 P) + H 2 CO (a 3 A") products after traversing through both CoIn regions. Here, the complementary experimental and theoretical investigation provides insight on the photodissociation of CH 2 OO via multiple dissociation pathways through two regions of CoIn that control the branching and energy distributions of products.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Physics guided machine learning using simplified theories
Recent applications of machine learning, in particular deep learning, motivate the need to address the generalizability of the statistical inference approaches in physical sciences. In this Letter, we introduce a modular physics guided machine learning framework to improve the accuracy of such data-driven predictive engines. The chief idea in our approach is to augment the knowledge of the simplified theories with the underlying learning process. To emphasize their physical importance, our architecture consists of adding certain features at intermediate layers rather than in the input layer. To demonstrate our approach, we select a canonical airfoil aerodynamic problem with the enhancement of the potential flow theory. We include the features obtained by a panel method that can be computed efficiently for an unseen configuration in our training procedure. By addressing the generalizability concerns, our results suggest that the proposed feature enhancement approach can be effectively used in many scientific machine learning applications, especially for the systems where we can use a theoretical, empirical, or simplified model to guide the learning module.
42 ENGINEERING↗
Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials
Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Active sampling for neural network potentials: Accelerated simulations of shear-induced deformation in Cu–Ni multilayers
Neural network potentials (NNPs) can greatly accelerate atomistic simulations relative to ab initio methods, allowing one to sample a broader range of structural outcomes and transformation pathways. In this work, we demonstrate an active sampling algorithm that trains an NNP that is able to produce microstructural evolutions with accuracy comparable to those obtained by density functional theory, exemplified during structure optimizations for a model Cu–Ni multilayer system. We then use the NNP, in conjunction with a perturbation scheme, to stochastically sample structural and energetic changes caused by shear-induced deformation, demonstrating the range of possible intermixing and vacancy migration pathways that can be obtained as a result of the speedups provided by the NNP. The code to implement our active learning strategy and NNP-driven stochastic shear simulations is openly available at https://github.com/pnnl/Active-Sampling-for-Atomistic-Potentials .
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗