Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimization analysis”

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,531 records · Page 85

Quantifying and reducing model-form uncertainties in Reynolds-averaged Navier–Stokes simulations: A data-driven, physics-informed Bayesian approach

Despite their well-known limitations, Reynolds-Averaged Navier–Stokes (RANS) models are still the workhorse tools for turbulent flow simulations in today's engineering analysis, design and optimization. While the predictive capability of RANS models depends on many factors, for many practical flows the turbulence models are by far the largest source of uncertainty. As RANS models are used in the design and safety evaluation of many mission-critical systems such as airplanes and nuclear power plants, quantifying their model-form uncertainties has significant implications in enabling risk-informed decision-making. In this work we develop a data-driven, physics-informed Bayesian framework for quantifying model-form uncertainties in RANS simulations. Uncertainties are introduced directly to the Reynolds stresses and are represented with compact parameterization accounting for empirical prior knowledge and physical constraints (e.g., realizability, smoothness, and symmetry). An iterative ensemble Kalman method is used to assimilate the prior knowledge and observation data in a Bayesian framework, and to propagate them to posterior distributions of velocities and other Quantities of Interest (QoIs). We use two representative cases, the flow over periodic hills and the flow in a square duct, to evaluate the performance of the proposed framework. Both cases are challenging for standard RANS turbulence models. Simulation results suggest that, even with very sparse observations, the obtained posterior mean velocities and other QoIs have significantly better agreement with the benchmark data compared to the baseline results. At most locations the posterior distribution adequately captures the true model error within the developed model form uncertainty bounds. The framework is a major improvement over existing black-box, physics-neutral methods for model-form uncertainty quantification, where prior knowledge and details of the models are not exploited. This approach has potential implications in many fields in which the governing equations are well understood but the model uncertainty comes from unresolved physical processes. - Highlights: • Proposed a physics–informed framework to quantify uncertainty in RANS simulations. • Framework incorporates physical prior knowledge and observation data. • Based on a rigorous Bayesian framework yet fully utilizes physical model. • Applicable for many complex physical systems beyond turbulent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Learning occupants’ indoor comfort temperature through a Bayesian inference approach for office buildings in United States

A carefully chosen indoor comfort temperature as the thermostat set-point is the key to optimizing building energy use and occupants’ comfort and well-being. ASHRAE Standard 55 or ISO Standard 7730 uses the PMV-PPD model or the adaptive comfort model that is based on small-sized or outdated sample data, which raises questions on whether and how ranges of occupant thermal comfort temperature should be revised using more recent larger-sized dataset. In this paper, a Bayesian inference approach has been used to derive new occupant comfort temperature ranges for U.S. office buildings using the ASHRAE Global Thermal Comfort Database. Bayesian inference can express uncertainty and incorporate prior knowledge. The comfort temperatures were found to be higher and less variable at cooling mode than at heating mode, and with significant overlapped variation ranges between the two modes. The comfort operative temperature of occupants varies between 21.9 and 25.4 °C for the cooling mode with a median of 23.7 °C, and between 20.5 and 24.9 °C for the heating mode with a median of 22.7 °C. These comfort temperature ranges are similar to the current ASHRAE standard 55 in the heating mode but 2–3 °C lower in the cooling mode. Finally, the results of this study could be adopted as more realistic thermostat set-points in building design, operation, control optimization, energy performance analysis, and policymaking.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Ligand Many-Body Expansion as a General Approach for Accelerating Transition Metal Complex Discovery

Methods that accelerate the evaluation of molecular properties are essential for chemical discovery. While some degree of ligand additivity has been established for transition metal complexes, it is underutilized in asymmetric complexes, such as the square pyramidal coordination geometries highly relevant to catalysis. To develop predictive methods beyond simple additivity, we apply a many-body expansion to octahedral and square pyramidal complexes and introduce a correction based on adjacent ligands (i.e., the cis interaction model). We first test the cis interaction model on adiabatic spin-splitting energies of octahedral Fe(II) complexes, predicting DFT-calculated values of unseen binary complexes to within an average of 1.4 kcal/mol. Uncertainty analysis reveals the optimal basis, comprising the homoleptic and mer symmetric complexes. We next show that the cis model (i.e., the cis interaction model solved for the optimal basis) infers both DFT- and CCSD(T)-calculated model catalytic reaction energies to within 1 kcal/mol on average. The cis model predicts low-symmetry complexes with reaction energies outside the range of binary complex reaction energies. We observe that trans interactions are unnecessary for most monodentate systems but can be important for some combinations of ligands, such as complexes containing a mixture of bidentate and monodentate ligands. Lastly, we demonstrate that the cis model may be combined with Δ-learning to predict CCSD(T) reaction energies from exhaustively calculated DFT reaction energies and the same fraction of CCSD(T) reaction energies needed for the cis model, achieving around 30% of the error from using the CCSD(T) reaction energies in the cis model alone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hardware-accelerated inference for real-time gravitational-wave astronomy

The field of transient astronomy has seen a revolution with the first gravitational-wave detections and the arrival of multi-messenger observations they enabled. Transformed by the first detection of binary black hole and binary neutron star mergers, computational demands in gravitational-wave astronomy are expected to grow by at least a factor of two over the next five years as the global network of kilometer-scale interferometers are brought to design sensitivity. With the increase in detector sensitivity, real-time delivery of gravitational-wave alerts will become increasingly important as an enabler of multi-messenger followup. In this work, we report a novel implementation and deployment of deep learning inference for real-time gravitational-wave data denoising and astrophysical source identification. This is accomplished using a generic Inference-as-a-Service model that is capable of adapting to the future needs of gravitational-wave data analysis. Overall, our implementation allows seamless incorporation of hardware accelerators and also enables the use of commercial or private (dedicated) as-a-service computing. Based on our results, we propose a paradigm shift in low-latency and offline computing in gravitational-wave astronomy. Such a shift can address key challenges in peak-usage, scalability and reliability, and provide a data analysis platform particularly optimized for deep learning applications. The achieved sub-millisecond scale latency will also be relevant for any machine learning-based real-time control systems that may be invoked in the operation of near-future and next generation ground-based laser interferometers, as well as the front-end collection, distribution and processing of data from such instruments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Benchmarking second and third-generation sequencing platforms for microbial metagenomics

Shotgun metagenomic sequencing is a common approach for studying the taxonomic diversity and metabolic potential of complex microbial communities. Current methods primarily use second generation short read sequencing, yet advances in third generation long read technologies provide opportunities to overcome some of the limitations of short read sequencing. Here, we compared seven platforms, encompassing second generation sequencers (Illumina HiSeq 300, MGI DNBSEQ-G400 and DNBSEQ-T7, ThermoFisher Ion GeneStudio S5 and Ion Proton P1) and third generation sequencers (Oxford Nanopore Technologies MinION R9 and Pacific Biosciences Sequel II). We constructed three uneven synthetic microbial communities composed of up to 87 genomic microbial strains DNAs per mock, spanning 29 bacterial and archaeal phyla, and representing the most complex and diverse synthetic communities used for sequencing technology comparisons. Our results demonstrate that third generation sequencing have advantages over second generation platforms in analyzing complex microbial communities, but require careful sequencing library preparation for optimal quantitative metagenomic analysis. Our sequencing data also provides a valuable resource for testing and benchmarking bioinformatics software for metagenomics.

59 BASIC BIOLOGICAL SCIENCES↗

In situ self-assembled organic interface layers for the controlled growth of oligothiophene thin films on ferroelectric Pb(Zr{sub 0.2}Ti{sub 0.8})O{sub 3}

We introduce an in situ vacuum procedure for the optimal preparation and analysis of self-assembled monolayers (SAMs) as used in organic molecular electronics on ferroelectric lead zirconate titanate (PZT) substrates. Excellent ordering of oligothiophene semiconductor layers is heavily promoted through the presence of an interfacial bi-functional SAM layer that binds to both the oxidic PZT surface and the organic semiconductor molecules. The described method can be extended to other material combinations, featuring a variety of substrate materials and molecular functionalities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Periodic boundary conditions for QM/MM calculations: Ewald summation for extended Gaussian basis sets

An implementation of Ewald summation for use in mixed quantum mechanics/molecular mechanics (QM/MM) calculations is presented, which builds upon previous work by others that was limited to semi-empirical electronic structure for the QM region. Unlike previous work, our implementation describes the wave function's periodic images using “ChElPG” atomic charges, which are determined by fitting to the QM electrostatic potential evaluated on a real-space grid. This implementation is stable even for large Gaussian basis sets with diffuse exponents, and is thus appropriate when the QM region is described by a correlated wave function. Derivatives of the ChElPG charges with respect to the QM density matrix are a potentially serious bottleneck in this approach, so we introduce a ChElPG algorithm based on atom-centered Lebedev grids. The ChElPG charges thus obtained exhibit good rotational invariance even for sparse grids, enabling significant cost savings. Detailed analysis of the optimal choice of user-selected Ewald parameters, as well as timing breakdowns, is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The electric dipole moment of cobalt monoxide, CoO

A number of low-rotational lines of the E{sup 4}Δ{sub 7/2} ← X{sup 4}Δ{sub 7/2} (1,0) band system of cobalt monoxide, CoO, were recorded field free and in the presence of a static electric field. The magnetic hyperfine parameter, h{sub 7/2}, and the electron quadrupole parameter, eQq{sub 0}, for the E{sup 4}Δ{sub 7/2}(υ = 1) state were optimized from the analysis of the field-free spectrum. The permanent electric dipole moment, μ{sup -vector}{sub el}, for the X{sup 4}Δ{sub 7/2} (υ = 0) and E{sup 4}Δ{sub 7/2} (υ = 1) states were determined to be 4.18 ± 0.05 D and 3.28 ± 0.05 D, respectively, from the analysis of the observed Stark spectra of F′ = 7 ← F″ = 6 branch feature in the Q(7/2) line and the F′ = 8 ← F″ = 7 branch feature in the R(7/2) line. The measured dipole moments of CoO are compared to those from theoretical predictions and the trend across the 3d-metal monoxide series discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gigahertz-gated InGaAs/InP single-photon detector with detection efficiency exceeding 55% at 1550 nm

We report on a gated single-photon detector based on InGaAs/InP avalanche photodiodes (APDs) with a single-photon detection efficiency exceeding 55% at 1550 nm. Our detector is gated at 1 GHz and employs the self-differencing technique for gate transient suppression. It can operate nearly dead time free, except for the one clock cycle dead time intrinsic to self-differencing, and we demonstrate a count rate of 500 Mcps. We present a careful analysis of the optimal driving conditions of the APD measured with a dead time free detector characterization setup. It is found that a shortened gate width of 360 ps together with an increased driving signal amplitude and operation at higher temperatures leads to improved performance of the detector. We achieve an afterpulse probability of 7% at 50% detection efficiency with dead time free measurement and a record efficiency for InGaAs/InP APDs of 55% at an afterpulse probability of only 10.2% with a moderate dead time of 10 ns.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

HIGH-RESOLUTION FOURIER-TRANSFORM MICROWAVE SPECTROSCOPY OF METHYL- AND DIMETHYLNAPTHALENES

High-resolution pure rotational spectra of four alkylnaphthalenes were measured in the range of 6–15 GHz using a molecular-beam Fourier-transform microwave spectrometer. Both a- and b-type transitions were observed for 1-methylnaphthalene (1-MN), 1,2-dimethylnaphthalene (1,2-DMN), and 1,3-dimethylnaphthalene (1,3-DMN); only a-type transitions were observed for 2-methylnaphthalene (2-MN). Geometry optimization and vibrational analysis calculations at the B3LYP/6-311++G(d,p) level of theory aided in the assignments of the spectra and the characterization of the structures. Differences between the experimental and predicted rotational constants are small, and they can be attributed in part to low-lying out-of-plane vibrations, which distort the alkylnaphthalenes out of their equilibrium geometries. Splittings of rotational lines due to methyl internal rotation were observed in the spectra of 2-MN, 1,2-DMN, and 1,3-DMN, and allowed for the determination of the barriers to methyl internal rotation, which are compared to values from density functional theory calculations. All four species are moderately polar, so they are candidate species for detection by radio astronomy, by targeting the transition frequencies reported here.

79 ASTRONOMY AND ASTROPHYSICS↗

The fraction of muon tracks in cosmic neutrinos

The study of the distintive signatures of the ultra high energy events recently seen by IceCube [1-4] can allow to single the neutrino origin out. The detection of tau neutrinos would be a clear way to prove that they come from cosmic distances, but at the highest energies currently seen, about 1 PeV, an experimental characterization of tau events is difficult. The study of the fraction of the muon tracks seems more promising. In fact, for any initial composition, because of the occurrence of flavor oscillations and despite their uncertainties, the fraction of muon tracks in the cosmic neutrinos is smaller than the one of atmospheric neutrinos, even hypothesizing an arbitrarily large contribution from charmed mesons. A good understanding of the detection efficiencies and the optimization of the analysis cuts, along with a reasonable increase in the statistics, should provide us with a significant test of the cosmic origin of these events.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Search for the weak decay ψ (3686) → $Λ^{+}_{c}$ $\overline{Σ}$ - + c.c. *

Using (448.1 ± 2.9) x 10 6 ψ (3686) events collected with the BESIII detector, we perform the first search for the weak baryonic decay ψ (3686) → $Λ^{+}_{c}$ $\overline{Σ}$ - + c.c. . The analysis procedure is optimized using a blinded method. No significant signal is observed, and the upper limit on the branching fraction ($\mathcal{B}$) of ψ (3686) → $Λ^{+}_{c}$ $\overline{Σ}$ - + c.c. is set as 1.4 x 10 -5 at the 90% confidence level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Understanding effects of printhead geometry in aerosol jet printing

Aerosol jet printing offers a versatile, high-resolution digital patterning capability broadly relevant to flexible and printed electronic systems. Despite its promise and numerous demonstrations, the theoretical principles driving process outputs have not been thoroughly explored. In this study, a custom-built, modular printing system is developed to provide a head-to-head comparison of two print nozzle geometries to better understand the technology. Print resolution data from a range of process parameters are analyzed using a support vector machine framework. The linear deposition rate is identified as a key variable, which can confound careful studies of printing performance. Taking this into account, a clear difference is observed between the printheads, corresponding to a difference in resolution of 57% ± 11% under typical conditions. Models to understand differences in aerodynamic and mass transport effects identify enhanced drying within the NanoJet printhead as a likely cause of this difference. Overall, this study provides improved understanding of the aerosol jet printing process, including valuable insight to inform process optimization, robust data analysis, ink formulation, and printer geometric design.

42 ENGINEERING↗

Two-nucleon 𝑆-wave interactions at the SU(3) flavor-symmetric point with 𝑚 𝑢⁢𝑑 ≈ 𝑚$^{phys}_𝑠$: A first lattice QCD calculation with the stochastic Laplacian Heaviside method

We report on the first application of the stochastic Laplacian Heaviside method for computing multiparticle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a set of positive-definite two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the SU(3)-flavor-symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of 𝑎 ≈0.086 fm, lattice volume of 𝐿 = 48⁢𝑎 ≈ 4.1 fm and pion mass 𝑚 𝜋 ≈ 714 MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Novel technique for the study of pileup events in cryogenic bolometers

Precise characterization of detector time resolution is of crucial importance for next-generation cryogenic-bolometer experiments searching for neutrinoless double-beta decay, such as CUPID, in order to reject background due to pile-up of two-neutrino double-beta decay events. In this paper, we describe a technique developed to study the pile-up rejection capability of cryogenic bolometers. Our approach, which consists of producing controlled pile-up events with a programmable waveform generator, has the benefit that we can reliably and reproducibly control the time separation and relative energy of the individual components of the generated pile-up events. The resulting data allow us to optimize and benchmark analysis strategies to discriminate between individual and pile-up pulses. We describe a test of this technique performed with a small array of detectors at the Laboratori Nazionali del Gran Sasso, in Italy; we obtain a 90% rejection efficiency against pulser-generated pile-up events with rise time of ~15 ms down to time separation between the individual events of 2 ms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Search for the nonresonant production of Higgs boson pairs via gluon fusion and vector-boson fusion in the $b\bar{b}$τ + τ – final state in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for the nonresonant production of Higgs boson pairs in the HH → $b\overline{b}$$\tau^{+}\tau^{-}$ channel is performed using 140 fb -1 of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. The analysis strategy is optimized to probe anomalous values of the Higgs boson self-coupling modifier $\kappa_{\tau}$ and of the quartic HHVV (v = W,Z) coupling modifier $\kappa$ 2⁢V . No significant excess above the expected background from Standard Model processes is observed. An observed (expected) upper limit μ HH < 5.9⁢(3.3) is set at 95% confidence-level on the Higgs boson pair production cross section normalized to its Standard Model prediction. The coupling modifiers are constrained to an observed (expected) 95% confidence interval of -3.1 <$\kappa_{\tau}$ < 9.0 (-2.5 < $\kappa_{\tau}$ < 9.3) and -0.5 <$\kappa$ 2⁢V < 2.7 (-0.2 <$\kappa$ 2⁢V < 2.4), assuming all other Higgs boson couplings are fixed to the Standard Model prediction. The results are also interpreted in the context of effective field theories via constraints on anomalous Higgs boson couplings and Higgs boson pair production cross sections assuming different kinematic benchmark scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the release of cppxfel for processing X-ray free-electron laser images

As serial femtosecond crystallography expands towards a variety of delivery methods, including chip-based methods, and smaller collected data sets, the requirement to optimize the data analysis to produce maximum structure quality is becoming increasingly pressing. Here cppxfel , a software package primarily written in C++, which showcases several data analysis techniques, is released. This software package presently indexes images using DIALS (diffraction integration for advanced light sources) and performs an initial orientation matrix refinement, followed by post-refinement of individual images against a reference data set. Cppxfel is released with the hope that the unique and useful elements of this package can be repurposed for existing software packages. However, as released, it produces high-quality crystal structures and is therefore likely to be also useful to experienced users of X-ray free-electron laser (XFEL) software who wish to maximize the information extracted from a limited number of XFEL images.

97 MATHEMATICS AND COMPUTING↗