SEARCH · Engineering Papers

Results for “asymptotic analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Analysis of the interaction of a weak normal shock wave with a turbulent boundary layer

The method of matched asymptotic expansions is used to analyze the interaction of a normal shock wave with an unseparated turbulent boundary layer on a flat surface at transonic speeds. The theory leads to a three-layer description of the interaction in the double limit of Reynolds number approaching infinity and Mach number approaching unity. The interaction involves an outer, inviscid rotational layer, a constant shear-stress wall layer, and a blending region between them. The pressure distribution is obtained from a numerical solution of the outer-layer equations by a mixed-flow relaxation procedure. An analytic solution for the skin friction is determined from the inner-layer equations. The significance of the mathematical model is discussed with reference to existing experimental data.

Melnik, R. E.↗

DOE OSTI2022

Bohm Criterion of Plasma Sheaths away from Asymptotic Limits

The plasma exit flow speed at the sheath entrance is constrained by the Bohm criterion. The so-called Bohm speed regulates the plasma particle and power exhaust fluxes to the wall, and it is commonly deployed as a boundary condition to exclude the sheath region in quasineutral plasma modeling. In this work, the Bohm criterion analysis is performed in the intermediate plasma regime away from the previously known limiting cases of adiabatic laws and the asymptotic limit of infinitesimal Debye length in a finite-size system, using the transport equations of an anisotropic plasma. The resulting Bohm speed has explicit dependence on local plasma heat flux, temperature isotropization, and thermal force. Comparison with kinetic simulations demonstrates its accuracy over the plasma-sheath transition region in which quasineutrality is weakly perturbed and the Bohm criterion applies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1994

Multigrid one shot methods for optimal control problems: Infinite dimensional control

The multigrid one shot method for optimal control problems, governed by elliptic systems, is introduced for the infinite dimensional control space. ln this case, the control variable is a function whose discrete representation involves_an increasing number of variables with grid refinement. The minimization algorithm uses Lagrange multipliers to calculate sensitivity gradients. A preconditioned gradient descent algorithm is accelerated by a set of coarse grids. It optimizes for different scales in the representation of the control variable on different discretization levels. An analysis which reduces the problem to the boundary is introduced. It is used to approximate the two level asymptotic convergence rate, to determine the amplitude of the minimization steps, and the choice of a high pass filter to be used when necessary. The effectiveness of the method is demonstrated on a series of test problems. The new method enables the solutions of optimal control problems at the same cost of solving the corresponding analysis problems just a few times.

Arian, Eyal↗

DOE OSTI2021

Theoretical study of intrinsic defects in cubic silicon carbide 3C-SiC

Using the local moment counter charge (LMCC) method to accurately represent the asymptotic electrostatic boundary conditions within density functional theory supercell calculations, we present a comprehensive analysis of the atomic structure and energy levels of point defects in cubic silicon carbide (3C-SiC). Finding that the classical long-range dielectric screening outside the supercell induced by a charged defect is a significant contributor to the total energy. we describe and validate a modified Jost screening model to evaluate this polarization energy. This leads to bulk-converged defect levels in finite size supercells. With the LMCC boundary conditions and a standard Perdew-Burke-Ernzerhof (PBE) exchange correlation functional, the computed defect level spectrum exhibits no band gap problem: the range of defect levels spans ~2.4 eV, an effective defect band gap that agrees with the experimental band gap. Comparing with previous literature, our LMCC-PBE defect results are in consistent agreement with the hybrid-exchange functional results of Oda et al. [J. Chem. Phys. 139, 124707 (2013)] rather than their PBE results. The difference with their PBE results is attributed to their use of a conventional jellium approximation rather than the more rigorous LMCC approach for handling charged supercell boundary conditions. The difference between standard dft and hybrid functional results for defect levels lies not in a band gap problem but rather in solving a boundary condition problem. The LMCC-PBE entirely mitigates the effect of the band gap problem on defect levels. The more computationally economical PBE enables a systematic exploration of 3C-SiC defects, where, most notably, we find that the silicon vacancy undergoes Jahn-Teller-induced distortions from the previously assumed T d symmetry, and that the divacancy, like the silicon vacancy, exhibits a site-shift bistability in p -type conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Preparation of excited states for nuclear dynamics on a quantum computer

We study two different methods to prepare excited states on a quantum computer, a key initial step to study dynamics within linear response theory. The first method uses unitary evolution for a short time T = O(√1 - F) to approximate the action of an excitation operator Ô with fidelity F and success probability P ≈ 1 – F. The second method probabilistically applies the excitation operator using the Linear Combination of Unitaries (LCU) algorithm. We benchmark these techniques on emulated and real quantum devices, using a toy model for thermal neutron-proton capture. Despite its larger memory footprint, the LCU-based method is efficient even on current generation noisy devices and can be implemented at a lower gate cost than a naive analysis would suggest. Here, these findings show that quantum techniques designed to achieve good asymptotic scaling on fault tolerant quantum devices might also provide practical benefits on devices with limited connectivity and gate fidelity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1967

A Two-Variable Asymptotic Solution for Three - Dimensional, Solar-Powered, Low-Thrust Trajectories in the Vicinity of the Ecliptic Plane

Two-variable asymptotic solution for three-dimensional solar powered low thrust trajectories in vicinity of ecliptic plane

ELECTRIC PROPULSION↗

NASA NTRS1977

Theory of viscous transonic flow over airfoils at high Reynolds number

This paper considers viscous flows with unseparated turbulent boundary layers over two-dimensional airfoils at transonic speeds. Conventional theoretical methods are based on boundary layer formulations which do not account for the effect of the curved wake and static pressure variations across the boundary layer in the trailing edge region. In this investigation an extended viscous theory is developed that accounts for both effects. The theory is based on a rational analysis of the strong turbulent interaction at airfoil trailing edges. The method of matched asymptotic expansions is employed to develop formal series solutions of the full Reynolds equations in the limit of Reynolds numbers tending to infinity. Procedures are developed for combining the local trailing edge solution with numerical methods for solving the full potential flow and boundary layer equations. Theoretical results indicate that conventional boundary layer methods account for only about 50% of the viscous effect on lift, the remaining contribution arising from wake curvature and normal pressure gradient effects.

Melnik, R. E.↗

NASA NTRS1985

Vortex-dominated flow with viscous core structure

Recent theoretical studies of vortex-dominated flows are reviewed with special emphasis on those for which the viscous core structures play an important role. The problems to be described are: The interaction and merging of two-dimensional vortices and of curved vortex filaments, the roll-up and decay of trailing far wakes, and the initiation of vortex breakdown. The analysis utilizes finite-difference solutions of the Navier-Stokes equations complemented by asymptotic expansion techniques.

Liu, C. H.↗

NASA NTRS1985

An asymptotic investigation of the stationary modes of instability of the boundary layer on a rotating disc

An investigation of high Reynolds number stationary instabilities in the boundary layer on a rotating disc is given. It is shown that in addition to the inviscid mode at high Reynolds numbers, there is a stationary short wavelength mode. This mode has its structure fixed by a balance between viscous and Coriolis forces and cannot be described by an inviscid theory. The asymptotic structure of the wavenumber and orientation of this mode is obtained. A similar analysis is given for the inviscid mode, the expansion procedure used is capable of taking nonparallel effects into account in a self consistent manner. The results are compared to numerical calculations and experimental observations.

Hall, P.↗

NASA NTRS1986

An asymptotic investigation of the stationary modes of instability of the boundary layer on a rotating disc

Hall, P.↗

NASA NTRS1992

Lyapunov function-based control laws for revolute robot arms - Tracking control, robustness, and adaptive control

A new class of joint level control laws for all-revolute robot arms is introduced. The analysis is similar to a recently proposed energy-like Liapunov function approach, except that the closed-loop potential function is shaped in accordance with the underlying joint space topology. This approach gives way to a much simpler analysis and leads to a new class of control designs which guarantee both global asymptotic stability and local exponential stability. When Coulomb and viscous friction and parameter uncertainty are present as model perturbations, a sliding mode-like modification of the control law results in a robustness-enhancing outer loop. Adaptive control is formulated within the same framework. A linear-in-the-parameters formulation is adopted and globally asymptotically stable adaptive control laws are derived by simply replacing unknown model parameters by their estimates (i.e., certainty equivalence adaptation).

Wen, John T.↗

NASA NTRS2013

Investigation of Control Inceptor Dynamics and Effect on Human Subject Performance

The control inceptor used in a vehicle simulation is an important part of adequately representing the dynamics of the vehicle. The inceptor characteristics are typically based on a second order spring mass damper system with damping, force gradient, breakout force, and natural frequency parameters. Changing these parameters can have a great effect on pilot control of the vehicle. A quasi transfer of training experiment was performed employing a high fidelity and a low fidelity control inceptor. A disturbance compensatory task was employed which involved a simple horizon line disturbed in roll by a sum of sinusoids presented in an out-the-window display. Vehicle dynamics were modeled as 1/s and 1/s2. The task was to maintain level flight. Twenty subjects were divided between the high and the low fidelity training groups. Each group was trained to a performance asymptote, and then transferred to the high fidelity simulation. RMS tracking error, a PSD analysis, and a workload analysis were performed to quantify the transfer of training effect. Quantitative results of the experiments show that there is no significant difference between the high and low fidelity training groups for 1/s plant dynamics. For 1/s2 plant dynamics there is a greater difference in tracking performance and PSD; and the subjects are less correlated with the input disturbance function

Stanco, Anthony A.↗

DOE OSTI2023

Finite-Time Analysis of Whittle Index based Q-Learning for Restless Multi-Armed Bandits with Neural Network Function Approximation

Whittle index policy is a heuristic to the intractable restless multi-armed bandits (RMAB) problem. Although it is provably asymptotically optimal, finding Whittle indices remains difficult. In this paper, we present Neural-Q-Whittle, a Whittle index based Q-learning algorithm for RMAB with neural network function approximation, which is an example of nonlinear two-timescale stochastic approximation with Q-function values updated on a faster timescale and Whittle indices on a slower timescale. Despite the empirical success of deep Q-learning, the non-asymptotic convergence rate of Neural-Q-Whittle, which couples neural networks with two-timescale Q-learning largely remains unclear. This paper provides a finite-time analysis of Neural-Q-Whittle, where data are generated from a Markov chain, and Q-function is approximated by a ReLU neural network. Our analysis leverages a Lyapunov drift approach to capture the evolution of two coupled parameters, and the nonlinearity in value function approximation further requires us to characterize the approximation error. Combing these provide Neural-Q-Whittle with convergence rate, where is the number of iterations.

reinforcement learning, structured learning, conve↗

DOE OSTI2021

Excitation spectra of quantum matter without quasiparticles. I. Sachdev-Ye-Kitaev models

Here, we study the low-frequency spectra of complex Sachdev-Ye-Kitaev models at general densities. The analysis applies also to SU(M) magnets with random exchange at large M. The spectral densities are computed by numerical analysis of the saddle-point equations on the real frequency ω axis at zero temperature T. The asymptotic low-ω behaviors are found to be in excellent agreement with the scaling dimensions of irrelevant operators which perturb the conformally invariant critical states. Of possible experimental interest is our computation of the universal spin spectral weight of the SU(M) magnets at low ω and T; this includes a contribution from the time reparametrization mode, which is the boundary graviton of the holographic dual. This analysis is extended to a random t-J model in the following paper [M. Tikhanovskaya et al., following paper, Phys. Rev. B 103, 075142 (2021)].

36 MATERIALS SCIENCE↗

DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

DOE OSTI2023

Stopping to reflect: Asymptotic static moving mirrors as quantum analogs of classical radiation

Radiation from an accelerating charge is a basic process that can serve as an intersection between classical and quantum physics. We present two exactly soluble electron trajectories that permit analysis of the radiation emitted, exploring its time evolution and spectrum by analogy with the moving mirror model of the dynamic Casimir effect. These classical solutions are finite energy, rectilinear (nonperiodic), asymptotically zero velocity worldlines with corresponding quantum analog beta Bogolyubov coefficients. One of them has an interesting connection to uniform acceleration and Leonardo da Vinci's water pitcher experiment.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1992

Pulsatile instability in rapid directional solidification - Strongly-nonlinear analysis

In the rapid directional solidification of a dilute binary alloy, analysis reveals that, in addition to the cellular mode of Mullins and Sekerka (1964), there is an oscillatory instability. For the model analyzed by Merchant and Davis (1990), the preferred wavenumber is zero; the mode is one of pulsation. Two strongly nonlinear analyses are performed that describe this pulsatile mode. In the first case, nonequilibrium effects that alter solute rejection at the interface are taken asymptotically small. A nonlinear oscillator equation governs the position of the solid-liquid interface at leading order, and amplitude and phase evolution equations are derived for the uniformly pulsating interface. The analysis provides a uniform description of both subcritical and supercritical bifurcation and the transition between the two. In the second case, nonequilibrium effects that alter solute rejection are taken asymptotically large, and a different nonlinear oscillator equation governs the location of the interface to leading order. A similar analysis allows for the derivation of an amplitude evolution equation for the uniformly pulsating interface. In this case, the bifurcation is always supercritical. The results are used to make predictions about the characteristics of solute bands that would be frozen into the solid.

Merchant, G. J.↗

NASA NTRS1974

Geophysical subsurface probing with radio-frequency interferometry

The radio-frequency interferometry method can be used to probe interiors of celestial bodies and terrestrial areas with low conductivity. In order to interpret the interference patterns, a theoretical study is made of the electromagnetic fields due to a dipole antenna on the surface of a horizontally stratified n-layered medium. Three approaches are used to calculate the interference patterns: direct numerical integration, asymptotic evaluation by the saddle point method, and a residue series approach. The asymptotic approach leads to the geometrical-optics interpretation. The residue approach leads to modal analysis. The validity of the formulation is checked by comparisons with analog model tank experiments and actual field data obtained from glaciers.

Kong, J. A.↗