Engineering PapersSearch

SEARCH · Engineering Papers

Results for “VARIATION METHOD”

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 127 records · Page 7

Design of a Variational Multiscale Method for Turbulent Compressible Flows

A spectral-element framework is presented for the simulation of subsonic compressible high-Reynolds-number flows. The focus of the work is maximizing the efficiency of the computational schemes to enable unsteady simulations with a large number of spatial and temporal degrees of freedom. A collocation scheme is combined with optimized computational kernels to provide a residual evaluation with computational cost independent of order of accuracy up to 16th order. The optimized residual routines are used to develop a low-memory implicit scheme based on a matrix-free Newton-Krylov method. A preconditioner based on the finite-difference diagonalized ADI scheme is developed which maintains the low memory of the matrix-free implicit solver, while providing improved convergence properties. Emphasis on low memory usage throughout the solver development is leveraged to implement a coupled space-time DG solver which may offer further efficiency gains through adaptivity in both space and time.

Design

An historical survey of computational methods in optimal control.

Review of some of the salient theoretical developments in the specific area of optimal control algorithms. The first algorithms for optimal control were aimed at unconstrained problems and were derived by using first- and second-variation methods of the calculus of variations. These methods have subsequently been recognized as gradient, Newton-Raphson, or Gauss-Newton methods in function space. A much more recent addition to the arsenal of unconstrained optimal control algorithms are several variations of conjugate-gradient methods. At first, constrained optimal control problems could only be solved by exterior penalty function methods. Later algorithms specifically designed for constrained problems have appeared. Among these are methods for solving the unconstrained linear quadratic regulator problem, as well as certain constrained minimum-time and minimum-energy problems. Differential-dynamic programming was developed from dynamic programming considerations. The conditional-gradient method, the gradient-projection method, and a couple of feasible directions methods were obtained as extensions or adaptations of related algorithms for finite-dimensional problems. Finally, the so-called epsilon-methods combine the Ritz method with penalty function techniques.

Polak, E.

Coupling finite and boundary element methods for two-dimensional potential problems

A finite element-boundary element (FE-BE) coupling method based on a weighted residual variational method is presented for potential problems, governed by either the Laplace or the Poisson equations. In this method, a portion of the domain of interest is modeled by finite elements (FE) and the remainder of the region by boundary elements (BE). Because the BE fundamental solutions are valid for infinite domains, a procedure that limits the effects of the BE fundamental solution to a small region adjacent to the FE region, called the transition region (TR), is developed. This procedure involves a judicious choice of functions called the transition (T) functions that have unit values on the BE-TR interface and zero values on the FE-TR interface. The present FE-BE coupling algorithm is shown to be independent of the extent of the transition region and the choice of the transition functions. Therefore, transition regions that extend to only one layer of elements between FE and BE regions and the use of simple linear transition functions work well.

Krishnamurthy, T.

Three-dimensional inversion of satellite observed radiances A proposal

Fully three-dimensional temperature fields are obtained from observed satellite radiances through the use of variational methods which generalize the method of Wahba and Wendelberger (1980). This variational formalism allows a variety of different information types to contribute to the inversion process. Full use can be made of the horizontal smoothness of the temperature field in the calculation of three-dimensional retrievals. The present method facilitates the acquisition of maximum vertical resolution, accomplishes cloud effect filtering, and minimizes the effects of instrument noise, to yield horizontally coherent retrievals.

Hoffman, R. N.

A Multifunctional Interface Method for Coupling Finite Element and Finite Difference Methods: Two-Dimensional Scalar-Field Problems

A multifunctional interface method with capabilities for variable-fidelity modeling and multiple method analysis is presented. The methodology provides an effective capability by which domains with diverse idealizations can be modeled independently to exploit the advantages of one approach over another. The multifunctional method is used to couple independently discretized subdomains, and it is used to couple the finite element and the finite difference methods. The method is based on a weighted residual variational method and is presented for two-dimensional scalar-field problems. A verification test problem and a benchmark application are presented, and the computational implications are discussed.

Ransom, Jonathan B.

Variational neural network approach to QFT in the field basis

We present a variational neural network approach for solving quantum field theories in the field basis, focusing on the free Klein-Gordon model formulated in momentum space. While recent studies have explored neural-network-based variational methods for scalar field theory in position space, a systematic benchmark of the analytically solvable Klein-Gordon ground state—particularly in the momentum-space field basis—has been lacking. In this work, we represent the ground-state wavefunctional as a neural network defined on a discretized set of field configurations and train it by minimizing the Hamiltonian expectation value. This framework enables direct comparison to exact analytic results for a range of key observables, including the ground-state energy, two-point correlators, expectation value of the field, and the structure of the learned wavefunctional itself. Our results provide quantitative diagnostics of accuracy and establish a validated foundation for extending neural-network wavefunctional methods to interacting field theories and position-space formulations.

Klein-Gordon model

Constraints on the finite volume two-nucleon spectrum at 𝑚𝜋 ≈806 MeV

The low-energy, finite-volume spectrum of the two-nucleon system at a quark mass corresponding to a pion mass of 𝑚𝜋≈806 MeV is studied with lattice quantum chromodynamics (LQCD) using variational methods. The interpolating-operator sets used in [Variational study of two-nucleon systems with lattice QCD, Phys. Rev. D 107, 094508 (2023).] are extended by including a complete basis of local hexaquark operators, as well as plane-wave dibaryon operators built from products of both positive- and negative-parity nucleon operators. Results are presented for the isosinglet and isotriplet two-nucleon channels. In both channels, noticeably weaker variational bounds on the lowest few energy eigenvalues are obtained from operator sets which contain only hexaquark operators or operators constructed from the product of two negative-parity nucleons, while other operator sets produce low-energy variational bounds which are consistent within statistical uncertainties. The consequences of these studies for the LQCD understanding of the two-nucleon spectrum are investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Absence of diurnal variation of C-reactive protein concentrations in healthy human subjects

BACKGROUND: The concentration of C-reactive protein (CRP) in otherwise healthy subjects has been shown to predict future risk of myocardial infarction and stroke. CRP is synthesized by the liver in response to interleukin-6, the serum concentration of which is subject to diurnal variation. METHODS: To examine the existence of a time-of-day effect for baseline CRP values, we determined CRP concentrations in hourly blood samples drawn from healthy subjects (10 males, 3 females; age range, 21-35 years) during a baseline day in a controlled environment (8 h of nighttime sleep). RESULTS: Overall CRP concentrations were low, with only three subjects having CRP concentrations >2 mg/L. Comparison of raw data showed stability of CRP concentrations throughout the 24 h studied. When compared with cutoff values of CRP quintile derived from population-based studies, misclassification of greater than one quintile did not occur as a result of diurnal variation in any of the subjects studied. Nonparametric ANOVA comparing different time points showed no significant differences for both raw and z-transformed data. Analysis for rhythmic diurnal variation using a method fitting a cosine curve to the group data was negative. CONCLUSIONS: Our data show that baseline CRP concentrations are not subject to time-of-day variation and thus help to explain why CRP concentrations are a better predictor of vascular risk than interleukin-6. Determination of CRP for cardiovascular risk prediction may be performed without concern for diurnal variation.

Non-NASA Center

A design method for minimizing sensitivity to plant parameter variations

A method is described for minimizing the sensitivity of multivariable systems to parameter variations. The variable parameters are considered as random variables and their effect is included in a quadratic performance index. The performance index is a weighted sum of the state and control covariances that stem from both the random system disturbances and the parameter uncertainties. The numerical solution of the problem is described and application of the method to several initially sensitive tracking systems is discussed. The sensitivity factor of reduction was typically 2 or 3 over a system based on random system noise only, and yet resulted in state RMS increases of only about a factor of two.

Hadass, Z.

Partial spline models for the inclusion of tropopause and frontal boundary information in otherwise smooth two- and three-dimensional objective analysis

A new method, based on partial spline models, is developed for including specified discontinuities in otherwise smooth two- and three-dimensional objective analyses. The method is appropriate for including tropopause height information in two- and three-dimensinal temperature analyses, using the O'Sullivan-Wahba physical variational method for analysis of satellite radiance data, and may in principle be used in a combined variational analysis of observed, forecast, and climate information. A numerical method for its implementation is described and a prototype two-dimensional analysis based on simulated radiosonde and tropopause height data is shown. The method may also be appropriate for other geophysical problems, such as modeling the ocean thermocline, fronts, discontinuities, etc.

Shiau, Jyh-Jen