Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling method”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗
DOE OSTI2026

Nonmonotonic-potential description of polarization effects, fusion, and nuclear rainbows in elastic scattering of 6 Li + 12 C at 4.5–600 MeV

The experimental differential cross-section (CS) and analyzing power (AP) data of the 6 Li + 12 C elastic scattering over a wide laboratory energy scale (4.5MeV≤𝐸 lab ≤600MeV) are analyzed within the framework of the optical model (OPM) using nonmonotonic (NM) nucleus–nucleus potentials. The real part of the NM potentials is derived from the Pauli-embodied energy density-functional (EDF) formalism with the sudden approximation. The real part of the noncentral spin-orbit and tensor terms, as well as the imaginary parts, are treated phenomenologically. The effect of the radius of sensitivity on the CS and AP data is found to be more important at lower energies. The diffractive and refractive scattering with Airy structures in the whole angular region of the elastic scattering across the studied energy range is successfully described within the OPM using the NM 6 Li + 12 C potential. The near- and far-side (N and F) decomposition of the total elastic-scattering amplitudes has also been studied using our NM potentials. The evolution of the Airy minima in the angular distributions, coupled with the fitting of the AP data, provides an accurate description of Airy minima of different orders. The OPM calculations with the NM potentials describe exceptionally well the CS, vector analyzing power (VAP), and tensor analyzing power data at 𝐸 lab =9.0,19.24,20,30, and 50MeV. In agreement with our past successful descriptions of CS and the opposite signs of the VAP data for the 6 Li and 7 Li elastic scattering using NM potentials in OPM, the present results appear to provide a better fit, so far, than those obtained from the coupled-channels method. The fusion cross sections of 6 Li + 12 C have been predicted in the energy range (4.5MeV≤𝐸 lab ≤20MeV), fitting the experimental data well in the range 𝐸 lab =2.97–11.77MeV. The EDF potential without any energy dependence and renormalization is also found to describe satisfactorily the experimental CS and AP data at energies up to several hundreds of MeV.

6 ≤ A ≤ 19↗
DOE OSTI2022

Convergent close-coupling calculations of electron scattering on HeH +

Here, we use the molecular convergent close-coupling (MCCC) method to perform calculations of 10–1000 eV electron scattering on the ground state of HeH + . Cross sections are presented for excitation of the n = 2 – 3 singlet and triplet states (where n is the united-atoms-limit principle quantum number), as well as ionization. We also present cross sections for He + and H + ion production following dissociative excitation and ionization. The He + production cross section is compared with the measurements of Lecointre, Jureta, Urbain, and Defrance [J. Lecointre, J. J. Jureta, X. Urbain, and P. Defrance, J. Phys. B At. Mol. Opt. Phys. 47, 015203 (2014)]. We find that the MCCC results are up to 30% higher than experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Elastic scattering and rotational excitation of H 2 by electron impact: Convergent close-coupling calculations

We apply the adiabatic-nuclei molecular convergent close-coupling (MCCC) method to the study of elastic scattering and rotational excitation of H 2 by 0.01–20-eV electrons. Integral cross sections are presented for all rotational transitions with |Δ⁢N| ≤ 2 and N i = 0–31 within the v = 0 vibrational level and differential cross sections for a selection of transitions. Agreement with the available measurements and previous calculations is mostly excellent, depending on the transition and incident energy. Finally, we suggest possible reasons, and argue for the accuracy of the MCCC data.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Charge radii of exotic neon and magnesium isotopes

We compute the charge radii and ground-state energies of even-mass neon and magnesium isotopes from neutron number N = 8 to the dripline. Our calculations are based on nucleon-nucleon and three-nucleon potentials from chiral effective field theory that include Δ isobars. These potentials yield an accurate saturation point and symmetry energy of nuclear matter. We use the coupled-cluster method and start from an axially symmetric reference state. Binding energies and two-neutron separation energies largely agree with data, and the dripline in neon is accurate. The computed charge radii are accurate for many isotopes where data exist. Finer details, such as isotope shifts, however, are not accurately reproduced. These chiral potentials indicate a subshell closure at N = 14 for the radii (but not for two-neutron separation energies) and a decrease in charge radii at N = 8 (observed in neon and predicted for magnesium). Furthermore, they yield a continued increase of charge radii as neutrons are added beyond N = 14 yet underestimate the large increase at N = 20 in magnesium.

20 ≤ A ≤ 38↗
DOE OSTI2023

Entanglement entropy of nuclear systems

We study entanglement entropies between the single-particle states of the hole space and its complement in nuclear systems. Analytical results based on the coupled-cluster method show that entanglement entropies are proportional to the particle number fluctuation and the depletion number of the hole space for sufficiently weak interactions. General arguments also suggest that the entanglement entropy in nuclear systems fulfills a volume instead of an area law. Here, we test and confirm these results by computing entanglement entropies of the pairing model and neutron matter, and the depletion number of finite nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Theoretical methods to design and test quantum simulators for the compact Abelian Higgs model

The lattice compact Abelian Higgs model is a nonperturbative regularized formulation of low-energy scalar quantum electrodynamics. In 1+1 dimensions, this model can be quantum simulated using a ladder-shaped optical lattice with Rydberg-dressed atoms. In this setup, one spatial dimension is used to carry the angular momentum of the quantum rotors. One can use truncations corresponding to spin-2 and spin-1 to build local Hilbert spaces associated with the links of the lattice. We argue that ladder-shaped configurable arrays of Rydberg atoms can be used for the same purpose. We make concrete proposals involving two and three Rydberg atoms to build one local spin-1 space (a qutrit). We show that the building blocks of the Hamiltonian calculations are models with one and two spins. We compare target and simulators using perturbative and numerical methods. The two-atom setup provides an easily controllable simulator of the one-spin model while the three-atom setup involves solving nonlinear equations. We discuss approximate methods to couple two spin-1 spaces. The article provides analytical and numerical tools necessary to design and build the proposed simulators with current technology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Dynamics of Density Fluctuations in Atomic Nuclei

We study the spatiotemporal patterns of density fluctuations in 16,24O and 48Ca using nuclear interactions from chiral effective field theory and the time-dependent coupled-cluster method. We find that two-particle–two-hole excitations generate small-amplitude fluctuations that are fast, short-ranged, and of stochastic character.

Bonaiti, Francesca [ORNL] (ORCID:0000000239261609)↗
DOE OSTI2022

Start-up Optimization Considering Integrated Power Electronic Systems

The startup of systems composed of power electronic systems is becoming more complex. As the systems grow, standard pre-configuration sequencing methods will struggle. This work proposes a linear programming optimization method to couple power electronic systems in a series of startup commands. This formulation is demonstrated in a real-time implementation on hardware.

Starke, Michael↗
DOE OSTI2023

Large Eddy Simulation of Gasoline Sprays in a Lagrangian–Eulerian Framework Using the High-Order Spectral Element Method

Predicting the spray evolution using simulations requires accurate modeling of the turbulent gas-phase flow field. Here, in this study, the high-order spectral-element method (SEM), implemented in the code Nek5000, was used to provide highly resolved solutions to the turbulent flow field. Spray modeling capabilities were implemented into the Nek5000 code. The spray is modeled in a Lagrangian–Eulerian (LE) framework, where the liquid is represented by discrete parcels of droplets. The method for coupling liquid and gas in the context of SEM is described, which allows for very fine meshes to be used without affecting the stability of the solution. Large-eddy simulations (LES) of the eight-hole ECN Spray G gasoline injector were conducted. Numerical results are compared against experimental data for liquid penetration, droplet size and gas velocity. The morphology of the multiplume spray is compared against experimental data. The effect of different spray injection inputs is analyzed. It was found that using a plume direction of 33 deg and an injection cone angle of 30 deg produced the best results overall. This work shows the applicability of SEM for spray modeling applications, where use of a high-order flow solver can help us understand the multiplume spray aerodynamics and how it leads to plume collapse under certain conditions. Results also highlight the need for tuning spray input parameters in the LE framework, even when high-fidelity gas flow solutions are possible.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

NuLattice: Ab initio computations of atomic nuclei on lattices

Here, we introduce NuLattice, a Python software package for ab initio computations of atomic nuclei on lattices. The computational tools consist of Hartree Fock, the coupled cluster method, the in-medium similarity renormalization group, and full configuration interaction. At present, the employed interactions are from pion-less effective field theory at leading order and consist of two-body and three-body contacts. We present results for light nuclei 2 H, 3,4 He, 8 Be, 12 C, and 16 O. NuLattice algorithms exploit the sparsity and locality of lattice interactions, and as a result computations can be run on laptops.

Rothman, Maxwell [Univ. of Tennessee, Knoxville, T↗
DOE OSTI2020

Modeling synergistic fuel cell membrane degradation with mitigating effects of cerium

During operation, polymer-electrolyte fuel cells undergo mechanical and chemical degradation, which behave synergistically and lead to accelerated membrane deterioration over time. This study builds upon previous modeling work on mechanical degradation, as described by a pinhole in the membrane, and chemical degradation mitigation by using cerium. By combining these two models, analysis can be carried out on the coupled degradation methods and how the mitigation effects of cerium disrupt the degradation cycle. The model results show how the presence of a pinhole in the membrane changes the distribution of cerium in the cell and modifies chemical degradation rates. In addition, the model shows how cerium affects the rate of pinhole growth by slowing down the rate of change in mechanical properties due to chemical degradation. Finally, the model shows how cerium modifies the mechanical and chemical degradation rates of the membrane under humidity- and voltage-cycling conditions.

Ehlinger, VM↗
DOE OSTI2020

Modeling Synergistic Fuel Cell Membrane Degradation with Mitigating Effects of Cerium

During operation, polymer-electrolyte-membrane (PEM) fuel cells undergo mechanical and chemical degradation mechanisms, which behave synergistically and lead to accelerated membrane degradation over time. This study builds upon previous modeling work on mechanical degradation as described by a pinhole in the membrane and the effects of cerium on chemical degradation.1 By combining these two models, analysis can be carried out on the coupled degradation methods and how the mitigation effects of cerium disrupt the degradation cycle. The mechanical model represents a pinhole in the membrane using an effective void fraction, which allows for increased gas crossover through the membrane and impacts the membrane transport and mechanical properties. A microkinetic model for the chemical degradation is included in the model, including attack of the membrane polymer by hydroxyl radicals as well as the mitigation reaction for quenching of hydroxyl radicals with cerium. A concentrated solution theory approach is used to model the transport of cerium ions throughout the cell.2 The model results show how the location of cerium in the cell can be used to prevent chemical degradation. In addition, the model shows how cerium slows down the rate of pinhole growth by reducing the gas crossover and membrane thinning rate. Finally, the model can be used to optimize the distribution of cerium in the membrane and catalyst layers by balancing trade-offs between lowering degradation rates and decreasing fuel cell performance.

25 ENERGY STORAGE↗
DOE OSTI2022

Xyce(™) Parallel Electronic Simulator v.7.5

The Xyce Parallel Electronic Simulator simulates electronic circuit behavior in DC, AC, HB, MPDE and transient mode using standard analog (DAE) and/or device (PDE) device models including several age and radiation aware devices. It supports a variety of computing platforms (both serial and parallel) computers. Lastly, it uses a variety of modern solution algorithms dynamic parallel load-balancing and iterative solvers.! ! Xyce is primarily used to simulate the voltage and current behavior of a circuit network (a network of electronic devices connected via a conductive network). As a tool, it is mainly used for the design and analysis of electronic circuits.! ! Kirchoff's conservation laws are enforced over a network using modified nodal analysis. This results in a set of differential algebraic equations (DAEs). The resulting nonlinear problem is solved iteratively using a fully coupled Newton method, which in turn results in a linear system that is solved by either a standard sparse-direct solver or iteratively using Trilinos linear solver packages, also developed at Sandia National Laboratories.

Source record↗
DOE OSTI2023

A Novel Method for Diagnosing Land–Atmosphere Coupling Sensitivity in a Single-Column Model

The response of boundary layer properties and cloudiness to changes in surface evaporative fraction (EF) is investigated in a single-column model to quantify the locally coupled impact of subgrid surface variations on the atmosphere during summer. Sensitive coupling days are defined when the model atmosphere exhibits large variations across a range of EFs centered on the analyzed value. Coupling sensitivity exists as both positive feedback (cloudiness increases with EF) and negative feedback (clouds increase with decreasing EF) regimes. The positive regime manifests in shallow convection situations, which are capped by a strengthened inversion and subsidence, restricting the vertical extent of convection to just above the boundary layer. Surfaces with larger EF (greater surface latent heat flux) can inject more moisture into the vertically confined system, lowering the cloud base and an increasing cloud liquid water path (LWP). Negative feedback regimes tend to manifest when large-scale deep convection, such as from mesoscale convective systems and fronts, is advected through the domain, where convection strengthens over surfaces with a lower EF (greater surface sensible heat flux). The invigoration of these systems by the land surface leads to an increase in LWP through strengthened updrafts and stronger coupling between the boundary layer and the free atmosphere. These results apply in the absence of heterogeneity-induced mesoscale circulations, providing a one-dimensional dynamical perspective on the effect of surface heterogeneity. Finally, this study provides a framework of intermediate complexity, lying between parcel theory and high-resolution coupled land–atmosphere modeling, and therefore isolates the relevant first-order processes in land–atmosphere interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗