SEARCH · Search NASA
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.
With a view to applying the generator coordinate method to large configuration spaces, we propose a simple approximate formula to compute diabatic many-body matrix elements without having to evaluate two-body interaction matrix elements. Here, the method is illustrated with two analytically solvable Hamiltonians based on the harmonic oscillator.
We adapt the proton-neutron finite-amplitude method, which in its original form is an efficient implementation of the Skyrme quasiparticle random phase approximation, to include the coupling of quasiparticles to like-particle phonons. The approach allows us to add beyond-quasiparticle random-phase approximation correlations to computations of Gamow-Teller strength and 𝛽-decay rates in deformed nuclei for the first time. We test the approach in several deformed isotopes for which measured strength distributions are available. Furthermore, the additional correlations dramatically improve agreement with the data, and will lead to improved global 𝛽-decay rates.
It is well known that the neutrino flavor in extreme astrophysical environments changes under the effect of three contributions: the vacuum oscillation, the interaction with the surrounding matter, and the collective oscillations due to interactions between different neutrinos. The latter adds a nonlinear contribution to the equations of motion, making the description of their dynamics complex. In this work we study various strategies to simulate the coherent collective oscillations of a system of N neutrinos in the two-flavor approximation using quantum computation. This was achieved by using a pair-neutrino decomposition designed to account for the fact that the flavor Hamiltonian, in the presence of the neutrino-neutrino term, presents an all-to-all interaction that makes the implementation of the evolution dependent on the qubit topology. We analyze the Trotter error caused by the decomposition demonstrating that the complexity of the implementation of time evolution scales polynomially with the number of neutrinos and that the noise from near-term quantum device simulation can be reduced by optimizing the quantum circuit decomposition and exploiting a full-qubit connectivity. We find that the gate complexity using second order Trotter-Suzuki formulas scales better with system size than with other decomposition methods such as quantum signal processing. In conclusion, we finally present the application and the results of our algorithm on a real quantum device based on trapped-ion qubits.
Imaginary-time response functions of finite-temperature quantum systems are often obtained with methods that exhibit stochastic or systematic errors. Reducing these errors comes at a large computational cost—in quantum Monte Carlo simulations, the reduction of noise by a factor of two incurs a simulation cost of a factor of four. In this paper, we relate certain imaginary-time response functions to an inner product on the space of linear operators on Fock space. We then show that data with noise typically does not respect the positive definiteness of its associated Gramian. The Gramian has the structure of a Hankel matrix. As a method for denoising noisy data, we introduce an alternating projection algorithm that finds the closest positive definite Hankel matrix consistent with noisy data. We test our methodology at the example of fermion Green's functions for continuous-time quantum Monte Carlo data and show remarkable improvements of the error, reducing noise by a factor of up to 20 in practical examples. We argue that Hankel projections should be used whenever finite-temperature imaginary-time data of response functions with errors is analyzed, be it in the context of quantum Monte Carlo, quantum computing, or in approximate semianalytic methodologies. Published by the American Physical Society 2024
Abstract We consider the problem of including $$\varLambda $$ Λ hyperons into the ab initio framework of nuclear lattice effective field theory. In order to avoid large sign oscillations in Monte Carlo simulations, we make use of the fact that the number of hyperons is typically small compared to the number of nucleons in the hypernuclei of interest. This allows us to use the impurity lattice Monte Carlo method, where the minority species of fermions in the full nuclear Hamiltonian is integrated out and treated as a worldline in Euclidean projection time. The majority fermions (nucleons) are treated as explicit degrees of freedom, with their mutual interactions described by auxiliary fields. This is the first application of the impurity lattice Monte Carlo method to systems where the majority particles are interacting. Here, we show how the impurity Monte Carlo method can be applied to compute the binding energies of the light hypernuclei. In this exploratory work we use spin-independent nucleon–nucleon and hyperon–nucleon interactions to test the computational power of the method. We find that the computational effort scales approximately linearly in the number of nucleons. The results are very promising for future studies of larger hypernuclear systems using chiral effective field theory and realistic hyperon–nucleon interactions, as well as applications to other quantum many-body systems.
Doppler and Raman lidar observations of vertical velocity and water vapor mixing ratio are used to probe the physics and statistics of subcloud and cloud-base latent heat fluxes during cumulus convection at the ARM Southern Great Plains (SGP) site in Oklahoma, United States. The statistical results show that latent heat fluxes increase with height from the surface up to ~0.8 Z i (where Z i is the convective boundary layer depth) and then decrease to ~0 at Z i . Peak fluxes aloft exceeding 500 W m −2 are associated with periods of increased cumulus cloud cover and stronger jumps in the mean humidity profile. These entrainment fluxes are much larger than the surface fluxes, indicating substantial drying over the 0–0.8 Z i layer accompanied by moistening aloft as the CBL deepens over the diurnal cycle. We also show that the boundary layer humidity budget is approximately closed by computing the flux divergence across the 0–0.8 Z i layer. Composite subcloud velocity and water vapor anomalies show that clouds are linked to coherent updraft and moisture plumes. The moisture anomaly is Gaussian, most pronounced above 0.8 Z i and systematically wider than the velocity anomaly, which has a narrow central updraft flanked by downdrafts. This size and shape disparity results in downdrafts characterized by a high water vapor mixing ratio and thus a broad joint probability density function (JPDF) of velocity and mixing ratio in the upper CBL. We also show that cloud-base latent heat fluxes can be both positive and negative and that the instantaneous positive fluxes can be very large (~10 000 W m −2 ). However, since cloud fraction tends to be small, the net impact of these fluxes remains modest.
An efficient semi-Lagrangian (SL) numerical scheme for the advection of positive-definite microphysical variables in cloud-resolving models is proposed. The scheme is a linear combination of SL schemes of the first and second order. The proposed scheme is compared with two high-order, monotonic, nonlinear Eulerian schemes, in two 2D tests. In the first test, advection of a positive scalar field with large gradients is performed within complicated idealized steering flows. In the second test, a hail storm is simulated using the 2D Hebrew University Cloud Model, describing the interactions between 15 size distributions of different microphysical quantities (each distribution is described by 43 mass bins). The proposed scheme is computationally efficient, has a low numerical diffusion and a high level of mass conservation accuracy, and preserves the sum of multiple advected variables as well as the linear relationships between the advected scalar variables. Although the proposed SL scheme is not exactly positive definite, the negative values and their number, which may appear only in the case of very strong gradients, are negligible. The scheme produces results similar to those of nonlinear Eulerian advection schemes while reducing computation time by approximately a factor of 10.
Abstract We present the formulation and optimization of a Runge–Kutta-type time-stepping scheme for solving the shallow-water equations, aimed at substantially increasing the effective allowable time step over that of comparable methods. This scheme, called FB-RK(3,2), uses weighted forward–backward averaging of thickness data to advance the momentum equation. The weights for this averaging are chosen with an optimization process that employs a von Neumann–type analysis, ensuring that the weights maximize the admittable Courant number. Through a simplified local truncation error analysis and numerical experiments, we show that the method is at least second-order in time for any choice of weights and exhibits low dispersion and dissipation errors for well-resolved waves. Further, we show that an optimized FB-RK(3,2) can take time steps up to 2.8 times as large as a popular three-stage, third-order strong stability-preserving Runge–Kutta method in a quasi-linear test case. In fully nonlinear shallow-water test cases relevant to oceanic and atmospheric flows, FB-RK(3,2) outperforms SSPRK3 in admittable time step by factors roughly between 1.6 and 2.2, making the scheme approximately twice as computationally efficient with little to no effect on solution quality. Significance Statement The purpose of this work is to develop and optimize time-stepping schemes for models relevant to oceanic and atmospheric flows. Specifically, for the shallow-water equations we optimize for schemes that can take time steps as large as possible while retaining solution quality. We find that our optimized schemes can take time steps between 1.6 and 2.2 times larger than schemes that cost the same number of floating point operations, translating directly to a corresponding speedup. Our ultimate goal is to use these schemes in climate-scale simulations.
In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.
Sudden approximation applied to computation of rotational transition probability and inelastic total cross sections for scattering of polar and nonpolar diatomic molecules by atoms
Approximate equations for computing enthalpy and velocity derivatives at surface of flat plate, cone, or blunt axisymmetric body
Direct summation of series involving higher transcendental functions, integrals of confluent hypergeometric functions, and computer methods for approximating continuous functions
Algorithmic and digital computer method for approximating irrational transfer functions by rational functions in control systems
In considering transfer phenomena in reacting gas mixtures, the following topics are discussed: (1) dynamics of paired collisions; (2) Enskog method for solving the Boltzmann equation, taking into account inelastic collisions; (3) use of the variation principle for approximate solutions; (4) computations of inelastic bracketed expressions; (5) calculation of collision integrals; and (6) formulas for transfer coefficients.
Review of recent work in recasting the energy type of approximation to aircraft flight in terms of singular perturbation theory and in extending this theory to three-dimensional maneuvers. Singular perturbations for differential equations arising in optimal control are first examined for a system of fairly general form but low order. The attitude dynamics for optimal flight of a rocket in vacuum is then studied as an introductory, and fairly transparent, example. The question of the choice of variables is then discussed. Finally, optimal aircraft flight in various reduced-order approximations is investigated. In particular, the problem of three-dimensional aircraft flight is formulated for singular perturbation treatment, and possibilities for decoupling into several lower-order problems are illustrated. The use of reduced-order approximation facilitates numerical computations by reducing the number of multiplier initial values that must be determined simultaneously and by improving the conditioning of the differential equations.
A fast direct (noniterative) 'Cauchy-Riemann Solver' is developed for solving the finite-difference equations representing systems of first-order elliptic partial differential equations in the form of the nonhomogeneous Cauchy-Riemann equations. The method is second-order accurate and requires approximately the same computer time as a fast cyclic-reduction Poisson solver. The accuracy and efficiency of the direct solver are demonstrated in an application to solving an example problem in aerodynamics: subsonic inviscid flow over a biconvex airfoil. The analytical small-perturbation solution contains singularities, which are captured well by the computational technique. The algorithm is expected to be useful in nonlinear subsonic and transonic aerodynamics.
Description of a model developed to predict the interactions between multiple electrostatic probes separated by distances that are shorter than or comparable to the characteristic plasma shielding length. The computation involves crude approximations accounting for space charge effects, currents due to the emission and exchange of photoelectrons, and currents from the ambient plasma. Floating potentials are calculated by requiring current balance for each probe and solving by iterating the resulting coupled nonlinear equations. Although only qualitative accuracy is claimed, predictions using the model are found to be in good agreement with observations from the Explorer 45 satellite. The Explorer 45 satellite measures the potential difference between a pair of symmetrical spherical floating probes. The Explorer 45 data show that with increasing altitude, first, the sunward sphere will typically assume a more negative potential than the antisunward sphere; second, the polarities of the probes in the region L = 2 to 4 may switch, sometimes several times; and, finally, near the plasmapause the detector output may saturate, the sunward sphere being more negative.
The damping method was applied successfully to real structures as well as analytical models. It depends on the ability to determine an appropriate modal damping matrix for each substructure. In the past, modal damping matrices were assumed diagonal for lack of being able to determine the coupling terms which are significant in the general case of nonproportional damping. This problem was overcome by formulating the damped equations of motion as a linear perturbation of the undamped equations for light structural damping. Damped modes are defined as complex vectors derived from the complex frequency response vectors of each substructure and are obtained directly from sinusoidal vibration tests. The damped modes are used to compute first order approximations to the modal damping matrices. The perturbation approach avoids ever having to solve a complex eigenvalue problem.