Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed computation”

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,189 records · Page 66

Temporal breakdown and Borel resummation in the complex Langevin method

We reexamine the Parisi-Klauder conjecture for complex e{sup i{theta}/2}{phi}{sup 4} measures with a Wick rotation angle 0{<=}{theta}/2{<=}{pi}/2 interpolating between Euclidean signature and Lorentzian signature. Our main result is that the asymptotics for short stochastic times t encapsulates information also about the equilibrium aspects. The moments evaluated with the complex measure and with the real measure defined by the stochastic Langevin equation have the same t{yields}0 asymptotic expansion which is shown to be Borel summable. The Borel transform correctly reproduces the time dependent moments of the complex measure for all t, including their t{yields}{infinity} equilibrium values. On the other hand the results of a direct numerical simulation of the Langevin moments are found to disagree from the 'correct' result for t larger than a finite t{sub c}. The breakdown time t{sub c} increases powerlike for decreasing strength of the noise's imaginary part but cannot be excluded to be finite for purely real noise. To ascertain the discrepancy we also compute the real equilibrium distribution for complex noise explicitly and verify that its moments differ from those obtained with the complex measure. - Highlights: Black-Right-Pointing-Pointer The Parisi-Klauder conjecture is reexamined for complex e{sup i{theta}/2}{phi}{sup 4} measures. Black-Right-Pointing-Pointer The time dependent moments are evaluated by temporal Borel resummation. Black-Right-Pointing-Pointer The results disagree with the Langevin simulations beyond a critical time t{sub c}. Black-Right-Pointing-Pointer t{sub c} increases with decreasing strength of the noise's imaginary part. Black-Right-Pointing-Pointer The technical reason for the breakdown is identified.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Estimating the reproductive number, total outbreak size, and reporting rates for Zika epidemics in South and Central America

As South and Central American countries prepare for increased birth defects from Zika virus outbreaks and plan for mitigation strategies to minimize ongoing and future outbreaks, understanding important characteristics of Zika outbreaks and how they vary across regions is a challenging and important problem. We developed a mathematical model for the 2015/2016 Zika virus outbreak dynamics in Colombia, El Salvador, and Suriname. We fit the model to publicly available data provided by the Pan American Health Organization, using Approximate Bayesian Computation to estimate parameter distributions and provide uncertainty quantification. The model indicated that a country-level analysis was not appropriate for Colombia. We then estimated the basic reproduction number to range between 4 and 6 for El Salvador and Suriname with a median of 4.3 and 5.3, respectively. We estimated the reporting rate to be around 16% in El Salvador and 18% in Suriname with estimated total outbreak sizes of 73,395 and 21,647 people, respectively. The uncertainty in parameter estimates highlights a need for research and data collection that will better constrain parameter ranges.

60 APPLIED LIFE SCIENCES↗

Novel artificial neural network model for instantaneous power losses and operational efficiency mapping of MW-scale vanadium redox flow battery for improved technoeconomic analysis

A novel data-driven, machine-learning-based method for modeling the instantaneous power losses of a distribution-sited 2 MW/8MWh vanadium redox flow battery (VRFB), a grid-scale electrochemical storage technology, is introduced and compared against benchmark empirical modeling approaches, including symmetric and asymmetric models, as well as a recent convex hull modeling approach. The novel loss modeling method introduces several advantages over the benchmark models and over simplistic efficiency estimates, the most significant of which is that the model can accurately reflect the stepwise and non-linear parasitic losses associated with the duty cycles of mechanical auxiliary systems like pump motor drives and blower fans. Residuals of the models are compared; the proposed data driven model features significantly improved accuracy over the benchmark models. The model's coefficient of determination is also improved relative to that of the benchmark models. Furthermore, a novel method for visualization of operational efficiency of the grid-scale storage technology is introduced. To demonstrate the benefits of the novel data-driven method for modeling the VRFB, the benchmark models and the proposed models are embedded into an Open DSS distribution network model to study two applications of the grid-scale electrical storage system: load leveling for grid support and energy arbitrage. This article demonstrates that the accuracy of the instantaneous power loss model significantly impacts the understanding of the state of charge of the VRFB. In turn, the accuracy of the efficiency modeling of the VRFB impacts the understanding of the potential economic value and technical benefits to the distribution network operators. In conclusion, the presented power loss modeling approach is, therefore, highly relevant for utility-stakeholders, battery asset owners, system engineers, system designers, and financial planners interested in evaluating or optimizing the operation of grid-scale VRFBs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An experimental investigation of gas fuel injection with X-ray radiography

In this paper, an outward-opening compressed natural gas, direct injection fuel injector has been studied with single-shot x-ray radiography. Three dimensional simulations have also been performed to compliment the x-ray data. Argon was used as a surrogate gas for experimental and safety reasons. This technique allows the acquisition of a quantitative mapping of the ensemble-average and standard deviation of the projected density throughout the injection event. Two dimensional, ensemble average and standard deviation data are presented to investigate the quasi-steady-state behavior of the jet. Upstream of the stagnation zone, minimal shot-to-shot variation is observed. Downstream of the stagnation zone, bulk mixing is observed as the jet transitions to a subsonic turbulent jet. From the time averaged data, individual slices at all downstream locations are extracted and an Abel inversion was performed to compute the radial density distribution, which was interpolated to create three dimensional visualizations. The Abel reconstructions reveal that upstream of the stagnation zone, the gas forms an annulus with high argon density and large density gradients. Inside this annulus, a recirculation region with low argon density exists. Downstream, the jet transitions to a fully turbulent jet with Gaussian argon density distributions. This experimental data is intended to serve as a quantitative benchmark for simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Hierarchical optimization for neutron scattering problems

We present a scalable optimization method for neutron scattering problems that determines confidence regions of simulation parameters in lattice dynamics models used to fit neutron scattering data for crystalline solids. The method uses physics-based hierarchical dimension reduction in both the computational simulation domain and the parameter space. We demonstrate for silicon that after a few iterations the method converges to parameters values (interatomic force-constants) computed with density functional theory simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Formulation and characterization of a continuous crystal lattice orientation finite element method (LOFEM) and its application to dislocation fields

We report that since the 1950s, a large body of work has been published on connecting the curvature of a crystal lattice to geometrically necessary dislocations densities of a crystal lattice. Studying dislocation transmission through grains and across their boundaries requires the lattice curvature to be preserved. However, traditional crystal plasticity models and their numerical implementations do not formally preserve lattice curvature. In this paper, a continuous crystal lattice orientation finite element method (LOFEM) is proposed to rectify this impediment to the inclusion of dislocation-based constitutive models. The methodology is first presented, and then it is demonstrated for tension and compression deformations of a copper polycrystal. It is shown that under the same deformation histories, the lattice continuity constraint alters the evolving state in comparison to the traditional approach, including retarding the rate at which the crystallographic texture strengthens under monotonic deformation. Lastly, taking advantage of the finite element representation of the lattice orientation, the Nye tensor is computed on lattices misoriented by deformation and is subsequently used to compute evolving dislocation density distributions.

36 MATERIALS SCIENCE↗

Sampling two-dimensional isometric tensor network states

Sampling a quantum system’s underlying probability distributions is an important computational task, e.g., for quantum advantage experiments and quantum Monte Carlo algorithms. Tensor networks are an invaluable tool for efficiently representing states of large quantum systems with limited entanglement. Algorithms for sampling one-dimensional (1D) tensor networks are well-established and utilized in several 1D tensor network methods. In this paper we introduce two novel sampling algorithms for two-dimensional (2D) isometric tensor network states (isoTNS) that generalize existing 1D tensor network sampling algorithms. Our first proposed algorithm performs independent sampling and yields a single configuration together with its associated probability. The second algorithm employs a greedy search strategy to identify high-probability configurations and their corresponding probabilities. Numerical results demonstrate the effectiveness of these algorithms across quantum states with varying entanglement and system size.

Dumitrescu, Eugene [ORNL] (ORCID:0000000158519567)↗

Mechanistic Pathways for N 2 O Elimination from trans -R 3 Sn-O-N=N-O-SnR 3 and for Reversible Binding of CO 2 to R 3 Sn-O-SnR 3 (R = Ph, Cy)

The rate and mechanism of the elimination of N 2 O from trans-R 3 Sn-O-N=N-O-SnR 3 (R = Ph ( 1 Ph ) and R = Cy ( 1 Cy )) to form R 3 Sn-O-SnR 3 (R = Ph ( 2 Ph ) and R = Cy ( 2 Cy )) have been studied using both NMR and IR techniques to monitor the reactions in the temperature range of 39–79 °C in C 6 D 6 . Activation parameters for this reaction are ΔH ‡ = 15.8 ± 2.0 kcal·mol –1 and ΔS ‡ = –28.5 ± 5 cal·mol –1 ·K –1 for 1 Ph and ΔH ‡ = 22.7 ± 2.5 kcal·mol –1 and ΔS ‡ = –12.4 ± 6 cal·mol –1 ·K –1 for 1 Cy . Addition of O 2 , CO 2 , N 2 O, or PPh 3 to sealed tube NMR experiments did not alter in a detectable way the rate or product distribution of the reactions. Computational DFT studies of elimination of hyponitrite from trans-Me 3 Sn-O-N=N-O-SnMe 3 ( 1 Me ) yield a mechanism involving initial migration of the R 3 Sn group from O to N passing through a marginally stable intermediate product and subsequent N 2 O elimination. Reactions of 1 Ph with protic acids HX are rapid and lead to formation of R 3 SnX and trans-H 2 N 2 O 2 . Reaction of 1 Ph with the metal radical •Cr(CO) 3 C 5 Me 5 at low concentrations results in rapid evolution of N 2 O. At higher •Cr(CO) 3 C 5 Me 5 concentrations, evolution of CO 2 rather than N 2 O is observed. Addition of 1 atm or less CO 2 to benzene or toluene solutions of 2 Ph and 2 Cy resulted in very rapid reaction to form the corresponding carbonates R 3 Sn-O-C(=O)-O-SnR 3 (R = Ph ( 3 Ph ) and R = Cy ( 3 Cy )) at room temperature. Evacuation results in fast loss of bound CO 2 and regeneration of 2 Ph and 2 Cy . Variable temperature data for formation of 3 Cy yield ΔH o = –8.7 ± 0.6 kcal·mol –1 , ΔS o = –17.1 ± 2.0 cal·mol –1 ·K –1 , and ΔG o 298K = –3.6 ± 1.2 kcal·mol –1 . Furthermore, DFT studies were performed and provide additional insight into the energetics and mechanisms for the reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MARTINI-Compatible Coarse-Grained Model for the Mesoscale Simulation of Peptoids

Peptoids (poly-N-substituted glycines) are a class of synthetic polymers that are regioisomers of peptides (poly-C-substituted glycines), in which the point of side-chain connectivity is shifted from the backbone C to the N atom. Peptoids have found diverse applications as peptidomimetic drugs, protein mimetic polymers, surfactants, and catalysts. Computational modeling is valuable in the understanding and design of peptoid-based nanomaterials. In this work, we report the bottom-up parameterization of coarse-grained peptoid force fields based on the MARTINI peptide force field against all-atom peptoid simulation data. Our parameterization pipeline iteratively refits coarse-grained bonded interactions using iterative Boltzmann inversion and nonbonded interactions by matching the potential of mean force for chain extension. We assure good sampling of the amide bond cis/trans isomerizations in the all-atom simulation data using parallel bias metadynamics. We develop coarse-grained models for two representative peptoids—polysarcosine (poly(N-methyl glycine)) and poly(N-((4-bromophenyl)ethyl)glycine)—and show their structural and thermodynamic properties to be in excellent accord with all-atom calculations but up to 25-fold more efficient and compatible with MARTINI force fields. Here, this work establishes a new rigorously parameterized coarse-grained peptoid force field for the understanding and design of peptoid nanomaterials at length and time scales inaccessible to all-atom calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Capacitance of Carbon Nanotube/Graphene Composite Electrodes with [BMIM + ][BF$_4^–$]/Acetonitrile: Fixed Voltage Molecular Dynamics Simulations

The capacitance of nanoporous carbon electrode materials is dictated by the physical interactions with electrolyte ions and molecules at the accessible interior electrode surface. While significant progress has been made in designing and synthesizing carbon materials with well-defined and relatively homogeneous nanoporosity, the majority of materials remain heterogeneous, and the resulting properties reflect an average over this distribution. In this regard, computer simulations can be a valuable tool to predict or validate structure/property relationships by systematically investigating well-defined, model electrode morphologies. In this work, we utilize fixed-voltage molecular dynamics simulations to predict structure/capacitance relationships for five different model morphologies of carbon nanotube/graphene (CNT/G) composite electrodes with 1-butyl-3-methylimidazolium tetrafluoroborate/acetonitrile electrolyte. The CNT/G electrode models are inspired by experimental “layer-by-layer” syntheses and strike a balance between realism and computational tractability. Comparison between different model electrode architectures elucidates important structure/property relationships. We find that CNT/graphene contact points serve as “hot spots” with significantly enhanced charge separation relative to the rest of the electrode. Furthermore, we demonstrate a specific nanoconfinement motif that provides substantial 3–4$\times$ enhancement of local capacitance, resulting in a ~40% increase of the total electrode differential capacitance. Because the accessible surface area of the model CNT/G electrodes is precisely determined, a comparison of per-area capacitance across systems is unambiguous. Our results thus complement a prior computational demonstration of capacitance enhancement in nanoconfinement while elucidating additional interaction motifs at CNT/G electrochemical interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Light-Promoted C(sp 3 )–C(sp 3 ) Reductive Elimination from Dialkyl NiII Complexes

Ni-catalyzed cross-coupling is a powerful strategy to forge C(sp 3 )–C(sp 3 ) bonds. Typically, to do so requires overcoming a challenging C–C bond-forming reductive elimination, often enabled by the intermediacy of highly oxidized Ni species or outer-sphere processes. While direct C(sp 3 )–C(sp 3 ) reductive elimination from the Ni II base oxidation state is normally thermally inaccessible, light-activation provides an avenue to affect such transformations. Here, we investigate the mechanism of light-induced C(sp 3 )–C(sp 3 ) bond formation from dialkyl bipyridine Ni II complexes through a variety of organometallic, spectroscopic, and computational studies. Wavelength-dependent quantum yields, ligand electronics–reactivity relationships, excited-state lifetimes, computed barriers, and product distributions from crossover studies support a photolysis/radical rebound mechanism. Furthermore, this reactivity paradigm complements existing strategies in the literature to promote reductive elimination from Ni II , such as the use of destabilizing, sterically hindered ligands and reduction of electron density at Ni through the binding of electron-deficient olefins. Hence, we envision that light-induced reductive elimination may enable the development of challenging C(sp 3 )–C(sp 3 ) couplings.

Alkyls↗

Down-selection of biomolecules to assemble “reverse micelle” with perovskites

Abstract Biological molecule-semiconductor interfacing has triggered numerous opportunities in applied physics such as bio-assisted data storage and computation, brain-computer interface, and advanced distributed bio-sensing. The introduction of electronics into biological embodiment is being quickly developed as it has great potential in providing adaptivity and improving functionality. Reciprocally, introducing biomaterials into semiconductors to manifest bio-mimetic functionality is impactful in triggering new enhanced mechanisms. In this study, we utilize the vulnerable perovskite semiconductors as a platform to understand if certain types of biomolecules can regulate the lattice and endow a unique mechanism for stabilizing the metastable perovskite lattice. Three tiers of biomolecules have been systematically tested and the results reveal a fundamental mechanism for the formation of a “reverse-micelle” structure. Systematic exploration of a large set of biomolecules led to the discovery of guiding principle for down-selection of biomolecules which extends the classic emulsion theory to this hybrid systems. Results demonstrate that by introducing biomaterials into semiconductors, natural phenomena typically observed in biological systems can also be incorporated into semiconducting crystals, providing a new perspective to engineer existing synthetic materials.

59 BASIC BIOLOGICAL SCIENCES↗

New Flexible Channels for Room Temperature Tunneling Field Effect Transistors

Tunneling field effect transistors (TFETs) have been proposed to overcome the fundamental issues of Si based transistors, such as short channel effect, finite leakage current, and high contact resistance. Unfortunately, most if not all TFETs are operational only at cryogenic temperatures. Here we report that iron (Fe) quantum dots functionalized boron nitride nanotubes (QDs-BNNTs) can be used as the flexible tunneling channels of TFETs at room temperatures. The electrical insulating BNNTs are used as the one-dimensional (1D) substrates to confine the uniform formation of Fe QDs on their surface as the flexible tunneling channel. Consistent semiconductor-like transport behaviors under various bending conditions are detected by scanning tunneling spectroscopy in a transmission electron microscopy system (insitu STM-TEM). Ultimately, as suggested by computer simulation, the uniform distribution of Fe QDs enable an averaging effect on the possible electron tunneling pathways, which is responsible for the consistent transport properties that are not sensitive to bending.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Planck 2013 results. XXI. Power spectrum and high-order statistics of the Planck all-sky Compton parameter map

We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. This map shows an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales (&ell; < 60), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales (&ell; > 500) the clustered cosmic infrared background and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We thus measure the tSZ power spectrum over angular scales 0.17° ≲ θ ≲ 3.0° that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with clear evidence of additional signal from unresolved clusters and, potentially, diffuse warm baryons. Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. The measured tSZ power spectrum and high order statistics are used to place constraints on σ 8 .

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic Analysis of Levee Stability Subject to Variable Seepage Conditions

© ASCE. Levee stability is highly influenced by seepage. Specifically, hydraulic conductivity distribution within a levee influences pore pressure distribution and controls the drained strength of the soil. In this study the influence of hydraulic conductivity and blanket layer thickness on failure probability is evaluated within the context of reliability analyses that also include soil strength and unit weight as random variables. First-Order Reliability Method (FORM) is used to evaluate reliability, rank random variables by importance and to obtain sensitivity of the solution to each random variable and its distribution parameters. Stability is computed using Spencer's method of slices coupled to a finite element seepage code to directly evaluate pore pressure. In addition, response surface method solutions are compared to the direct reliability solution to assess response surface accuracy. The results show that blanket layer thickness is more important than most strength parameters and that the uncertainty in hydraulic conductivity is less important to the variance in safety factor for the case of a low-permeability blanket layer. Numerical challenges caused by implementing a finite element limit-state function are discussed. Response surface methods are found to give a reasonable approximation to the direct reliability solution when the design point is between response surface fitting points.

Lanzafame, Robert↗