Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient methods”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Unsteady turbulent boundary layers in two-dimensional, incompressible flow.

The paper presents the results of some calculations aimed at studying the characteristics of turbulent boundary layers in unsteady flow. The calculations were performed using the method of Singleton and Nash, which is based on the turbulent kinetic energy equation. Emphasis is placed on the approach to zero wall stress, and the results show the extent to which onset of this condition can be delayed by the effects of time dependence. The delay is not associated merely with alleviation of the pressure gradients, and quasi-steady methods cannot readily be used to calculate this type of flow.

Nash, J. F.↗
DOE OSTI2023

Design Criteria for High-Gradient Radio-Frequency Linacs

This article will review methods used at the SLAC National Accelerator Laboratory and other world accelerator laboratories to design high-gradient normal conducting accelerating structures. A quest for compact radio-frequency linacs fueled decades of studies toward a higher accelerating gradient. A major phenomena limiting the increase of the gradient is vacuum radio-frequency breakdown; therefore, this paper will address the breakdown physics and discuss approaches that reduce the breakdown probability. This discussion will cover both the electrical design and fabrication technology of the accelerating structures to achieve practical operating accelerating gradients in excess of 100 MV/m. Most of the data described here were obtained during the development of 11 GHz linacs for electron–positron linear colliders, so extrapolation of the results to other frequencies should be performed cautiously.

Dolgashev, Valery A. (ORCID:0000000166715402)↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
NASA NTRS1992

Iterative solution of linear systems

Recent advances in the field of iterative methods for solving large linear systems are reviewed. The main focus is on developments in the area of conjugate gradient-type algorithms and Krylov subspace methods for nonHermitian matrices.

Freund, Roland W.↗
NASA NTRS1996

Computation of Kinetics for the Hydrogen/Oxygen System Using the Thermodynamic Method

A new method for predicting chemical rate constants using thermodynamics has been applied to the hydrogen/oxygen system. This method is based on using the gradient of the Gibbs free energy and a single proportionality constant D to determine the kinetic rate constants. Using this method the rate constants for any gas phase reaction can be computed from thermodynamic properties. A modified reaction set for the H/O system is determined. A11 of the third body efficiencies M are taken to be unity. Good agreement was obtained between the thermodynamic method and the experimental shock tube data. In addition, the hydrogen bromide experimental data presented in previous work is recomputed with M's of unity.

Marek, C. John↗
DOE OSTI2015

Efficient approach to obtain free energy gradient using QM/MM MD simulation

The efficient computational approach denoted as charge and atom dipole response kernel (CDRK) model to consider polarization effects of the quantum mechanical (QM) region is described using the charge response and the atom dipole response kernels for free energy gradient (FEG) calculations in the quantum mechanical/molecular mechanical (QM/MM) method. CDRK model can reasonably reproduce energies and also energy gradients of QM and MM atoms obtained by expensive QM/MM calculations in a drastically reduced computational time. This model is applied on the acylation reaction in hydrated trypsin-BPTI complex to optimize the reaction path on the free energy surface by means of FEG and the nudged elastic band (NEB) method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Tool employing homotopy-based approaches in finding the controlling unstable equilibrium point in the electric power grid

The analysis tool employs a computer-implemented algorithm that uses homotopy-based approaches to map the solution from the exit point to the controlling unstable equilibrium point (UEP). The computational time is reduced significantly by using an approximate exit point rather than computing an accurate exit point as it is required in finding the controlling UEPs using traditional transient stability direct methods. In addition, this method eliminates the necessity of computing the minimum gradient point (MGP) which is a key element in using Newton methods. These properties provide an advantage to homotopy-based approaches over traditional iterative methods in terms of both speed of computation and reliability of finding solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Gradient-Based Optimization of Wind Farms with Different Turbine Heights: Preprint

Turbine wakes reduce power production in a wind farm. Current wind farms are generally built with turbines that are all the same height, but if wind farms included turbines with different tower heights, the cost of energy (COE) may be reduced. We used gradient-based optimization to demonstrate a method to optimize wind farms with varied hub heights. Our study includes a modified version of the FLORIS wake model that accommodates three-dimensional wakes integrated with a tower structural model. Our purpose was to design a process to minimize the COE of a wind farm through layout optimization and varying turbine hub heights. Results indicate that when a farm is optimized for layout and height with two separate height groups, COE can be lowered by as much as 5%-9%, compared to a similar layout and height optimization where all the towers are the same. The COE has the best improvement in farms with high turbine density and a low wind shear exponent.

17 WIND ENERGY↗
NASA NTRS1975

Airship stresses due to vertical velocity gradients and atmospheric turbulence

Munk's potential flow method is used to calculate the resultant moment experienced by an ellipsoidal airship. This method is first used to calculate the moment arising from basic maneuvers considered by early designers, and then expended to calculate the moment arising from vertical velocity gradients and atmospheric turbulence. This resultant moment must be neutralized by the transverse force of the fins. The results show that vertical velocity gradients at a height of 6000 feet in thunderstorms produce a resultant moment approximately three to four times greater than the moment produced in still air by realistic values of pitch angle or steady turning. Realistic values of atmospheric turbulence produce a moment which is significantly less than the moment produced by maneuvers in still air.

Sheldon, D.↗
DOE OSTI2020

DPM: A deep learning PDE augmentation method with application to large-eddy simulation

A framework is introduced that leverages known physics to reduce overfitting in machine learning for scientific applications. The partial differential equation (PDE) that expresses the physics is augmented with a neural network that uses available data to learn a description of the corresponding unknown or unrepresented physics. Training within this combined system corrects for missing, unknown, or erroneously represented physics, including discretization errors associated with the PDE's numerical solution. For optimization of the network within the PDE, an adjoint PDE is solved to provide high-dimensional gradients, and a stochastic adjoint method (SAM) further accelerates training. Additionally, the approach is demonstrated for large-eddy simulation (LES) of turbulence. High-fidelity direct numerical simulations (DNS) of decaying isotropic turbulence provide the training data used to learn sub-filter-scale closures for the filtered Navier–Stokes equations. Out-of-sample comparisons show that the deep learning PDE method outperforms widely-used models, even for filter sizes so large that they become qualitatively incorrect. It also significantly outperforms the same neural network when a priori trained based on simple data mismatch, not accounting for the full PDE. Measures of discretization errors, which are well-known to be consequential in LES, point to the importance of the unified training formulation's design, which without modification corrects for them. For comparable accuracy, simulation runtime is significantly reduced. A relaxation of the typical discrete enforcement of the divergence-free constraint in the solver is also successful, instead allowing the DPM to approximately enforce incompressibility physics. Since the training loss function is not restricted to correspond directly to the closure to be learned, training can incorporate diverse data, including experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Pseudo-Z symmetric space-times

In this paper, we investigate Pseudo-Z symmetric space-time manifolds. First, we deal with elementary properties showing that the associated form A{sub k} is closed: in the case the Ricci tensor results to be Weyl compatible. This notion was recently introduced by one of the present authors. The consequences of the Weyl compatibility on the magnetic part of the Weyl tensor are pointed out. This determines the Petrov types of such space times. Finally, we investigate some interesting properties of (PZS){sub 4} space-time; in particular, we take into consideration perfect fluid and scalar field space-time, and interesting properties are pointed out, including the Petrov classification. In the case of scalar field space-time, it is shown that the scalar field satisfies a generalized eikonal equation. Further, it is shown that the integral curves of the gradient field are geodesics. A classical method to find a general integral is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Scalable Simulation of Pressure Gradient-Driven Transport of Rarefied Gases in Complex Permeable Media Using Lattice Boltzmann Method

Accurate representations of slip and transitional flow regimes present a challenge in the simulation of rarefied gas flow in confined systems with complex geometries. In these regimes, continuum-based formulations may not capture the physics correctly. This work considers a regularized multi-relaxation time lattice Boltzmann (LB) method with mixed Maxwellian diffusive and halfway bounce-back wall boundary treatments to capture flow at high Kn. The simulation results are validated against atomistic simulation results from the literature. We examine the convergence behavior of LB for confined systems as a function of inlet and outlet treatments, complexity of the geometry, and magnitude of pressure gradient and show that convergence is sensitive to all three. The inlet and outlet boundary treatments considered in this work include periodic, pressure, and a generalized periodic boundary condition. Compared to periodic and pressure treatments, simulations of complex domains using a generalized boundary treatment conserve mass but require more iterations to converge. Convergence behavior in complex domains improves at higher magnitudes of pressure gradient across the computational domain, and lowering the porosity deteriorates the convergence behavior for complex domains.

42 ENGINEERING↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2022

General-Purpose Bayesian Tensor Learning With Automatic Rank Determination and Uncertainty Quantification

A major challenge in many machine learning tasks is that the model expressive power depends on model size. Low-rank tensor methods are an efficient tool for handling the curse of dimensionality in many large-scale machine learning models. The major challenges in training a tensor learning model include how to process the high-volume data, how to determine the tensor rank automatically, and how to estimate the uncertainty of the results. While existing tensor learning focuses on a specific task, this paper proposes a generic Bayesian framework that can be employed to solve a broad class of tensor learning problems such as tensor completion, tensor regression, and tensorized neural networks. We develop a low-rank tensor prior for automatic rank determination in nonlinear problems. Our method is implemented with both stochastic gradient Hamiltonian Monte Carlo (SGHMC) and Stein Variational Gradient Descent (SVGD). We compare the automatic rank determination and uncertainty quantification of these two solvers. We demonstrate that our proposed method can determine the tensor rank automatically and can quantify the uncertainty of the obtained results. We validate our framework on tensor completion tasks and tensorized neural network training tasks.

Bayesian inference↗
NASA NTRS1983

Vertical gradients of sunspot magnetic fields

The results of a Solar Maximum Mission (SMM) guest investigation to determine the vertical gradients of sunspot magnetic fields for the first time from coordinated observations of photospheric and transition-region fields are described. Descriptions are given of both the photospheric vector field of a sunspot, derived from observations using the NASA Marshall Space Flight Center vector magnetograph, and of the line-of-sight component in the transition region, obtained from the SMM Ultraviolet Spectrometer and Polarimeter instrument. On the basis of these data, vertical gradients of the line-of-sight magnetic field component are calculated using three methods. It is found that the vertical gradient of Bz is lower than values from previous studies and that the transition-region field occurs at a height of approximately 4000-6000 km above the photosphere.

Hagyard, M. J.↗
DOE OSTI2026

Renormalized quark masses using gradient flow

We propose a new and simple method for determining the renormalized quark masses from lattice simulations. Renormalized quark masses are an important input to many phenomenological applications, including searching and modeling physics beyond the Standard Model. The nonperturbative renormalization is performed using gradient flow combined with the short-flow-time expansion that is improved by renormalization group (RG) running to match to the $\overline{MS}$ scheme. Implementing the RG running perturbatively, we demonstrate this method works reliably at least up to the charm-quark mass and exhibits an easily attainable “windowing condition.” Using RBC/UKQCD’s (2+1)-flavor Shamir domain-wall fermion ensembles with Iwasaki gauge action, we find 𝑚 $\overline{MS}$ 𝑠⁡ (𝜇 = 2 GeV) = 89⁢(3) MeV and 𝑚 $\overline{MS}$ 𝑐 ⁡(𝜇 = 3 GeV) = 972⁢(16) MeV. These results predict the scale-independent ratio 𝑚 𝑐 /𝑚 𝑠 = 12.1⁢(4). Generalization to other observables is possible, providing an efficient approach to determine nonperturbatively renormalized fermionic observables like form factors or bag parameters from lattice simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Stability of compressible reduced magnetohydrodynamic equilibria-Analogy with magnetorotational instability

Stability analyses for equilibria of the compressible reduced magnetohydrodynamics (CRMHD) model are carried out by means of the Energy-Casimir (EC) method. Stability results are compared with those obtained for ideal magnetohydrodynamics (MHD) from the classical {delta}W criterion. An identification of the terms in the second variation of the free energy functional for CRMHD with those of {delta}W is made: two destabilizing effects present for CRMHD turn out to correspond to the kink and interchange instabilities in usual MHD, while the stabilizing roles of field line bending and compressibility are also identified in the reduced model. Also, using the EC method, stability conditions in the presence of toroidal flow are obtained. A formal analogy between CRMHD and a reduced incompressible model for magnetized rotating disks, due to Julien and Knobloch [EAS Pub. Series, 21, 81 (2006)], is discovered. In light of this analogy, energy stability analysis shows that the condition for magnetorotational instability (MRI) for the latter model corresponds to the condition for interchange instability in CRMHD, with the Coriolis term and shear velocity playing the roles of the curvature term and pressure gradient, respectively. Using the EC method, stability conditions for the rotating disk model, for a large class of equilibria with possible non-uniform magnetic fields, are obtained. In particular, this shows it is possible for the MRI system to undergo, in addition to the MRI, another instability that is analogous to the kink instability. For vanishing magnetic field, the Rayleigh hydrodynamical stability condition is recovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Peak-Seeking Control Using Gradient and Hessian Estimates

A peak-seeking control method is presented which utilizes a linear time-varying Kalman filter. Performance function coordinate and magnitude measurements are used by the Kalman filter to estimate the gradient and Hessian of the performance function. The gradient and Hessian are used to command the system toward a local extremum. The method is naturally applied to multiple-input multiple-output systems. Applications of this technique to a single-input single-output example and a two-input one-output example are presented.

Ryan, John J.↗