Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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 919 records · Page 51

Competition between Exceptionally Long‐Range Alkyl Sidechain Ordering and Backbone Ordering in Semiconducting Polymers and Its Impact on Electronic and Optoelectronic Properties

Abstract Intra‐ and intermolecular ordering greatly impacts the electronic and optoelectronic properties of semiconducting polymers. The interrelationship between ordering of alkyl sidechains and conjugated backbones has yet to be fully detailed, despite much prior effort. Here, the discovery of a highly ordered alkyl sidechain phase in six representative semiconducting polymers, determined from distinct spectroscopic and diffraction signatures, is reported. The sidechain ordering exhibits unusually large coherence lengths (≥70 nm), induces torsional/twisting backbone disorder, and results in a vertically multilayered nanostructure with ordered sidechain layers alternating with disordered backbone layers. Calorimetry and in situ variable temperature scattering measurements in a model system poly{4‐(5‐(4,8‐bis(3‐butylnonyl)‐6‐methylbenzo[1,2‐b:4,5‐b′]dithiophen‐2‐yl)thiophen‐2‐yl)‐2‐(2‐butyloctyl)‐5,6‐difluoro‐7‐(5‐methylthiophen‐2‐yl)‐2H‐benzo[d][1,2,3]triazole} (PBnDT‐FTAZ) clearly delineate this competition of ordering that prevents simultaneous long‐range order of both moieties. The long‐range sidechain ordering can be exploited as a transient state to fabricate PBnDT‐FTAZ films with an atypical edge‐on texture and 2.5× improved field‐effect transistor mobility. The observed influence of ordering between the moieties implies that improved molecular design can produce synergistic rather than destructive ordering effects. Given the large sidechain coherence lengths observed, such synergistic ordering should greatly improve the coherence length of backbone ordering and thereby improve electronic and optoelectronic properties such as charge transport and exciton diffusion lengths.

Carpenter, Joshua H.↗

A Two Colorable Fourth Order Compact Difference Scheme and Parallel Iterative Solution of the 3D Convection Diffusion Equation

A new fourth order compact difference scheme for the three dimensional convection diffusion equation with variable coefficients is presented. The novelty of this new difference scheme is that it Only requires 15 grid points and that it can be decoupled with two colors. The entire computational grid can be updated in two parallel subsweeps with the Gauss-Seidel type iterative method. This is compared with the known 19 point fourth order compact differenCe scheme which requires four colors to decouple the computational grid. Numerical results, with multigrid methods implemented on a shared memory parallel computer, are presented to compare the 15 point and the 19 point fourth order compact schemes.

Zhang, Jun↗

A reconstructed discontinuous Galerkin method for compressible flows in Lagrangian formulation

In this work, a high-order accurate reconstructed discontinuous Galerkin (rDG) method is developed for solving two-dimensional hydrodynamic problems in cell-centered updated Lagrangian formulation. This method is the Lagrangian limit of the unsplit rDG-ALE formulation, and is obtained by assuming the equality of the grid velocity to the fluid velocity only at cell boundaries. The conservative variables and the Taylor basis defined on the time-dependent moving mesh, provide the piece-wise polynomial expansion in the updated Lagrangian formulation. A multi-directional nodal Riemann solver is implemented for computing the grid velocity at the vertices and the numerical flux at the cell boundaries. A characteristic limiting procedure is extended from the primitive variable version to the conservative variable version, and its performance is compared with the limiter on physical variables. A number of benchmark test cases are conducted to assess the accuracy, robustness, and non-oscillatory property of the DG(P0), DG(P1) and rDG(P1P2) methods. The numerical experiments demonstrate that the developed rDG method is able to attain the designed order of accuracy and the characteristic limiting procedure outperforms the limiter on physical variables in terms of the monotonicity and symmetry preservation for shock problems.

42 ENGINEERING↗

High-Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

It has been known for some time that Taylor series (TS) integration is among the most efficient and accurate numerical methods in solving differential equations. However, the full benefit of the method has yet to be realized in calculating spacecraft trajectories, for two main reasons. First, most applications of Taylor series to trajectory propagation have focused on relatively simple problems of orbital motion or on specific problems and have not provided general applicability. Second, applications that have been more general have required use of a preprocessor, which inevitably imposes constraints on computational efficiency. The latter approach includes the work of Berryman et al., who solved the planetary n-body problem with relativistic effects. Their work specifically noted the computational inefficiencies arising from use of a preprocessor and pointed out the potential benefit of manually coding derivative routines. In this Engineering Note, we report on a systematic effort to directly implement Taylor series integration in an operational trajectory propagation code: the Spacecraft N-Body Analysis Program (SNAP). The present Taylor series implementation is unique in that it applies to spacecraft virtually anywhere in the solar system and can be used interchangeably with another integration method. SNAP is a high-fidelity trajectory propagator that includes force models for central body gravitation with N X N harmonics, other body gravitation with N X N harmonics, solar radiation pressure, atmospheric drag (for Earth orbits), and spacecraft thrusting (including shadowing). The governing equations are solved using an eighth-order Runge-Kutta Fehlberg (RKF) single-step method with variable step size control. In the present effort, TS is implemented by way of highly integrated subroutines that can be used interchangeably with RKF. This makes it possible to turn TS on or off during various phases of a mission. Current TS force models include central body gravitation with the J2 spherical harmonic, other body gravitation, thrust, constant atmospheric drag from Earth's atmosphere, and solar radiation pressure for a sphere under constant illumination. The purpose of this Engineering Note is to demonstrate the performance of TS integration in an operational trajectory analysis code and to compare it with a standard method, eighth-order RKF. Results show that TS is 16.6 times faster on average and is more accurate in 87.5% of the cases presented.

Scott, James R.↗

Interacting nuclei in distant galaxies

The N-galaxy 3C 390.3 has been monitored spectroscopically since 1974 (Osterbrock, Koski and Phillips 1975; Oke 1988). From various archives and literature, it is found that the Balmer lines change their intensities and profiles in a dramatic manner. The H alpha profile is very broad and peculiar, and the relative intensities of its two humps changes consistently with time, possibly periodically. Before 1980, the blue hump was significantly stronger than the one in the red. From 1980 to 1983 the blue hump became stronger (see Oke 1988). After 1983 the H alpha profile has returned to its early shape and seems to have completed a full circle. Unlike the rapid (on the order of a month or even less) and aperiodic variation in the continuum and integrated line intensities, the change in broad profile seems slow and consistent. Taking the analogy of cataclysmic variables, the double-horn profiles have been observed in cases of interacting stars. For example, the emission lines, both in He II and hydrogen Balmer lines in GD 552 (Stover 1985) show double-horn profiles and periodical changes in their line profiles, including the change in ratios of two humps. It is understood that the D-wave components (Smak 1976) are the signature of an emitting disk and the S-wave component is from the emission at a hot spot which rotates and results in a moving component in the velocity space. The mass flow from the nearby interacting star provides the stream toward the core of a neutron star or white dwarf. Therefore, it is proposed that the variation of broad line profiles observed in 3C 390.3 may be the result of a pair of interacting massive cores. The rotational velocity dominates and produces a variable double-horn profile. However, the line widths observed in broad line radio galaxies are one order larger than that in interacting stars. The Balmer decrements imply a much smaller density (10(exp 10-12) cm(exp-3)) than that in the cataclysmic variables. The much larger velocity and much thinner density make it unlikely that the broad line emission is simply formed in an accretion disk. The authors postulate that a significant rotational motion is involved. If the observed squared profiles are indeed due to the rotational velocity field, one can naturally explain their rare occurrence.

Zheng, Wei↗

Heptagon functions and seven-gluon amplitudes in multi-Regge kinematics

We compute all 2 → 5 gluon scattering amplitudes in planar $ \mathcal{N} $ N = 4 super-Yang-Mills theory in the multi-Regge limit that is sensitive to the non-trivial (“long”) Regge cut. We provide the amplitudes through four loops and to all logarithmic accuracy at leading power, in terms of single-valued multiple polylogarithms of two variables. To obtain these results, we leverage the function-level results for the amplitudes in the Steinmann cluster bootstrap. To high powers in the series expansion in the two variables, our results agree with the recently conjectured all-order central emission vertex used in the Fourier-Mellin representation of amplitudes in multi-Regge kinematics. Our results therefore provide a resummation of the Fourier-Mellin residues into single-valued polylogarithms, and constitute an important cross-check between the bootstrap approach and the all-orders multi-Regge proposal.

1/N expansion↗

Differentiating Structural Changes of Glycoproteins in Solution using Small Angle Scattering Analysis

Alpha 1-Antitrypsin (A1AT) and Fibrinogen are plasma glycoproteins with different, but specific biological functions. A1AT (which forms part of the serpin family) has been shown to have protective roles of lung cells against emphysema, a disease characterized by lung tissue destruction, while fibrinogen is a major factor in the blood clotting process. While biologically distinct, the two proteins present similarities in biophysical studies due to their size and polymerization tendency. Previous approximations, using radius of gyration analysis (shape independent), have shown that there is a phase transition for A1AT under osmotic stress. The goals of this study are to better understand: the stability of A1AT's native conformation when its environment deviates from physiological conditions; and the effects on fibrin's final structure as fibrinogen undergoes polymerization under different conditions. Environmental variables include osmotic stress, crowding, and presence of other molecules (e.g. lipids, catalytic enzyme thrombin, or cholesterol) found in plasma. The method used to assess the stability of the native conformation was to observe changes in the structural envelope of the protein. Static and dynamic small angle x-ray scattering (SAXS) and small angle neutron scattering (SANS) experiments of the proteins in solution were conducted for this purpose. Scattering data was fitted using various shape dependent models in SASView. Pair distance functions (PDF) were used to validate any changes identified by the model fits. Multi-variable phase diagrams were constructed from the fitting results in order to quantify structural and conformational changes of the proteins.

Segally, Taylor N.↗

Scalable Preconditioning of Block-Structured Linear Algebra Systems using ADMM

In this paper, we study the solution of block-structured linear algebra systems arising in optimization by using iterative solution techniques. These systems are the core computational bottleneck of many problems of interest such as parameter estimation, optimal control, network optimization, and stochastic programming. Our approach uses a Krylov solver (GMRES) that is preconditioned with an alternating method of multipliers (ADMM). We show that this ADMM-GMRES approach overcomes well-known scalability issues of Schur complement decomposition in problems that exhibit a high degree of coupling. The effectiveness of the approach is demonstrated using linear systems that arise in stochastic optimal power flow problems and that contain up to 2 million total variables and 4,000 coupling variables. We find that ADMM-GMRES is nearly an order of magnitude faster than Schur complement decomposition. Moreover, we demonstrate that the approach is robust to the selection of the augmented Lagrangian penalty parameter, which is a key advantage over the direct use of ADMM.

97 MATHEMATICS AND COMPUTING↗

The small-angle neutron scattering extension in MCNPX and the SANS cross section for nanodiamonds

Detonation nanodiamonds provide an intriguing prospect for the development of a very cold neutron (VCN) source or as a high albedo material for VCN reflection helping to control leakage of cold and very cold neutrons in a moderator. The investigation of source and reflector designs necessitates the development of new simulation tools for existing Monte Carlo codes that include small angle neutron scattering (SANS). The SANS extension in MCNPX 2.7 provides an analytical model of hard-sphere scattering with variable scattering particle sizes, size distributions, and packing fractions in order to supplement currently existing scattering kernels as well as the capability to read an experimentally determined small angle neutron scattering profile from a data file. The SANS extension is discussed and compared to nanodiamond SANS measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Traleika Glacier: A hardware-software co-designed approach to exascale computing

The move from current petascale machines to future exascale machines will need both hardware improvements and software changes. Hardware will need to evolve to focus primarily on features that lower energy consumption: near-threshold voltage operation, fine-grained power and clock management and heterogeneity. Software will also need to evolve and be able to express more parallelism, become more dynamic and adaptable in order to be able to operate on a much more variable hardware. Here, we present Traleika Glacier, an effort that seeks to evaluate radical design changes to meet the constraints, both in terms of power and cost, of exascale computing. The salient features of the hardware design presented in the work include a) a use of heterogeneous cores, b) a redesign of the memory system that centers around hierarchical scratchpads and a global address space, c) the hardware acceleration of certain memory and network operations through specialized engines and, d) very fine-grained control and monitoring capabilities. On the software side, we describe a task-based runtime system, the Open Community Runtime (OCR) which aims to express a wide range of higher-level programming models with a very limited set of core concepts: event-driven tasks for computation, events for synchronization and relocatable data-blocks for data management.

97 MATHEMATICS AND COMPUTING↗

Use of the rVV10 Nonlocal Correlation Functional in the B97M-V Density Functional: Defining B97M-rV and Related Functionals [On the Use of the rVV10 Nonlocal Correlation Functional in the B97M-V Density Functional: Defining B97M-rV and Related Functionals]

The VV10 and rVV10 nonlocal correlation functionals are consistently implemented and assessed, with the goal of determining if the rVV10 nonlocal correlation functional can replace the VV10 nonlocal correlation functional in the recently developed B97M-V density functional, to give the B97M-rV density functional. Along the way, four density functionals are simultaneously tested: VV10, rVV10, B97M-V, and B97M-rV. An initial assessment is carried out across the S22 data set, and the short-range damping variable, b, is varied for all four density functionals in order to determine the sensitivity of the functionals to the empirical parameter. The results of this test indicate that a value of b = 6 (fortuitously the same as that in B97M-V) is suitable for B97M-rV. The functionals are then compared across an extensive database of interaction energies, and it is demonstrated that B97M-rV either matches or outperforms B97M-V for all of the tests considered. Finally, the optimization of b across the S22 data set is extended to two range-separated hybrid density functionals, ωB97X-V and ωB97M-V, and a value of b = 6.2 is recommended for both ωB97X-rV and ωB97M-rV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of Grafting Density on Block Polymer Self-Assembly: From Linear to Bottlebrush

Grafting density is an important structural parameter that imparts significant influences over the physical properties of architecturally complex polymers. In this paper, the physical consequences of varying the grafting density (z) were studied in the context of block polymer self-assembly. Well-defined block polymers spanning the linear, comb, and bottlebrush regimes (0 ≤ z ≤ 1) were prepared via grafting-through ring-opening-metathesis polymerization (ROMP). ω-norbornenyl poly(D,L-lactide) (PLA) and polystyrene (PS) macromonomers were copolymerized with discrete co-monomers in different feed ratios, enabling precise control over the grafting density. Small-angle X-ray scattering (SAXS) experiments demonstrate that these graft block polymers can self-assemble into long-range-ordered lamellar structures. For seventeen series of block polymers with variable z, the scaling of the lamellar period with the total backbone degree of polymerization (d* ~ N bb α ) was studied. The scaling exponent α monotonically decreases with decreasing z and exhibits an apparent transition at z ≈ 0.2, suggesting significant changes in the chain conformations. Comparison of two block polymer systems, one that is strongly segregated for all z (System I) and one that experiences weak segregation at low z (System II), indicates that the observed trends are primarily caused by the polymer architectures, instead of segregation strengths. A model is pro-posed in which the characteristic ratio (C ∞ ), a proxy for the backbone stiffness, scales with N bb as a function of the grafting density: C ∞ ~ N bb f(z) . To the best of our knowledge, this report represents the first study of scaling behavior for the self-assembly of block polymers with variable grafting density. Lastly, the relationships disclosed herein provide valuable insights into conformational changes with grafting density, thus introducing new opportunities for future block polymer design.

36 MATERIALS SCIENCE↗

Influence of temporal variations in plasma conditions on the electric potential near self-organized dust chains

Self-organization of dust grains into stable filamentary dust structures (or “chains”) largely depends on dynamic interactions between individual charged dust grains and complex electric potential arising from the distribution of charges within a local plasma environment. Recent studies have shown that the positive column of the gas discharge plasma in the Plasmakristall-4 (PK-4) experiment at the International Space Station supports the presence of fast-moving ionization waves, which lead to variations of plasma parameters by up to an order of magnitude from the average background values. The highly variable environment resulting from ionization waves may have interesting implications for the dynamics and self-organization of dust particles, particularly concerning the formation and stability of dust chains. Here, we investigate the electric potential surrounding dust chains in the PK-4 experiment by employing a molecular dynamics model of the dust and ions with boundary conditions supplied by a particle-in-cell with Monte Carlo collision simulation of the ionization waves. The model is used to examine the effects of the plasma conditions within different regions of the ionization wave and compare the resulting dust structure to that obtained by employing the time-averaged plasma conditions. Here, the comparison between simulated dust chains and experimental data from the PK-4 experiment shows that the time-averaged plasma conditions do not accurately reproduce observed results for dust behavior, indicating that more careful treatment of plasma conditions in the presence of ionization waves is required. It is further shown that commonly used analytic forms of the electric potential do not accurately describe the electric potential near charged dust grains under these plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Grad–Shafranov equilibria via data-free physics informed neural networks

A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad–Shafranov equation by using physics-informed neural networks (PINNs). For PINNs, we optimize neural networks by directly minimizing the residual of the partial differential equation as a loss function. We show that PINNs can accurately and effectively solve the Grad–Shafranov equation with several different boundary conditions, making it more flexible than traditional solvers. This method is flexible as it does not require any mesh and basis choice, thereby streamlining the computational process. We also explore the parameter space by varying the size of the model, the learning rate, and boundary conditions to map various tradeoffs such as between reconstruction error and computational speed. Additionally, we introduce a parameterized PINN framework, expanding the input space to include variables such as pressure, aspect ratio, elongation, and triangularity in order to handle a broader range of plasma scenarios within a single network. Parameterized PINNs could be used in future work to solve inverse problems such as shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural basis for aminoacylation of cellular modified tRNALys3 by human lysyl-tRNA synthetase

Abstract The average eukaryotic transfer ribonucleic acid (tRNA) contains 13 post-transcriptional modifications; however, their functional impact is largely unknown. Our understanding of the complex tRNA aminoacylation machinery in metazoans also remains limited. Herein, using a series of high-resolution cryo-electron microscopy (cryo-EM) structures, we provide the mechanistic basis for recognition and aminoacylation of fully modified cellular tRNALys3 by human lysyl-tRNA synthetase (h-LysRS). The tRNALys3 anticodon loop modifications S34 (mcm5s2U) and R37 (ms2t6A) play an integral role in recognition by h-LysRS. Modifications in the T-, variable-, and D-loops of tRNALys3 are critical for ordering the metazoan-specific N-terminal domain of LysRS. The two catalytic steps of tRNALys3 aminoacylation are structurally ordered; docking of the 3′-CCA end in the active site cannot proceed until the lysyl–adenylate intermediate is formed and the pyrophosphate byproduct is released. Association of the h-LysRS–tRNALys3 complex with a multi-tRNA synthetase complex-derived peptide shifts the equilibrium toward the 3′-CCA end “docked” conformation and allosterically increases h-LysRS catalytic efficiency. The insights presented here have broad implications for understanding the role of tRNA modifications in protein synthesis, the human aminoacylation machinery, and the growing catalog of metabolic and neurological diseases linked to it.

Devarkar, Swapnil C. (ORCID:000000029271243X)↗

Time-dependent orbital-free density functional theory: Background and Pauli kernel approximations

Time-dependent orbital-free density functional theory (DFT) is an efficient method for calculating the dynamic properties of large-scale quantum systems due to the low computational cost compared to standard time-dependent DFT. In this work, we formalize this method by mapping the real system of interacting fermions onto a fictitious system of noninteracting bosons. The dynamic Pauli potential and associated kernel emerge as key ingredients of time-dependent orbital-free DFT. Using the uniform electron gas as a model system, we derive an approximate frequency-dependent Pauli kernel. Pilot calculations suggest that space nonlocality is a key feature for this kernel. Nonlocal terms arise already in the second-order expansion with respect to unitless frequency and reciprocal space variable ($\frac{ω}{qk_F}$ and $\frac{q}{2k_F}$, respectively). Given the encouraging performance of the proposed kernel, we expect it will lead to more accurate orbital-free DFT simulations of nanoscale systems out of equilibrium. Additionally, the proposed path to formulate nonadiabatic Pauli kernels presents several avenues for further improvements which can be exploited in future work to improve the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Collaborative Research: Cloudiness transitions within shallow marine clouds near the Azores

Marine stratocumulus clouds are low, persistent, liquid phase clouds that cover large areas and play a significant role in moderating the climate by reflecting large quantities of incoming solar radiation. The deficiencies in simulating these clouds in global climate models are widely recognized. Much of the uncertainty arises from sub-grid scale variability in the cloud albedo that is not accurately parameterized in climate models. The Clouds, Aerosol and Precipitation in the Marine Boundary Layer (CAP–MBL) observational campaign and the ongoing ARM site measurements on Graciosa Island in the Azores aim to sample the Northeast Atlantic low cloud regime. These data represent, the longest continuous research quality cloud radar/lidar/radiometer/aerosol data set of open-ocean shallow marine clouds in existence. Data coverage from CAP–MBL and the series of cruises to the southeast Pacific culminating in VOCALS will both be of sufficient length to contrast the two low cloud regimes and explore the joint variability of clouds in response to several environmental factors implicated in cloudiness transitions. Our research seeks to better understand cloud system processes in an underexplored but climatologically important maritime region. Our primary goal is an improved physical understanding of low marine clouds on temporal scales of hours to days. It is well understood that aerosols, synoptic-scale forcing, surface fluxes, mesoscale dynamics, and cloud microphysics all play a role in cloudiness transitions. However, the relative importance of each mechanism as a function of different environmental conditions is unknown. To better understand cloud forcing and response, we are documenting the joint variability of observed environmental factors and associated cloud characteristics. In order to narrow the realm of likely parameter ranges, we assess the relative importance of parameter conditions based primarily on two criteria: how often the condition occurs (frequency) and to what degree varying that condition within its typically observed range affects cloud characteristics (magnitude of impact given the condition). In this manner we will be able to address the relative importance of individual factors within a multivariate range of environmental conditions. We will determine the relative roles of the thermodynamic, aerosol, and synoptic environmental factors on low cloud and drizzle formation and lifetime.

54 ENVIRONMENTAL SCIENCES↗

Thermal field reconstruction and compressive sensing using proper orthogonal decomposition

Model order reduction allows critical information about sensor placement and experiment design to be distilled from raw fluid mechanics simulation data. In many cases, sensed information in conjunction with reduced order models can also be used to regenerate full field variables. In this paper, a proper orthogonal decomposition (POD) inferencing method is extended to the modeling and compressive sensing of temperature, a scalar field variable. The method is applied to a simulated, critically stable, incompressible flow over a heated cylinder (Re = 1000) with Prandtl number varying between 0.001 and 50. The model is trained on pressure and temperature data from simulations. Field reconstructions are then generated using data from selected sensors and the POD model. Finally, the reconstruction error is evaluated across all Prandtl numbers for different numbers of retained modes and sensors. The predicted trend of increasing reconstruction accuracy with decreasing Prandtl number is confirmed and a Prandtl number/sensor count error matrix is presented.

Matulis, John↗