Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 415 records · Page 23

Interactions of πK, ππK and KKπ systems at maximal isospin from lattice QCD

We study the interactions of systems of two and three nondegenerate mesons composed of pions and kaons at maximal isospin using lattice QCD, specifically $π^+K^+, π^+π^+K^+$ and $K^+K^+π^+$. Utilizing the stochastic LapH method, we determine the spectrum of these systems on two CLS $N_f$ = 2 + 1 ensembles with pion masses of 200 MeV and 340 MeV, and include many levels in different momentum frames. We constrain the K matrices describing two- and three-particle interactions by fitting the spectrum to the results predicted by the finite-volume formalism, including up to p waves. This requires also results for the $π^+π^+$ and $K^+K^+$ spectrum, which have been obtained previously on the same configurations. We explore different fitting strategies, comparing fits to energy shifts with fits to energies boosted to the rest frame, and also comparing simultaneous global fits to all relevant two- and three-particle channels to those where we first fit twoparticle channels and then add in the three-particle information. We provide the first determination of the three-particle K matrix in $π^+π^+K^+$ and $K^+K^+π^+$ systems, finding statistically significant nonzero results in most cases. We include s and $\textit{p}$ waves in the K matrix for $π^+K^+$ scattering, finding evidence for an attractive p-wave scattering length. We compare our results to Chiral Perturbation Theory, including an investigation of the impact of discretization errors, for which we provide the leading order predictions obtained using Wilson Chiral Perturbation Theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Real-time dynamics of the Schwinger model as an open quantum system with Neural Density Operators

Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle demonstration in a QCD-like theory, the approach is applied to solve the Lindblad master equation in the 1 + 1d lattice Schwinger Model as an open quantum system. Neural Density Operators enable the study of in-medium dynamics on large lattice volumes, where multiple-string interactions and their effects on string-breaking and recombination phenomena can be studied. Thermal properties of the system at equilibrium can also be probed with these methods by variationally constructing the steady state of the Lindblad master equation. Scaling of this approach with system size is studied, and numerical demonstrations on up to 32 spatial lattice sites and with up to 3 interacting strings are performed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Three-particle formalism for multiple channels: the ηππ + $ K\overline{K}\pi $ system in isosymmetric QCD

We generalize previous three-particle finite-volume formalisms to allow for multiple three-particle channels. For definiteness, we focus on the two-channel ηππ and $ K\overline{K}\pi $ system in isosymmetric QCD, considering the positive G parity sector of the latter channel, and neglecting the coupling to modes with four or more particles. The formalism we obtain is thus appropriate to study the b 1 (1235) and η(1295) resonances. The derivation is made in the generic relativistic field theory approach using the time-ordered perturbation theory method. We study how the resulting quantization condition reduces to that for a single three-particle channel when one drops below the upper ($ K\overline{K}\pi $) threshold. We also present parametrizations of the three-particle K matrices that enter into the formalism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Jet angularities in dijet production in proton-proton and heavy-ion collisions at RHIC

We study jet angularities for dijet production at the Relativistic Heavy Ion Collider (RHIC) in proton-proton (pp) and nucleus-nucleus (AA) collisions at 200 GeV nucleon-nucleon center-of-mass collision energy. In particular, we provide NLL resummed predictions for angularity observables of groomed and ungroomed jets produced in pp collisions matched to next-to-leading order QCD calculations resulting in NLO + NLL ' accuracy. Our parton-level predictions are corrected for non-perturbative effects, such as hadronization and underlying event, using parton-to-hadron level transfer matrices obtained with the Sherpa event generator. Furthermore, we use the Q-PYTHIA and JEWEL generators to estimate the impact of the interaction between quarks and gluons produced by the parton shower with the dense medium formed in heavy-ion collisions on the considered jet angularities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

S-matrix bootstrap in 3+1 dimensions: regularization and dual convex problem

The S-matrix bootstrap maps out the space of S-matrices allowed by analyticity, crossing, unitarity, and other constraints. For the 2 → 2 scattering matrix S 2→2 such space is an infinite dimensional convex space whose boundary can be determined by maximizing linear functionals. On the boundary interesting theories can be found, many times at vertices of the space. Here we consider 3 + 1 dimensional theories and focus on the equivalent dual convex minimization problem that provides strict upper bounds for the regularized primal problem and has interesting practical and physical advantages over the primal problem. Its variables are dual partial waves k ℓ (s) that are free variables, namely they do not have to obey any crossing, unitarity or other constraints. Nevertheless they are directly related to the partial waves f ℓ (s), for which all crossing, unitarity and symmetry properties result from the minimization. Numerically, it requires only a few dual partial waves, much as one wants to possibly match experimental results. We consider the case of scalar fields which is related to pion physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

T$ \overline{T} $ deformation in SCFTs and integrable supersymmetric theories

We calculate the \( \mathcal{S} \) -multiplets for two-dimensional Euclidean \( \mathcal{N} \) = (0 , 2) and \( \mathcal{N} \) = (2 , 2) superconformal field theories under the T \( \overline{T} \) deformation at leading order of perturbation theory in the deformation coupling. Then, from these \( \mathcal{N} \) = (0 , 2) deformed multiplets, we calculate two- and three-point correlators. We show the \( \mathcal{N} \) = (0 , 2) chiral ring’s elements do not flow under the T \( \overline{T} \) deformation. Specializing to integrable supersymmetric seed theories, such as \( \mathcal{N} \) = (2 , 2) Landau-Ginzburg models, we use the thermodynamic Bethe ansatz to study the S-matrices and ground state energies. From both an S-matrix perspective and Melzer’s folding prescription, we show that the deformed ground state energy obeys the inviscid Burgers’ equation. Finally, we show that several indices independent of D -term perturbations including the Witten index, Cecotti-Fendley-Intriligator-Vafa index and elliptic genus do not flow under the T \( \overline{T} \) deformation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interactions of two and three mesons including higher partial waves from lattice QCD

We study two- and three-meson systems composed either of pions or kaons at maximal isospin using Monte Carlo simulations of lattice QCD. Utilizing the stochastic LapH method, we are able to determine hundreds of two- and three-particle energy levels, in nine different momentum frames, with high precision. We fit these levels using the relativistic finite-volume formalism based on a generic effective field theory in order to determine the parameters of the two- and three-particle K-matrices. We find that the statistical precision of our spectra is sufficient to probe not only the dominant s-wave interactions, but also those in d waves. In particular, we determine for the first time a term in the three-particle K-matrix that contains two-particle d waves. We use three N f = 2 + 1 CLS ensembles with pion masses of 200, 280, and 340 MeV. This allows us to study the chiral dependence of the scattering observables, and compare to the expectations of chiral perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dispersive sum rules in AdS 2

Dispersion relations for S-matrices and CFT correlators translate UV consistency into bounds on IR observables. In this note, we construct dispersive sum rules for 1D CFTs. We use them to prove bounds on higher-derivative couplings in weakly-coupled non-gravitational EFTs in AdS2. At the leading order in the bulk-point limit, the bounds agree with the flat-space result. We compute the leading universal effect of finite AdS radius on the bounds. Along the way, we give an explicit formula for anomalous dimensions in general higher-derivative contact Witten diagrams in AdS 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Nonperturbative quantum gravity in a closed Lorentzian universe

We study how meaningful physical predictions can arise in nonperturbative quantum gravity in a closed Lorentzian universe. In such settings, recent developments suggest that the quantum gravitational Hilbert space is one-dimensional and real for each α-sector, as induced by spacetime wormholes. This appears to obstruct the conventional quantum-mechanical prescription of assigning probabilities via projection onto a basis of states. While previous approaches have introduced external observers or augmented the theory to resolve this issue, we argue that quantum gravity itself contains all the necessary ingredients to make physical predictions. We demonstrate that the emergence of classical observables and probabilistic outcomes can be understood as a consequence of partial observability: physical observers access only a subsystem of the universe. Tracing out the inaccessible degrees of freedom yields reduced density matrices that encode classical information, with uncertainties exponentially suppressed by the environment’s entropy. We develop this perspective using both the Lorentzian path integral and operator formalisms and support it with a simple microscopic model. Our results show that quantum gravity in a closed universe naturally gives rise to meaningful, robust predictions without recourse to external constructs.

AdS-CFT Correspondence↗

Extracting the asymptotic behavior of S-matrix elements from their phases

The asymptotic kinematic limits of S-matrices are dominated by large logarithms which, roughly speaking, fall into two categories: those which are controlled by a renormalization group (RG) scale, which we may think of as logs involving ratios of invariant mass scales, and those which are functions of ratios of rapidities, so called “rapidity logs”. It has been pointed out by Caron-Huot and Wilhlem [1] that RG anomalous dimension can be extracted from the phase of the S-matrix, which can greatly simplify calculations via unitarity methods. In this paper we generalize the results of [1] to show that the phase can be used to reconstruct rapidity anomalous dimensions, by performing a special type of complex boost. The methodology introduced allows one to calculate without the need for a rapidity regulator which can lead to significant simplifications. We demonstrate the use of this method to derive the rapidity anomalous dimensions in the Sudakov form factor and the two parton soft function at two loops order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The endpoint of partial deconfinement

We study the matrix quantum mechanics of two free hermitian N × N matrices subject to a singlet constraint in the microcanonical ensemble. This is the simplest example of a theory that at large N has a confinement/deconfinement transition. In the microcanonical ensemble, it also exhibits partial deconfinement with a Hagedorn density of states. We argue that the entropy of these configurations, based on a combinatorial counting of Young diagrams, are dominated by Young diagrams that have the VKLS shape. When the shape gets to the maximal depth allowed for a Young diagram of SU(N), namely N, we argue that the system stops exhibiting the Hagedorn behavior. The number of boxes (energy) at the transition is N 2 /4, independent of the charge of the state.

1/N Expansion↗

Acylation as a successful derivatization strategy for the analysis of pinacolyl alcohol in a glycerol-rich matrix by GC-MS: application during an OPCW Proficiency Test

A derivatization protocol based on the acylation of pinacolyl alcohol (PA), an important marker for the nerve agent soman, is presented. The procedure provides a convenient means of detecting, by gas chromatography-mass spectrometry (GC-MS), PA when present at a low concentration in a complex glycerol/alcohol-rich matrix. While there are only two reports describing the specific analysis of PA in matrices at low concentrations, the protocol described herein represents the first of its kind in the analysis of PA in a highly reactive matrix. Two alternative paths for the protocol’s execution are presented. The first involves the direct derivatization of the PA with either acetyl or benzoyl chloride; both reactions yield ester products with significantly different retention times than those of the interferences of the reactive glycerol-rich matrix and in areas of the GC-chromatogram featuring lower levels of matrix interferences. A second procedure involved an initial diethyl ether/aqueous extraction of the matrix; while the extraction was found to substantially remove many of the hydrophilic matrix components and improve the overall derivatization, it also led to some loss of PA available for the derivatization. Both protocols were applied to the successful derivatization and analysis of PA by GC-MS when present at a 5 μg.mL -1 concentration in a glycerol-rich matrix sample administered during the 48th Proficiency Test administered by the Organisation for the Prohibition of Chemical Weapons (OPCW).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mycorrhizal status and host genotype interact to shape plant nutrition in field grown maize ( Zea mays ssp. mays )

Arbuscular mycorrhizal fungi (AMF) establish symbioses with the major cereal crops, providing plants with increased access to nutrients while enhancing their tolerance to toxic heavy metals. However, not all plant varieties benefit equally from this association. In this study, we used quantitative trait loci (QTL) mapping to evaluate the combined effect of host genotypic variation (G) and AMF across 141 genotypes on the concentration of 20 mineral elements in the leaves and grain of field grown maize (Zea mays spp. mays). Our mapping design included selective incorporation of a castor AMF-incompatibility mutation, allowing estimation of AMF, QTL and QTLxAMF effects by comparison of mycorrhizal and non-mycorrhizal plants. Overall, AMF compatibility was associated with higher concentrations of boron (B), copper (Cu), molybdenum (Mo), phosphorus (P), selenium (Se) and zinc (Zn) and lower concentrations of arsenic (As), iron (Fe), magnesium (Mg), manganese (Mn), potassium (K) and strontium (Sr). In addition to effects on individual elements, pairwise correlation matrices for element concentration differed between mycorrhizal and non-mycorrhizal plants. We mapped 22 element QTLs, including 18 associated with QTLxAMF effects that indicate plant genotype-specific differences in the impact of AMF on the host ionome. Although there is considerable interest in AMF as biofertilizers, it remains challenging to estimate the impact of AMF in the field. Our design illustrates an effective approach for field evaluation of AMF effects. Furthermore, we demonstrate the capacity of the ionome to reveal host genotype-specific variation in the impact of AMF on plant nutrition.

59 BASIC BIOLOGICAL SCIENCES↗

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗

Exploiting chitosan to improve the interface of nanocellulose reinforced polymer composites

Cellulose nanofibrils (CNFs) have been widely used as a nanofiller for polymer composite reinforcement due to their excellent mechanical properties. However, CNF is produced in water and needs to be dried prior to use in composite materials. The presence of hydroxyl groups on the surface of CNF creates strong hydrogen bonding that makes it difficult and costly to dry. Additionally, the hydrophilicity at the fiber surface results in agglomeration of CNFs within many polymer matrices. In this study, chitosan (CS) was co-precipitated with CNF to produce a dual-bonding filler for use in poly (lactic acid) (PLA) composites. CS promotes improved interfacial interaction within the polymer matrix by forming strong hydrogen bonds with the CNF and potential covalent bonds with the PLA. Finally, the results confirmed that the addition of a small amount of CS significantly improved the mechanical properties compared to PLA + CNF composites and neat PLA. The detailed study of the PLA + CNF/CS composites reveals the synergetic effect of the hydrogen and covalent bonding mechanism for PLA reinforcement.

59 BASIC BIOLOGICAL SCIENCES↗

Effects of mercerization and fiber sizing of coir fiber for utilization in polypropylene composites

The use of natural fibers as an alternative to synthetic fibers for reinforcing composites is increasing. However, the poor interfacial adhesion between natural fibers and polymer matrices limits their applications. Several approaches have been considered to improve fiber-matrix adhesion via chemical and/or physical treatment. However, the effectiveness of these treatments varies based on the type of fiber, its source, and its composition. Thus, it is imperative to understand the effectiveness of treatment conditions. In this study, we investigated the influence of alkali treatment and fiber sizing on the chemical, thermal, morphological, and mechanical properties of coir fibers and the interface between coir fiber and polypropylene matrix. Further, it was found that using a 5 wt% sodium hydroxide solution for 6 h at room temperature was the optimal treatment condition that led to an improvement in tensile strength by 58%, tensile modulus by 71%, and elongation at break by 37% compared to untreated fibers, and an increment in interfacial shear strength (IFSS) between coir fibers and polypropylene matrix by 32%. The alkali treatment removed the fiber surface impurities, made the fiber surface rough, and enhanced the fiber crystallinity. Sizing of the alkali-treated fiber led to an improvement in IFSS by 87% with no modification of the fiber’s mechanical properties.

36 MATERIALS SCIENCE↗

Improving the Recyclability of Polymer Composites With Cellulose Nanofibrils

Cellulose nanofibers (CNFs) have been widely studied for their reinforcing potential in high-performance composites. While there are numerous publications on CNF-reinforced composites in a variety of polymer matrices, few have considered the recyclability of such thermoplastic composites and whether the incorporation of CNFs deteriorates or improves their performance upon reprocessing. In this study, two thermoplastic resins, poly(lactic acid) (PLA), and glycol-modified polyethylene terephthalate (PETg), were prepared with CNF reinforcement and thermomechanically recycled to investigate the effect of CNF inclusion on the composite properties after reprocessing as well as their effect on the composites’ number of useful life cycles. Additionally, changes in mechanical, thermal, rheological, molecular, and microstructural properties of the composites and/or base resins were monitored as a function of cycle numbers. As is typical, the polymers’ molecular weight and mechanical performance deteriorated with continued processing. However, the addition of spray dried CNF was found to better maintain the mechanical performance of both polymers throughout multiple recycling steps as compared to neat samples. For example, the tensile strength of PETg with 20 wt% CNF after 6 processing cycles was found to exceed that of virgin neat PETg, and higher loadings of CNF were found to preserve a higher yield strength during multiple rounds of reprocessing compared to PETg composites with lower CNF loadings. Ultimately this study indicates that the addition of CNF to some thermoplastic materials can increase both their sustainability by offsetting the use of high-embodied energy resins and their circularity by enabling performance retention over more use cycles.

42 ENGINEERING↗