Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative 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 541 records · Page 30

The equilibrium dayside magnetosphere

A method is presented of computing the dayside global earth magnetic field which is in equilibrium with the plasma pressure, based on satellite observations at a local region of the magnetosphere. The method, which utilizes a perturbation around a dipole magnetic field, involves computation of the global plasma pressure profile based on the equatorial (anisotropic) pressure data, derivation of the current profile which satisfies the equilibrium condition, and computation of the magnetic field using the current profile and the boundary current produced by the solar wind. The method is applied for the Active Magnetospheric Particle Tracer Explorers data, and the result of the computation is found to compare reasonably well with the observed magnetic field profile near the geomagnetic equator.

Zavriyev, Anton↗

Velocity perturbation distributions in the breakup of artificial satellites

A method is presented for calculating the three orthogonal components of the velocity perturbations of satellite fragments, with a view to ascertaining the nature and intensity of the satellite breakup. The method employs three simultaneous equations furnished by changes in fragment specific energy, specific angular momentum, and plane orientation. Velocity perturbations are thereby calculated for fragments from 20 major satellite breakup events; these results, in conjunction with a technique for determining fragment masses, yield a description of the breakup process.

Badhwar, Gautam D.↗

Sensitivity of magnetic islands in permanent magnet stellarators using the gradient and Hessian methods

Stellarator plasmas are known to be very sensitive to perturbations in the magnetic field. The permanent magnet stellarator was in part developed as a solution to high machining tolerances placed on the shape properties of electromagnetic coils in traditional stellarators. However, as a consequence of this high sensitivity to the field structure, sensitivities of permanent magnet stellarator plasmas to perturbations of permanent magnet properties must necessarily be well-understood. The gradient and Hessian matrix methods have been previously demonstrated to be useful sensitivity analysis methods for modular coils. We apply these two methods to the study of island width sensitivities in both the MUSE and PM4STELL permanent magnet stellarator projects. These sensitivity methods were used to determine the relative impacts of permanent magnet parameter perturbations on island widths in the vacuum field approximation of both stellarator equilibria. The square of resonant magnetic field perturbation is used here as a proxy for island width. In particular, gradients of magnetizations of individual magnets were examined in MUSE, as well as gradients of magnet group displacements informed by device design. Three different forms of permanent magnet magnetization perturbations are investigated for MUSE, and the flux surface response to perturbations is demonstrated. The Hessian matrix method is applied to PM4STELL, illustrating the sensitivity of dominant island widths to displacements of toroidal wedge structures. These methods allow for selective direction of experimental resources toward regions of heightened sensitivity, while constraints on less impactful permanent magnet parameters can be relaxed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics-informed Estimation of the Covariance Matrix for Various Neutron Spectra

A method for estimating covariance matrices which capture the uncertainties in calculated reactor spectra has been developed. This method is based on perturbing the parameters of a physics-based analytic model fitted to a calculated spectrum. The covariance of the perturbed analytic spectra imposes energy-dependent correlations due to the physics of the neutron processes in the reactor, i.e., a fission component, a 1/E down-scatting component, and a thermal Maxwellian component. An analytic model is developed which is shown to produce good fits to several reactor environments. The covariance matrices produced via this method are then used as the prior spectrum in STAYSL least squares spectrum adjustment where it is combined with integral metrics, such as activation measurements, to produce a high-fidelity neutron spectrum characterization. It was concluded that the methodology showed agreeable results for the ACRR free-field spectrum adjustment in STAYSL resulting in a 𝜒 2 value of 2.21 (per degree of freedom), but further work is needed to describe scattering and interface regions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Detailed comparison of renormalization scale-setting procedures based on the principle of maximum conformality

The Principle of Maximum Conformality (PMC), which generalizes the conventional Gell-Mann-Low method for scale-setting in perturbative QED to non-Abelian QCD, provides a rigorous method for achieving unambiguous scheme-independent, fixed-order predictions for physical observables consistent with the principles of the renormalization group. In addition to the original multi-scale-setting approach (PMCm), two variations of the PMC have been proposed to deal with ambiguities associated with the uncalculated higher order terms in the pQCD series, i.e. the single-scale-setting approach (PMCs) and the procedures based on ``intrinsic conformality" (PMC ∞ ). In this paper, we will give a detailed comparison of these PMC approaches by comparing their predictions for three important quantities R e+e– , R τ , and $Γ(H→b\bar{b}$) up to four-loop pQCD corrections. The PMCs approach determines an overall effective running coupling α s (Q) by the recursive use of the renormalization group equation, whose argument Q represents the actual momentum flow of the process. Our numerical results show that the PMCs method, which involves a somewhat simpler analysis, can serve as a reliable substitute for the full multi-scale PMCm method, and that it leads to more precise pQCD predictions with small residual scale dependence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Recent Developments in the General Atomic and Molecular Electronic Structure System

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods like the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite like GAMESS is maintained and developed is considered. Future developments are briefly summarized.

74 ATOMIC AND MOLECULAR PHYSICS↗

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↗

Solution of the two-dimensional Euler equations with generalized coordinate transformation using flux vector splitting

An implicit finite difference code using flux vector splitting has been developed for solving the two-dimensional inviscid gas dynamics equations. The method is spatially second-order acurate, fully conservative, and uses body-conforming generalized coordinates for treating complex geometries. Numerical results have been obtained for transonic flow over a circular cylinder and airfoils. Steady results for a half cylinder (top and bottom symmetry-imposed) range from critical flow to a strong shock case with rotationally induced flow separation. Full cylinder solutions at freestream Mach number values of 0.5, however, show unsteady oscillation. A perturbation form of the method has also been developed and used to compute both fore and aft inviscid flow separation about a cylinder for a nonuniform incoming stream.

Buning, P. G.↗

Solutions of the motion of synchronous satellites with arbitrary eccentricity and inclination

A first order, semianalytical theory for the long term motion of resonant satellites is presented. The theory is valid for all eccentricities and inclinations and for all commensurability ratios. The method allows the inclusion of all the zonal and tesseral harmonics as well as luni solar perturbations and radiation pressure. The method is applied to a synchronous satellite including only the J sub 2 and J sub 22 harmonics. Global, long term solutions for this problem, eccentricity, argument of perigee, and inclination are obtained.

Nacozy, P. E.↗

Boundary layer transition studies

A small-scale wind tunnel previously used for turbulent boundary layer experiments was modified for two sets of boundary layer transition studies. The first study concerns a laminar separation/turbulent reattachment. The pressure gradient and unit Reynolds number are the same as the fully turbulent flow of Spalart and Watmuff. Without the trip wire, a laminar layer asymptotes to a Falkner & Skan similarity solution in the FPG. Application of the APG causes the layer to separate and a highly turbulent and approximately 2D mean flow reattachment occurs downstream. In an effort to gain some physical insight into the flow processes a small impulsive disturbance was introduced at the C(sub p) minimum. The facility is totally automated and phase-averaged data are measured on a point-by-point basis using unprecedently large grids. The evolution of the disturbance has been tracked all the way into the reattachment region and beyond into the fully turbulent boundary layer. At first, the amplitude decays exponentially with streamwise distance in the APG region, where the layer remains attached, i.e. the layer is viscously stable. After separation, the rate of decay slows, and a point of minimum amplitude is reached where the contours of the wave packet exhibit dispersive characteristics. From this point, exponential growth of the amplitude of the disturbance is observed in the detached shear layer, i.e. the dominant instability mechanism is inviscid. A group of large-scale 3D vortex loops emerges in the vicinity of the reattachment. Remarkably, the second loop retains its identify far downstream in the turbulent boundary layer. The results provide a level of detail usually associated with CFD. Substantial modifications were made to the facility for the second study concerning disturbances generated by Suction Holes for laminar flow Control (LFC). The test section incorporates suction through interchangeable porous test surfaces. Detailed studies have been made using isolated holes in the impervious test plate that used to establish the Blasius base flow. The suction is perturbed harmonically and data are averaged on the basis of the phase of the disturbance, for conditions corresponding to strong suction and without suction. The technique was enhanced by using up to nine multiple probes to reduce the experimental run-time. In both cases, 3D contour surfaces in the vicinity of the hole show highly 3D TS waves which fan out in the spanwise direction forming bow-shaped waves downstream. The case without suction has proved useful for evaluating calculation methods. With suction, the perturbations on the centerline are much stronger and decay less rapidly, while the TS waves in the far field are similar to the case without suction. Downstream, the contour surfaces of the TS waves develop spanwise irregularities which eventually form into clumps. The spanwise clumping is evidence of a secondary instability that could be associated with suction vortices. Designers of porous surfaces use Goldsmith's Criterion to minimize cross-stream interactions. It is shown that partial TS wave cancellation is possible, depending on the hole spacing, disturbance frequency and free-stream velocity. New high-performance Constant Temperature Hot-Wire Anemometers were designed and built, based on a linear system theory analysis that can be extended to arbitrary order. The motivation was to achieve the highest possible frequency reponse while ensuring overall system stability. The performance is equal to or superior to commercially available instruments at about 10% of the cost. Details, such as fabrication drawings and a parts list, have been published to enable the instrument to be construced by others.

Watmuff, Jonathan H.↗

Computing Geopotential Perturbations

Improved formulation reduces error and number of computational steps. Report describes method for calculating partial derivative of geopotential perturbation with respect to time. Useful in computations of orbit of lightweight satellite in which total orbital energy or related parameter one of dependent variables.

Bond, Victor R.↗

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗

Control strategies

Techniques for improving the performance of robotic-manipulator path-control systems comprising independent SISO feedback controllers for each joint are discussed and illustrated with block diagrams, reviewing the results of recent analytical investigations. Topics examined include the servo design for a single link, the equations of motion for manipulators, SISO servo design for multiple links, inverse methods, pole placement with compensation of the gravity terms, linear state-feedback control based on the perturbation equations, and adaptive control methods. Consideration is given to variable-structure systems, suboptimal controllers, and the optimal-design-strategy approach.

Wang, J. C.↗

Aerodynamic sensitivities from subsonic, sonic and supersonic unsteady, nonplanar lifting-surface theory

The technique of implicit differentiation has been used in combination with linearized lifting-surface theory to derive analytical expressions for aerodynamic sensitivities (i.e., rates of change of lifting pressures with respect to general changes in aircraft geometry, including planform variations) for steady or oscillating planar or nonplanar lifting surfaces in subsonic, sonic, or supersonic flow. The geometric perturbation is defined in terms of a single variable, and the user need only provide simple expressions or similar means for defining the continuous or discontinuous global or local perturbation of interest. Example expressions are given for perturbations of the sweep, taper, and aspect ratio of a wing with trapezoidal semispan planform. In addition to direct computational use, the analytical method presented here should provide benchmark criteria for assessing the accuracy of aerodynamic sensitivities obtained by approximate methods such as finite geometry perturbation and differencing. The present process appears to be readily adaptable to more general surface-panel methods.

Yates, E. Carson, Jr.↗

Using GPUs and the Parameterization Method for Rapid Search and Refinement of Connections between Tori in Periodically Perturbed Planar Circular Restricted 3-Body Problems

When the planar circular restricted 3-body problem is periodically perturbed, most unstable periodic orbits become invariant tori. However, 2D Poincare ́ sections no longer work to find their manifolds’ intersections; new methods are needed. In this study, we first review a method of restricting the intersection search to only certain manifold subsets. We then implement this search using Julia and OpenCL, representing the manifolds as triangular meshes and gaining a 30x speedup using GPUs. We finally show how to use manifold parametrizations to refine the ap- proximate connections found in the mesh search. We demonstrate the tools on the planar elliptic RTBP.

de la Llave, Rafael↗

Microstrata

Abstract Microstrata are the non-extremal analogues of superstrata: they are smooth, non-extremal (non-BPS) solitonic solutions to IIB supergravity whose deep-throat limits approximate black holes. Using perturbation theory and numerical methods, we construct families of solutions using a consistent truncation to three-dimensional supergravity. The most general families presented here involve two continuous parameters, or amplitudes, and four quantized parameters that set the angular momenta and energy levels. Our solutions are asymptotic to the vacuum of the D1-D5 system: AdS 3 ×S 3 × 𝕋 4 . Using holography, we show that the they are dual to multi-particle states in the D1-D5 CFT involving a large number of mutually non-BPS supergravitons and we determine the anomalous dimensions of these states from the binding energies in supergravity. These binding energies are uniformly negative and depend non-linearly on the amplitudes of the states. In one family of solutions, smoothness restricts some of the fields to lie on a special locus of the parameter space. Using precision holography we show that this special locus can be identified with the multi-particle states constructed via the standard OPE of the single-particle constituents. Our numerical analysis shows that microstrata are robust at large amplitudes and the solutions can be obtained to very high precision.

Physics↗

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗

Optimization of nonlinear turbulence in stellarators

We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC . Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasi-symmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗