Engineering Papers⌕ Search

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.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Some spectral approximation of one-dimensional fourth-order problems

Some spectral type collocation method well suited for the approximation of fourth-order systems are proposed. The model problem is the biharmonic equation, in one and two dimensions when the boundary conditions are periodic in one direction. It is proved that the standard Gauss-Lobatto nodes are not the best choice for the collocation points. Then, a new set of nodes related to some generalized Gauss type quadrature formulas is proposed. Also provided is a complete analysis of these formulas including some new issues about the asymptotic behavior of the weights and we apply these results to the analysis of the collocation method.

Bernardi, Christine↗
NASA NTRS1999

Directional Solidification of a Binary Alloy into a Cellular Convective Flow: Localized Morphologies

A steady, two dimensional cellular convection modifies the morphological instability of a binary alloy that undergoes directional solidification. When the convection wavelength is far longer than that of the morphological cells, the behavior of the moving front is described by a slow, spatial-temporal dynamics obtained through a multiple-scale analysis. The resulting system has a "parametric-excitation" structure in space, with complex parameters characterizing the interactions between flow, solute diffusion, and rejection. The convection stabilizes two dimensional disturbances oriented with the flow, but destabilizes three dimensional disturbances in general. When the flow is weak, the morphological instability behaves incommensurably to the flow wavelength, but becomes quantized and forced to fit into the flow-box as the flow gets stronger. At large flow magnitudes the instability is localized, confined in narrow envelopes with cells traveling with the flow. In this case the solutions are discrete eigenstates in an unbounded space. Their stability boundary and asymptotics are obtained by the WKB analysis.

Chen, Y.- J.↗
NASA NTRS2019

Multilevel Monte Carlo Estimation of Unbiased Expectation via Sample Reuse and the Low Variance Estimation of Asymptotic Rates

A new variant of the multilevel Monte Carlo estimator [5, 3, 9, 12] is presented for the estimation of expectation statistics that utilizes sample reuse in specified levels, explicitly removes approximation error bias associated with numerically computed output quantities of interest that have an asymptotic limit behavior, and permits a low variance estimate of the asymptotic rate of convergence to that limit. In addition, it is shown that this new multilevel Monte Carlo variant can yield a computational cost savings. A review of Monte Carlo and multilevel Monte Carlo estimators is presented that includes analysis of expected value, expected mean squared error, and the calculation of optimized multilevel sample size parameters. The multilevel Monte Carlo estimator produces estimates of expectation for numerically approximated output quantities of interest that are biased by approximation error. When the quantity of interest can be modeled as the asymptotic limit of numerically approximated output quantities of interest, it is theoretically possible to remove this approximation error bias in the multilevel Monte Carlo estimator. In actual implementations, however, this procedure is unreliable due to statistical variability and inaccuracy in estimating the needed asymptotic limit. Analysis and numerical experiment show that the proposed variant of the multilevel Monte Carlo method greatly reduces (in some cases eliminates) the statistical variability in this limit estimation.

Barth, Timothy↗
DOE OSTI2022

Generalized Analysis of Dynamic Flow Fouling on Heat Transfer Surfaces

Several energy applications involve flow of untreated heat transfer fluids that contain trace amounts of mineral salts, such as calcium carbonate and calcium sulfate, which show reverse solubility with temperature. The gradual accretion of these salts onto the heat transfer surfaces causes fouling that degrades thermal performance over time. While fouling of bare surfaces has been extensively studied, fouling of nonwetting superhydrophobic or liquid-infused surfaces is less understood, especially under dynamic flow conditions. This study presents a systematic experimental analysis to elucidate dynamic flow fouling of calcium carbonate and calcium sulfate on smooth, superhydrophobic and liquid-infused surfaces in a generalized manner. Here, copper tubes with modified inner surface wettability are fabricated and subjected to forced convection at different flow rates and fouling salt concentrations. Fouling characteristics of the two salts on the different surfaces are examined through microstructural analysis, and quantified in terms of asymptotic fouling resistance that measures the total extent of fouling, and a threshold time, which correlates to the onset of significant fouling. Using Taguchi design of experiments, the study reports, for the first time, closed-form analytical relationships for the asymptotic fouling thermal resistance and threshold time on the Reynolds number, dimensionless concentration of fouling agent and dimensionless infused liquid viscosity that represents the different surface types in a unified manner. It is shown that the analytical model accurately predicts the fouling characteristics in up to about 97% of the data, based on which contour maps of optimum surface designs for minimizing fouling resistance are presented. The results provide information on designing surfaces for reduced fouling or for estimating the fouling characteristics for a given surface in applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Semiglobal Safety-Filtered Extremum Seeking With Unknown CBFs

We introduce a safe extremum-seeking (Safe ES) algorithm which achieves the minimization of an unknown objective function while ensuring that an unknown, yet measured, control barrier function (CBF) remains above an arbitrarily small negative value for all time. In other words, “practical safety” is maintained during the entire period of convergence to the constrained extremum. Our design is based on quadratic program (QP) CBF style filters for safety, which is applied in an average and estimated sense. Using nonsmooth analysis tools, we guarantee semiglobal practical asymptotic (SPA) stability of the global constrained optimum, practical convergence to the safe set if starting in a condition violating the CBF, and practical safety for all time—semiglobally—if starting in safe set. The safety result of the paper is analogous with modern notions of SPA stability, guaranteeing that, for any small violation of safety, there exist design coefficients which guarantee that such a small violation is not exceeded. The paper outlines a set of sufficient conditions on the barrier and objective functions, and by way of a Lyapunov argument, we demonstrate that nonconvex constrained optimization problems can be solved. We present these results in the setting of a static map and a dynamical system. A simulation example illustrates the results.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

Viscous, radiating hypersonic flow about a blunt body

The viscous, radiating hypersonic flow past an axisymmetric blunt body is analyzed based on the Navier-Stokes equations, plus a radiative equation of transfer derived from the Milne-Eddington differential approximation. The fluid is assumed to be a perfect gas with constant specific heats, a constant Prandtl number of order unity, a viscosity coefficient varying as a power of the temperature, and an absorption coefficient varying as the first power of the density and as a power of the temperature. The gray gas assumption is invoked, thereby making the absorption coefficient independent of the spectral frequency. Limiting forms of the solutions are studied as the freestream Mach number freestream Reynolds number and the temperature ratio across the shock wave, go to infinity, and as the Bouguer number and the density ratio across the shock wave go to zero. The method of matched asymptotic expansions is used in the analysis, and it is shown that there is a far-field precursor, composed of two regions, in which the fluid mechanics can be neglected for all practical purposes but included for completeness.

Passamaneck, R. S.↗
NASA NTRS1993

An augmented Young-Laplace model of an evaporating meniscus in a micro-channel with high heat flux

High flux evaporations from a steady meniscus formed in a 2 micron channel is modeled using the augmented Young-Laplace equation. The heat flux is found to be a function of the long range van der Waals dispersion force which represents interfacial conditions between heptane and various substrates. Heat fluxes of (1.3-1.6) x 10(exp 6) W/m(exp 2) based on the width of the channel are obtained for heptane completely wetting the substrate at 100 C. Small channels are used to obtain these large fluxes. Even though the real contact angle is 0 deg, the apparent contact angle is found to vary between 24.8 deg and 25.6 deg. The apparent contact angle, which represents viscous losses near the contact line, has a large effect on the heat flow rate because of its effect on capillary suction and the area of the meniscus. The interfacial heat flux is modeled using kinetic theory for the evaporation rate. The superheated state depends on the temperature and the pressure of the liquid phase. The liquid pressure differs from the pressure of the vapor phase due to capillarity and long range van der Waals dispersion forces which are relevant in the ultra think film formed at the leading edge of the meniscus. Important pressure gradients in the thin film cause a substantial apparent contact angle for a complete wetting system. The temperature of the liquid is related to the evaporation rate and to the substrate temperature through the steady heat conduction equation. Conduction in the liquid phase is calculated using finite element analysis except in the vicinity of the thin film. A lubrication theory solution for the thin film is combined with the finite element analysis by the method of matched asymptotic expansions.

Wayner, P. C., Jr.↗
NASA NTRS1999

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

An optimal control approach using variable-structure (sliding-mode) tracking for large angle spacecraft maneuvers is presented. The approach expands upon a previously derived regulation result using a quaternion parameterization for the kinematic equations of motion. This parameterization is used since it is free of singularities. The main contribution of this paper is the utilization of a simple term in the control law that produces a maneuver to the reference attitude trajectory in the shortest distance. Also, a multiplicative error quaternion between the desired and actual attitude is used to derive the control law. Sliding-mode switching surfaces are derived using an optimal-control analysis. Control laws are given using either external torque commands or reaction wheel commands. Global asymptotic stability is shown for both cases using a Lyapunov analysis. Simulation results are shown which use the new control strategy to stabilize the motion of the Microwave Anisotropy Probe spacecraft.

John L. Crassidis↗
DOE OSTI2026

Control Parameter Sensitivity Study for Inverter-Based-Resource Dominated Grids: A Small Signal Stability Approach and Framework

The growing adoption of renewable energy is driving the prevalence of inverter-based resources (IBRs) within power grids. Future power grids will integrate both grid-following IBRs (GFM-IBRs) and grid-forming IBRs (GFL-IBRs) alongside synchronous generators. Therefore, it is crucial to perform stability studies that account for all components and especially control interactions related to IBRs. Extensive research has performed to study the IBR-related stability, however, the sensitivity study of IBRs' control parameters on system stability has not been adequately studied yet, especially from a systematic way. Therefore, this paper conducts a small signal stability analysis for a generic grid with multiple types of resources and develops an analytical framework for assessing the sensitivity of control parameters affecting stability margins. To achieve that, the non-autonomous reduced-order non-linear dynamic model is developed for a generic power system with multiple synchronous generator-based resources (SGBRs), GFM-IBRs, and GFL-IBRs. Based on the analytic model, a systematic framework for parametric sensitivity on systems' asymptotic stability is developed. A parameter sensitivity analysis based on eigenvalue methods is proposed. The impact of the droop controllers of GFM-IBRs, PQ-dispatch and the PLL controller of GFL-IBR on the system asymptotic stability is discussed. This sensitivity study is aiming to provide deep insights on control parameters' impact on system stability, and gives direction for parameter tuning in case of instability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1991

The asymptotic spectra of banded Toeplitz and quasi-Toeplitz matrices

Toeplitz matrices occur in many mathematical, as well as, scientific and engineering investigations. This paper considers the spectra of banded Toeplitz and quasi-Toeplitz matrices with emphasis on non-normal matrices of arbitrarily large order and relatively small bandwidth. These are the type of matrices that appear in the investigation of stability and convergence of difference approximations to partial differential equations. Quasi-Toeplitz matrices are the result of non-Dirichlet boundary conditions for the difference approximations. The eigenvalue problem for a banded Toeplitz or quasi-Toeplitz matrix of large order is, in general, analytically intractable and (for non-normal matrices) numerically unreliable. An asymptotic (matrix order approaches infinity) approach partitions the eigenvalue analysis of a quasi-Toeplitz matrix into two parts, namely the analysis for the boundary condition independent spectrum and the analysis for the boundary condition dependent spectrum. The boundary condition independent spectrum is the same as the pure Toeplitz matrix spectrum. Algorithms for computing both parts of the spectrum are presented. Examples are used to demonstrate the utility of the algorithms, to present some interesting spectra, and to point out some of the numerical difficulties encountered when conventional matrix eigenvalue routines are employed for non-normal matrices of large order. The analysis for the Toeplitz spectrum also leads to a diagonal similarity transformation that improves conventional numerical eigenvalue computations. Finally, the algorithm for the asymptotic spectrum is extended to the Toeplitz generalized eigenvalue problem which occurs, for example, in the stability of Pade type difference approximations to differential equations.

Beam, Richard M.↗
DOE OSTI2023

Resurgence, conformal blocks, and the sum over geometries in quantum gravity

In two dimensional conformal field theories the limit of large central charge plays the role of a semi-classical limit. Certain universal observables, such as conformal blocks involving the exchange of the identity operator, can be expanded around this classical limit in powers of the central charge c. This expansion is an asymptotic series, so — via the same resurgence analysis familiar from quantum mechanics — necessitates the existence of non-perturbative effects. In the case of identity conformal blocks, these new effects have a simple interpretation: the CFT must possess new primary operators with dimension of order the central charge. This constrains the data of CFTs with large central charge in a way that is similar to (but distinct from) the conformal bootstrap. We study this phenomenon in three ways: numerically, analytically using Zamolodchikov’s recursion relations, and by considering non-unitary minimal models with large (negative) central charge. In the holographic dual to a CFT2, the expansion in powers of c is the perturbative loop expansion in powers of ћ. So our results imply that the graviton loop expansion is an asymptotic series, whose cure requires the inclusion of new saddle points in the gravitational path integral. In certain cases these saddle points have a simple interpretation: they are conical excesses, particle-like states with negative mass which are not in the physical spectrum but nevertheless appear as non-manifold saddle points that control the asymptotic behaviour of the loop expansion. This phenomenon also has an interpretation in SL(2, R) Chern-Simons theory, where the non-perturbative effects are associated with the non-Teichmüller component of the moduli space of flat connections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

Estimating the change in asymptotic direction due to secular changes in the geomagnetic field

The concept of geomagnetic optics, as described by the asymptotic directions of approach, is extremely useful in the analysis of cosmic radiation data. However, when changes in cutoff occur as a result of evolution in the geomagnetic field, there are corresponding changes in the asymptotic cones of acceptance. A method is introduced of estimating the change in the asymptotic direction of approach for vertically incident cosmic ray particles from a reference set of directions at a specific epoch by considering the change in the geomagnetic cutoff.

Flueckiger, E. O.↗
DOE OSTI2021

Static response and Love numbers of Schwarzschild black holes

We derive the quadratic action for the physical degrees of freedom of massless spin-0, spin-1, and spin-2 perturbations on a Schwarzschild-(A)dS background in arbitrary dimensions. Here, we then use these results to compute the static response of asymptotically flat Schwarzschild black holes to external fields. Our analysis reproduces known facts about black hole Love numbers — in particular that they vanish for all types of perturbation in four spacetime dimensions — but also leads to new results. For instance, we find that neutral Schwarzschild black holes polarize in the presence of an electromagnetic background in any number of spacetime dimensions except four. Moreover, we calculate for the first time black hole Love numbers for vector-type gravitational perturbations in higher dimensions and find that they generically do not vanish. Along the way, we shed some light on an apparent discrepancy between previous results in the literature, and clarify some aspects of the matching between perturbative calculations of static response on a Schwarzschild background and the point-particle effective theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1975

Analysis of transonic flow about lifting wing-body configurations

An analytical solution was obtained for the perturbation velocity potential for transonic flow about lifting wing-body configurations with order-one span-length ratios and small reduced-span-length ratios and equivalent-thickness-length ratios. The analysis is performed with the method of matched asymptotic expansions. The angles of attack which are considered are small but are large enough to insure that the effects of lift in the region far from the configuration are either dominant or comparable with the effects of thickness. The modification to the equivalence rule which accounts for these lift effects is determined. An analysis of transonic flow about lifting wings with large aspect ratios is also presented.

Barnwell, R. W.↗