Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Based Procedures”

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 199 records · Page 11

A mass conserving level set method for detailed numerical simulation of liquid atomization

An improved mass conserving level set method for detailed numerical simulations of liquid atomization is developed to address the issue of mass loss in the existing level set method. This method introduces a mass remedy procedure based on the local curvature at the interface, and in principle, can ensure the absolute mass conservation of the liquid phase in the computational domain. Three benchmark cases, including Zalesak's disk, a drop deforming in a vortex field, and the binary drop head-on collision, are simulated to validate the present method, and the excellent agreement with exact solutions or experimental results is achieved. It is shown that the present method is able to capture the complex interface with second-order accuracy and negligible additional computational cost. The present method is then applied to study more complex flows, such as a drop impacting on a liquid film and the swirling liquid sheet atomization, which again, demonstrates the advantages of mass conservation and the capability to represent the interface accurately.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adiabatic elimination-based coupling control in densely packed subwavelength waveguides

The ability to control light propagation in photonic integrated circuits is at the foundation of modern light-based communication. However, the inherent crosstalk in densely packed waveguides and the lack of robust control of the coupling are a major roadblock toward ultrahigh density photonic integrated circuits. As a result, the diffraction limit is often considered as the lower bound for ultra-dense silicon photonics circuits. Here we experimentally demonstrate an active control of the coupling between two closely packed waveguides via the interaction with a decoupled waveguide. This control scheme is analogous to the adiabatic elimination, a well-known procedure in atomic physics. This approach offers an attractive solution for ultra-dense integrated nanophotonics for light-based communications and integrated quantum computing

36 MATERIALS SCIENCE↗

MAVTgsa: An R Package for Gene Set (Enrichment) Analysis

Gene set analysis methods aim to determine whether an a priori defined set of genes shows statistically significant difference in expression on either categorical or continuous outcomes. Although many methods for gene set analysis have been proposed, a systematic analysis tool for identification of different types of gene set significance modules has not been developed previously. This work presents an R package, called MAVTgsa, which includes three different methods for integrated gene set enrichment analysis. (1) The one-sided OLS (ordinary least squares) test detects coordinated changes of genes in gene set in one direction, either up- or downregulation. (2) The two-sided MANOVA (multivariate analysis variance) detects changes both up- and downregulation for studying two or more experimental conditions. (3) A random forests-based procedure is to identify gene sets that can accurately predict samples from different experimental conditions or are associated with the continuous phenotypes. MAVTgsa computes the P values and FDR (false discovery rate) q -value for all gene sets in the study. Furthermore, MAVTgsa provides several visualization outputs to support and interpret the enrichment results. This package is available online.

Chien, Chih-Yi↗

An Approach to Dynamic Human Reliability Analysis Using the EMRALD Dynamic Risk Assessment Tool

There is a need for research into dynamic human reliability analysis (HRA) (i.e., simulation- or computation-based HRA), as many researchers have emphasized the importance of applying dynamic approaches to probabilistic safety assessments (PSAs). This study proposes the Procedure-based Investigation Method of EMRALD Risk Assessment – Human Reliability Analysis (PRIMERA-HRA), a dynamic HRA approach based on the Event Modeling Risk Assessment Using Linked Diagram (EMRALD) software developed by Idaho National Laboratory (INL). This study suggests how this software can be used to model human actions and evaluate error probabilities. The applicability of this approach is also be investigated in light of an extended loss of AC power (ELAP) scenario. This paper then discusses the major insights derived by comparing this dynamic HRA approach against the static one.

99 GENERAL AND MISCELLANEOUS↗

Approximation of refrigerant thermophysical properties using neural networks to speed up transient thermofluid simulations

Accurate and efficient evaluations of refrigerant thermophysical properties and their partial derivatives are essential for transient simulations of thermofluid systems, where several computations need to be executed at each integration time step. Since the utilization of an Equation of State for retrieving properties based on a pair of independent inputs typically involves numerical iterations in solution procedures, when the input variables differ from the refrigerant state variables employed in dynamic models, a variety of approaches including lookup table interpolation and curve fitting have been developed to explicitly approximate these properties based on the state variables, and consequently eliminate internal iterations. This paper presents an alternative method that exploits derivative-informed neural networks to model refrigerant properties explicitly from inputs of pressure and enthalpy, while ensuring consistent partial derivatives generated by differentiating the neural networks. Computational speed and accuracy of the proposed approach are demonstrated via transient simulations of a discretized heat exchanger model in Modelica, and comparisons against other property evaluation routines. Simulation results indicate that the proposed approach can realize a significant speedup with negligible discrepancies in predicted transients. The method is implemented in an open-source Modelica library.

Ma, Jiacheng↗

Field Programmable Gate Array-Based Reactor Protection Systems and Potential for Inclusion of Secure Elements to Improve Cybersecurity

For acceptable implementations of technologies like wireless communications, remote monitoring, etc., strong mitigations must be developed and evaluated to ensure that new attack pathways do not increase risk for Advanced Reactors. Secure Elements can be adopted and adapted for this purpose based on tamper resistance and cryptographic abilities, but research must be done to properly integrate into critical components such as FPGA-based Important to Safety systems in conjunction with current and future regulations on cyber security features in Advanced Reactors. Typically, the integration of a Secure Element happens during the POST and UEFI boot of a computing platform, performed by the Operating System, which is not possible with FPGAs because they do not include these firmware components. Work must be done to identify a reliable and secure method for integration in FPGA-based systems which lack Operating Systems and therefore complex boot procedures, system calls, etc.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ANALYTIC H i-to-H{sub 2} PHOTODISSOCIATION TRANSITION PROFILES

We present a simple analytic procedure for generating atomic (H i) to molecular (H{sub 2}) density profiles for optically thick hydrogen gas clouds illuminated by far-ultraviolet radiation fields. Our procedure is based on the analytic theory for the structure of one-dimensional H i/H{sub 2} photon-dominated regions, presented by Sternberg et al. Depth-dependent atomic and molecular density fractions may be computed for arbitrary gas density, far-ultraviolet field intensity, and the metallicity-dependent H{sub 2} formation rate coefficient, and dust absorption cross section in the Lyman–Werner photodissociation band. We use our procedure to generate a set of H I-to-H{sub 2} transition profiles for a wide range of conditions, from the weak- to strong-field limits, and from super-solar down to low metallicities. We show that if presented as functions of dust optical depth, the H I and H{sub 2} density profiles depend primarily on the Sternberg “αG parameter” (dimensionless) that determines the dust optical depth associated with the total photodissociated H I column. We derive a universal analytic formula for the H I-to-H{sub 2} transition points as a function of just αG. Our formula will be useful for interpreting emission-line observations of H i/H{sub 2} interfaces, for estimating star formation thresholds, and for sub-grid components in hydrodynamics simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Learning to Optimize Variational Quantum Circuits to Solve Combinatorial Problems

Quantum computing is a computational paradigm with the potential to outperform classical methods for a variety of problems. Proposed recently, the Quantum Approximate Optimization Algorithm (QAOA) is considered as one of the leading candidates for demonstrating quantum advantage in the near term. QAOA is a variational hybrid quantum-classical algorithm for approximately solving combinatorial optimization problems. The quality of the solution obtained by QAOA for a given problem instance depends on the performance of the classical optimizer used to optimize the variational parameters. In this paper, we formulate the problem of finding optimal QAOA parameters as a learning task in which the knowledge gained from solving training instances can be leveraged to find high-quality solutions for unseen test instances. To this end, we develop two machine-learning-based approaches. Our first approach adopts a reinforcement learning (RL) framework to learn a policy network to optimize QAOA circuits. Our second approach adopts a kernel density estimation (KDE) technique to learn a generative model of optimal QAOA parameters. In both approaches, the training procedure is performed on small-sized problem instances that can be simulated on a classical computer; yet the learned RL policy and the generative model can be used to efficiently solve larger problems. Furthermore, extensive simulations using the IBM Qiskit Aer quantum circuit simulator demonstrate that our proposed RL- and KDE-based approaches reduce the optimality gap by factors up to 30.15 when compared with other commonly used off-the-shelf optimizers.

97 MATHEMATICS AND COMPUTING↗

Fast, clash-free RNA conformational morphing using molecular junctions

Non-coding ribonucleic acids (ncRNA) are functional RNA molecules that are not translated into protein. They are extremely dynamic, adopting diverse conformational substates, which enables them to modulate their interaction with a large number of other molecules. The flexibility of ncRNA provides a challenge for probing their complex 3D conformational landscape, both experimentally and computationally. As a result, despite their conformational diversity, ncRNAs mostly preserve their secondary structure throughout the dynamic ensemble. Here we present a kinematics-based procedure to morph an RNA molecule between conformational substates, while avoiding inter-atomic clashes. We represent an RNA as a kinematic linkage, with fixed groups of atoms as rigid bodies and rotatable bonds as degrees of freedom. Our procedure maintains RNA secondary structure by treating hydrogen bonds between base pairs as constraints. The constraints define a lower-dimensional, secondary-structure constraint manifold in conformation space, where motions are largely governed by molecular junctions of unpaired nucleotides. On a large benchmark set, we show that our morphing procedure compares favorably to peer algorithms, and can approach goal conformations to within a low all-atom RMSD by directing fewer than 1% of its atoms. Furthermore, our results suggest that molecular junctions can modulate 3D structural rearrangements, while secondary structure elements guide large parts of the molecule along the transition to the correct final conformation.

59 BASIC BIOLOGICAL SCIENCES↗

Wide Bandgap Generation (WBGen): Developing the Future Wide Bandgap Power Electronics Engineering Workforce

This Final Technical Report (FTR) summarizes the work conducted at the Center for Power Electronics Systems (CPES) under the Wide-Bandgap Generation (WBGen) fellowship and traineeship program established by the Advanced Manufacturing and Technologies Office (AMMTO) of the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy (EERE) at Virginia Tech, which had as main objective to train the next generation of U.S. citizen power engineers with wide-bandgap (WBG) power semiconductor expertise, with the intent to aid in fulfilling the future workforce needs in this field. The latter was deemed of strategic importance given the fast-paced growth observed—and predicted—in the demand of this technical expertise, whose practitioners have become the enablers and executioners of the electrification transformation process that not just the U.S., but the whole world, is currently undergoing as it seeks for more effective and efficient ways to use energy. As such, the WBGen program set forth to achieve its educational goals, which in addition sought to broaden the range of WBG-based power electronics by conducting research and development on high-efficiency grid apparatus and high-efficiency electrical power systems, and to also enhance the power engineering curriculum by formalizing WBG-oriented design procedures replacing existent yet now obsolete design procedures developed for Silicon (Si) based power electronics. This effort led CPES to spearhead the development of a new major within The Bradley Electrical and Computer Engineering (ECE) Department at Virginia Tech, namely Electronic Power and Energy Systems (EPES), which coalesced power electronics and power systems courses to provide undergraduate students with a strong formation in the power engineering field, while creating a pipeline of graduate students that could pursue the WBG-based curriculum and conduct research at CPES. In all, in what is considered a true success, eight of the twenty WBGen fellows that graduated program were recruited from the ECE department undergraduate cohort. The traineeship emphasized as well, from its beginning, the partnership with industry and national laboratories, which took advantage of the successful industry consortium at CPES that has historically been formed by 80–90 power and energy companies working in close collaboration with the center. This gave WBGen fellows the accessibility and possibility to conduct internships at partner facilities during the summer months, focused solely on the evaluation, testing, and adoption of WBG devices, which were many times tightly related to their respective research work and plans. In addition, WBGen fellows conducted their main research work within the confines of research programs at CPES conducted with these industry partners, providing them with a unique opportunity to develop not just their technical expertise—while advancing their knowledge, but to also learn and practice a slew of skills needed for their professional growth. As such, the fellows tackled a variety of WBG-related research topics, from device reliability and capability aspects as well as packaging and integration, encompassing the use of advanced materials and new structures, current sharing challenges, and insulation systems, to advanced gate-drivers with integrated sensors and protection mechanisms and active current- and voltage-based control, to optimized layouts seeking to maximize the switching and power processing performance of these devices, to power processing solutions adopting Gallium-Nitride (GaN) and Silicon-Carbide (SiC) power semiconductors for a variety of applications; including radiation-hardened converters for space dc distribution systems, direct three-phase ac-to-ac power converters for aerospace systems, dc-ac inverters for automotive traction drives, high-frequency isolated dc-dc battery chargers also for heavy transportation systems, and medium-voltage dc-dc and dc-ac converters for distribution systems and future power grids. The WBGen program ultimately graduated a total of 20 power engineers, all experts on WBG-based power electronics, awarding 18 M.S. and 2 PhD degrees in the process. These fellows, all U.S. citizens, allowed CPES to increase the number of citizens students to 33 % at the peak of the program, as the traineeship made possible the recruitment of talent with more attractive graduate research assistantship (GRA) contracts. Unfortunately, the present U.S. citizen enrollment at CPES has declined back to historic levels—approximately 10 %, as the regular GRA rates are not competitive enough when compared with entry-level industry jobs. In all, WBGen fellows published a total of 7 peer-reviewed journal articles, 44 papers at international technical conferences, made 42 presentations at international technical conferences, and filed 4 invention disclosures and patent applications, which have since then been granted by the U.S. Patent and Trademark Office (USPTO). Their contribution to CPES, Virginia Tech, the United States, and the world, has been significant, and continues to yield results thanks to the exemplary career that the fellows have initiated at many of the partners of the program, which include Wolfspeed, Raytheon Technologies, Infineon, Lockheed Martin, Dominion Energy, Northrop Grumman, Rivian, Aerospace Corporation, Sandia National Laboratory, National Renewable Energy Laboratory, John Hopkins University, and Virginia Tech.

14 SOLAR ENERGY↗

Achieving stress-constrained topological design via length scale control

A new suite of computational procedures for stress-constrained continuum topology optimization is presented. In contrast to common approaches for imposing stress constraints, herein it is proposed to limit the maximum stress by controlling the length scale of the optimized design. Several procedures are formulated based on the treatment of the filter radius as a design variable. This enables to automatically manipulate the minimum length scale such that stresses are constrained to the allowable value, while the optimization is driven to minimizing compliance under a volume constraint – without any direct constraints on stresses. Numerical experiments are presented that incorporate the following : 1) Global control over the filter radius that leads to a uniform minimum length scale throughout the design; 2) Spatial variation of the filter radius that leads to local manipulation of the minimum length according to stress concentrations; and 3) Combinations of the two above. The optimized designs provide high-quality trade-offs between compliance, stress and volume. From a computational perspective, the proposed procedures are efficient and simple to implement: essentially, stress-constrained topology optimization is posed as a minimum compliance problem with additional treatment of the length scale.

42 ENGINEERING↗

Learning compact physics‐aware delayed photocurrent models using dynamic mode decomposition

Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.

Hanson, Joshua↗

A bidirectional coupling procedure applied to multiscale respiratory modeling

In this study, we present a novel multiscale computational framework for efficiently linking multiple lower-dimensional models describing the distal lung mechanics to imaging-based 3D computational fluid dynamics (CFDs) models of the upper pulmonary airways in order to incorporate physiologically appropriate outlet boundary conditions. The framework is an extension of the modified Newton’s method with nonlinear Krylov accelerator developed by Carlson and Miller [1], Miller [2] and Scott and Fenves [3]. Our extensions include the retention of subspace information over multiple timesteps, and a special correction at the end of a timestep that allows for corrections to be accepted with verified low residual with as little as a single residual evaluation per timestep on average. In the case of a single residual evaluation per timestep, the method has zero additional computational cost compared to uncoupled or unidirectionally coupled simulations. We expect these enhancements to be generally applicable to other multiscale coupling applications where timestepping occurs. In addition we have developed a “pressure-drop” residual which allows for stable coupling of flows between a 3D incompressible CFD application and another (lower-dimensional) fluid system. We expect this residual to also be useful for coupling non-respiratory incompressible fluid applications, such as multiscale simulations involving blood flow. The lower-dimensional models that are considered in this study are sets of simple ordinary differential equations (ODEs) representing the compliant mechanics of symmetric human pulmonary airway trees. To validate the method, we compare the predictions of hybrid CFD–ODE models against an ODE-only model of pulmonary airflow in an idealized geometry. Subsequently, we couple multiple sets of ODEs describing the distal lung to an imaging-based human lung geometry. Boundary conditions in these models consist of atmospheric pressure at the mouth and intrapleural pressure applied to the multiple sets of ODEs. In both the simplified geometry and in the imaging-based geometry, the performance of the method was comparable to that of monolithic schemes, in most cases requiring only a single CFD evaluation per time step. Thus, this new accelerator allows us to begin combining pulmonary CFD models with lower-dimensional models of pulmonary mechanics with little computational overhead. Moreover, because the CFD and lower-dimensional models are totally separate, this framework affords great flexibility in terms of the type and breadth of the adopted lower-dimensional model, allowing the biomedical researcher to appropriately focus on model design. Research funded by the National Heart and Blood Institute Award 1RO1HL073598.

60 APPLIED LIFE SCIENCES↗

Extreme ultraviolet microscope characterization using photomask surface roughness

We demonstrate a method for characterizing the field-dependent aberrations of a full-field synchrotron-based extreme ultraviolet microscope. The statistical uniformity of the inherent, atomic-scale roughness of readily-available photomask blanks enables a self-calibrating computational procedure using images acquired under standard operation. We characterize the aberrations across a 30-um field-of-view, demonstrating a minimum aberration magnitude of smaller than λ/21 rms averaged over the center 5-um area, with a measurement accuracy better than λ/180 rms. The measured field variation of aberrations is consistent with system geometry and agrees with prior characterizations of the same system. In certain cases, it may be possible to additionally recover the illumination wavefront from the same images. Our method is general and is easily applied to coherent imaging systems with steerable illumination without requiring invasive hardware or custom test objects; hence, it provides substantial benefits when characterizing microscopes and high-resolution imaging systems in situ.

47 OTHER INSTRUMENTATION↗

Fluid dynamic modeling of nano-thermite reactions

This paper presents a direct numerical method based on gas dynamic equations to predict pressure evolution during the discharge of nanoenergetic materials. The direct numerical method provides for modeling reflections of the shock waves from the reactor walls that generates pressure-time fluctuations. The results of gas pressure prediction are consistent with the experimental evidence and estimates based on the self-similar solution. Artificial viscosity provides sufficient smoothing of shock wave discontinuity for the numerical procedure. The direct numerical method is more computationally demanding and flexible than self-similar solution, in particular it allows study of a shock wave in its early stage of reaction and allows the investigation of “slower” reactions, which may produce weaker shock waves. Moreover, numerical results indicate that peak pressure is not very sensitive to initial density and reaction time, providing that all the material reacts well before the shock wave arrives at the end of the reactor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluation of phase shifts for nonrelativistic elastic scattering using quantum computers

Simulations of scattering processes are essential in understanding the physics of our universe. Computing relevant scattering quantities from ab initio methods is extremely difficult on classical devices because of the substantial computational resources needed. Here, this work reports the development of an algorithm that makes it possible to obtain phase shifts for generic nonrelativistic elastic scattering processes on a quantum computer. This algorithm is based on extracting phase shifts from the direct implementation of the real-time evolution. The algorithm is improved by a variational procedure making it more accurate and resistant to the quantum noise. The reliability of the algorithm is first demonstrated by means of classical numerical simulations for different potentials and later tested on existing quantum hardware, specifically on IBM quantum processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Step-Down Test Procedure for Wavelet Shrinkage Using Bootstrapping

Wavelet thresholding (or shrinkage) attempts to remove the noises existing in the signals while preserving inherent pattern characteristics in the reconstruction of true signals. For data-denoising purpose, we present a new wavelet thresholding procedure which employs the step-down testing idea of identifying active contrasts in unreplicated fractional factorial experiments. The proposed method employs bootstrapping methods to a step-down test for thresholding wavelet coefficients. By introducing the concept of a false discovery error rate in testing wavelet coefficients, we shrink the wavelet coefficients with p -values higher than the error rate. The error rate controls the expected proportion of wrongly accepted coefficients among chosen wavelet coefficients. Bootstrap samples are used to approximate the p -value for computational efficiency. We also present some guidelines for selecting the values of hyper-parameters which affect the performance in the step-down thresholding procedure. Based on some common testing signals and an air-conditioner sounds example, the comparison of our proposed procedure with other thresholding methods in the literature is performed. The analytical results show that the proposed procedure has a potential in data-denoising and data-reduction in a variety of signal reconstruction applications.

42 ENGINEERING↗

An Approach to Dynamic Human Reliability Analysis using EMRALD Dynamic Risk Assessment Tool

Research for dynamic human reliability analysis (HRA) (a.k.a., simulation-based or computation-based HRA) has been required, as many researchers have emphasized the importance of dynamic approaches to probabilistic safety assessment (PSA). This study proposes the PRocedure-based Investigation Method of EMRALD Risk Assessment – Human Reliability Analysis (PRIMERA-HRA) method, which is an approach to dynamic HRA using the Event Modeling Risk Assessment using Linked Diagram (EMRALD) software developed by Idaho National Laboratory (INL). In this study, how to model human actions using the software and evaluate the error probabilities is suggested. The applicability of this approach will be also investigated through an extended loss of AC power (ELAP) scenario. Then, major insights from this dynamic HRA approach against the static one will be discussed in the paper.

99 GENERAL AND MISCELLANEOUS↗