Search NASA⌕ Search

SEARCH · Search NASA

Results for “classical simulation”

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 469 records · Page 26

Efficient Multi-Dimensional Simulation of Quantum Confinement Effects in Advanced MOS Devices

We investigate the density-gradient (DG) transport model for efficient multi-dimensional simulation of quantum confinement effects in advanced MOS devices. The formulation of the DG model is described as a quantum correction to the classical drift-diffusion model. Quantum confinement effects are shown to be significant in sub-100nm MOSFETs. In thin-oxide MOS capacitors, quantum effects may reduce gate capacitance by 25% or more. As a result, the inclusion or quantum effects in simulations dramatically improves the match between C-V simulations and measurements for oxide thickness down to 2 nm. Significant quantum corrections also occur in the I-V characteristics of short-channel (30 to 100 nm) n-MOSFETs, with current drive reduced by up to 70%. This effect is shown to result from reduced inversion charge due to quantum confinement of electrons in the channel. Also, subthreshold slope is degraded by 15 to 20 mV/decade with the inclusion of quantum effects via the density-gradient model, and short channel effects (in particular, drain-induced barrier lowering) are noticeably increased.

Biegel, Bryan A.↗

Transonic Flutter Suppression Control Law Design Using Classical and Optimal Techniques with Wind-Tunnel Results

The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using (1) classical, (2) linear quadratic Gaussian (LQG), and (3) minimax techniques are described. A unified general formulation and solution for the LQG and minimax approaches, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf. The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Mukhopadhyay, Vivek↗

Transonic Flutter Suppression Control Law Design Using Classical and Optimal Techniques with Wind-Tunnel Results

The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using (1) classical, (2) linear quadratic Gaussian (LQG), and (3) minimax techniques are described. A unified general formulation and solution for the LQG and minimax approaches, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf. The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Mukhopadhyay, Vivek↗

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗

Evidence of short chains in liquid sulfur

High energy x-ray pair distribution function measurements show the average coordination number of the first shell in liquid sulfur is 1.86 ± 0.04 across the λ-transition, not precisely 2.0 as widely accepted. This indicates that upon melting, liquid sulfur does not comprise solely of S 8 rings but also possesses a significant number of short chains. Intensities of the pre-peak and first diffraction peak of the x-ray structure factor and third peak height of the pair distribution function all show deviations at the λ-transition temperature T λ , associated with the break-up of S 8 rings and the start of oligomer polymerization. A significant number of non-bonded or loosely bonded “interstitial atoms,” with an average coordination number of 0.20 ± 0.005, are also observed in the so-called “forbidden zone” between the first and second shells upon melting. The number of interstitial atoms is found to decrease to a minimum at the λ-transition, but the majority persist into the high temperature polymerized liquid. Furthermore, the existence of short chains and nearby interstitial atoms represent the two main factors required to initiate the S 8 -ring to chain transition, as proposed by recent molecular dynamics simulations.

Chemical bonding↗

Validity test for linear error analysis

To determine whether estimation process simulated by linear error analysis will converge, criterion has been developed based on extension of classical observability. Particular application of technique is with groups of batched navigation data where statistics of estimation errors are derived with classical minimum-variance methods.

Diamant, L. S.↗

Nonlinear structural and life analyses of a combustor liner

Three-dimensional, nonlinear, finite element structural analyses were performed for a simulated aircraft combustor liner specimen in order to assess the capability of nonlinear analyses using classical inelastic material models to represent the thermoplastic-creep response of the component. In addition, the computed stress-strain history at the critical location was input into life prediction methods in order to evaluate the ability of these procedures to predict crack initiation life. It is concluded that: (1) elastic analysis is adequate for obtaining strain range and critical location; (2) inelastic analyses did not accurately represent cyclic behavior of materials; and (3) none of the crack initiation life prediction methods were satisfactory.

Kaufman, A.↗

Towards the Identification of the Keeper Erosion Cause(s): Numerical Simulations of the Plasma and Neutral Gas Using the Global Cathode Model OrCa2D-II

Numerical simulations with the time-dependent Orificed Cathode (OrCa2D-II) computer code show that classical enhancements of the plasma resistivity can not account for the elevated electron temperatures and steep plasma potential gradients measured in the plume of a 25-27.5 A discharge hollow cathode. The cathode, which employs a 0.11-in diameter orifice, was operated at 5.5 sccm without an applied magnetic field using two different anode geometries. It is found that anomalous resistivity based on electron-driven instabilities improves the comparison between theory and experiment. It is also estimated that other effects such as the Hall-effect from the self-induced magnetic field, not presently included in OrCa2D-II, may contribute to the constriction of the current density streamlines thus explaining the higher plasma densities observed along the centerline.

plasma modeling↗

Solar Contamination in Extreme-precision Radial-velocity Measurements: Deleterious Effects and Prospects for Mitigation

Solar contamination, due to moonlight and atmospheric scattering of sunlight, can cause systematic errors in stellar radial velocity (RV) measurements that significantly detract from the ~10 cm/s sensitivity required for the detection and characterization of terrestrial exoplanets in or near habitable zones of Sun-like stars. The addition of low-level spectral contamination at variable effective velocity offsets introduces systematic noise when measuring velocities using classical mask-based or template-based cross-correlation techniques. Here we present simulations estimating the range of RV measurement error induced by uncorrected scattered sunlight contamination. We explore potential correction techniques, using both simultaneous spectrometer sky fibers and broadband imaging via coherent fiber imaging bundles, that could reliably reduce this source of error to below the photon-noise limit of typical stellar observations. We discuss the limitations of these simulations, the underlying assumptions, and mitigation mechanisms. We also present and discuss the components designed and built into the NEID (NN-EXPLORE Exoplanet Investigations with Doppler spectroscopy) precision RV instrument for the WIYN 3.5m telescope, to serve as an ongoing resource for the community to explore and evaluate correction techniques. We emphasize that while “bright time” has been traditionally adequate for RV science, the goal of 10 cm/s precision on the most interesting exoplanetary systems may necessitate access to darker skies for these next generation instruments.

Arpita Roy↗

High-Fidelity Computational Aerodynamics of Multi-Rotor Unmanned Aerial Vehicles

High-fidelity Computational Fluid Dynamics (CFD) simulations have been carried out for several multi-rotor Unmanned Aerial Vehicles (UAVs). Three vehicles have been studied: the classic quadcopter DJI Phantom 3, an unconventional quadcopter specialized for forward flight, the SUI Endurance, and an innovative concept for Urban Air Mobility (UAM), the Elytron 4S UAV. The three-dimensional unsteady Navier-Stokes equations are solved on overset grids using high-order accurate schemes, dual-time stepping, and a hybrid turbulence model. The DJI Phantom 3 is simulated with different rotors and with both a simplified airframe and the real airframe including landing gear and a camera. The effects of weather are studied for the DJI Phantom 3 quadcopter in hover. The SUI En- durance original design is compared in forward flight to a new configuration conceived by the authors, the hybrid configuration, which gives a large improvement in forward thrust. The Elytron 4S UAV is simulated in helicopter mode and in airplane mode. Understanding the complex flows in multi-rotor vehicles will help design quieter, safer, and more efficient future drones and UAM vehicles.

High-Fidelity↗

Quantum Vision Transformers for Quark–Gluon Classification

We introduce a hybrid quantum-classical vision transformer architecture, notable for its integration of variational quantum circuits within both the attention mechanism and the multi-layer perceptrons. The research addresses the critical challenge of computational efficiency and resource constraints in analyzing data from the upcoming High Luminosity Large Hadron Collider, presenting the architecture as a potential solution. In particular, we evaluate our method by applying the model to multi-detector jet images from CMS Open Data. The goal is to distinguish quark-initiated from gluon-initiated jets. We successfully train the quantum model and evaluate it via numerical simulations. Using this approach, we achieve classification performance almost on par with the one obtained with the completely classical architecture, considering a similar number of parameters.

Comajoan Cara, Marçal (ORCID:0009000126263752)↗

Planar shock-induced bubble collapse and jetting in water captured via X-ray phase contrast imaging

Shock wave–bubble interactions in water manifest rich dynamics driven by a combination of strong pressure and density mismatches. They have a wide variety of applications, including the injection of pharmaceuticals, and through scaling, enable the exploration of various aspects of high-energy-density systems such as inertial confinement fusion. In this work, the interaction between a micrometric nitrogen bubble and a planar shock wave, characterized by a Mach number of M = 1.24 and a peak pressure of p max = 0.57 , is experimentally recorded using ultra-high-speed x-ray phase contrast imaging. Highly resolved radiographs provide access to all phase discontinuities along the beam path, offering quantities such as the time-varying bubble size, the speed of a jet produced during the bubble collapse, and the time evolution of the shock wave front, which are critical benchmark data for numerical scheme validation. This study addresses the lack of well-characterized, repeatable, and high spatiotemporal resolution experiments at negative Atwood numbers by providing shock–bubble visualization and corresponding numerical simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adaptive variational quantum dynamics simulations with compressed circuits and fewer measurements

The adaptive variational quantum dynamics simulation (AVQDS) method performs real-time evolution of quantum states using automatically generated parametrized quantum circuits that often contain substantially fewer gates than Trotter circuits. Here we report an improved version of the method, which we call AVQDS(T), by porting the tiling efficient trial circuits with rotations implemented simultaneously technique. The algorithm adaptively adds layers of disjoint unitary gates to the ansatz circuit so as to keep the McLachlan distance, a measure of the accuracy of the variational dynamics, below a fixed threshold. Here we perform benchmark noiseless AVQDS(T) simulations of quench dynamics in local spin models and compare with an alternative adaptive variational approach on quantum resource requirement. Quantum dynamical simulations implementing realistic noise channels are also reported. Finally, we propose a way to substantially alleviate the measurement overhead of AVQDS(T) while maintaining high accuracy by synergistically integrating quantum circuit calculations on quantum processing units with classical calculations using, e.g., tensor networks to evaluate the quantum geometric tensor. We showcase that this approach enables AVQDS(T) to deliver more accurate results than simulations using a fixed ansatz of comparable final depth for a significant time duration with fewer quantum resources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Report for the ASCR Workshop on Basic Research Needs in Quantum Computing and Networking - 2023

Employing quantum mechanical resources in computing and networking opens the door to new computation and communication models and potential disruptive advantages over classical counterparts. However, quantifying and realizing such advantages face extensive scientific and engineering challenges. Investments by the Department of Energy (DOE) have driven progress toward addressing such challenges. Quantum algorithms have been recently developed, in some cases offering asymptotic exponential advantages in speed or accuracy, for fundamental scientific problems such as simulating physical systems, solving systems of linear equations, or solving differential equations. Empirical demonstrations on nascent quantum hardware suggest better performance than classical analogs on specialized computational tasks favorable to the quantum computing systems. However, demonstration of an end-to-end, substantial and rigorously quantifiable quantum performance advantage over classical analogs remains a grand challenge, especially for problems of practical value. The definition of requirements for quantum technologies to exhibit scalable, rigorous, and transformative performance advantages for practical applications also remains an outstanding open question, namely, what will be required to ultimately demonstrate practical quantum advantage?

97 MATHEMATICS AND COMPUTING↗

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗

Interacting scales and energy transfer in isotropic turbulence

The dependence of the energy transfer process on the disparity of the interacting scales is investigated in the inertial and far-dissipation ranges of isotropic turbulence. The strategy for generating the simulated flow fields and the choice of a disparity parameter to characterize the scaling of the interactions is discussed. The inertial range is found to be dominated by relatively local interactions, in agreement with the Kolmogorov assumption. The far-dissipation is found to be dominated by relatively non-local interactions, supporting the classical notion that the far-dissipation range is slaved to the Kolmogorov scales. The measured energy transfer is compared with the classical models of Heisenberg, Obukhov, and the more detailed analysis of Tennekes and Lumley. The energy transfer statistics measured in the numerically simulated flows are found to be nearly self-similar for wave numbers in the inertial range. Using the self-similar form measured within the limited scale range of the simulation, an 'ideal' energy transfer function and the corresponding energy flux rate for an inertial range of infinite extent are constructed. From this flux rate, the Kolmogorov constant is calculated to be 1.5, in excellent agreement with experiments.

Zhou, YE↗

Improved honeycomb and hyperhoneycomb lattice Hamiltonians for quantum simulations of non-Abelian gauge theories

Improved Kogut-Susskind Hamiltonians for quantum simulations of non-Abelian Yang-Mills gauge theories are developed for honeycomb (2+1⁢D) and hyperhoneycomb (3+1⁢D) spatial tessellations. This is motivated by the desire to identify lattices for quantum simulations that involve only 3-link vertices among the gauge field group spaces in order to reduce the complexity in applications of the plaquette operator. For the honeycomb lattice, we derive a classically 𝒪⁡(𝑏 2 )-improved Hamiltonian, with 𝑏 being the lattice spacing. Tadpole improvement via the mean-field value of the plaquette operator is used to provide the corresponding quantum improvements. We have identified the (nonchiral) hyperhoneycomb as a candidate spatial tessellation for 3+1⁢D quantum simulations of gauge theories, and determined the associated 𝒪⁡(𝑏)-improved Hamiltonian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling of an Atmospheric-Boundary-Layer Profile in Support of Experiments in the NASA Langley Transonic Dynamics Tunnel

This paper presents the results of the Computational Fluid Dynamics (CFD) analyses to model two atmospheric-boundary-layer (ABL) profiles inside the NASA Langley Research Center Transonic Dynamics Tunnel (TDT). The CFD models include tunnel walls and all additional hardware needed to simulate the ABL profiles. This hardware includes Irwin spires and floor-mounted roughness elements. The numerical simulations show that the application of higher-fidelity numerical methods is necessary to compute boundary-layer and turbulenceintensity profiles that match experimental data. The ABL profiles are computed both inside the numerical model of the TDT and in a classical free-air model. Both Unsteady Reynolds Averaged Navier-Stokes (URANS) with Spalart-Allmaras (SA) turbulence model and Modified Delayed Detached Eddy (MDDES) simulation methods are used. The results show that the MDDES-simulation results match the experimental data very well while URANS-SA does not. The results also show that the wind-tunnel walls have a significant effect on ABL prediction.

Chwalowski, Pawel↗