Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Number Theory”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Complex Dynamics of Air Traffic Flow

Air traffic in the United States has continued to grow at a steady pace since 1980, except for a dip immediately after the tragic events of September 11, 2001. There are different growth scenarios associated both with the magnitude and the composition of the future air traffic. The Terminal Area Forecast (TAF), prepared every year by the FAA, projects the growth of traffic in the United States. Both Boeing and Airbus publish market outlooks for air travel annually. Although predicting the future growth of traffic is difficult, there are two significant trends: heavily congested major airports continue to see an increase in traffic, and the emergence of regional jets and other smaller aircraft with fewer passengers operating directly between non-major airports. The interaction between air traffic demand and the ability of the system to provide the necessary airport and airspace resources can be modeled as a network. The size of the resulting network varies depending on the choice of its nodes. It would be useful to understand the properties of this network to guide future design and development. Many questions, such as the growth of delay with increasing traffic demand and impact of the en route weather on future air traffic, require a systematic understanding of the properties of the air traffic network. There has been a major advance in the understanding of the behavior of networks with a large number of components. Several theories have been advanced about the evolution of large biological and engineering networks by authors in diversified disciplines like physics, mathematics, biology and computer science. Several networks exhibit a scale-free property in the sense that the probabilistic distribution of their nodes as a function of connections decreases slower than an exponential. These networks are characterized by the fact that a small number of components have a disproportionate influence on the performance of the network. Scale-free networks are tolerant to random failure of components, but are vulnerable to selective attack on components. This paper examines two network representations for the baseline air traffic system. A network defined with the 40 major airports as nodes and with standard flight routes as links has a characteristic scale: all nodes have 60 or more links and no node has more than 460 links. Another network is defined with baseline aircraft routing structure exhibits an exponentially truncated scale-free behavior. Its degree ranges from 2 connections to 2900 connections, and 225 nodes have more than 250 connections. Furthermore, those high-degree nodes are homogeneously distributed in the airspace. A consequence of this scale-free behavior is that the random loss of a single node has little impact, but the loss of multiple high-degree nodes (such as occurs during major storms in busy airspace) can adversely impact the system. Two future scenarios of air traffic growth are used to predict the growth of air traffic in the United States. It is shown that a three-times growth in the overall traffic may result in a ten-times impact on the density of traffic in certain parts of the United States.

Scale-free Networks↗
DOE OSTI2024

A collinear perspective on the Regge limit

The high energy (Regge) limit provides a playground for understanding all loop structures of scattering amplitudes, and plays an important role in the description of many phenomenologically relevant cross-sections. While well understood in the planar limit, the structure of non-planar corrections introduces many fascinating complexities, for which a general organizing principle is still lacking. We study the structure of multi-reggeon exchanges in the context of the effective field theory for forward scattering, and derive their factorization into collinear operators (impact factors) and soft operators. We derive the structure of the renormalization group consistency equations in the effective theory, showing how the anomalous dimensions of the soft operators are related to those of the collinear operators, allowing us to derive renormalization group equations in the Regge limit purely from a collinear perspective. The rigidity of the consistency equations provides considerable insight into the all orders organization of Regge amplitudes in the effective theory, as well as its relation to other approaches. Along the way we derive a number of technical results that improve the understanding of the effective theory. We illustrate this collinear perspective by re-deriving all the standard BFKL equations for two-Glauber exchange from purely collinear calculations, and we show that this perspective provides a number of conceptual and computational advantages as compared to the standard view from soft or Glauber physics. We anticipate that this formulation in terms of collinear operators will enable a better understanding of the relation between BFKL and DGLAP in gauge theories, and facilitate the analysis of renormalization group evolution equations describing Reggeization beyond next-to-leading order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Machine learned potential for high-throughput phonon calculations of metal—organic frameworks

Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.

Elena, Alin Marin↗
DOE OSTI2025

Efficient Hamiltonian encoding algorithms for extracting quantum control mechanism as interfering pathway amplitudes in the Dyson series

Hamiltonian encoding is a methodology for revealing the mechanism behind the dynamics governing controlled quantum systems. In this paper, following Mitra and Rabitz \cite{abhra_1}, we define mechanism via pathways of eigenstates that describe the evolution of the system, where each pathway is associated with a complex-valued amplitude corresponding to a term in the Dyson series. The evolution of the system is determined by the constructive and destructive interference of these pathway amplitudes. Pathways with similar attributes can be grouped together into pathway classes. The amplitudes of pathway classes are computed by modulating the Hamiltonian matrix elements and decoding the subsequent evolution of the system rather than by direct computation of the individual terms in the Dyson series. The original implementation of Hamiltonian encoding was computationally intensive and became prohibitively expensive in large quantum systems. This paper presents two new encoding algorithms that calculate the amplitudes of pathway classes by using techniques from graph theory and algebraic topology to exploit patterns in the set of allowed transitions, greatly reducing the number of matrix elements that need to be modulated. These new algorithms provide an exponential decrease in both computation time and memory utilization with respect to the Hilbert space dimension of the system. To demonstrate the use of these techniques, they are applied to two illustrative state-to-state transition problems.

Abrams, Erez [Princeton University, Massachusetts ↗
NASA NTRS1979

Inverse boundary-layer technique for airfoil design

A description is presented of a technique for the optimization of airfoil pressure distributions using an interactive inverse boundary-layer program. This program allows the user to determine quickly a near-optimum subsonic pressure distribution which meets his requirements for lift, drag, and pitching moment at the desired flow conditions. The method employs an inverse turbulent boundary-layer scheme for definition of the turbulent recovery portion of the pressure distribution. Two levels of pressure-distribution architecture are used - a simple roof top for preliminary studies and a more complex four-region architecture for a more refined design. A technique is employed to avoid the specification of pressure distributions which result in unrealistic airfoils, that is, those with negative thickness. The program allows rapid evaluation of a designed pressure distribution off-design in Reynolds number, transition location, and angle of attack, and will compute an airfoil contour for the designed pressure distribution using linear theory.

Henderson, M. L.↗
NASA NTRS1992

Renormalization group methods for the Reynolds stress transport equations

The Yakhot-Orszag renormalization group is used to analyze the pressure gradient-velocity correlation and return to isotropy terms in the Reynolds stress transport equations. The perturbation series for the relevant correlations, evaluated to lowest order in the epsilon-expansion of the Yakhot-Orszag theory, are infinite series in tensor product powers of the mean velocity gradient and its transpose. Formal lowest order Pade approximations to the sums of these series produce a rapid pressure strain model of the form proposed by Launder, Reece, and Rodi, and a return to isotropy model of the form proposed by Rotta. In both cases, the model constants are computed theoretically. The predicted Reynolds stress ratios in simple shear flows are evaluated and compared with experimental data. The possibility is discussed of deriving higher order nonlinear models by approximating the sums more accurately. The Yakhot-Orszag renormalization group provides a systematic procedure for deriving turbulence models. Typical applications have included theoretical derivation of the universal constants of isotropic turbulence theory, such as the Kolmogorov constant, and derivation of two equation models, again with theoretically computed constants and low Reynolds number forms of the equations. Recent work has applied this formalism to Reynolds stress modeling, previously in the form of a nonlinear eddy viscosity representation of the Reynolds stresses, which can be used to model the simplest normal stress effects. The present work attempts to apply the Yakhot-Orszag formalism to Reynolds stress transport modeling.

Rubinstein, R.↗
NASA NTRS1993

Multigrid methods for flow transition in three-dimensional boundary layers with surface roughness

The efficient multilevel adaptive method has been successfully applied to perform direct numerical simulations (DNS) of flow transition in 3-D channels and 3-D boundary layers with 2-D and 3-D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semi-coarsening multigrid method associated with line distributive relaxation scheme, and an improved outflow boundary-condition treatment, which needs only a very short buffer domain to damp all order-one wave reflections, are developed. These approaches make the multigrid DNS code very accurate and efficient. This allows us not only to be able to do spatial DNS for the 3-D channel and flat plate at low computational costs, but also to do spatial DNS for transition in the 3-D boundary layer with 3-D single and multiple roughness elements, which would have extremely high computational costs with conventional methods. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments. The contribution of isolated and distributed roughness to transition is analyzed.

Liu, Chaoqun↗
NASA NTRS1994

Ion properties of the high-latitude middle-atmosphere

Investigations of high-latitude, middle-atmosphere ion properties have been conducted using rocket probe techniques. Electrodynamics studies have addressed: mesospheric dynamics and small-scale structure (the MAC/EPSILON Campaign during October-November 1977); ionization effects associated with highly energetic greater than or equal to 1 MeV) precipitating electrons (the MAC/REP Program during May 1990); and the characteristics of noctilucent cloud (NLC) and polar mesosphere summer echo (PMSE) regions (the NLC-91 Campaign during July-August 1991). Recent developments in Gerdien condenser probe theory, using computational fluid dynamics to model ion collection, have resulted in more accurate evaluation of the ion mobilities and number densities. Measurements of very small ion mobilities are thought to indicate an aerosol ion presence in the upper mesosphere.

Li, C.↗
DOE OSTI2026

Degenerate coupled-cluster theory

A size-extensive, converging, black-box, ab initio coupled-cluster (ΔCC) ansatz is introduced that computes the energies and wave functions of states from any degenerate or nondegenerate Slater-determinant references with any numbers of α- and β-spin electrons, any patterns of orbital occupancy, any spin multiplicities, and any spatial symmetries. For a nondegenerate reference, it reduces to the single-reference coupled-cluster ansatz. For a degenerate multireference, it is a natural coupled-cluster extension of degenerate Møller–Plesset perturbation (ΔMP) theory. For ionized and electron-attached references, it is a coupled-cluster Green’s function, although the present theory is convergent toward the full-configuration-interaction limits, while the Feynman–Dyson many-body Green’s function (MBGF) theory generally is not. Its single-excitation instance is a projection Hartree–Fock theory as per the Thouless theorem, which may be useful for core ionizations, high-spin states, and possibly electron affinities. Additionally, a new multireference coupled-cluster theory for a general model space is developed. This quasidegenerate coupled-cluster (QCC) theory is exactly converging, but not black-box, and intended for strong correlation. Determinant-based, general-order algorithms of ΔCC and QCC theories are implemented and compared with configuration-interaction (CI) and equation-of-motion coupled-cluster (EOM-CC) theories through octuple excitations and with ΔMP and MBGF theories up to the nineteenth order. An algebraic, optimal-scaling algorithm of the ΔCC theory is computer-synthesized at the levels of single excitations (ΔCCS) and of single and double excitations (ΔCCSD). As a result, the order of performance is QCC ≈ ΔCC > EOM-CC > CI at the same order or QCC ≈ ΔCC > ΔMP > MBGF at the same cost scaling.

Hirata, So [University of Illinois at Urbana-Champ↗
NASA NTRS1981

Progress on interacting boundary-layer computations at high Reynolds number

The purpose of this paper is to review progress made in the solution of the interacting boundary-layer equations for subsonic flow. The interrelationship of triple deck theory and the interacting boundary-layer approach is discussed with emphasis placed on the development of efficient and reliable algorithms for the solution of the interacting boundary-layer equations. Example studies are presented for laminar and turbulent finite flat plate flow, laminar flow past a flat plate with a separation causing depression, and laminar and turbulent flow past a blunt based trailing edge.

Davis, R. T.↗
DOE OSTI2026

Clifford transformations for fermionic quantum systems: From Pauli and Majorana operators to Dirac fermions

Clifford gates and transformations, which map products of elementary Pauli or Majorana operators to other such products, are foundational in quantum computing, underpinning the stabilizer formalism, error-correcting codes, magic state distillation, quantum communication and cryptography, and qubit tapering. Moreover, circuits composed entirely of Clifford gates are classically simulatable, highlighting their computational significance. In this article we extend the concept of Clifford transformations to Dirac fermions. We demonstrate that discrete Clifford transformations are generated by half-body and pair operators while continuous Clifford transformations are generated by number operators, providing a systematic framework for their characterization. Additionally, we establish connections with fermionic mean-field theories and applications in qubit tapering, offering insights into their broader implications in quantum computing.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1993

Stability and instability of thermocapillary convection in models of the float-zone crystal-growth process

This project was concerned with the determination of conditions of guaranteed stability and instability for thermocapillary convection in a model of the float-zone crystal-growth process. This model, referred to as the half-zone, was studied extensively, both experimentally and theoretically. Our own earlier research determined, using energy-stability theory, sufficient conditions for stability to axisymmetric disturbances. Nearly all results computed were for the case of a liquid with Prandtl Number Pr = 1. Attempts to compute cases for higher Prandtl numbers to allow comparison with the experimental results of other researchers were unsuccessful, but indicated that the condition guaranteeing stability against axisymmetric disturbances would be a value of the Marangoni number (Ma), significantly higher than that at which oscillatory convection was observed experimentally. Thus, additional results were needed to round out the stability picture for this model problem. The research performed under this grant consisted of the following: (1) computation of energy-stability limits for non-axisymmetric disturbances; (2) computation of linear-stability limits for axisymmetric and non-axisymmetric disturbances; (3) numerical simulation of the basic state for half- and full-zones with a deformable free surface; and (4) incorporation of radiation heat transfer into a model energy-stability problem. Each of these is summarized briefly below.

Neitzel, G. P.↗
NASA NTRS1956

Second-order subsonic airfoil theory including edge effects

Several recent advances in plane subsonic flow theory are combined into a unified second-order theory for airfoil sections of arbitrary shape. The solution is reached in three steps: the incompressible result is found by integration, it is converted into the corresponding subsonic compressible result by means of the second-order compressibility rule, and it is rendered uniformly valid near stagnation points by further rules. Solutions for a number of airfoils are given and are compared with the results of other theories and of experiment. A straight-forward computing scheme is outlined for calculating the surface velocities and pressures on any airfoil at any angle of attack

Van Dyke, Milton D↗
NASA NTRS2023

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow↗
NASA NTRS1989

On the Lagrangian description of unsteady boundary layer separation. Part 1: General theory

Although unsteady, high-Reynolds number, laminar boundary layers have conventionally been studied in terms of Eulerian coordinates, a Lagrangian approach may have significant analytical and computational advantages. In Lagrangian coordinates the classical boundary layer equations decouple into a momentum equation for the motion parallel to the boundary, and a hyperbolic continuity equation (essentially a conserved Jacobian) for the motion normal to the boundary. The momentum equations, plus the energy equation if the flow is compressible, can be solved independently of the continuity equation. Unsteady separation occurs when the continuity equation becomes singular as a result of touching characteristics, the condition for which can be expressed in terms of the solution of the momentum equations. The solutions to the momentum and energy equations remain regular. Asymptotic structures for a number of unsteady 3-D separating flows follow and depend on the symmetry properties of the flow. In the absence of any symmetry, the singularity structure just prior to separation is found to be quasi 2-D with a displacement thickness in the form of a crescent shaped ridge. Physically the singularities can be understood in terms of the behavior of a fluid element inside the boundary layer which contracts in a direction parallel to the boundary and expands normal to it, thus forcing the fluid above it to be ejected from the boundary layer.

Vandommelen, Leon L.↗
NASA NTRS1990

On the Lagrangian description of unsteady boundary-layer separation. I - General theory

Although unsteady, high-Reynolds number, laminar boundary layers have conventionally been studied in terms of Eulerian coordinates, a Lagrangian approach may have significant analytical and computational advantages. In Lagrangian coordinates the classical boundary layer equations decouple into a momentum equation for the motion parallel to the boundary, and a hyperbolic continuity equation (essentially a conserved Jacobian) for the motion normal to the boundary. The momentum equations, plus the energy equation if the flow is compressible, can be solved independently of the continuity equation. Unsteady separation occurs when the continuity equation becomes singular as a result of touching characteristics, the condition for which can be expressed in terms of the solution of the momentum equations. The solutions to the momentum and energy equations remain regular. Asymptotic structures for a number of unsteady 3-D separating flows follow and depend on the symmetry properties of the flow. In the absence of any symmetry, the singularity structure just prior to separation is found to be quasi 2-D with a displacement thickness in the form of a crescent shaped ridge. Physically the singularities can be understood in terms of the behavior of a fluid element inside the boundary layer which contracts in a direction parallel to the boundary and expands normal to it, thus forcing the fluid above it to be ejected from the boundary layer.

Van Dommelen, Leon L.↗
NASA NTRS2009

Roughness Based Crossflow Transition Control: A Computational Assessment

A combination of parabolized stability equations and secondary instability theory has been applied to a low-speed swept airfoil model with a chord Reynolds number of 7.15 million, with the goals of (i) evaluating this methodology in the context of transition prediction for a known configuration for which roughness based crossflow transition control has been demonstrated under flight conditions and (ii) of analyzing the mechanism of transition delay via the introduction of discrete roughness elements (DRE). Roughness based transition control involves controlled seeding of suitable, subdominant crossflow modes, so as to weaken the growth of naturally occurring, linearly more unstable crossflow modes. Therefore, a synthesis of receptivity, linear and nonlinear growth of stationary crossflow disturbances, and the ensuing development of high frequency secondary instabilities is desirable to understand the experimentally observed transition behavior. With further validation, such higher fidelity prediction methodology could be utilized to assess the potential for crossflow transition control at even higher Reynolds numbers, where experimental data is currently unavailable.

Li, Fei↗
NASA NTRS2015

Impedance Eduction in Sound Fields With Peripherally Varying Liners and Flow

A two-dimensional impedance eduction theory is extended to three-dimensional sound fields and peripherally varying duct liners. The approach is to first measure the acoustic pressure field at a series of flush-mounted wall microphones located around the periphery of the flow duct. The numerical solution for the acoustic pressure field at these microphones is also obtained by solving the three-dimensional convected Helmholtz equation using the finite element method. A quadratic objective function based on the difference between the measured and finite element solution is constructed and the unknown impedance function is obtained by minimizing this objective function. Impedance spectra educed for two uniform-structure liners (a wire-mesh and a conventional liner) and a hard-soft-hard peripherally varying liner (for which the soft segment is that of the conventional liner) are presented. Results are presented at three mean flow Mach numbers and fourteen sound source frequencies. The impedance spectra of the uniform-structure liners are also computed using a two-dimensional impedance eduction theory. The primary conclusions of the study are: 1) when measured data is used with the uniform-structure liners, the three-dimensional theory reproduces the same impedance spectra as the two-dimensional theory except for frequencies corresponding to very low or very high liner attenuation; and 2) good agreement between the educed impedance spectra of the uniform structure conventional liner and the soft segment of the peripherally varying liner is obtained.

Watson, W. R.↗