Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

Normal-mode-based analysis of electron plasma waves with second-order Hermitian formalism

The classic problem of the dynamic evolution and Landau damping of linear Langmuir electron waves in a collisionless plasma with Maxwellian background is cast as a second-order, self-adjoint problem with a continuum spectrum of real and positive squared frequencies. The corresponding complete basis of singular normal modes is obtained, along with their orthogonality relation. This yields easily the general expression of the time-reversal-invariant solution for any initial-value problem. Examples are then given for specific initial conditions that illustrate different behaviors of the Landau-damped macroscopic moments of the perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using Machine Learning in Adversarial Environments.

Intrusion/anomaly detection systems are among the first lines of cyber defense. Commonly, they either use signatures or machine learning (ML) to identify threats, but fail to account for sophisticated attackers trying to circumvent them. We propose to embed machine learning within a game theoretic framework that performs adversarial modeling, develops methods for optimizing operational response based on ML, and integrates the resulting optimization codebase into the existing ML infrastructure developed by the Hybrid LDRD. Our approach addresses three key shortcomings of ML in adversarial settings: 1) resulting classifiers are typically deterministic and, therefore, easy to reverse engineer; 2) ML approaches only address the prediction problem, but do not prescribe how one should operationalize predictions, nor account for operational costs and constraints; and 3) ML approaches do not model attackers’ response and can be circumvented by sophisticated adversaries. The principal novelty of our approach is to construct an optimization framework that blends ML, operational considerations, and a model predicting attackers reaction, with the goal of computing optimal moving target defense. One important challenge is to construct a realistic model of an adversary that is tractable, yet realistic. We aim to advance the science of attacker modeling by considering game-theoretic methods, and by engaging experimental subjects with red teaming experience in trying to actively circumvent an intrusion detection system, and learning a predictive model of such circumvention activities. In addition, we will generate metrics to test that a particular model of an adversary is consistent with available data.

97 MATHEMATICS AND COMPUTING↗

Solid rocket motor aft field joint flow field analysis

An efficient Navier-Stokes analysis was successfully applied to simulate the complex flow field in the vicinity of a slot in a solid rocket motor with segment joints. The capability of the computer code to resolve the flow near solid surfaces without using a wall function assumption was demonstrated. In view of the complex nature of the flow field in the vicinity of the slot, this approach is considered essential. The results obtained from these calculations provide valuable design information, which would otherwise be extremely difficult to obtain. The results of the axisymmetric calculations indicate the presence of a region of reversed axial flow at the aft-edge of the slot and show the over-pressure in the slot to be only about 10 psi. The results of the asymmetric calculations indicate that a pressure asymmetry more than two diameters downstream of the slot has no noticeable effect on the flow field in the slot. They also indicate that the circumferential pressure differential caused in the slot due to failure of a 15 deg section of the castable inhibitor will be approximately 1 psi.

Sabnis, Jayant S.↗

Ablative Rocket Deflector Testing and Computational Modeling

A deflector risk mitigation program was recently conducted at the NASA Stennis Space Center. The primary objective was to develop a database that characterizes the behavior of industry-grade refractory materials subjected to rocket plume impingement conditions commonly experienced on static test stands. The program consisted of short and long duration engine tests where the supersonic exhaust flow from the engine impinged on an ablative panel. Quasi time-dependent erosion depths and patterns generated by the plume impingement were recorded for a variety of different ablative materials. The erosion behavior was found to be highly dependent on the material s composition and corresponding thermal properties. For example, in the case of the HP CAST 93Z ablative material, the erosion rate actually decreased under continued thermal heating conditions due to the formation of a low thermal conductivity "crystallization" layer. The "crystallization" layer produced near the surface of the material provided an effective insulation from the hot rocket exhaust plume. To gain further insight into the complex interaction of the plume with the ablative deflector, computational fluid dynamic modeling was performed in parallel to the ablative panel testing. The results from the current study demonstrated that locally high heating occurred due to shock reflections. These localized regions of shock-induced heat flux resulted in non-uniform erosion of the ablative panels. In turn, it was observed that the non-uniform erosion exacerbated the localized shock heating causing eventual plume separation and reversed flow for long duration tests under certain conditions. Overall, the flow simulations compared very well with the available experimental data obtained during this project.

Allgood, Daniel C.↗

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗

Enhancing robustness of multiparty quantum correlations using weak measurement

Multipartite quantum correlations are important resources for the development of quantum information and computation protocols. However, the resourcefulness of multipartite quantum correlations in practical settings is limited by its fragility under decoherence due to environmental interactions. Though there exist protocols to protect bipartite entanglement under decoherence, the implementation of such protocols for multipartite quantum correlations has not been sufficiently explored. Here, we study the effect of local amplitude damping channel on the generalized Greenberger–Horne–Zeilinger state, and use a protocol of optimal reversal quantum weak measurement to protect the multipartite quantum correlations. We observe that the weak measurement reversal protocol enhances the robustness of multipartite quantum correlations. Further it increases the critical damping value that corresponds to entanglement sudden death. To emphasize the efficacy of the technique in protection of multipartite quantum correlation, we investigate two proximately related quantum communication tasks, namely, quantum teleportation in a one sender, many receivers setting and multiparty quantum information splitting, through a local amplitude damping channel. We observe an increase in the average fidelity of both the quantum communication tasks under the weak measurement reversal protocol. The method may prove beneficial, for combating external interactions, in other quantum information tasks using multipartite resources. - Highlights: • Extension of weak measurement reversal scheme to protect multiparty quantum correlations. • Protection of multiparty quantum correlation under local amplitude damping noise. • Enhanced fidelity of quantum teleportation in one sender and many receivers setting. • Enhanced fidelity of quantum information splitting protocol.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Edge-Enhanced Oxygen Evolution Reactivity at Ultrathin, Au-Supported Fe 2 O 3 Electrocatalysts

Transition metal oxides have gained attention as promising oxygen evolution reaction (OER) electrocatalysts in alkaline electrolytes, but heterogeneities in typical catalyst samples often obscure key structure–property relationships that are essential for developing higher performance materials. Here, we have combined ultrahigh vacuum surface science techniques, electrochemical measurements, and density functional theory (DFT) to quantify structure-dependent OER activity in a series of well-defined electrocatalysts. We describe a direct correlation between the population of Fe edge-site atoms and the OER activity of ultrathin Fe 2 O 3 nanostructures (~0.5 nm apparent height) grown on Au(111) substrates. Hydroxylated Fe atoms residing at edge-sites along the catalyst/support interface were spectroscopically identified as key reaction centers, and these Fe edge-site atoms were estimated to produce OER turnover frequencies approximately 150 times higher than that of Fe atoms on the catalyst surface at an applied potential of 1.8 V vs the reversible hydrogen electrode. Impressively, ultrathin Fe 2 O 3 /Au nanostructures with a high density of catalytically active Fe edge-site atoms outperformed an ultrathin IrOx/Au catalyst at moderate overpotentials. DFT calculations revealed more favorable OER at edge sites along the Fe 2 O 3 /Au interface, with lower predicted overpotentials due to beneficial modification of intermediate binding. Our results demonstrate how a combination of surface science, electrochemistry, and computational modeling can be used to identify key structure–property relationships in a well-defined electrocatalytic system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energetic particle marginal stability profile for HL-2M integrated simulation based on neural network module

Abstract A critical gradient model is employed to develop a module of energetic particle (EP) marginal stability profiles in OMFIT integrated simulations for studying EP transport. Currently, each iteration of transport evolution is approximately 10 min in the integrated simulation, whereas, the EP marginal stability profile, which serves as an input in the integrated simulation could take much longer; the reason being a combination of the TGLFEP and EPtran codes is employed in our previous investigation. To reduce the simulation time, the critical gradient is predicted by a neural network instead of the TGLFEP code, and the EPtran code is revised with parallel computing, so that the running time of this module can be controlled to within 5 min. The predictions are in good agreement with previous approaches. The integrated simulation of HL-2M with Alfven eigenmodes transported by neutral beam EP profiles indicates that EP transport reduces the total pressure and current as expected, but could also under some conditions raise the safety factor in the core, which is favorable for reversed magnetic shear and high-performance plasmas.

Physics↗

The numerical calculation of laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.

Klineberg, J. M.↗

On laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for a number of separated flows and the accuracy of the method is verified.

Klineberg, J. M.↗

Thickening oscillation of a delta wing using Navier-Stokes and Navier-displacement equations

The problem of unsteady, supersonic, locally conical, vortical flow around a delta wing undergoing thickening oscillation is solved using the unsteady, thin-layer Navier-Stokes equations and the unsteady linearized Navier-displacement equations. The unsteady, thin-layer Navier-Stokes equations are solved using the implicit approximate-factorization finite-volume scheme to compute the conservative components of the flow vector field. With the conservative components known at any time step, the linearized, Navier-displacement equations are solved using the alternating, direction-implicit scheme to obtain the grid points displacements due to known displacement boundary conditions. A grid-displacements limiter, in the form of a low mesh Reynolds number, is used to limit grid-folding in regions of highly reversed flow.

Chuang, Hsin-Kung A.↗

In situ studies of reversible solid–gas reactions of ethylene responsive silver pyrazolates

Solid–gas reactions and in situ powder X-ray diffraction investigations of trinuclear silver complexes {[3,4,5-(CF 3 ) 3 Pz]Ag} 3 and {[4-Br-3,5-(CF 3 ) 2 Pz]Ag} 3 supported by highly fluorinated pyrazolates reveal that they undergo intricate ethylene-triggered structural transformations in the solid-state producing dinuclear silver–ethylene adducts. Despite the complexity, the chemistry is reversible producing precursor trimers with the loss of ethylene. Less reactive {[3,5-(CF 3 ) 2 Pz]Ag} 3 under ethylene pressure and low-temperature conditions stops at an unusual silver–ethylene complex in the trinuclear state, which could serve as a model for intermediates likely present in more common trimer–dimer reorganizations described above. Complete structural data of three novel silver–ethylene complexes are presented together with a thorough computational analysis of the mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward full simulations for a liquid metal blanket: part 2. Computations of MHD flows with volumetric heating for a PbLi blanket prototype at Ha ~10 4 and Gr ~10 12

On the pathway toward full simulations for a liquid metal (LM) blanket, this part 2 extends a previous study of purely magnetohydrodynamic (MHD) flows in a DCLL blanket in reference Chen et al (2020 Nucl. Fusion 60 076003) to more general conditions when the MHD flow is coupled with heat transfer. The simulated prototypic blanket module includes all components of a real LM blanket system, such as supply ducts, inlet and outlet manifolds, multiple poloidal ducts and a U-turn zone. Volumetric heating generated by fusion neutrons is added to simulate thermal effects in the flowing lead–lithium (PbLi) breeder. The MHD flow equations and the energy equation are solved with a DNS-type finite-volume code ‘MHD-UCAS’ on a very fine mesh of 470 × 10 6 cells. The applied magnetic field is 5 T (Hartmann number Ha ~ 10 4 ), the PbLi velocity in the poloidal ducts is 10 cm s –1 (Reynolds number Re ~10 5 ), whereas the maximum volumetric heating is 30 MW m –3 (Grashof number Gr ~ 10 12 ). Four cases have been simulated, including forced- and mixed-convection flows, and either an electrically conducting or insulating blanket structure. Various comparisons are made between the four computed cases and also against the purely MHD flows computed earlier in reference Chen et al (2020 Nucl. Fusion 60 076003) with regards to the (1) MHD pressure drop, (2) flow balancing, (3) temperature field, (4) flows in particular blanket components, and (5) 3D and turbulent flow effects. The strongest buoyancy effects were found in the poloidal ducts. In the electrically non-conducting blanket, the buoyancy forces lead to significant modifications of the flow structure, such as formation of reverse flows, whereas their effect on the MHD pressure drop is relatively small. In the electrically conducting blanket case, the buoyancy effects on the flow and MHD pressure drop are almost negligible.

Physics↗

Investigation of two-dimensional radio-frequency sheath properties using a microscale fluid model

Abstract In previous work (Kohno and Myra 2023 Comput. Phys. Commun. 291 108841), we developed a numerical scheme based on a two-dimensional microscale radio-frequency (RF) sheath model with periodically curved wall boundaries. Here, we expand the capability of this scheme through modification of the boundary conditions (BCs) on the conducting walls, which allows the ion flow to turn back to the plasma at locations on the walls where the electromagnetic force on the ions is reversed from its usual direction. Numerical simulations are carried out to investigate the dependences of the surface-integrated admittances on the wall bump height, ion magnetization, ion mobility, and the magnetic field angle, and to visualize the sheath structures in several cases. One of the main results is the ion cyclotron admittance resonance observed under the condition of low ion mobility (high normalized frequency). It is shown that the amplitude of the resonance peak depends on the wall bump height and the ion velocity is reversed on the sides of the bump in an RF cycle for the resonance cases. Furthermore, the differences in the admittances between the one- and two-dimensional microscale models are assessed for the purpose of understanding non-locality of the sheath near the wall surface for the parameters considered in this study. This information will be essential for improving the sheath BC for macroscale calculations in the future.

Kohno, H. (ORCID:0000000332690010)↗

Wake Velocimetry of a Sphere-Cone Model in a Mach 10 Air Freestream using FLEET

Velocity measurements in the wake behind a sphere-cone-shaped vehicle wind tunnel model in a Mach 10 hypersonic air flow using femtosecond laser electronic excitation tagging (FLEET) are reported in this work. The FLEET measurements were performed at 1 kHz using a femtosecond laser centered at 800 nm and an ultrahigh-speed intensified camera system. The results reported here are for a single freestream condition corresponding to approximately Re ∞ /L = 3.6∙10 6 m -1 . The performance of the velocimetry instrument was evaluated in the quiescent test section at conditions relevant to the wake measurements. For velocity measurements in the wake, the FLEET measurement volume was translated to several radial locations from the model centerline in the vertical direction and at a streamwise location corresponding to approximately 12 mm from the payload. Reverse flow with mean velocities ranging from -12 m/s to -48 m/s were observed in the recirculating region of the wake, whereas hypersonic velocities in the range of 1079 m/s to 1183 m/s were observed in the free shear layer. A velocity uncertainty methodology is also outlined and applied for the reported measurements. The velocity data reported in this work is expected to play a significant role in the validation of computational codes modeling the hard-to-predict separated hypersonic wake.

FLEET↗

Wake Velocimetry of a Sphere-Cone Model in a Mach 10 Air Freestream using FLEET

Velocity measurements in the wake behind a sphere-cone-shaped vehicle wind tunnel model in a Mach 10 hypersonic air flow using femtosecond laser electronic excitation tagging (FLEET) are reported in this work. The FLEET measurements were performed at 1 kHz using a femtosecond laser centered at 800 nm and an ultrahigh-speed intensified camera system. The results reported here are for a single freestream condition corresponding to approximately Re ∞ /L = 3.6∙10 6 m -1 . The performance of the velocimetry instrument was evaluated in the quiescent test section at conditions relevant to the wake measurements. For velocity measurements in the wake, the FLEET measurement volume was translated to several radial locations from the model centerline in the vertical direction and at a streamwise location corresponding to approximately 12 mm from the payload. Reverse flow with mean velocities ranging from -12 m/s to -48 m/s were observed in the recirculating region of the wake, whereas hypersonic velocities in the range of 1079 m/s to 1183 m/s were observed in the free shear layer. A velocity uncertainty methodology is also outlined and applied for the reported measurements. The velocity data reported in this work is expected to play a significant role in the validation of computational codes modeling the hard-to-predict separated hypersonic wake.

FLEET↗

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Correlating dynamic microstructure to observed color in electrophoretic displays via in situ small-angle x-ray scattering

Electrophoretic deposition (EPD) is an industrially relevant and scalable technique used to form particle deposits from colloidal suspensions. Highly concentrated particle suspensions generally prevent real-time in situ microscopy observations which limit the characterization of EPD films to ex situ, or postprocessed, laboratory techniques. For dynamic systems, such as tunable amorphous photonic crystals (APCs), only reversible deposits are formed during the EPD process. Since reversible deposits cannot be characterized with standard ex situ methods, the particle-particle and particle-field interactions that govern the displayed color and crystallinity of these systems are not well understood. In this paper, we present in situ small-angle x-ray scattering and UV-Vis techniques for measuring both the structural and optical response of an APC under applied electric fields. Furthermore, we also develop a computational model based on colloidal interactions to explain the observed change in the interparticle spacing of APCs due to the applied electric field which correlates to displayed color. Ultimately, this work provides a new in situ characterization method that could be expanded for other dynamic, tunable colloidal systems.

36 MATERIALS SCIENCE↗