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 343 records · Page 19

Principal normal curvature of surfaces

Certain principal normal curvatures of differential geometry were developed for use in curvature matrices associated with the asymptotic solution of electromagnetic diffraction problems. The effort is directed toward microwave antenna simulations and high speed digital computer analysis of radiometric instruments used to obtain soil moisture, sea state, salinity and temperature data. It is shown that the methods used to develop the principal normal curvatures for paraboloid, hyperboloid, ellipsoid, sphere, and cone can be applied to other radiometer geometries such as the parabolic torus, even though the surface parameterizations are different. It is concluded that deployable offset geometries, distorted by rotational forces and solar loads may be analyzed by similar means given a suitable surface description.

Schmidt, R. F.↗

Self-similar perturbations of a Friedmann universe

The present analysis of spherically symmetric self-similar solutions to the Einstein equations gives attention to those solutions that are asymptotically k = 0 Friedmann at large z values, and possess finite but perturbed density at the origin. Such solutions represent nonlinear density fluctuations which grow at the same rate as the universe's particle horizon. The overdense solutions span only a narrow range of parameters, and resemble static isothermal gas spheres just within the sonic point; the underdense solutions may have arbitrarily low density at the origin while exhibiting a unique relationship between amplitude and scale. Their relevance to large-scale void formation is considered.

Carr, Bernard J.↗

Simple robust control laws for robot manipulators. Part 2: Adaptive case

A new class of asymptotically stable adaptive control laws is introduced for application to the robotic manipulator. Unlike most applications of adaptive control theory to robotic manipulators, this analysis addresses the nonlinear dynamics directly without approximation, linearization, or ad hoc assumptions, and utilizes a parameterization based on physical (time-invariant) quantities. This approach is made possible by using energy-like Lyapunov functions which retain the nonlinear character and structure of the dynamics, rather than simple quadratic forms which are ubiquitous to the adaptive control literature, and which have bound the theory tightly to linear systems with unknown parameters. It is a unique feature of these results that the adaptive forms arise by straightforward certainty equivalence adaptation of their nonadaptive counterparts found in the companion to this paper (i.e., by replacing unknown quantities by their estimates) and that this simple approach leads to asymptotically stable closed-loop adaptive systems. Furthermore, it is emphasized that this approach does not require convergence of the parameter estimates (i.e., via persistent excitation), invertibility of the mass matrix estimate, or measurement of the joint accelerations.

Bayard, D. S.↗

Base pressure in laminar supersonic flow.

An asymptotic description is proposed for supersonic laminar flow over a wedge or a backward-facing step, for large Reynolds number and for a base or step height which is small compared with the boundary-layer length. The analysis is carried out for adiabatic wall conditions and a viscosity coefficient proportional to temperature. In a particular limit corresponding to a very thick boundary layer, a similarity law is obtained for the base pressure. For a thinner boundary layer an asymptotic form for the base pressure is obtained which shows the dependence on the parameters explicitly and which permits good agreement with experiment. This latter result is based on an inviscid-flow approximation for the corner expansion and for reattachment with viscous forces important primarily in a thin sublayer about the dividing streamline. A prediction of the pressure distribution at reattachment is given and the result is compared with experimental pressure distributions.

Messiter, A. F.↗

Eigenfrequency constrained topology optimization of finite strain hyperelastic structures

This paper incorporates hyperelastic materials, nonlinear kinematics, and preloads in eigenfrequency constrained density–based topology optimization. The formulation allows for initial finite deformations and subsequent small harmonic oscillations. The optimization problem is solved by the method of moving asymptotes, and the gradients are calculated using the adjoint method. Both simple and degenerate eigenfrequencies are considered in the sensitivity analysis. A well-posed topology optimization problem is formulated by filtering the volume fraction field. Numerical issues associated with excessive distortion and spurious eigenmodes in void regions are reduced by removing low volume fraction elements. The optimization objective is to maximize stiffness subject to a lower bound on the fundamental eigenfrequency. Numerical examples show that the eigenfrequencies drastically change with the load magnitude, and that the optimization is able to produce designs with the desired fundamental eigenfrequency.

42 ENGINEERING↗

Reciprocal swimming at intermediate Reynolds number

In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers Re≈1–1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, in this work, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate Re. In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.

42 ENGINEERING↗

Thermocapillary migration of bubbles and droplets

A discussion is presented concerning the motion of droplets in a surrounding fluid induced by a temperature gradient. Such 'thermocapillary migration' may be characterized as the tangential stress at the interface that is induced by the temperature gradient's variation of interfacial tension, causing motion in the neighboring fluid both inside and outside the droplet. As a result, the droplet moves in the direction of the pole with the lower interfacial tension. An analysis is conducted in the creeping flow limit, for the case of small Marangoni number values. The method of matched asymptotic expansions is used for the solution of the field equations.

Subramanian, R. S.↗

On the application of transonic similarity rules to wings of finite span

The transonic aerodynamic characteristics of wings of finite span are discussed from the point of view of a unified small perturbation theory for subsonic, transonic, and supersonic flows about thin wings. This approach avoids certain ambiguities which appear if one studies transonic flows by means of equations derived under the more restrictive assumption that the local velocities are everywhere close to sonic velocity. The relation between the two methods of analysis of transonic flow is examined, the similarity rules and known solutions of transonic flow theory are reviewed, and the asymptotic behavior of the lift, drag, and pitching-moment characteristics of wings of large and small aspect ratio is discussed. It is shown that certain methods of data presentation are advantageous for the effective display of these characteristics.

Spreiter, John R↗

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗

The stability of motion of satellites with flexible appendages

The mathematical formulation associated with the problem of stability of motion of a satellite consisting of a main rigid body and three (or less) pairs of flexible rods is presented. The rods are capable of flexure in two orthogonal directions. Whereas the rotational motion of the body is described by generalized coordinates depending on time alone, the elastic displacements of the rods depend both on spatial position and time. Assuming no external torques, there exist motion integrals in the form of momentum integrals. These integrals can be regarded as constraint equations relating the system velocities, and used to reduce the number of variables describing the motion. The stability analysis has been carried out by means of an extension of the Liapunov direct method. Since the elastic vibrations result in energy dissipation, it is shown that the equilibrium position is asymptotically stable if the Hamiltonian is positive definite and unstable if it can take negative values in the neighborhood of the equilibrium. Determining the sign definiteness of the Hamiltonian is complicated by the fact that it contains spatial derivatives of the elastic displacements. Two methods are presented to cope with this problem. The first, the standard modal analysis in conjunction with series truncation, develops criteria in terms of infinite series associated with the natural modes and frequencies of the elastic rods. The second, the method of integral coordinates, yields closed-form stability criteria involving the system parameters, such as the body moments of inertia, the length and mass distribution of the rods, the lowest natural frequencies of the rods, and the satellite spin velocity.

Meirovitch, L.↗

Quasisteady evolution of fast neutrino-flavor conversions

In astrophysical environments such as core-collapse supernovae (CCSNe) and binary neutron star mergers (BNSMs), neutrinos potentially experience substantial flavor mixing due to the refractive effects of neutrino self-interactions. Determining the survival probability of neutrinos in asymptotic states is paramount to incorporating flavor conversions’ effects in the theoretical modeling of CCSN and BNSM. Some phenomenological schemes have shown good performance in approximating asymptotic states of fast neutrino-flavor conversions (FFCs), known as one of the collective neutrino oscillation modes induced by neutrino self-interactions. However, a recent study showed that they would yield qualitatively different asymptotic states of FFC if the neutrino number is forced to evolve. It is not yet fully understood why the canonical phenomenological models fail to predict asymptotic states. In this paper, we perform detailed investigations through numerical simulations and then provide an intuitive explanation with a quasihomogeneous analysis. Based on the analysis, we propose a new phenomenological model, in which the quasisteady evolution of FFCs is analytically determined. The model also allows us to express the convolution term of spatial wave number as a concise form, which corresponds to useful information on analyses for the nonlinear feedback from small-scale flavor conversions to large-scale ones. Furthermore, our model yields excellent agreement with numerical simulations, which lends support to our interpretation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Elliptic Capture Orbits for Missions to the Near Planets

Elliptic capture orbits around Mars and Venus have often been considered as means for reducing arrival and departure energy requirements for two-way missions. It had also generally been feared that the energy savings obtained by capturing a spacecraft into a highly elliptical orbit (rather than a near circular orbit of the same periapsis) would largely be offset by the penalties incurred in aligning the semi-major axis of the ellipse in such a way as to obtain the proper orientation of the departure hyperbola. This paper, presents the results of an analysis which takes into consideration the penalties arising from the requirement to match the orientation of the elliptical orbit with the asymptote of the departure hyperbola. The scientific aspects of elliptical orbits around the target planet are discussed, and it is shown that such orbits exhibit characteristics which may be considered advantageous or disadvantageous depending on the purpose of the mission. Alignment of ' the semi-major axis of the capture, ellipse relative to the, asymptote of the escape hyperbola was found not to be a critical requirement since the kinetic energy remains high over a substantial portion of the elliptical capture orbit. This 'means that the escape stage can operate efficiently even when ignited at some angle from the true periapsis point. Considerable freedom in choosing this angle is available at little propulsive cost. The resulting latitude in the choice of angles between arrival and escape asymptotes makes it possible to consider a wide variety of interplanetary transfers and planetary staytimes without the need for separate propulsive maneuvers to realign the capture ellipse before departure., Special consideration has also been g1ven to plane change maneuvers around the planet. These may be required for reasons of orbit dynamics or scientific experimentation and are not uniquely tied to elliptical captures. The sensitivity of the mass of the excursion module to the eccentricity of the capture orbit is discussed and mass-penalty diagrams are presented. It is shown that these penalties do not materially offset the large gains obtained through the use of the elliptical capture mode.

Casal, Federico G.↗

Small-scale yielding analysis of mixed mode plane-strain crack problems

The small-scale yielding analysis of an elastic-plastic body with a line crack under plane-strain conditions subject to combinations of Mode I and II loadings is examined. The analysis of the near-tip field follows the works of Hutchinson and Rice and Rosengren. Dominant singularity solutions governing the asymptotic behavior of the stresses and strains at the crack tip are obtained for the complete range of loadings between Mode I and II. The results of an accurate finite element technique, which imbeds the dominant singularity solutions, directly relates the near-tip behavior to the elastic stress intensity factors K(I) and K(II). Implications of this study to mixed mode fracture mechanics is also discussed, particularly with respect to the direction of crack initiation and the relation of fracture toughness under mixed modes to that in Mode I. Details of the mixed mode plastic zone sizes and shapes are also given.

Shih, C. F.↗

Asymptotic behavior of vector recurrences with applications

An investigation is conducted of the behavior of a vector recurrence relation under very weak assumptions. The obtained results are used in an analysis of the iterations for nonlinear equations and a study concerning the composition of such iterations. The results can be employed in studies of one-point iterations involving the use of memory elements. Iterations for solving nonlinear equations on parallel computers can be investigated.

Feldstein, A.↗

Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations

This paper addresses the problem of detecting and characterizing local variability in time series and other forms of sequential data. The goal is to identify and characterize statistically significant variations, at the same time suppressing the inevitable corrupting observational errors. We present a simple nonparametric modeling technique and an algorithm implementing it-an improved and generalized version of Bayesian Blocks [Scargle 1998]-that finds the optimal segmentation of the data in the observation interval. The structure of the algorithm allows it to be used in either a real-time trigger mode, or a retrospective mode. Maximum likelihood or marginal posterior functions to measure model fitness are presented for events, binned counts, and measurements at arbitrary times with known error distributions. Problems addressed include those connected with data gaps, variable exposure, extension to piece- wise linear and piecewise exponential representations, multivariate time series data, analysis of variance, data on the circle, other data modes, and dispersed data. Simulations provide evidence that the detection efficiency for weak signals is close to a theoretical asymptotic limit derived by [Arias-Castro, Donoho and Huo 2003]. In the spirit of Reproducible Research [Donoho et al. (2008)] all of the code and data necessary to reproduce all of the figures in this paper are included as auxiliary material.

signal detection↗

Resolution requirements for numerical modeling of neutrino quantum kinetics

Neutrino quantum kinetics is a rapidly evolving field in computational astrophysics, with a primary focus on collective neutrino oscillations in core-collapse supernovae and postmerger phases of binary neutron star mergers. In recent years, there has been considerable debate concerning resolution dependence in numerical simulations. In this paper, we conduct a comprehensive resolution study in both angular- and spatial directions by using two independent schemes of quantum kinetic neutrino transport: finite volume and pseudospectral methods. We complement our discussion by linear stability analysis including inhomogeneous modes. Our result suggests that decreasing spatial resolutions underestimates the growth of flavor instability, and then leads to wrong asymptotic states of flavor conversions, which potentially has a critical impact on astrophysical consequences. We further delve into numerical results of low resolution simulations, that reveals the underlying mechanism responsible for numerical artifacts caused by insufficient resolutions. Finally, this study settles the debate on requirements of resolutions and serves as a guideline for numerical modeling of quantum kinetic neutrino transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Influence of conical wire array geometry on plasma flow and temperature profiles of radiatively cooled jets

The influence of conical wire array geometry on the formation and dynamics of pulsed-power driven plasma jets is investigated. In the experiments, the jet becomes isolated from the inflows as it passes through an aperture, allowing the study of its intrinsic evolution for different conical angles. Here, our results show that, regardless of the array opening angle, the jets are supersonic, highly collisional, and exhibit an exponential axial density decay with a characteristic scale length of 𝐿 𝑛 ≈ 3 mm, significantly shorter than the overall length of the jet. In contrast, axial velocity systematically increases with larger array opening angles. The near invariance of the density profile is attributed to a compensating mechanism between geometric divergence and axial acceleration, consistent with an asymptotically steady flow regime. Additionally, temperature measurements reveal ion-electron thermal decoupling near the base of the jet, with equilibration downstream. Analysis of radiative cooling and collisional energy exchange timescales indicates that the plasma evolves in a radiatively stable regime and that the two-temperature profile is fully described by the competition between both methods. These findings provide valuable information into the internal structure and evolution of pulsed-power plasma jets in regimes dominated by geometric shaping and radiative cooling, demonstrating experimental control over flow acceleration and collimation relevant to both laboratory and astrophysical contexts.

Physics - Plasma physics↗