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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Energy Models for One-Carrier Transport in Semiconductor Devices

Moment models of carrier transport, derived from the Boltzmann equation, made possible the simulation of certain key effects through such realistic assumptions as energy dependent mobility functions. This type of global dependence permits the observation of velocity overshoot in the vicinity of device junctions, not discerned via classical drift-diffusion models, which are primarily local in nature. It was found that a critical role is played in the hydrodynamic model by the heat conduction term. When ignored, the overshoot is inappropriately damped. When the standard choice of the Wiedemann-Franz law is made for the conductivity, spurious overshoot is observed. Agreement with Monte-Carlo simulation in this regime required empirical modification of this law, or nonstandard choices. Simulations of the hydrodynamic model in one and two dimensions, as well as simulations of a newly developed energy model, the RT model, are presented. The RT model, intermediate between the hydrodynamic and drift-diffusion model, was developed to eliminate the parabolic energy band and Maxwellian distribution assumptions, and to reduce the spurious overshoot with physically consistent assumptions. The algorithms employed for both models are the essentially non-oscillatory shock capturing algorithms. Some mathematical results are presented and contrasted with the highly developed state of the drift-diffusion model.

Jerome, Joseph W.↗
DOE OSTI2024

Designing observables for measurements with deep learning

Many analyses in particle and nuclear physics use simulations to infer fundamental, effective, or phenomenological parameters of the underlying physics models. When the inference is performed with unfolded cross sections, the observables are designed using physics intuition and heuristics. We propose to design targeted observables with machine learning. Unfolded, differential cross sections in a neural network output contain the most information about parameters of interest and can be well-measured by construction. The networks are trained using a custom loss function that rewards outputs that are sensitive to the parameter(s) of interest while simultaneously penalizing outputs that are different between particle-level and detector-level (to minimize detector distortions). We demonstrate this idea in simulation using two physics models for inclusive measurements in deep inelastic scattering. We find that the new approach is more sensitive than classical observables at distinguishing the two models and also has a reduced unfolding uncertainty due to the reduced detector distortions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1989

Consequences of nonclassical measurement for the algorithmic description of continuous dynamical systems

Continuous dynamical systems intuitively seem capable of more complex behavior than discrete systems. If analyzed in the framework of the traditional theory of computation, a continuous dynamical system with countably many quasistable states has at least the computational power of a universal Turing machine. Such an analysis assumes, however, the classical notion of measurement. If measurement is viewed nonclassically, a continuous dynamical system cannot, even in principle, exhibit behavior that cannot be simulated by a universal Turing machine.

Fields, Chris↗
NASA NTRS1989

Consequences of nonclassical measurement for the algorithmic description of continuous dynamical systems

Continuous dynamical systems intuitively seem capable of more complex behavior than discrete systems. If analyzed in the framework of the traditional theory of computation, a continuous dynamical system with countablely many quasistable states has at least the computational power of a universal Turing machine. Such an analyses assumes, however, the classical notion of measurement. If measurement is viewed nonclassically, a continuous dynamical system cannot, even in principle, exhibit behavior that cannot be simulated by a universal Turing machine.

Fields, Chris↗
NASA NTRS2015

Microstructural Influence on Deformation and Fatigue Life of Composites Using the Generalized Method of Cells

A fully coupled deformation and damage approach to modeling the response of composite materials and composite laminates is presented. It is based on the semi-­‐analytical generalized method of cells (GMC) micromechanics model as well as its higher fidelity counterpart, HFGMC, both of which provide closed-form constitutive equations for composite materials as well as the micro scale stress and strain fields in the composite phases. The provided constitutive equations allow GMC and HFGMC to function within a higher scale structural analysis (e.g., finite element analysis or lamination theory) to represent a composite material point, while the availability of the micro fields allow the incorporation of lower scale sub­‐models to represent local phenomena in the fiber and matrix. Further, GMC's formulation performs averaging when applying certain governing equations such that some degree of microscale field accuracy is surrendered in favor of extreme computational efficiency, rendering the method quite attractive as the centerpiece in a integrated computational material engineering (ICME) structural analysis; whereas HFGMC retains this microscale field accuracy, but at the price of significantly slower computational speed. Herein, the sensitivity of deformation and the fatigue life of graphite/epoxy PMC composites, with both ordered and disordered microstructures, has been investigated using this coupled deformation and damage micromechanics based approach. The local effects of fiber breakage and fatigue damage are included as sub‐models that operate on the microscale for the individual composite phases. For analysis of laminates, classical lamination theory is employed as the global or structural scale model, while GMC/HFGMC is embedded to operate on the microscale to simulate the behavior of the composite material within each laminate layer. A key outcome of this study is the statistical influence of microstructure and micromechanics idealization (GMC or HFGMC) on the overall accuracy of unidirectional and laminated composite deformation and fatigue response.

Microscale Models↗
NASA NTRS2017

Hybrid Guidance Control for a Hypervelocity Small Size Asteroid Interceptor Vehicle

Near-Earth Objects (NEOs) are comets and/or asteroids that have orbits in proximity with Earth's own orbit. NEOs have collided with the Earth in the past, which can be seen at such places as Chicxulub crater, Barringer crater, and Manson crater, and will continue in the future with potentially significant and devastating results. Fortunately such NEO collisions with Earth are infrequent, but can happen at any time. Therefore it is necessary to develop and validate techniques as well as technologies necessary to prevent them. One approach to mitigate future NEO impacts is the concept of high-speed interceptor. This concept is to alter the NEO's trajectory via momentum exchange by using kinetic impactors as well as nuclear penetration devices. The interceptor has to hit a target NEO at relative velocity which imparts a sufficient change in NEO velocity. NASA's Deep Impact mission has demonstrated this scenario by intercepting Comet Temple 1, 5 km in diameter, with an impact relative speed of approximately 10 km/s. This paper focuses on the development of hybrid guidance navigation and control (GNC) algorithms for precision hypervelocity intercept of small sized NEOs. The spacecraft's hypervelocity and the NEO's small size are critical challenges for a successful mission as the NEO will not fill the field of view until a few seconds before intercept. The investigation needs to consider the error sources modeled in the navigation simulation such as spacecraft initial state uncertainties in position and velocity. Furthermore, the paper presents three selected spacecraft guidance algorithms for asteroid intercept and rendezvous missions. The selected algorithms are classical Proportional Navigation (PN) based guidance that use a first order difference to compute the derivatives, Three Plane Proportional Navigation (TPPN), and the Kinematic Impulse (KI). A manipulated Bennu orbit that has been changed to impact Earth will be used as a demonstrative example to compare the three methods. In addition, a hybrid approach that is a combination between proportional navigation and kinematic impulse will be investigated to find an effective, error tolerant, and power saving approach. A 3-dimension mission scenario for both the asteroid and the interceptor spacecraft software simulator is developed for testing of the controllers. The current result demonstrates that a miss distance magnitude of less than 10m is found using the PN and TPPN guidance laws for small asteroid in the presence of error in the spacecraft states. Moreover, the paper presents these results and also the hybrid control approach simulation results.

asteroid↗
DOE OSTI2026

Quantum Solver Using Singular Value Decomposition for Computational Fluid Dynamics

Numerical solutions for fluid flow problems are challenging and have been focus of Computational Fluid Dynamics (CFD) research for past several decades. The advent of quantum computing promises exponential speedup in comparison to existing classical methods and alleviate computational constraints posed by CFD problems. Although solutions for most problems of interest in fluid dynamics using quantum computing are distant, recent advances in algorithms, software and hardware provide a path towards realizing this goal. Quantum linear solver algorithms (QLSA) such as Harrow–Hassidim–Lloyd (HHL) and Variational Quantum Linear Solver (VQLS) have been successfully implemented to solve for canonical problems such as Hele-Shaw flow. However, these algorithms still suffer to scale and address problems with ill-conditioned Jacobians. In the current paper, we alleviate these restrictions with a new quantum solver based on Singular Value Decomposition (SVD) and simulate flow past a 2D cylinder. The fidelity of the SVD based quantum solver in predicting the flow past 2D cylinder is computed along with an assessment of errors. Classical and quantum solutions for the flow are compared for different resolutions. Finally, we discuss variation in the solutions based on number of shots used.

Gottiparthi, Kalyan [ORNL] (ORCID:0000000213540255↗
NASA NTRS1994

The effect of wave-particle interactions on the polar winds O(+)

The escape of the polar wind plasma is an important element in the ionosphere-magnetosphere coupling. Both theory and observations indicate that the wave-particle interactions (WPI) play a significant role in the dynamics of ion outflow along open geomagnetic field lines. A Monte Carlo simulation was developed in order to include the effect of the WPI in addition to the factors that are traditionally included in the 'classical' polar wind (i.e. gravity, electrostatic field, and divergence of geomagnetic field lines). The ion distribution function (f(sub j)), as well as the profiles of its moments (density, drift velocity, temperature, etc.) were found for different levels of WPI, that is, for different values of the normalized diffusion rate in the velocity space (D(tilde) (sub j perpendicular). Although the model included O(=), H(+) and eletrons, we presented only the results related to the O(+) ion. We found that (1) both the density and drift velocity of O(+) increased with the WPI strength, and consequently, the O(+) escape flux was enhanced by a factor of up to 10(exp 5), (2) The O(+) ions could be energized up to a few electron volts; (3) for moderate and high levels of WPI D(tilde) (sub perpendicular) (O(+) greater than (tilde) 1, the distribution function f(O(+)) displayed very pronounced conic features at altitudes around 3 R(sub e). Finally, the interplay between the downward body force, the upward mirror force, and the perpendicular heating resulted in the formation of the 'pressure cooker' effect. This phenomena explained some interesting features of our solution, such as, the peak in the O(=) temperature, and the formation of 'ears' and conics for f(O(+)) around 2.5 R(sub e).

Barakat, A. R.↗
NASA NTRS2000

2D Quantum Simulation of MOSFET Using the Non Equilibrium Green's Function Method

The objectives this viewgraph presentation summarizes include: (1) the development of a quantum mechanical simulator for ultra short channel MOSFET simulation, including theory, physical approximations, and computer code; (2) explore physics that is not accessible by semiclassical methods; (3) benchmarking of semiclassical and classical methods; and (4) study other two-dimensional devices and molecular structure, from discretized Hamiltonian to tight-binding Hamiltonian.

Svizhenko, Alexel↗
NASA NTRS1999

Transonic Flutter Suppression Control Law Design, Analysis and 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 classical, and minimax techniques are described. A unified general formulation and solution for the minimax approach, 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↗
NASA NTRS2018

Demonstration of the Space Launch System Augmenting Adaptive Control Algorithm on Pole-Cart Platform

NASA's baseline Space Launch System (SLS) ight control system (FCS) includes an adaptive augmenting control (AAC) portion in addition to the ight-heritage nominal classical controller. The AAC algorithm is intended to improve the robustness and performance of the classical controller. Over the past several years, the AAC algorithm developed at NASA Marshall Space Flight Center (MSFC) has matured significantly through extensive simulation, rigorous analytical proofs, and a series of successful ight tests on a F18 aircraft. This study was part of a SLS program and NASA Engineering and Safety Center (NESC) joint e ort to further increase the confidence level of the AAC algorithm by demonstrating its key functionalities on a classroom type of example, the pole-cart sys- tem, at the NASA Langley Research Center (LaRC) dynamics and control laboratory. The fundamental dynamics behind balancing an inverted pendulum is similar to controlling an aerodynamically unstable rocket. Both systems are inherently open-loop unstable and requires feedback control for attitude stabilization. The principles behind the AAC algorithm is applicable to a wide range of conditionally stable dynamical systems. Hence, the outcomes from this simple and inexpensive exercise has provided the SLS program with additional confidence into the AAC design, operation, robustness, and application.

Pei, Jing↗
NASA NTRS1999

Quantum Decision-Maker

A quantum device simulating human decision making process is introduced. It consists of quantum recurrent nets generating stochastic processes which represent the motor dynamics, and a classical neural nets describing evolution of probabilities of these processes which represent the mental dynamics.

quantum↗
NASA NTRS1983

Integrated Application of Active Controls (IAAC) technology to an advanced subsonic transport project-longitudinal handling qualities study of a relaxed-stability airplane

The results of a piloted simulation of longitudinal handling qualities of an airplane with relaxed static stability are described. This task was performed under the Integrated Application of Active Controls (IAAC) Technology Project within the NASA Energy Efficient Transport Program. A representative medium range transport airplane, the Boeing Model 757, was simulated. Evaluations were made of the unaugmented airplane and of the airplane with an Essential Pitch Augmented Stability (PAS) System and with a Primary PAS System at various center of gravity (cg) conditions. Level 2 pilot ratings were attained with cg locations aft to about 57% mean aerodynamic chord (MAC) or 6% aft of the neutral point for unaugmented landing approach. For Mach = 0.80, unaugmented cruise Level 2 ratings were attained to 47% MAC or 5% forward of the maneuver point. The augmented airplane model provided handling qualities close to or within the Level 1 boundary at all cg locations for both Essential and Primary PAS. Analyses of the test conditions when compared with existing handling qualities criteria based on classical unaugmented airplane characteristics agreed well with the pilot ratings. The unaugmented results are comparable to those reported by both the Douglas Aircraft Company and Lockheed California Company from simulation investigations of transport configurations with roughly similar dimensional and mass characteristics.

Source record↗
DOE OSTI2026

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗
NASA NTRS2008

Quantum-Inspired Maximizer

A report discusses an algorithm for a new kind of dynamics based on a quantum- classical hybrid-quantum-inspired maximizer. The model is represented by a modified Madelung equation in which the quantum potential is replaced by different, specially chosen 'computational' potential. As a result, the dynamics attains both quantum and classical properties: it preserves superposition and entanglement of random solutions, while allowing one to measure its state variables, using classical methods. Such optimal combination of characteristics is a perfect match for quantum-inspired computing. As an application, an algorithm for global maximum of an arbitrary integrable function is proposed. The idea of the proposed algorithm is very simple: based upon the Quantum-inspired Maximizer (QIM), introduce a positive function to be maximized as the probability density to which the solution is attracted. Then the larger value of this function will have the higher probability to appear. Special attention is paid to simulation of integer programming and NP-complete problems. It is demonstrated that the problem of global maximum of an integrable function can be found in polynomial time by using the proposed quantum- classical hybrid. The result is extended to a constrained maximum with applications to integer programming and TSP (Traveling Salesman Problem).

Zak, Michail↗
NASA NTRS1993

Liquid propellant rocket engine combustion simulation with a time-accurate CFD method

Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.

Chen, Y. S.↗
DOE OSTI2025

Growth of Hexagonal Boron Nitride from Molten Nickel Solutions: A Reactive Molecular Dynamics Study

Metal flux methods are excellent for synthesizing high-quality hexagonal boron nitride (hBN) crystals, but the atomic mechanisms of hBN nucleation and growth in these systems are poorly understood and difficult to probe experimentally. Here, we harness classical reactive molecular dynamics (ReaxFF) to unravel the mechanisms of hBN synthesis from liquid nickel solvent over time scales up to 30 ns. These simulations mimic experimental conditions by including relatively large liquid nickel slabs containing dissolved boron and a molecular nitrogen gas phase. Overall, the reaction takes place almost exclusively on the surface of the liquid nickel, owing to the low solubility of nitrogen in bulk nickel and the intermediate species’ preference for the metal–gas interface. The formation of hBN invariably begins by reaction of dinitrogen with nickel-solvated boron atoms at the surface, forming intermediate N–N–B species, which typically evolve into B–N–B units through a short-lived intermediate where a single nitrogen atom is coordinated by one nitrogen and two boron atoms. The resulting B–N–B units, in turn, coalesce with growing hBN nuclei and carry nitrogen between hBN nanocrystals in an Ostwald ripening process. The amount of hBN produced on the tens of nanosecond time scale depends critically on the boron concentration, while having a much weaker dependence on the N 2 pressure for the regime considered (N 2 pressures of 2.5–10 MPa, Ni–B solutions with 6–12% boron by atom fraction). The highest rate of hBN formation occurs at the lowest temperature considered (1750 K, just above the melting point of nickel), while no hBN sheets are formed at 2000 K or above. An analysis of the transition pathways for nitrogen atoms shows that the final step, incorporation of small B–N motifs into larger hBN sheets, is the rate-limiting step in the regimes considered. While raising the temperature from 1750 to 2000 K has little effect on the formation of intermediates (N–N–B, B–N–B, etc.), the lack of large hBN sheets at temperatures >1900 K is explained by decreased probability of the final step and increased probability of breakup of hBN into B–N motifs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Exponential concentration in quantum kernel methods

Kernel methods in Quantum Machine Learning (QML) have recently gained significant attention as a potential candidate for achieving a quantum advantage in data analysis. Among other attractive properties, when training a kernel-based model one is guaranteed to find the optimal model’s parameters due to the convexity of the training landscape. However, this is based on the assumption that the quantum kernel can be efficiently obtained from quantum hardware. In this work we study the performance of quantum kernel models from the perspective of the resources needed to accurately estimate kernel values. We show that, under certain conditions, values of quantum kernels over different input data can be exponentially concentrated (in the number of qubits) towards some fixed value. Thus on training with a polynomial number of measurements, one ends up with a trivial model where the predictions on unseen inputs are independent of the input data. We identify four sources that can lead to concentration including expressivity of data embedding, global measurements, entanglement and noise. For each source, an associated concentration bound of quantum kernels is analytically derived. Lastly, we show that when dealing with classical data, training a parametrized data embedding with a kernel alignment method is also susceptible to exponential concentration. Our results are verified through numerical simulations for several QML tasks. Altogether, we provide guidelines indicating that certain features should be avoided to ensure the efficient evaluation of quantum kernels and so the performance of quantum kernel methods.

97 MATHEMATICS AND COMPUTING↗