Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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 1,351 records · Page 75

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗

ESPPU INPUT: C$^3$ within the "Linear Collider Vision"

The Linear Collider Vision calls for a Linear Collider Facility with a physics reach from a Higgs Factory to the TeV-scale with $e^+e^{-}$ collisions. One of the technologies under consideration for the accelerator is a cold-copper distributed-coupling linac capable of achieving high gradient. This technology is being pursued by the C$^3$ collaboration to understand its applicability to future colliders and broader scientific applications. In this input we share the baseline parameters for a C$^3$ Higgs-factory and the energy reach of up to 3 TeV in the 33 km tunnel foreseen under the Linear Collider Vision. Recent results, near-term plans and future R&D needs are highlighted.

43 PARTICLE ACCELERATORS↗

A pseudo-sound constitutive relationship for the dilatational covariances in compressible turbulence: An analytical theory

The mathematical consequences of a few simple scaling assumptions about the effects of compressibility are explored using a simple singular perturbation idea and the methods of statistical fluid mechanics. Representations for the pressure-dilation and dilatational dissipation covariances appearing in single-point moment closures for compressible turbulence are obtained. While the results are expressed in the context of a second-order statistical closure they provide some interesting and very clear physical metaphors for the effects of compressibility that have not been seen using more traditional linear stability methods. In the limit of homogeneous turbulence with quasi-normal large-scales the expressions derived are - in the low turbulent Mach number limit - asymptotically exact. The expressions obtained are functions of the rate of change of the turbulence energy, its correlation length scale, and the relative time scale of the cascade rate. The expressions for the dilatational covariances contain constants which have a precise and definite physical significance; they are related to various integrals of the longitudinal velocity correlation. The pressure-dilation covariance is found to be a nonequilibrium phenomena related to the time rate of change of the internal energy and the kinetic energy of the turbulence. Also of interest is the fact that the representation for the dilatational dissipation in turbulence, with or without shear, features a dependence on the Reynolds number. This article is a documentation of an analytical investigation of the implications of a pseudo-sound theory for the effects of compressibility.

Ristorcelli, J. R.↗

Linear in temperature resistivity in the limit of zero temperature from the time reparameterization soft mode

The most puzzling aspect of the ‘strange metal’ behavior of correlated electron compounds is that the linear in temperature resistivity often extends down to low temperatures, lower than natural microscopic energy scales. We consider recently proposed deconfined critical points (or phases) in models of electrons in large dimension lattices with random nearest-neighbor exchange interactions. The criticality is in the class of Sachdev–Ye–Kitaev models, and exhibits a time reparameterization soft mode representing gravity in dual holographic theories. We compute the low temperature resistivity in a large M limit of models with SU(M) spin symmetry, and find that the dominant temperature dependence arises from this soft mode. The resistivity is linear in temperature down to zero temperature at the critical point, with a co-efficient universally proportional to the product of the residual resistivity and the co-efficient of the linear in temperature specific heat. We argue that the time reparameterization soft mode offers a promising and generic mechanism for resolving the strange metal puzzle.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

EpitopeVec: linear epitope prediction using deep protein sequence embeddings

Motivation: B-cell epitopes (BCEs) play a pivotal role in the development of peptide vaccines, immuno-diagnostic reagents and antibody production, and thus in infectious disease prevention and diagnostics in general. Experimental methods used to determine BCEs are costly and time-consuming. Therefore, it is essential to develop computational methods for the rapid identification of BCEs. Although several computational methods have been developed for this task, generalizability is still a major concern, where cross-testing of the classifiers trained and tested on different datasets has revealed accuracies of 51–53%. Results: We describe a new method called EpitopeVec, which uses a combination of residue properties, modified antigenicity scales, and protein language model-based representations (protein vectors) as features of peptides for linear BCE predictions. Extensive benchmarking of EpitopeVec and other state-of-the-art methods for linear BCE prediction on several large and small datasets, as well as cross-testing, demonstrated an improvement in the performance of EpitopeVec over other methods in terms of accuracy and area under the curve. As the predictive performance depended on the species origin of the respective antigens (viral, bacterial and eukaryotic), we also trained our method on a large viral dataset to create a dedicated linear viral BCE predictor with improved cross-testing performance.

59 BASIC BIOLOGICAL SCIENCES↗

Streaming Instability with Multiple Dust Species – I. Favourable Conditions for the Linear Growth

A recent study suggests that the streaming instability, one of the leading mechanisms for driving the formation of planetesimals, may not be as efficient as previously thought. Under some disc conditions, the growth time-scale of the instability can be longer than the disc lifetime when multiple dust species are considered. To further explore this finding, we use both linear analysis and direct numerical simulations with gas fluid and dust particles tomutually validate and study the unstablemodes of the instability in more detail.We extend the previously studied parameter space by one order of magnitude in both the range of the dust-size distribution [Ts, min, Ts, max] and the total solid-to-gas mass ratio ε and introduce a third dimension with the slope q of the size distribution. We find that the fast-growth regime and the slow-growth regime are distinctly separated in the ε–Ts, max space, while this boundary is not appreciably sensitive to q or Ts, min.With a wide range of dust sizes present in the disc (e.g. Ts, min  10−3), the growth rate in the slow-growth regime decreases as more dust species are considered.With a narrow range of dust sizes (e.g. Ts, max/Ts, min = 5), on the other hand, the growth rate in most of the ε–Ts, max space is converged with increasing dust species, but the fast and the slow growth regimes remain clearly separated. Moreover, it is not necessary that the largest dust species dominate the growth of the unstable modes, and the smaller dust species can affect the growth rate in a complicated way. In any case, we find that the fast-growth regime is bounded by ε  1 or Ts, max  1, which may represent the favourable conditions for planetesimal formation.

Zhaohuan Zhu↗

Coronal heating by the resonant absorption of Alfven waves: Wavenumber scaling laws.

The importance of global modes in coronal loop heating is well established. In the present work the scaling of the global-mode resonant heating rate with the perturbation wavenumbers is studied with the numerical solution of the linearized time-dependent MHD equations for a full compressible, low-beta, resistive plasma using an implicit integration scheme. The numerical simulations demonstrate that the dissipation on inhomogeneties in the background Alfven speed occurs in narrow resonant layer with the highest heating rate at the global-mode frequency. The global-mode heating rate H (sub r) was found to scale as H (Sub r) approximately k (sub y) (exp 1.03) when k (sub z) = 0.1, and as H (sub r) approximately k (sub y) (exp -1.93) when k (sub z) = 0.75, where k (sub y) and k (sub z) are the wavenumbers in the perpendicular and parallel to the magnetic field directions, respectively, while the dependence of H (sub r) on k (sub z) is more complex. The quality factor Q of the MHD resonance cavity scales as Q approximately k (sub y) (exp -1.8) for k (sub z) = 0.75 and as Q approximately k (sub y) (exp -1.46) for k (sub z) = 0.1. The numerically determined heating rate scaling, the global-mode fequency, and the quality factor are in good agreement with the analytical linear theory. The magnitude of the perturbed velocities was found to decrease with k (sub y). Assuming typical coronal loop parameters (B (sub 0) = 100-200 G, upsilon (sub A) = 2000-4000 km/s), the Alfven waves can supply the required heating to a low-Q loops.

Ofman, L.↗

A Gaussian Process Enhancement to Linear Parameter Varying Models

Simulation and analysis for modern engineering systems now routinely requires the merging of multiple disciplines, physical-domains, time-scales, and data sets — all at ever increasing levels. These capabilities are especially needed in the domain of Advanced Air Mobility, where rapidly emerging vehicle designs are significantly more complex, while having to be both cost-effective and safe. To meet these engineering challenges, machine learning methods are an attractive option for merging models and data across multiple areas while providing uncertainty quantification and maintaining computational efficiency. This paper examines the use of Gaussian process machine learning to generalize and enhance the commonly used class of quasi-Linear Parameter Varying models for fast full-envelope simulation while also supporting control system design and analysis with model uncertainty. Gaussian process machine learning is selected because it: can fuse multiple data sets, enables an easy trade-off between data fitting and smoothing, provides model uncertainty quantification, scales well with increasing complexity, and does not generally require starting from a large training data set. To demonstrate the benefits of the approach, a robust stability analysis with Gaussian process uncertainty is shown for a NASA reference design of an electric quad-rotor air-taxi concept vehicle with motor parameter uncertainty.

Gaussian Process↗

Coronal Heating and the Magnetic Flux Content of the Network

Previously, from analysis of SOHO coronal images in combination with Kitt Peak magnetograms, we found that the quiet corona is the sum of two components: the large-scale corona and the coronal network. The large-scale corona consists of all coronal-temperature (T approximately 10(exp 6) K) structures larger than supergranules (greater than approximately 30,000 kilometers). The coronal network (1) consists of all coronal-temperature structures smaller than supergranules, (2) is rooted in and loosely traces the photospheric magnetic network, (3) has its brightest features seated on polarity dividing lines (neutral lines) in the network magnetic flux, and (4) produces only about 5% of the total coronal emission in quiet regions. The heating of the coronal network is apparently magnetic in origin. Here, from analysis of EIT coronal images of quiet regions in combination with magnetograms of the same quiet regions from SOHO/MDI and from Kitt Peak, we examine the other 95% of the quiet corona and its relation to the underlying magnetic network. We find: (1) Dividing the large-scale corona into its bright and dim halves divides the area into bright "continents" and dark "oceans" having spans of 2-4 supergranules. (2) These patterns are also present in the photospheric magnetograms: the network is stronger under the bright half and weaker under the dim half. (3) The radiation from the large-scale corona increases roughly as the cube root of the magnetic flux content of the underlying magnetic network. In contrast, the coronal radiation from an active region increases roughly linearly with the magnetic flux content of the active region. We assume, as is widely held, that nearly all of the large-scale corona is magnetically rooted in the network. Our results suggest that either the coronal heating in quiet regions has a large non-magnetic component, or, if the heating is predominantly produced via the magnetic field, the mechanism is significantly different than in active regions.

Falconer, D. A.↗

Quasi-linear regime of gravitational instability: Implication to density-velocity relation

The well known linear relation between density and peculiar velocity distributions is a powerful tool for studying the large-scale structure in the Universe. Potentially it can test the gravitational instability theory and measure Omega. At present it is used in both ways: the velocity is reconstructed, provided the density is given, and vice versa. Reconstructing the density from the velocity field usually makes use of the Zel'dovich approximation. However, the standard linear approximation in Eulerian space is used when the velocity is reconstructed from the density distribution. I show that the linearized Zel'dovich approximation, in other words the linear approximation in the Lagrangian space, is more accurate for reconstructing velocity. In principle, a simple iteration technique can recover both the density and velocity distributions in Lagrangian space, but its practical application may need an additional study.

Shandarin, Sergei F.↗

Spectral discrimination of lithologic facies in the granite of the Pedra Branca Goias using LANDSAT 1 digital imagery

The applicability of LANDSAT MSS imagery for discriminating geobotanical associations observed in zones of cassiterite-rich metasomatic alterations in the granitic body of Serra da Pedra Branca was investigated. Computer compatible tapes of dry and rainy season imagery were analyzed. Image enlargement, corrections, linear contrast stretch, and ratioing of noncorrelated spectral bands were performed using the Image 100 with a grey scale of 256 levels between zero and 255. Only bands 5 and 7 were considered. Band ratioing of noncorrelated channels (5 and 7) of rainy season imagery permits distinction of areas with different vegetation coverage percentage, which corresponds to geobotanial associations in the area studied. The linear contrast stretch of channel 5, especially of the dry season image is very unsatisfactory in this area.

Parada, N. D. J.↗

A simplified integrated framework for predicting the economic impacts of feedstock variations in a catalytic fast pyrolysis conversion process

Feedstock attributes of lignocellulosic biomass, such as particle size, compositional makeup, and moisture content, can vary substantially even within pre-processed materials and have a significant effect on conversion in fast pyrolysis-based processes. However, the economic impacts of these attributes are not well understood. To address this, biomass deconstruction phenomena captured with a versatile particle-scale simulation were linked to techno-economic impacts via reduced-order models. Parametric analysis of the particle-scale model, which was validated using literature data, was used in combination with multiple linear regression models to develop correlations between feedstock attributes and yields of pyrolysis oil, gas, and char. Yields were then correlated with the minimum fuel selling price (MFSP) using a techno-economic model, bridging the gap between physics-based biomass conversion simulations and predictions of MFSP for a catalytic fast-pyrolysis process. Empirical correlations derived from the literature regarding the impact of mineral matter (ash) on oil yield were also considered. The model correlations deployed in the integrated framework capture the impacts of variation in feedstock attributes on the MFSP. Variations in ash were shown to have the biggest impact, varying MFSP by -13%/+22% due to catalytic effects and lower relative amounts of convertible lignocellulosic material. It was also found that, if ash can be controlled to low levels, the increased extractives in forest residues can help compensate for some yield losses associated with increased ash. As a result, other inputs considered (particle size, moisture content, and reactor temperature) had relatively negligible effects on process economics within the ranges analyzed considering particle-scale effects alone.

BIOMASS FUELS↗

DNA Reaction–Diffusion Attractor Patterns

Abstract Living systems can form and recover complex chemical patterns with precisely sized features in the ranges of tens or hundreds of microns. We show how designed reaction–diffusion processes can likewise produce precise patterns, termed attractor patterns, that reform their precise shape after being perturbed. We use oligonucleotide reaction networks, photolithography, and microfluidic delivery to form precisely controlled attractor patterns and study the responses of these patterns to different localized perturbations. Linear and “hill”‐shaped patterns formed and stabilized into shapes and at time scales consistent with reaction–diffusion models. When patterns were perturbed in particular locations with UV light, they reliably reformed their steady‐state profiles. Recovery also occurred after repeated perturbations. By designing the far‐from‐equilibrium dynamics of a chemical system, this study shows how it is possible to design spatial patterns of molecules that are sustained and regenerated by continually evolving towards a specific steady state configuration.

Dorsey, Phillip James↗

DNA Reaction–Diffusion Attractor Patterns

Living systems can form and recover complex chemical patterns with precisely sized features in the ranges of tens or hundreds of microns. We show how designed reaction-diffusion processes can likewise produce precise patterns, termed attractor patterns, that reform their precise shape after being perturbed. We use oligonucleotide reaction networks, photolithography and microfluidic delivery to form precisely controlled attractor patterns and study the responses of these patterns to different localized perturbations. Linear and ‘hill’-shaped patterns formed and stabilized into shapes and at time scales consistent with reaction-diffusion models. When patterns were perturbed in particular locations with UV light, they reliably reformed their steady state profiles. Recovery also occurred after repeated perturbations. As a result, by designing the far-from-equilibrium dynamics of a chemical system, this study shows how it is possible to design spatial patterns of molecules that are sustained and regenerated by continually evolving towards a specific steady state configuration

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Droplet breakup by multimodal nonlinear Rayleigh Taylor instability

A droplet impacted by a shock wave will undergo a process of fragmentation due to the development of interfacial hydrodynamic instabilities. The interface experiences variable acceleration and shear that result in the development of both inertial (Rayleigh–Taylor) and shear (Kelvin–Helmholtz) instabilities. These perturbations grow in time and drive the fragmentation and breakup of the deformed droplet. Experiments are performed on nominally 0.86 mm water droplet subjected to a Mach 7.6 detonation wave, resulting in a high Weber number (~36,000) breakup event. Perturbation growth is measured from a series of high-speed (> MHz) shadowgraph images. It is proposed that, given the size of the large-scale perturbations observed in experiments, these instabilities are growing in the non-linear regime and can be described by bubble-merger models for nonlinear mixing. Calculations are performed for the growth rates and size of these instabilities using deformation and external flow models to establish the time-dependent boundary conditions. The concurrence of the measured perturbation widths and the predictions of the simple model lend credence to the theory. This novel approach serves to open a new avenue in the characterization of droplet breakup via hydrodynamic instabilities.

Droplet breakup↗

A scalable matrix-free spectral element approach for unsteady PDE constrained optimization using PETSc/TAO

In this work, we provide a new approach for the efficient matrix-free application of the transpose of the Jacobian for the spectral element method for the adjoint-based solution of partial differential equation (PDE) constrained optimization. This results in optimizations of nonlinear PDEs using explicit integrators where the integration of the adjoint problem is not more expensive than the forward simulation. Solving PDE constrained optimization problems entails combining expertise from multiple areas, including simulation, computation of derivatives, and optimization. The Portable, Extensible Toolkit for Scientific computation (PETSc) together with its companion package, the Toolkit for Advanced Optimization (TAO), is an integrated numerical software library that contains an algorithmic/software stack for solving linear systems, nonlinear systems, ordinary differential equations, differential algebraic equations, and large-scale optimization problems and, as such, is an ideal tool for performing PDE-constrained optimization. This paper describes an efficient approach in which the software stack provided by PETSc/TAO can be used for large-scale nonlinear time-dependent problems. Time integration can involve a range of high-order methods, both implicit and explicit. The PDE-constrained optimization algorithm used is gradient-based and seamlessly integrated with the simulation of the physical problem.

97 MATHEMATICS AND COMPUTING↗

Patterned anodes with sub-millimeter spatial resolution for large-area MCP-based photodetector systems

Micro-channel-plate-based photo-detectors are unique in being capable of covering very large areas such as those required in elementary particle and nuclear physics, while providing sub-millimeter space resolution, time resolutions of less than 10 picoseconds for charged particles, and time resolutions of 30 psec–50 psec for single photons. In such systems the electronic channel count is a major cost driver. Incorporating a capacitively-coupled anode allows the use of external pickup electrodes with patterns of individual channels optimized for occupancy, rate, and time/space resolution. The signal pickup antenna can be economically implemented as a printed circuit card with a 2-dimensional array of pads for high-occupancy/high-rate applications such as in particle colliders and medical imaging, or a 1-dimensional array of strips for a lower channel count in low-occupancy/low-rate applications such as large neutrino detectors. In this paper we present pad patterns that enhance signal-sharing between pads to lower the channel count per unit area in large-area systems by factors up to 4, while maintaining spatial resolutions of approximately 100 to 200 μ m for charged particles and 400 to 1000 μ m for single photons. Patterns that use multiple signal layers in the signal-pickup board can lower the channel count even further, moving the scaling behavior in the number of pads versus total area from quadratic to linear.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gyrokinetic simulation of turbulence and transport in the SPARC tokamak

The turbulence and transport expected in the SPARC tokamak Primary Reference Discharge (PRD) [P. Rodriguez-Fernandez et al., J. Plasma Phys. 86, 865860503 (2020)] have been investigated with the gyrokinetic code CGYRO [J. Candy et al., J. Comput. Phys. 324, 73–93 (2016)]. Linear and nonlinear simulations that focus on ion (kθρs<1.0) and electron-scale (kθρs>1.0) turbulence were used to probe the nature of the turbulence and the resulting transport in the fusion core. It is found that in the SPARC PRD, ion temperature gradient (ITG) turbulence is expected to dominate transport over most of the profile with some potential trapped electron mode impact in the near edge. Stiff turbulence is observed over a part of the plasma core such that SPARC's ion temperature profile will likely be pinned to just above the critical gradient for ITG. The role of electromagnetic turbulence, rotation, and electron-scale turbulence was investigated to provide some insight into the physics required to accurately predict SPARC performance via gyrokinetics. Additionally, predictions of impurity peaking for potential low- and high-Z SPARC first-wall materials are probed using ion-scale simulation. The dominance of low-k turbulence in SPARC provides a potential opportunity for more tractable prediction of plasma profiles using nonlinear gyrokinetics. This work is the first step toward full gyrokinetic profile prediction of SPARC kinetic profiles and the resulting fusion power and plasma gain.

Howard, N. T. (ORCID:0000000287876309)↗