Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational physics applications”

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 577 records · Page 32

Physics-constrained Gaussian process model for prediction of hydrodynamic interactions between wave energy converters in an array

To improve the efficiency of wave farms and achieve maximum power generation, the layout of wave energy converters (WECs) in an array needs to be carefully designed so that the hydrodynamic interactions can be positively exploited. For this, the hydrodynamic characteristics of the WEC array in different layouts need to be calculated. However, such calculations using numerical models usually entail significant computational cost, especially for large arrays of WECs. To address the computational challenge, a physics-constrained Gaussian process (GP) model is proposed to replace the original expensive numerical model and predict the hydrodynamic characteristics of the WECs for any array layout. By exploring the relationship between the WEC array (i.e., the input) and different hydrodynamic characteristics (i.e., the output), here we summarize a set of physical constraints/features, including invariance, symmetry, and additivity. This prior knowledge about the input-output relationship is then directly embedded in the constructed GP model through the design of physics-constrained kernels. In particular, a double-sum invariant kernel is first developed to incorporate the invariance and symmetry features, and then an additive kernel is developed to incorporate the additive feature of the problem. The invariant kernel and the additive kernel are then integrated to construct the physics-constrained GP model. Compared to the standard GP model, the proposed physics-constrained GP models require less training data to achieve the desired accuracy in predicting the hydrodynamic characteristics and are also less vulnerable to the curse of dimensionality (i.e., good scalability for large arrays) due to the use of an additive kernel. The efficiency, accuracy, and scalability of the proposed approach are demonstrated through an application to predict the hydrodynamic characteristics for WEC arrays of different sizes and layouts.

16 TIDAL AND WAVE POWER↗

Nondestructive Evaluation (NDE) of Cable Anomalies using Frequency Domain Reflectometry (FDR) and Spread Spectrum Time Domain Reflectometry (SSTDR)

This report presents a comparative assessment of the performance of frequency domain reflectometry (FDR) and spread spectrum time domain reflectometry (SSTDR) in detecting a wide range of electrical cable anomalies. All tests and results reported herein were performed at the PNNL Accelerated and Real-Time Environmental Nodal Assessment (ARENA) cable and motor test bed. The primary objective of this work was to evaluate the effectiveness of SSTDR, a fledgling cable monitoring technique that shows promise for application in online monitoring of energized cable systems, against FDR, an offline technique widely employed in the nuclear power plant (NPP) industry. FDR tests are becoming more widely used in nuclear power plant cable aging management and test programs – particularly for low voltage cables. FDR capabilities for these kinds of tests have been reported by PNNL and others. The FDR test is performed on de-energized cables by connecting the FDR instrument to two of the cable conductors, or one conductor and the shield. A broad band low voltage (< 5 V) chirp is introduced in the cable, and any reflected response is captured in the frequency domain. The captured reflection is then processed by performing an inverse Fourier transform to a time domain response which can then be converted to a distance response based on the cable velocity of propagation (VoP). SSTDR measurements are functionally similar to FDR measurements in that a broad-band voltage signal composed of a square or sine wave modulated pseudo-random sequence of chips (< 5 volts), is injected onto one of the cable conductors. The injected signal will experience partial energy reflection and transmission at each impedance discontinuity along the transmission line. Any reflected response is detected by computing a cross-correlation between the reflected signals and a delayed copy of the incident SSTDR signal. the time delay for the reflected signal to experience the best matched correlation with the incident signal, indicates the travel time for the signal to reach a change in impedance. By knowing this time delay and velocity of propagation (VoP) of the signal, one can compute the physical distance. A big advantage that SSTDR measurements have over other methods is the ability to be connected to energized or live wires (currently up to 1kV) thereby enabling online monitoring of cables. SSTDR has been used successfully in several applications, e.g., aircraft, rail, and photovoltaic systems. In this work FDR and SSTDR cable assessment techniques were used to characterize a variety of cable anomalies and faults including: (1) Presence or absence of a motor; (2) Ground faults and short circuit faults; (3) Moist environments and water ingress faults; (4) Accelerated thermal aging. Both shielded and non-shielded cables were evaluated in this report. Offline measurements were made using FDR and online measurements were made by SSTDR for a range of test scenarios. Based on the results across all cable anomalies evaluated in this study, FDR displayed high sensitivity towards cable condition assessment, while SSTDR showed promise for future application in monitoring NPP cable systems. However, further developments are suggested to improve the resolution and sensitivity of SSTDR towards faults and anomalies in low voltage cables.rt presents a comparative

42 ENGINEERING↗

Toward High-Throughput Deposition of III-V Materials and Devices Using Halide Vapor Phase Epitaxy: Preprint

III-V devices are used in countless applications due to their excellent physical properties. They could become more prevalent, especially in area-intensive applications such as solar power, if they can achieve significant cost decreases through increasing scale. The development of high-throughput growth systems can help to achieve this scale, leading to the use of III-V devices in areas where they are not currently economically feasible. Here, we describe a pilot production, pseudo inline HVPE reactor with the potential to greatly increase the throughput of III-V devices. We show computational modeling results that both informed system design and the understanding of the impact of different process parameters on the deposition. We show the throughput possibilities of this reactor with an example solar cell device design but note that this system is agnostic to the device structure and can be used to increase the throughput of lasers, LEDs, transistors, and more III-V devices.

high throughput↗

US Department of Energy, Office of Science High Performance Computing Facility Operational Assessment 2019 Oak Ridge Leadership Computing Facility

Oak Ridge National Laboratory's (ORNL's) Leadership Computing Facility (OLCF) continues to surpass its operational target goals: supporting users; delivering fast, reliable computational ecosystems; creating innovative solutions for high performance computing (HPC) needs; and managing risks, safety, and security associated with operating some of the most powerful computers in the world. The results can be seen in the cutting-edge science conducted by users and the praise from the research community. Calendar year (CY) 2019 was a big year as OLCF staff ran five world-class resources (the leadershipclass computers Titan and Summit, the large analysis cluster called Eos, and the massive parallel filesystems called Atlas and Alpine)) and also began power and cooling upgrades for a 2021 exascale system called Frontier. While continuing exceptional operation of Titan, Eos, and Rhea, the OLCF released the Summit supercomputer for production on January 1, 2019. Summit debuted as the most capable and efficient system in its class and has been recognized as the most powerful system in the world for its performance on both the high performance linpack (HPL) and conjugate gradient (HPCG) benchmark applications since June 2018 according to TOP500. Summit represents the culmination of a multiyear effort between the OLCF, IBM, NVIDIA, and Mellanox to deliver a system that is unmatched for modeling, simulation, data analysis, and learning. To hit the ground running with science-ready applications on day one, application teams worked closely with the OLCF through the Center for Accelerated Application Readiness (CAAR) program for years in advance of the Summit deployment. CY 2019 was filled with outstanding results and accomplishments: a very high rating from users on overall satisfaction for the sixth year in a row; a tremendous amount of core-hours delivered to researchers from two leadership-class systems; and success in delivering on the allocation split of roughly 60%, 30%, and 10% of core-hours offered for the Innovative and Novel Computational Impact on Theory and Experiment (INCITE), Advanced Scientific Computing Research Leadership Computing Challenge (ALCC), and Director's Discretionary (DD) programs, respectively (see Operational Performance section). These accomplishments, coupled with the high utilization rates (overall and capability usage), represent the fulfillment of the promise of both leadership-class machines: efficient facilitation of leadership-class computational applications. Table ES.1 presents a summary of the 2019 OLCF metric targets and the associated results. More information can be found in the Operational Performance section for each OLCF resource. The scientific accomplishments of OLCF users are a strong indication of long-term operational success, with publications this year in such notable journals and publications as Nature, Nature Physics, Nature Plants, Physical Review X, Journal of the American Physical Society, Cell, Nano Letters, and Trends in Biotechnology. Crucial domain-specific discoveries facilitated by resources at the OLCF are described in the High Performance Computing Facility Operational Assessment 2019 Oak Ridge Leadership Computing Facility (OAR) Strategic Results section. For example, researchers used Summit to pinpoint and understand the production of proteins from genetic information, including mutations and the functional expression of disease (Section 8.2).

97 MATHEMATICS AND COMPUTING↗

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

Probability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude-based model for describing advective-dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub-continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a manyto-one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub-continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Additionally, differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude-based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

54 ENVIRONMENTAL SCIENCES↗

Integrating magnetic Co-nanopillars in a NbN-based VAN thin film as a multifunctional hybrid metamaterial

NbN thin films are considered as a promising alternative candidate to conventional Al superconducting electrodes for Josephson junction devices, superconducting qubits and quantum logic circuits. Herein, by co-growing ferromagnetic Co with superconducting NbN, self-assembled NbN–Co vertically aligned nanocomposite (VAN) thin films have been successfully deposited on a MgO substrates, with NbN in its cubic superconducting phase. The obtained VAN metamaterials exhibit epitaxial growth of NbN with uniformly dispersed Co nanoplates. Further physical property characterization reveals that the NbN–Co VAN presents ferromagnet properties with strong out-of-plane magnetic anisotropy and optical anisotropy. This hybrid metamaterial system could find future applications in superconducting spintronic devices and quantum computing devices.

Zhang, Yizhi [Purdue Univ., West Lafayette, IN (Un↗

Reconfigurable plasma-dielectric hybrid photonic crystal as a platform for electromagnetic wave manipulation and computing

The introduction of 3D printing has enabled fabrication of photonic crystal devices with complex crystal structures that would be challenging to construct using other fabrication methods. In this paper, we construct and characterize a photonic crystal consisting of two 4 x 8 x 8 cubic lattices composed of spherical silicon nitride elements straddling a layer of 8 plasma discharge tubes, creating a 3D hybrid plasma photonic crystal device. Integrating under-dense gaseous plasma elements provides a unique coupling dynamic between the dielectric spheres and the cylindrical plasma discharges, creating a monolithic hybrid photonic crystal with solid state and reconfigurable elements. The device has resonant modes that have attenuation peaks that are either switchable, tunable in amplitude, or tunable in frequency with variations in plasma density. The response of these bands seen with varying plasma density is confirmed through simulations when effects due to the heating of the photonic crystal from the gaseous plasma elements are accounted for in the experiments. We discuss how this reconfigurable device may be used and expanded upon for applications in photonic artificial neural networks and optical computing systems.

3D printing↗

The “good,” the “bad,” and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl 2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca 2+ –Cl - interactions were primarily “hidden” in the Ca-NDIS signal due to overlapping with Ca 2+ –O w and Ca 2+ –H w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models “on average,” the question “which model is the best” has not been answered but replaced by the question “which model is better for a given research.” An overall “good” model can still be inappropriate or, in some instances, “bad” and, unfortunately, produce erroneous results. Additionally, the accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SAGIPS: a physics-inspired scalable asynchronous generative inverse-problem solver

Abstract Solving large-scale inverse problems using deep-learning algorithms have become an essential part of modern research and industrial applications. The complexity of the underlying inverse problem may require the utilization of high performance computing systems which poses a challenge on the algorithmic design of the inverse problem solver. Most deep learning algorithms require, due to their design, custom parallelization techniques in order to be resource efficient while showing a reasonable convergence. In this paper we introduce a S calable A synchronous G enerative I nverse P roblem S olver (SAGIPS) on high-performance computing systems. We present a workflow that utilizes an asynchronous ring-allreduce algorithm to transfer the gradients of the generator network across multiple GPUs. Experiments with a scientific proxy application demonstrate that SAGIPS shows near linear weak scaling, together with a convergence quality that is comparable to traditional methods. The approach presented here allows leveraging Generative Adverserial Network across multiple GPUs, promising advancements in solving complex inverse problems at scale.

97 MATHEMATICS AND COMPUTING↗

Electric field control of magnetism

Electric field control of magnetism is an extremely exciting area of research, from both a fundamental science and an applications perspective and has the potential to revolutionize the world of computing. To realize this will require numerous further innovations, both in the fundamental science arena as well as translating these scientific discoveries into real applications. Thus, this article will attempt to bridge the gap between condensed matter physics and the actual manifestations of the physical concepts into applications. We have attempted to paint a broad-stroke picture of the field, from the macroscale all the way down to the fundamentals of spin–orbit coupling that is a key enabler of the physics discussed. Here, we hope it will help spur more translational research within the broad materials physics community. Needless to say, this article is written on behalf of a large number of colleagues, collaborators and researchers in the field of complex oxides as well as current and former students and postdocs who continue to pursue cutting-edge research in this field.

36 MATERIALS SCIENCE↗

Exascale applications: skin in the game

As noted in Wikipedia, skin in the game refers to having ‘incurred risk by being involved in achieving a goal’, where ‘ skin is a synecdoche for the person involved, and game is the metaphor for actions on the field of play under discussion’. For exascale applications under development in the US Department of Energy Exascale Computing Project, nothing could be more apt, with the skin being exascale applications and the game being delivering comprehensive science-based computational applications that effectively exploit exascale high-performance computing technologies to provide breakthrough modelling and simulation and data science solutions. These solutions will yield high-confidence insights and answers to the most critical problems and challenges for the USA in scientific discovery, national security, energy assurance, economic competitiveness and advanced healthcare. This article is part of a discussion meeting issue ‘Numerical algorithms for high-performance computational science’.

97 MATHEMATICS AND COMPUTING↗

Four lectures on Euler integrals

These lecture notes provide a self-contained introduction to Euler integrals, which are frequently encountered in applications. In particle physics, they arise as Feynman integrals or string amplitudes. Our four selected topics demonstrate the diverse mathematical techniques involved in the study of Euler integrals, including polyhedral geometry, very affine varieties, differential equations, and computational algebra.

Matsubara-Heo, Saiei-Jaeyeong↗

An open-source hybrid unstructured mesh - CAD fusion multiphysics analysis workflow in SALAMANDER

Plasma facing components in fusion devices will endure extreme neutron and heat fluxes. To facilitate their design using simulation tools, the open-source Fusion Module, Fusion ENergy Integrated multiphys-X (FENIX) framework is being developed to model these components with a high-fidelity multi-physics multi-dimensional approach. It can iteratively resolve couplings between all the physics at play, from neutron radiation, to thermomechanics, to near-wall plasma dynamics. This framework is based on the Multiphysics Object Oriented Simulation Environment (MOOSE), which is developed by a collaboration of US National Laboratories since 2008, for advanced nuclear, geomechanics simulations and other applications. FENIX couples numerous simulation tools, including OpenMC, the Tritium Migration Analysis Program v8, the NekRS CFD software, and most MOOSE modules. For the coupling of radiation transport and other physics, FENIX supports a hybrid workflow between Computer Assisted Design (CAD) and unstructured mesh geometries. The CAD can be generated from skinning the unstructured mesh, to enable a coarse geometry for efficient particle transport, but still resolving the local material compositions and temperature gradients. Neutron transport is performed using DAGMC on the CAD, and Cardinal, integrated in FENIX, maps tallied quantities, such as the heat deposition or the tritium generation rates, from a tally volumetric mesh to the other physics’ unstructured mesh. This coupling was exercised on a simplified tokamak geometry, coupling neutron transport with the heat conduction equation, and on a monoblock divertor problem, coupling additionally with tritium migration. Mesh convergence studies highlight the importance of the mapping conservativeness. Coupling with thermo-mechanics is further enabled by the generalization of the approach to moving meshes. The presentation will include these coupled analysis as well as an update on status of the FENIX framework.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electronic properties and optical spectra of donor–acceptor conjugated organic polymers

Organic semiconductors based on conjugated donor-acceptor (D–A) polymers are a unique platform for electronic, spintronic, and energy-harvesting devices. Understanding the electronic structure of D–A polymers with a small band gap is essential for developing next-generation technologies. Here, we investigate the electronic structure and optical spectra of cyclopentadithiophene-based closed/open-shell D–A polymers using density functional theory and the Bethe–Salpeter equation based on G 0 W 0 approximation. We explored the role of different acceptor units and chemical substitutions on the structural changes and, more importantly, electronic, optical, and dielectric behavior. We found that the computed first exciton peak of the polymers agreed well with the available experimentally measured optical gap. Furthermore, D–A polymers with open-shell character display higher dielectric constant than the closed-shell polymers. We show that the exceptional performance of polycyclopentadithiophene-thiophenylthiadiazoloquinoxaline (PCPDT-TTQ) as a scalable n-type material for Faradaic supercapacitors can be partly ascribed to its elevated dielectric constant. Consequently, these D–A polymers, characterized by their high dielectric constants, exhibit significant potential for various applications, including energy storage, organic electronics, and the production of dielectric films.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular NMR shieldings, J -couplings, and magnetizabilities from numeric atom-centered orbital based density-functional calculations

This paper reports and benchmarks a new implementation of nuclear magnetic resonance shieldings, magnetizabilities, and J-couplings for molecules within semilocal density functional theory, based on numeric atom-centered orbital (NAO) basis sets. NAO basis sets are attractive for the calculation of these nuclear magnetic resonance (NMR) parameters because NAOs provide accurate atomic orbital representations especially near the nucleus, enabling high-quality results at modest computational cost. Moreover, NAOs are readily adaptable for linear scaling methods, enabling efficient calculations of large systems. Here, the paper has five main parts: (1) It reviews the formalism of density functional calculations of NMR parameters in one comprehensive text to make the mathematical background available in a self-contained way. (2) The paper quantifies the attainable precision of NAO basis sets for shieldings in comparison to specialized Gaussian basis sets, showing similar performance for similar basis set size. (3) The paper quantifies the precision of calculated magnetizabilities, where the NAO basis sets appear to outperform several established Gaussian basis sets of similar size. (4) The paper quantifies the precision of computed J-couplings, for which a group of customized NAO basis sets achieves precision of ~Hz for smaller basis set sizes than some established Gaussian basis sets. (5) The paper demonstrates that the implementation is applicable to systems beyond 1000 atoms in size.

74 ATOMIC AND MOLECULAR PHYSICS↗

Highly accelerated life testing (HALT): A review from a statistical perspective

Despite its use in one form or another for at least four decades, HALT and related techniques [e.g., highly accelerated-stress screening (HASS) and stress audits (HASA)] are not well understood within the statistical community and remain controversial. This largely reflects a conflict in motivation between engineers, testing under harsh conditions to discover and eliminate failure modes, and statisticians, taking a more cautious approach to develop quantitative estimates of parameters such as mean time between failures (MTBF). Here, this review article will clarify HALT concepts and methods and explain where it fits within the universe of methods that involve the application of accelerating factors to compress the time required to evaluate or enhance product reliability. A major distinction is between methods such as HALT, a high-stress test-analyze-fix-test iterative process directed at improving reliability by discovering and fixing weak points in a design, and quantitative accelerated life testing (QALT), whose goal is the estimation of product life for a fixed design. We discuss methods such as physics of failure that offer some hope of bridging the gap between the qualitative nature of HALT, and purely quantitative statistical methods. We present a variety of engineering applications of HALT including metal fatigue, piping and pressure vessels, structural damage, radiation damage, and rotating machinery. We also discuss potential synergies between HALT and QALT, such as rapid identification, through HALT, of failure modes requiring quantitative analysis. For further study, extensive references to the applicable literature are provided as well as an appendix that describes related methods.

97 MATHEMATICS AND COMPUTING↗

At-power subcritical multiplication in the Advanced Test Reactor during nuclear requalification testing

Power division information during nuclear requalification of the Advanced Test Reactor (ATR) is of considerable interest as an importance function for observed changes to core reactivity. The degree to which a given physical subdivision of a critical reactor acts as a neutron source for other lobes is not analytically characterized for general application. When ATR operates at power, individual power-producing lobes rely on each other as neutron sources in order to maintain constant power, which in general requires either exactly critical multiplication within a reactor or an external neutron source. Here, this work shows that fuel element and lobe powers in ATR can be related with subcritical multiplication theory. Subcritical multiplication factors are computed with a physically validated analytical method based on actual at-power operation, quantifying for each lobe its dependence on other lobes as an external neutron source. This explanation is significant for ATR due to the desire to irradiate a large variety of experiments simultaneously, each having its impact on the core neutron population. For any physical subdivision of any other critical reactor, it is likewise true that the subdivision undergoes only subcritical multiplication.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Reconstructing patchy reionization with deep learning

The precision anticipated from next-generation cosmic microwave background (CMB) surveys will create opportunities for characteristically new insights into cosmology. Secondary anisotropies of the CMB will have an increased importance in forthcoming surveys, due both to the cosmological information they encode and the role they play in obscuring our view of the primary fluctuations. Quadratic estimators have become the standard tools for reconstructing the fields that distort the primary CMB and produce secondary anisotropies. While successful for lensing reconstruction with current data, quadratic estimators will be suboptimal for the reconstruction of lensing and other effects at the expected sensitivity of the upcoming CMB surveys. In this paper we describe a convolutional neural network, ResUNet-CMB, that is capable of the simultaneous reconstruction of two sources of secondary CMB anisotropies, gravitational lensing and patchy reionization. Furthermore, we show that the ResUNet-CMB network significantly outperforms the quadratic estimator at low noise levels and is not subject to the lensing-induced bias on the patchy reionization reconstruction that would be present with a straightforward application of the quadratic estimator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗