SEARCH · Search NASA

Results for “interaction potential”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Unexpected Conformational Behavior of Poly(poly(ethylene glycol) methacrylate)-Poly(propylene carbonate)-Poly(poly(ethylene glycol) methacrylate) (PPEGMA-PPC-PPEGMA) Amphiphilic Block Copolymers in Micellar Solution and at the Air-Water Interface

Self-assemblies (bulk micelles and Langmuir monolayers) of a new amphiphilic block copolymer, poly(poly(ethylene glycol) methacrylate)-poly(propylene carbonate)-poly(poly(ethylene glycol) methacrylate) (PPEGMA-PPC-PPEGMA), were investigated in dilute aqueous solution and at the air-water interface. An important feature that distinguishes PPEGMA-PPC-PPEGMA from conventional linear PEG-based polymer surfactants is that in the hydrophilic PPEGMA block the PEG moieties exist as side chains attached to the poly(methacrylate) (PMA) backbone. This study was performed using three PPEGMA-PPC-PPEGMA samples having an identical PPC molecular weight (5.6 kDa) and different PPEGMA molecular weights (7.3, 2.8 and 2.1 kDa on either side) (named as “G7C6G7”, “G3C5G3”, and “G2C6G2”, respectively). Cryo-TEM images revealed that G2C6G2 micelles form large clusters, whereas G7C6G7 micelles exist as unimers; the behavior of G3C6G3 is intermediate. The average core surface distance between two adjacent G2C6G2 micelles in a cluster was estimated to be only about 3.7 nm, which suggests that in the micelle corona layer only the PEG side chains (but not the PMA backbone) exist as hydrated brush chains. A calculation showed that at the small micelle separation distance observed, the van der Waals attractive forces become significant, and are thus likely responsible for the cluster formation. Steady shear rheometry measurements confirmed that G2C6G2 micelles are indeed only weakly bound in clusters, and these micelle clusters are thus prone to breakup even under mild steady shear. The micelle breakup force estimated from the shear viscosity data is in reasonable agreement with an estimate obtained from the theoretical van der Waals interaction potential. Langmuir monolayers formed at the air-water interface by PPEGMA-PPC-PPEGMA were also characterized; surface pressure-area measurements were performed. Detailed features of the surface pressure-area isotherms further supported that only the PEG side chains are hydrated in the subphase, while the PMA backbone remains unhydrated and situated on the water surface.

Lee, Jaewon↗

DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗

DOE OSTI2020

Patchy particle model of hydrated amorphous silica

Glasses corroded in static aqueous conditions typically reach a slow, residual rate of dissolution as a result of the complex interplay between multiple processes. While intrinsically suited to address this problem, Monte Carlo (MC) models developed to date have relied on lattices to describe the amorphous structure of glass and are thus unable to reproduce the residual rate of glass corrosion. A recently developed MC approach based on amorphous structures derived from molecular dynamics simulations offers a solution [Kerisit and Du J. Non-Cryst. Solids 522 (2019) 119601], but it requires a simple model that can rapidly relax the glass-water interface upon each dissolution/condensation event to retain the ability of MC simulations to reach large spatial and temporal scales. To this end, a patchy particle model of hydrated amorphous silica was developed and evaluated against atomistic simulations in this work. Models of water and amorphous silica were first developed separately, and cross terms were then defined to enable simulations of hydrated amorphous silica structures with varying water content. While the nature of its interaction potential prevents the patchy particle model from reproducing simultaneously the structure and dynamics of the systems of interest with a single set of parameters, it described well the connectivity of hydrated amorphous silica structures and the local coordination geometry of individual species. Therefore, this model opens the door to MC simulations of the residual rate of glass corrosion and offers an alternative to atomistic models to investigate the evolution of silica gels over large spatial and temporal scales.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Neutron diffraction discriminates between models for the nanoarchitecture of graphene sheets in glassy carbon

Glassy carbon is a chemically inert, biocompatible, disordered material with the graphene sheet as its basic building block. Structural characterisation techniques have so far been unable to provide definitive distinction between the many proposed models for its structure. Computer based simulation methods have made a step forward by predicting structure from interaction potentials, but the results are sensitive to the choice of the potential. Here we use the white neutron beam of the Spallation Neutron Source at Oak Ridge National Laboratory coupled with the SNAP time-of-flight diffractometer to sample the reciprocal space of glassy carbon and calculate accurate radial distribution functions. From the radial distribution information, we determine graphene sheet dimensions, registration, curvature, and assess defect content to rank the proposed structures in terms of their agreement with experimental data. In this work, we find that the recent models derived from empirical potentials give the best agreement.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Kaon–proton strong interaction at low relative momentum via femtoscopy in Pb–Pb collisions at the LHC

In quantum scattering processes between two particles, aspects characterizing the strong and Coulomb forces can be observed in kinematic distributions of the particle pairs. The sensitivity to the interaction potential reaches a maximum at low relative momentum and vanishing distance between the two particles. Ultrarelativistic heavy-ion collisions at the LHC provide an abundant source of many hadron species and can be employed as a measurement method of scattering parameters that is complementary to scattering experiments. This study confirms that momentum correlations of particles produced in Pb–Pb collisions at the LHC provide an accurate measurement of kaon–proton scattering parameters at low relative momentum, allowing precise access to the K – p $\rightarrow$ K – p process. This work also validates the femtoscopic measurement in ultrarelativistic heavy-ion collisions as an alternative to scattering experiments and a complementary tool to the study of exotic atoms with comparable precision. In this work, the first femtoscopic measurement of momentum correlations of K – p (K + $\bar{p}$) and K + p (K – $\bar{p}$) pairs in Pb–Pb collisions at centre-of-mass energy per nucleon pair of $\sqrt{s_{NN}}$ =5.02 TeV registered by the ALICE experiment is reported. The components of the K – p complex scattering length are extracted and found to be $\mathfrak{R}$f 0 = -0.91±0.03(stat)$^{+0.17}_{-0.03}$(syst) and $\mathfrak{T}$f 0 =0.92±0.05(stat)$^{+0.12}_{-0.33}$(syst). The results are compared with chiral effective field theory predictions as well as with existing data from dedicated scattering and exotic kaonic atom experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Molecular Dynamics Driven by the Many-Body Expansion (MBE-MD)

Here, we present a protocol for classical and nuclear quantum dynamics, in which the energies and forces are generated by the many-body expansion (MBE), and apply it to water clusters using the TTM2.1-F and MB-Pol interaction potentials at various temperatures. We carry out MBE-MD classical and nuclear quantum dynamical simulations in which the energies and forces of the full system are approximated by the 2-, 3-, and 4-body terms of the MBE expansion and compare the average potential and the vibrational density of states (VDOS) with the full simulation, i.e., the one for which no MBE is used. Our results indicate that the thermally averaged potential energy from the MBE up to the 4-body term converges with near identical behavior to the one from the full simulation. The 3-body makes a substantial contribution (~20%) to the energy, whereas the 4- body is necessary for obtaining quantitatively accurate energetics and forces, albeit making a small contribution to each (~2%). We further show that the harmonic frequencies are reproduced to within a few wavenumbers (cm -1 ) at the 4-body level, and that the slowest modes to converge with the MBE rank are those involving the strongest hydrogen bonds. Anharmonicity exacerbates this effect, so that a 4-body description of the energies and forces is needed to achieve accurate anharmonic vibrational frequencies in the hydrogen-bonded OH-stretching region. We also discuss the asymptotic scaling of the MBE-MD protocol with respect to the cost of the underlying potential energy evaluation, suggesting that electronic structure methods that scale at least as N 4 ,N being the size of the system, are needed to result in savings over the traditional full MD simulation. We anticipate that MBE-MD protocol can evolve into a powerful and practical method, which will allow for highly accurate ab initio MD simulations on a much broader range of molecular systems than can be currently handled.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae

Here, we present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5–50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called “thermalized” rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e., the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [Panagiotopoulos, Mol. Phys. 1987, 61, 813−826]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multicomponent systems.

Azeotropes↗

DOE OSTI2021

Exploring the Ion Solvation Environments in Solid-State Polymer Electrolytes through Free-Energy Sampling

The success of poly(ethylene oxide) (PEO) in solid-state polymer electrolytes for lithium-ion batteries is well established. Recently, in order to understand this success and to explore possible alternatives, we studied polyacetal electrolytes to deepen the understanding of the effect of the local chemical structure on ion transport. Advanced molecular dynamics techniques using newly developed, tailored interaction potentials have helped elucidate the various coordination environments of ions in these systems. In particular, the competition between cation-Anion pairing and coordination by the polymer has been explored using free-energy sampling (metadynamics). At equivalent reduced temperatures, with respect to the polymer-specific glass-Transition temperature, two-dimensional free-energy plots reveal the existence of multiple coordination environments for the lithium (Li) ions in these systems and their relative stabilities. Furthermore, we observe that the Li-ion movement in PEO follows a serial, stepwise pathway when moving from one coordination state to another, whereas this happens in a more continuous and concerted fashion in a polyacetal such as poly(1,3-dioxalane) [P(EO-MO)]. The implication is that interconversion between coordination states of the Li ions may be easier in P(EO-MO). However, the overarching observation from our free-energy analysis is that Li-ion coordination is dominated by the polymer (in either case) and contact-ion pairs are rare. We rationalize the observed higher increase in glass-Transition temperature (Tg) with salt loading in polyacetals as due to intermolecular Li-ion coordination involving multiple polymer chains, rather than just one chain for PEO-based electrolytes. This interchain coupling in the polyacetals, resulting in the higher Tg, works against any gains due to variations in Li-ion coordination that might enhance transport processes over PEO. Further research is required to overcome the interdependence between local coordination and macroscopic properties to compete with PEO electrolytes at the same absolute working temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Segmental Dynamics and Vitrification in Associating Copolymer Melts: Role of Cluster Formation, Microdomains, and Cross-Linking

A microscopic statistical mechanical theory of the structure, self-assembly, and activated segmental relaxation is employed to study associating copolymer melts with high attractive sticker fractions, local clustering, and disordered microphase ordering. The stickers are dynamically pinned in a manner that does not affect equilibrium structure which mimics the much slower physical bond breaking process or postassembly cross-linking of sticky monomers. Local sticker clustering and microdomain spatial correlations significantly modify the activated relaxation of nonstickers and glass transition temperature, T g . A re-entrant glass-melting feature is predicted as sticker attraction strength is initially increased corresponding to a speed up of segmental relaxation, and hence reduction of T g relative to the cross-linked homopolymer network. A mechanistic analysis reveals three competing effects: a purely kinetic slowing down of nonstickers down due to cross-linking, disordering of the nonsticker local cage and weakening of effective forces they experience due to sticker physical clustering, and a longer range impact of microdomain scale correlations that results in nonmonotonic dynamical effects. At high enough attraction strength, a qualitative change emerges corresponding to a sticker fraction dependent elevation of T g , which eventually surpasses that of the cross-linked homopolymer network. Dynamically, the new physics arises from a complex evolution of the amplitude of the collective elastic field that dresses the large amplitude mobile segment hopping within a coupled local-nonlocal description of the alpha relaxation. Furthermore, the results are qualitatively consistent with recent experiments on associating PDMS and PPG telechelics of fixed sticker fraction but with chemically different end groups of variable attraction strengths. Finally, possible tests using simulation and the influence of material or model specific interaction potentials and other real world complications are discussed.

Copolymers↗

DOE OSTI2024

Assembly of respiratory syncytial virus matrix protein lattice and its coordination with fusion glycoprotein trimers

Respiratory syncytial virus (RSV) is an enveloped, filamentous, negative-strand RNA virus that causes significant respiratory illness worldwide. RSV vaccines are available, however there is still significant need for research to support the development of vaccines and therapeutics against RSV and related Mononegavirales viruses. Individual virions vary in size, with an average diameter of ~130 nm and ranging from ~500 nm to over 10 µm in length. Though the general arrangement of structural proteins in virions is known, we use cryo-electron tomography and sub-tomogram averaging to determine the molecular organization of RSV structural proteins. We show that the peripheral membrane-associated RSV matrix (M) protein is arranged in a packed helical-like lattice of M-dimers. We report that RSV F glycoprotein is frequently observed as pairs of trimers oriented in an anti-parallel conformation to support potential interactions between trimers. Our sub-tomogram averages indicate the positioning of F-trimer pairs is correlated with the underlying M lattice. These results provide insight into RSV virion organization and may aid in the development of RSV vaccines and anti-viral targets.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Nuclear spin conservation enables state-to-state control of ultracold molecular reactions

Quantum-state control of reactive systems has enabled microscopic probes of underlying interaction potentials and the alteration of reaction rates using quantum statistics. However, extending such control to the quantum states of reaction outcomes remains challenging. Here in this work, we realize this goal by utilizing the conservation of nuclear spins throughout the reaction. Using resonance-enhanced multiphoton ionization spectroscopy to investigate the products formed in bimolecular reactions between ultracold KRb molecules we find that the system retains a near-perfect memory of the reactants’ nuclear spins, manifested as a strong parity preference for the rotational states of the products. We leverage this effect to alter the occupation of these product states by changing the coherent superposition of initial nuclear spin states with an external magnetic field. In this way, we are able to control both the inputs and outputs of a reaction with quantum-state resolution. The techniques demonstrated here open up the possibilities to study quantum entanglement between reaction products, and ultracold reaction dynamics at the state-to-state level.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2021

Downscaling SSP-consistent global spatial urban land projections from 1/8-degree to 1-km resolution 2000–2100

Abstract Long-term, spatial urban land projections that simultaneously offer global coverage and local-scale empirical accuracy are rare. Recently a set of such projections was produced using data-science-based simulations and the Shared Socioeconomic Pathways (SSPs). These projections update at decadal time intervals from 2000 to 2100 with a spatial resolution of 1/8 degree, while many socio-environmental studies customarily run their analysis and modelling at finer spatial resolutions, e.g. 1-km. Here we develop and validate an algorithm to downscale the 1/8-degree spatial urban land projections to the 1-km resolution. The algorithm uses an iterative process to allocate the decadal amount of urban land expansion originally projected for each 1/8-degree grid to its constituent 1-km grids. The results are a set of global maps showing urban land fractions at the 1-km resolution, updated at decadal intervals from 2000 to 2100, under five different urban land expansion scenarios consistent with the SSPs. The data can support studies of potential interactions between future urbanization and environmental changes across spatial and temporal scales.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Metal-induced oxidative stress and human plasma protein oxidation after SARS-CoV-2 infection

Pathogenesis of COVID-19 by SARS-CoV-2 resulted in a global pandemic and public health emergency in 2020. Viral infection can induce oxidative stress through reactive oxygen species (ROS). Inflammation and environmental stress are major sources of oxidative stress after infection. Micronutrients such as iron, copper, zinc, and manganese play various roles in human tissues and their imbalance in blood can impact immune responses against pathogens including SARS CoV-2. We hypothesized that alteration of free metal ions during infection and metal-catalyzed oxidation plays a critical role towards pathogenesis after infection. We analyzed convalescent and hospitalized COVID-19 patient plasma using orthogonal analytical techniques to determine redox active metal concentrations, overall protein oxidation, oxidative modifications, and protein levels via proteomics to understand the consequences of metal-induced oxidative stress in COVID-19 plasma proteins. Metal analysis using ICP-MS showed significantly greater concentrations of copper in COVID-19 plasma compared to healthy controls. We demonstrate significantly greater total protein carbonylation, other oxidative modifications, and deamidation of plasma proteins in COVID-19 plasma compared to healthy controls. Proteomics analysis showed that levels of redox active proteins including hemoglobulin were elevated in COVID-19 plasma. Molecular modeling concurred with potential interactions between iron binding proteins and SARS CoV-2 surface proteins. Overall, increased levels of redox active metals and protein oxidation indicate that oxidative stress-induced protein oxidation in COVID-19 may be a consequence of the interactions of SARS-CoV-2 proteins with host cell metal binding proteins resulting in altered cellular homeostasis.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Magnetic-field-induced Wigner crystallization of charged interlayer excitons in van der Waals heterostructures

Abstract Ever since its inception, coherent excited states of semiconductors have been the focus of semiconductor materials research to evolve into a vibrant field of low-dimensional solid-state physics. The field is gaining new momentum these days due to emerging transdimensional semiconductors such as van der Waals bound layers of transition metal dichalcogenides (TMDs) of controlled thickness. Here, we develop the theory of magnetic-field-induced Wigner crystallization for charged interlayer excitons (CIE) discovered recently in TMD heterobilayers. We derive the ratio of the potential interaction energy to the kinetic energy for the many-particle CIE system in the perpendicular magnetostatic field of an arbitrary strength and predict the crystallization effect in the strong field regime. We show that magnetic-field-induced Wigner crystallization and melting of CIEs can be observed in magneto-photoluminescence experiments with TMD bilayers of systematically varied electron-hole doping concentrations. Our results advance the capabilities of this new generation of transdimensional quantum materials.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Unveiling Structural Heterogeneity and Imbalance of Gold Decahedral Nanoparticles using Four-dimensional Scanning Transmission Electron Microscopy

Multi-twinned structures have been observed in technologically important crystal systems, for example diamond cubic and face-centered cubic (FCC) lattices, that include materials such as diamond, silicon, a wide range of noble metals, and their nanoscale counterparts. Beyond atomic building blocks, the special arrangements also occur in the self-assembly of nanoparticles (NP) and μm-sized colloidal particles and occupy parts of their phase diagrams. Spanning a wide range of length scales, the universality of the structures arises when the systems attempt to achieve multitwinned structures by overcoming geometric misfits during minimizing surface energies with entirely {111} or close-packing facets. While it is fundamental to understand how strain is sustained upon twinned structures and symmetry breaking, the knowledge will be paramount in practical aspects such as guiding and controlling the thin film growth, anisotropic NP growth, and self-assembly of NPs. Au decahedral (Dh) NP, as the most prevalent multi-twinned model system, fits five tetrahedral motifs into a circle by sharing an axis resulting in a geometric misfit angle of 7.35°, or a disclination with power of -7.35°. Postulating how Au FCC lattice adopts the misfit, theoretical models have been developed to address the underlying lattice symmetry and inhomogeneous strain distribution separately. Yet, experimental reports regarding the former have been limited due to the relatively large X-ray beam sizes that do not fit the sizes of NPs. On the other hand, though the latter has been widely adapted in the thermodynamics of small (<10 nm) multi-twinned nanoparticles, previous literature has shown that, at edge length of 17 nm, the theory’s is invalidated by shear strain that is observed in a defect-free Au Dh NP by high-resolution transmission electron microscopy (HR-TEM) imaging. Though the advancement of aberrationcorrected scanning transmission electron microscopy (AC-STEM) imaging and ab initio calculation techniques brings new opportunities, along with challenges in complicated image analysis and limitation in particle size (usually below 10 nm), the gap between nanoscale and mesoscale has never been extended to gain insight from atomic system with straightforward interaction potentials.

4D-STEM↗

DOE OSTI2024

Spinodal enhancement of fluctuations in nucleus-nucleus collisions

Subensemble Acceptance Method (SAM) [1, 2] is an essential link between measured event-by-event fluctuations and their grand canonical theoretical predictions such as lattice QCD. The method allows quantifying the global conservation law effects in fluctuations. In its basic formulation, SAM requires a sufficiently large system such as created in central nucleus-nucleus collisions and sufficient space-momentum correlations. Directly in the spinodal region of the First Order Phase Transition (FOPT) different approximations should be used that account for finite size effects. Thus, we present the generalization of SAM applicable in both the pure phases, metastable and unstable regions of the phase diagram [3]. Obtained analytic formulas indicate the enhancement of fluctuations due to crossing the spinodal region of FOPT and are tested using molecular dynamics simulations. A rather good agreement is observed. Using transport model calculations with interaction potential we show that the spinodal enhancement of fluctuations survives till the later stages of collision via the memory effect [4]. However, at low collision energies the space-momentum correlation is not strong enough for this signal to be transferred to second and third order cumulants measured in momentum subspace. This result agrees well with recent HADES data on proton number fluctuations at $\sqrt{S_{NN}}$ = 2.4 GeV which are found to be consistent with the binomial momentum space acceptance [5].

Poberezhnyuk, Roman↗

DOE OSTI2023

First principles simulation of reacting hypersonic flow over a blunt wedge

This article presents molecular-level analysis of a reactive, near-continuum, Mach 21 nitrogen flow over a blunt wedge using the direct molecular simulation (DMS) method. The flow conditions lead to internal energy excitation and dissociation in the flow field, resulting in thermal and chemical nonequilibrium in the flow. Thermal nonequilibrium in the vibrational mode is observed to extend to the molecular level, where the vibrational energy distributions at various points in the flow field are observed to be non-Boltzmann. Furthermore, this is the first reactive DMS calculation where the wall is assumed to be isothermal and full momentum accommodation of the particles is enforced, hence incorporating viscous wall effects. Since the DMS method uses a quantum mechanically generated interaction potential as its only modeling input, all thermochemical and transport properties of the flow field can directly be attributed to the ab initio potential energy surface. Using the DMS solution as a benchmark, this article assesses the performance of Navier–Stokes computational fluid dynamics solutions using lower fidelity two-temperature models. Two models are chosen as points of comparison: the well-known Park two-temperature model and the recently developed modified Marrone and Treanor model.

Mechanics↗