Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied computing”

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 271 records · Page 15

Can exascale computing and explainable artificial intelligence applied to plant biology deliver on the United Nations sustainable development goals?

Human population growth and accelerated climate change necessitate agricultural improvements using designer crop ideotypes (idealized plants that can grow in niche environments). Diverse and highly skilled research groups must integrate efforts to bridge the gaps needed to achieve international goals toward sustainable agriculture. Given the scale of global agricultural needs and the breadth of multiple types of omics data needed to optimize these efforts, explainable artificial intelligence (AI with a decipherable decision making process that provides a meaningful explanation to humans) and exascale computing (computers that can perform 1018 floating-point operations per second, or exaflops) are crucial. Accurate phenotyping and daily-resolution climatype associations are equally important for refining ideotype production to specific environments at various levels of granularity. In this article, we review advances toward tackling technological hurdles to solve multiple United Nations Sustainable Development Goals and discuss a vision to overcome gaps between research and policy.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling betatron radiation using particle-in-cell codes for plasma wakefield accelerator diagnostics

The analysis of plasma wakefield acceleration experimental measurements, particularly in the characterization of photons emitted through the betatron radiation mechanism, requires the development of accurate numerical models. These computational models are crucial for supporting modern instrumentation designed to measure the single-shot, double-differential angular-energy radiation spectra resulting from interactions between beams and plasmas. Motivated by the needs of such applications, this paper presents detailed numerical models of betatron radiation generated in beam-plasma acceleration experiments. These models are based on the integration of the Liénard-Wiechert (LW) potentials, applied to computed particle trajectories. The particle trajectories are generated using three distinct methods: first, by tracking particles through idealized fields in the blowout regime of PWFA; second, by obtaining trajectories using the fast quasistatic particle-in-cell (PIC) code quickpic; and third, obtaining trajectories from the fully self-consistent PIC code osiris. To ensure the accuracy and reliability of these models, the paper includes various benchmark tests using analytical expressions, as well as employing the PIC code epoch, which takes an alternative approach by using a Monte Carlo quantum electrodynamics (QED)-based radiation model. Additionally, the paper presents simulations of the expected experimental betatron radiation spectra, taking into account parameters relevant to PWFA and plasma photocathode experiments at the SLAC FACET-II facility.

Yadav, M. [University of California, Los Angeles, ↗

FY 2021 Mid-Year HPC Utilization Report M3UF-21IN0701018

Idaho National Laboratory (INL), supported by the Department of Energy Office of Nuclear Energy (DOE-NE) through the Nuclear Science User Facilities, provides access to supercomputer systems and data storage along with support staff for system management, software installation, cybersecurity and user support to the broader DOE-NE user community. Users include individuals at universities, industry, and government laboratories enabling a wide range of research and development and mission-supporting activities. The availability of high-performance computing (HPC) capabilities is a key foundation of collaboration and innovation in nuclear energy systems research. HPC resources and INL directly support the mission and objectives of DOE-NE. From October 2020 through March 2021, INL HPC capabilities were utilized by a diverse set of computing and applied researchers, for a wide range of research and engineering activities. This report focuses on current INL HPC systems and utilization from October 1, 2020 through March 31, 2021.

99 GENERAL AND MISCELLANEOUS↗

The effects of self-generated and applied magnetic fields on the computation of flow over a Mars return aerobrake

A CFD technique is developed to calculate the electromagnetic phenomena simultaneously with the fluid flow in the shock layer over an axisymmetric blunt body in a thermal-equilibrium chemical-nonequilibrium environment. The flowfield is solved using an explicit time-marching, first-order spatially accurate scheme. The electromagnetic phenomena are coupled to the real-gas flow solver through an iterative procedure. The electromagnetic terms introduce a strong stiffness, which was overcome by using significantly smaller time steps for the electromagnetic conservation equation. The technique is applied in calculating the flow over a Mars return aerobrake vehicle entering the Earth's atmosphere. For the case where no external field is applied, the electromagnetic effects have little impact on the flowfield.

Palmer, Grant↗

A Simple, Scalable Large Deformation Solid Mechanics Implementation in the MOOSE Framework

This article describes a large deformation solid mechanics solver implemented as part of the freely available and open source MOOSE finite element simulation framework. The article documents the choices made in developing the solid mechanics framework and describes novel formulations for the gradient operator and constitutive modeling framework made to simplify implementations of different coordinate systems, stabilized gradient operators, and different constitutive model inputs and outputs. In the process, the article describes a new formulation that casts objective integration of the Cauchy stress as a linear transformation of the small stress rate. Finally, the article presents key implementation details and examines the parallel efficiency of the solid mechanics solver implemented in MOOSE. The implementation retains a good weak scaling efficiency beyond 1,000 parallel processes. The article includes a discussion of the factors limiting the parallel efficiency of implicit, large deformation solid mechanics codes on current high-performance computers, with the main current limitation being the scalability of the algebraic multigrid methods used to solve the linearized equilibrium equations.

Applied computing → Computer-aided design↗

The Control Point Library Building System

The Earth Resources Observation System (EROS) Data Center in Sioux Falls, South Dakota distributes precision corrected Landsat MSS and RBV data. These data are derived from master data tapes produced by the Master Data Processor (MDP), NASA's system for computing and applying corrections to the data. Included in the MDP is the Control Point Library Building System (CPLBS), an interactive, menu-driven system which permits a user to build and maintain libraries of control points. The control points are required to achieve the high geometric accuracy desired in the output MSS and RBV data. This paper describes the processing performed by CPLBS, the accuracy of the system, and the host computer and special image viewing equipment employed.

Niblack, W.↗

Disaster risk and artificial intelligence: A framework to characterize conceptual synergies and future opportunities

Artificial intelligence (AI) methods have revolutionized and redefined the landscape of data analysis in business, healthcare, and technology. These methods have innovated the applied mathematics, computer science, and engineering fields and are showing considerable potential for risk science, especially in the disaster risk domain. The disaster risk field has yet to define itself as a necessary application domain for AI implementation by defining how to responsibly balance AI and disaster risk. (1) How is AI being used for disaster risk applications; and how are these applications addressing the principles and assumptions of risk science, (2) What are the benefits of AI being used for risk applications; and what are the benefits of applying risk principles and assumptions for AI-based applications, (3) What are the synergies between AI and risk science applications, and (4) What are the characteristics of effective use of fundamental risk principles and assumptions for AI-based applications? This study develops and disseminates an online survey questionnaire that leverages expertise from risk and AI professionals to identify the most important characteristics related to AI and risk, then presents a framework for gauging how AI and disaster risk can be balanced. This study is the first to develop a classification system for applying risk principles for AI-based applications. This classification contributes to understanding of AI and risk by exploring how AI can be used to manage risk, how AI methods introduce new or additional risk, and whether fundamental risk principles and assumptions are sufficient for AI-based applications.

97 MATHEMATICS AND COMPUTING↗

Simulation-aided occupant-centric building design: A critical review of tools, methods, and applications

Occupants are active participants in their built environment, affecting its performance while simultaneously being affected by its design and indoor environmental conditions. With recent advances in computer modeling, simulation tools, and analysis techniques, topics such as human-building interactions and occupant behavior have gained significant interest in the literature given their premise of improving building design processes and operating strategies. In practice, the focus of occupant-centric literature has been mostly geared towards the latter (i.e., operation), leaving the implications on building design practices underexplored. This paper fills the gap by providing a critical review of existing studies applying computer-based modeling and simulation to guide occupant-centric building design. The reviewed papers are organized along four main themes, namely occupant-centric: (i) metrics of building performance, (ii) modeling and simulation approaches, (iii) design methods and applications, and (iv) supporting practices and mechanisms. Important barriers are identified for a more effective application of occupant-centric building design practices, including the limited consideration of metrics beyond energy efficiency (e.g., occupant well-being and space planning), the limited implementation and validation of the proposed methods, and the lack of integration of occupant behavior modeling in existing building performance simulation tools. Finally, future research directions are discussed, covering large-scale international data collection efforts to move from generic assumptions about occupant behavior to specific/localized knowledge, improved metrics of measuring building performance, and improved industry practices, such as building codes, to promote an occupant-in-the-loop approach to the building design process.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Analytical approximation of a distorted reflector surface defined by a discrete set of points

Reflector antennas on Earth orbiting spacecrafts generally cannot be described analytically. The reflector surface is subjected to a large temperature fluctuation and gradients, and is thus warped from its true geometrical shape. Aside from distortion by thermal stresses, reflector surfaces are often purposely shaped to minimize phase aberrations and scanning losses. To analyze distorted reflector antennas defined by discrete surface points, a numerical technique must be applied to compute an interpolatory surface passing through a grid of discrete points. In this paper, the distorted reflector surface points are approximated by two analytical components: an undistorted surface component and a surface error component. The undistorted surface component is a best fit paraboloid polynomial for the given set of points and the surface error component is a Fourier series expansion of the deviation of the actual surface points, from the best fit paraboloid. By applying the numerical technique to approximate the surface normals of the distorted reflector surface, the induced surface current can be obtained using physical optics technique. These surface currents are integrated to find the far field radiation pattern.

Acosta, Roberto J.↗

Self contamination effects in the TAUVEX UV Telescope: Ground testing and computer simulation

The contamination effects due to outgassing from construction materials of the TAUVEX (Tel Aviv University UV Telescope) were evaluated using a combination of ground testing and computer simulations. Tests were performed from the material level of the system level including: (1) High sensitivity CVCM(10(exp -3 percent) measurements of critical materials. (2) Optical degradation measurements of samples specially contaminated by outgassing products at different contamination levels. (3) FTIR studies of chemical composition of outgassed products on above samples. (4) High resolution AFM studies of surface morphology of contaminated surfaces. The expected degradation of TAUVEX performance in mission was evaluated applying a computer simulation code using input parameters determined experimentally in the above tests. The results have served as guidelines for the proper selection of materials, cleanliness requirements, determination of the thermal conditions of the system and bakeout processes.

Lifshitz, Y.↗

Readout rebalancing for near-term quantum computers

Readout errors are a significant source of noise for near-term intermediate-scale quantum computers. Mismeasuring a qubit as |1) when it should be |0) occurs much less often than mismeasuring a qubit as |0) when it should have been |1). We make the simple observation that one can improve the readout fidelity of quantum computers by applying targeted X gates prior to performing a measurement. These X gates are placed so that the expected number of qubits in the |1) state is minimized. Classical postprocessing can undo the effect of the X gates so that the expectation value of any observable is unchanged. This protocol is designed to have no effect when readout errors are symmetric. We show that the statistical uncertainty following readout error corrections is smaller when using readout rebalancing. The statistical advantage is circuit and computer dependent and is demonstrated for the W state, a Grover search, and for a Gaussian state. The benefit in statistical precision is most pronounced (and nearly a factor of 2 in some cases) when states with many qubits in the excited state have high probability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Advancing quantum simulations of the nuclear shell model with Gray-code–based resource-efficient protocols

Background: Some of the computational limitations in solving the nuclear many-body problem could be overcome by utilizing quantum computers. The nuclear shell-model calculations providing deeper insights into the properties of atomic nuclei are one such case with high demand for resources, as the size of the Hilbert space grows exponentially with the number of particles involved. Quantum algorithms are being developed to overcome these challenges and advance such calculations. Purpose: To develop quantum circuits for the nuclear shell-model, leveraging the capabilities of noisy intermediate-scale quantum (NISQ) devices. Here, we aim to minimize resource requirements (specifically in terms of qubits and gates) and strive to reduce the impact of noise by employing relevant mitigation techniques. Methods: We achieve noise resilience by designing an optimized Ansatz for the variational quantum eigensolver (VQE) based on Givens rotations and incorporating qubit-ADAPT-VQE in combination with variational quantum deflation (VQD) to compute ground and excited states, incorporating the zero-noise extrapolation mitigation technique. Furthermore, the qubit requirements are significantly reduced by mapping the basis states to qubits using Gray-code encoding and generalizing transformations of fermionic operators to efficiently represent many-body states. Results: By employing the resource-efficient protocols, we achieve the ground and excited state energy levels of 38 Ar and 6 Li with better accuracy. These energy levels are presented for noiseless simulations, noisy conditions, and after applying noise mitigation techniques. Results are compared for Jordan-Wigner and Gray-code encoding using VQE, qubit-ADAPT-VQE, and VQD. Conclusions: Our work highlights the potential of resource-efficient protocols to leverage the full potential of NISQ devices in scaling the nuclear shell model calculations, offering a pathway toward more complex quantum simulations in nuclear physics. This approach establishes a framework for studying other nuclear systems with improved quantum resource efficiency, marking a significant advancement in applying quantum computing to realistic nuclear physics applications.

Physics - Nuclear physics and radiation physics↗

Full-wave modeling of RF waves in fusion plasmas with finite element method: Progress in past decades and its future role

This paper reviews the progress in computing radio frequency (RF) wave fields in fusion plasmas, specifically focusing on simulations utilizing the finite element method (FEM) over the past few decades. Computing RF wave fields in fusion plasmas presents unique challenges due to large simulation domains, complex antenna structures, non-local dielectric properties, and wide ranges of spatial scales. It highlights key developments and outlines future directions, primarily addressing waves in the ion cyclotron (IC) to lower hybrid (LH) frequency range. We begin with briefly revisiting earlier developments before the widespread availability of modern computer-aided engineering (CAE) software based on FEM. This historical perspective illuminates early progress and the physics difficulties that motivated ongoing work within the community. Modern wave simulations for RF antennas based on FEM are characterized by the use of detailed 3D antenna model geometry generated from engineering CAD software and localized wave dielectric model. Significant advancements have also been made in improving physics models to include phenomena such as RF sheath rectification and wave scattering. FEM-based RF simulations have also been applied to compute wave propagation in the core region, where the inclusion of non-local dielectric response is crucial. This is a challenging goal, and several promising approaches have been proposed in this area. Additionally, RF simulation development initiatives based on open-source libraries have gained popularity, demonstrating scalability and flexibility in extending physics models. This paper will discuss the advantages and disadvantages of using such a publicly available FEM library.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-dimensional computation of mixing of transverse injector in a ducted supersonic airstream

Numerical solutions of the three-dimensional mass-averaged elliptic Navier-Stokes equations, including species transport, are obtained for nonreacting, turbulent, mixing flow fields for the case of transverse sonic injection of a secondary gas into a supersonic airstream through a circular orifice injector. Results are presented for flow through a constant area duct and through a duct with a rearward-facing step upstream of the injector. The equations are numerically integrated using MacCormack's explicit method and turbulence is included using the Baldwin-Lomax algebraic eddy viscosity model. In the species transport and energy equations, diffusion coefficients based on Fick's law and an assumption of unit Lewis number are applied. The computations were performed on a CDC-VPS-32 (extended version of Cyber-205) using a grid consisting of approximately 200,000 points. The computed results are compared with experimentally observed penetration and spreading boundaries for an injected gas at two dynamic pressure ratios. Three-dimensional flow field structures are dipicted in terms of static pressure, mass fractions of species and velocity vectors.

Uenishi, K.↗

Quantum sensing effect of electron tunneling in DQD/analyte complex

We investigate electron tunneling between quantum dots and molecules to propose a quantum sensor. This sensor consists of double quantum dots (DQD) with energy levels specifically tailored to mirror those of the target analyte. By analyzing the spectral distribution of electron localizations in the DQD system, we can delineate the analyte’s spectrum and deduce its composition by comparing it with a reference sample. To understand electron tunneling dynamics within the DQD/analyte complex, we performed three-dimensional computational modeling applying the effective potential approach to the InAs/GaAs heterostructure. In this modeling, we mimicked the analyte spectrum by utilizing a quantum well characterized by a quasi-discrete spectrum. Our calculations reveal the inherent potential of utilizing this method as a highly sensitive and selective sensor.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A Bilevel Approach for Identifying the Worst Contingencies for Nonconvex Alternating Current Power Systems

We address the bilevel optimization problem of identifying the most critical attacks to an alternating current (AC) power flow network. The upper-level binary maximization problem consists of choosing an attack that is treated as a parameter in the lower-level defender minimization problem. Instances of the lower-level global minimization problem by themselves are NP-hard due to the nonconvex AC power flow constraints, and bilevel solution approaches commonly apply a convex relaxation or approximation to allow for tractable bilevel reformulations at the cost of underestimating some power system vulnerabilities. Our main contribution is to provide an alternative branch-and-bound algorithm whose upper bounding mechanism (in a maximization context) is based on a reformulation that avoids relaxation of the AC power flow constraints in the lower-level defender problem. Lower bounding is provided with semidefinite programming (SDP) relaxed solutions to the lower-level problem. We establish finite termination with guarantees of either a globally optimal solution to the original bilevel problem, or a globally optimal solution to the SDP-relaxed bilevel problem which is included in a vetted list of upper-level attack solutions, at least one of which is a globally optimal solution to the bilevel problem. We demonstrate through computational experiments applied to IEEE case instances both the relevance of our contribution, and the effectiveness of our contributed algorithm for identifying power system vulnerabilities without resorting to convex relaxations of the lower-level problem. We conclude with a discussion of future extensions and improvements.

97 MATHEMATICS AND COMPUTING↗

Coupling of the FUN3D Unstructured Flow Solver and the LASTRAC Stability Code to Model Transition

We develop an iterative automated method to predict transition locations in boundary-layer flows by using the FUN3D solver to perform flow simulations and the LASTRAC code for linear stability computations. The coupling of FUN3D and LASTRAC allows for a robust physics-based approach to model boundary-layer transition by analyzing the growth of different instability waves and then using that information to iteratively update the resulting transition location. There is no user involvement during the iterative computations. We apply this automated method to subsonic flow over a flat plate with a sharp leading edge. The final solution has regions of laminar and turbulent flow with a transition onset location that agrees with experiments and stability-based correlations. This iterative automated method is also applied to an NLF(1)-0416 airfoil for conditions with and without a separation bubble. Along with predicting transition locations, we compare the streamwise distributions of surface-pressure and skin-friction coefficients to a transport-equation-based model. We consider a 6:1 prolate spheroid at three angles of attack, namely, five, ten, and fifteen degrees, where mixed-mode transition occurs due to both Tollmien-Schlichting and crossflow instabilities. The skin-friction contours and transition fronts at every angle of attack from our iterative automated method show good agreement with past experimental and computational results in the literature for the 6:1 prolate spheroid.

Transition↗