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 235 records · Page 13

Complete quantification of parametric uncertainties in ( d , p ) transfer reactions

Deuteron-induced transfer reactions are a popular probe in nuclear structure and nuclear astro physics studies. The interpretation of these transfer measurements relies on reaction theory that takes as input effective interactions between the nucleons and the target nucleus. Previous work quantified the uncertainty associated with the optical potentials between the nucleons and the target. In this study, we extend that work by also including the parameters of the mean field associated with the overlap function of the final bound state, thus obtaining the full parametric uncertainty on transfer observables. We use Bayesian Markov Chain Monte Carlo simulations to obtain parameter posterior distributions. We use elastic-scattering cross sections to constrain the optical potential parameters and use the asymptotic normalization coefficient of the final state to constrain the bound state interaction. We then propagate these posteriors to the transfer angular distributions and obtain confidence intervals for this observable. We study (d,p) reactions on 14 C, 16 O, and 48 Ca at energies in the range $E$ $d$ = 7 - 24 MeV. Our results show a strong reduction in uncertainty by using the asymptotic normalization coefficient as a constraint, particularly for those reactions most sensitive to ambiguities in the mean field. For those reactions, the importance of constraining the bound state interaction is equal to that of constraining the optical potentials. Furthermore, the case of 14 C is an outlier because the cross section is less sensitive to the nuclear interior. When minimal constraints are used on the parameters of the nucleon-target interaction, the 1σ uncertainties on the differential cross sections are large (~ 500-100%). However, if elastic-scattering data and the asymptotic normalization coefficient are used in the analysis, with an error of 10% (5%), this uncertainty reduces to ~ 30% (~ 15%).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The instability of the thin vortex ring of constant vorticity

A theoretical study indicating that vortex rings at moderate Reynolds numbers are unstable to azimuthal bending waves is presented. Only the case of a thin vortex ring with a core of constant vorticity in an inviscid flow is examined. The disturbance flow and the mean flow of the vortex ring are derived as asymptotic solutions near the core; the stability analysis is developed completely for a certain class of bending waves that are unstable on a line filament in the presence of strain. The vortex ring is found to be always unstable for at least two wavenumbers for which waves on a line filament of the same vorticity distribution would not rotate. Published experimental results are cited to support these conclusions.

Widnall, S. E.↗

Analytical prediction of EM scattering by complex obstacles within electrically large open cavities

An analytical technique is described for predicting the electromagnetic (EM) scattering by complex obstacles within electrically large open waveguide cavities. This topic is of relevance to EM scattering and EMC type applications. A general open ended waveguide cavity is illustrated and is assumed to be relatively arbitrarily shaped, and the open front end of the cavity is illuminated by an arbitrarily polarized external plane wave field. The interior obstacle is assumed to be located within the cavity near the back end, and it generally contributes significantly to the fields scattered into the exterior by the entire cavity configuration. The present scattering analysis is restricted to within approximately 70 degs from boresight with respect to the open front end; outside of this region, the scattering by the external features of the cavity are generally more significant than from the interior of the cavity. The scattering analysis is composed of two main parts. One analysis is of the coupling of the incident plane wave field, via the open front end using the high frequency ray or beam shooting methods. The other analysis deals with separately obtaining the plane wave scattering characteristics of the interior obstacle. Thus the present analysis represents a hybrid combination of an asymptotic high frequency ray/beam method for tracking the fields through the cavity, together with any suitable method for obtaining a plane wave characterization of the interior obstacle.

Pathak, Prabhakar H.↗

A note on the effects of viscosity on the stability of a trailing-line vortex

The linear stability of the Batchelor (1964) vortex is examined with emphasis on new viscous modes recently found numerically by Khorrami (1991). Unlike the previously reported inviscid modes of instability, these modes are destabilized by viscosity and exhibit small growth rates at large Reynolds numbers. The analysis presented here uses a combination of asymptotic and numerical techniques. The results confirm the existence of the additional modes of instability due to viscosity.

Duck, Peter W.↗

Slender Vortex Filament with Slowly Varying Core Structure

We give a brief review of the asymptotic theory of slender vortex filaments with emphases on the choices of scalings characterizing the physical problems and the corresponding assumptions and/or restrictions introduced in the formation of the asymptotic theory of Callegari and Ting (1978) and its extension by Klein and Ting (1992). In particular, the slender filaments considered are assumed to be forming loops or tori. Because of this restriction, the theory is not applicable to the trailing vortex system of a rotorcraft. We describe the multiple length scales characterizing the vortex system, formulate the expansion scheme, derive the governing equations and then identify the assumptions or restrictions inherent in the multi-scale analysis and needed for the validity of the asymptotic theory of the trailing vortex system.

Ting, Lu↗

Nonlinear motion analysis of flexible satellites

The article reviews research on pointing accuracies of flexible satellites with regard to derived equations of motion of individual elements and composite bodies. The solution of these equations may be achieved through the formulation of a set of singularly perturbed equations by eliminating the generalized position vector. These equations are generally solved by a method of asymptotic expansions if certain convergence criteria are satisfied. The analysis may be applied to dual-spin flexible satellites.

Huang, T. C.↗

The crack problem in a specially orthotropic shell with double curvature

The crack problem of a shallow shell with two nonzero curvatures is considered. It is assumed that the crack lies in one of the principal planes of curvature and the shell is under Mode I loading condition. The material is assumed to be specially orthotropic. After giving the general formulation of the problem the asymptotic behavior of the stress state around the crack tip is examined. The analysis is based on Reissner's transverse shear theory. Thus, as in the bending of cracked plates, the asymptotic results are shown to be consistent with that obtained from the plane elasticity solution of crack problems. Rather extensive numerical results are obtained which show the effect of material orthotropy on the stress intensity factors in cylindrical and spherical shells and in shells with double curvature. Other results include the stress intensity factors in isotropic toroidal shells with positive or negative curvature ratio, the distribution of the membrane stress resultant outside the crack, and the influence of the material orthotropy on the angular distribution of the stresses around the crack tip.

Delale, F.↗

The crack problem in a specially orthotropic shell with double curvature

The crack problem of a shallow shell with two nonzero curvatures is considered. It is assumed that the crack lies in one of the principal planes of curvature and the shell is under Mode I loading condition. The material is assumed to be specially orthotropic. After giving the general formulation of the problem the asymptotic behavior of the stress state around the crack tip is examined. The analysis is based on Reissner's transverse shear theory. Thus, as in the bending of cracked plates, the asymptotic results are shown to be consistent with that obtained from the plane elasticity solution of crack problems. Rather extensive numerical results are obtained which show the effect of material orthotropy on the stress intensity factors in cylindrical and spherical shells and in shells with double curvature. Other results include the stress intensity factors in isotropic toroidal shells with positive or negative curvature ratio, the distribution of the membrane stress resultant outside the crack, and the influence of the material orthotropy on the angular distribution of the stresses around the crack tip. Previously announced in STAR as N83-16782

Delale, F.↗

On the transition towards slow manifold in shallow-water and 3D Euler equations in a rotating frame

The long-time, asymptotic state of rotating homogeneous shallow-water equations is investigated. Our analysis is based on long-time averaged rotating shallow-water equations describing interactions of large-scale, horizontal, two-dimensional motions with surface inertial-gravity waves field for a shallow, uniformly rotating fluid layer. These equations are obtained in two steps: first by introducing a Poincare/Kelvin linear propagator directly into classical shallow-water equations, then by averaging. The averaged equations describe interaction of wave fields with large-scale motions on time scales long compared to the time scale 1/f(sub o) introduced by rotation (f(sub o)/2-angular velocity of background rotation). The present analysis is similar to the one presented by Waleffe (1991) for 3D Euler equations in a rotating frame. However, since three-wave interactions in rotating shallow-water equations are forbidden, the final equations describing the asymptotic state are simplified considerably. Special emphasis is given to a new conservation law found in the asymptotic state and decoupling of the dynamics of the divergence free part of the velocity field. The possible rising of a decoupled dynamics in the asymptotic state is also investigated for homogeneous turbulence subjected to a background rotation. In our analysis we use long-time expansion, where the velocity field is decomposed into the 'slow manifold' part (the manifold which is unaffected by the linear 'rapid' effects of rotation or the inertial waves) and a formal 3D disturbance. We derive the physical space version of the long-time averaged equations and consider an invariant, basis-free derivation. This formulation can be used to generalize Waleffe's (1991) helical decomposition to viscous inhomogeneous flows (e.g. problems in cylindrical geometry with no-slip boundary conditions on the cylinder surface and homogeneous in the vertical direction).

Mahalov, A.↗

Comparison of continuous and discrete adaptive identification algorithms

Discretization of a popular continuous-time control algorithm is effected and an equivalent discrete-time identification law developed and compared to a published discrete identification algorithm developed from Lyapunov Theory. Results are compared as regards asymptotic stability as insured using Lyapunov theory. Some analysis and design guidelines are proposed as regards implementation and practical utility.

Colburn, B. K.↗

Inferring spectral reflectances of plant elements by simple inversion of bidirectional reflectance measurements

Inverting previously developed explicit expressions for a vertical architecture, bidirectional reflectances measured over corn viewing from the solar quadrant at azimuths near the principal plane are used to determine the spectral reflectances of plant elements. The leaf reflectance values extracted in three visible bands at viewing zenith angles of 70 deg, 60 deg, and 45 deg agree closely with laboratory-measured reflectances of corn leaves. At viewing zenith angle of 30 deg, the inversion breaks down, inasmuch as the inferred plant element reflectances are too high. Satisfactory results are also achieved when the same approach is applied to bidirectional reflectances measured over potted balsam firs, but when applied to soybeans reflectances, the procedure yields unreasonably high leaf reflectances. The failure in this case is attributed to the nonvertical architecture of the soybean canopy; however, for this canopy, inversion based on horizontal architecture is possible. The bidirectional reflectances measured from the solar quadrant, at viewing angles appreciably far from 'hot spot' viewing, approximately equal in magnitude half of the leaf reflectance. The 0.5 ratio is predicted by a previous analysis of opaque horizontal Lambertian facets, as the asymptotic value for a dense canopy at any viewing angle. For soybeans, this ratio applies very closely at 15 deg viewing zenith angle. The results suggest that inversion based on simple architecture, applying explicit expressions, might be of value, either in itself or as a preliminary step before inversion applying complex models.

Otterman, J.↗

Analytic theory for the determination of velocity and stability of bubbles in a Hele-Shaw cell. I - Velocity selection. II - Stability

An asymptotic theory is presented for the determination of velocity and linear stability of a steady symmetric bubble in a Hele-Shaw cell for small surface tension. First the bubble velocity relative to the fluid velocity at infinity is determined for small surface tension by means of a transcendentally small correction to the asymptotic series solution. In addition, a linear stability analysis shows that only the solution branch corresponding to the largest possible bubble velocity for given surface tension is stable, while all the others are unstable.

Tanveer, S.↗

An Analysis of the Optimal Control Modification Method Applied to Flutter Suppression

Unlike basic Model Reference Adaptive Control (MRAC)l, Optimal Control Modification (OCM) has been shown to be a promising MRAC modification with robustness and analytical properties not present in other adaptive control methods. This paper presents an analysis of the OCM method, and how the asymptotic property of OCM is useful for analyzing and tuning the controller. We begin with a Lyapunov stability proof of an OCM controller having two adaptive gain terms, then the less conservative and easily analyzed OCM asymptotic property is presented. Two numerical examples are used to show how this property can accurately predict steady state stability and quantitative robustness in the presence of time delay, and relative to linear plant perturbations, and nominal Loop Transfer Recovery (LTR) tuning. The asymptotic property of the OCM controller is then used as an aid in tuning the controller applied to a large scale aeroservoelastic longitudinal aircraft model for flutter suppression. Control with OCM adaptive augmentation is shown to improve performance over that of the nominal non-adaptive controller when significant disparities exist between the controller/observer model and the true plant model.

aeroservoelasti↗

Influence of Incidence Angle on Sound Generation by Airfoils Interacting with High-Frequency Gusts

A theoretical model is developed for the sound generated when a convected vortical or entropic gust encounters an airfoil at nonzero angle of attack. The theory is based on a linearization of the Euler equations about the steady, subsonic flow past the airfoil. High-frequency gusts, whose wavelengths are short compared to the airfoil chord, but long compared to the displacement of the mean-flow stagnation point from the leading edge, are considered. The analysis utilizes singular perturbation techniques and involves four asymptotic regions. Local regions, which scale on the gust wavelength, are present at the airfoil leading and trailing edges. Behind the airfoil a "transition" region, which is similar to the transition zone between illuminated and shadow zones in optical problems, is present. In the outer region, far away from the airfoil edges and wake, the solution has a geometric-acoustics form. The primary sound generation is found to be concentrated in the local leading-edge region. The trailing edge plays a secondary role as a scatterer of the sound generated in the leading-edge region. Parametric calculations are presented which illustrate that moderate levels of airfoil steady loading can significantly affect the sound field produced by airfoil-gust interactions.

Matthew R. Myers↗

Optimal reusable-tug and expendable-kickstage trajectories for high-energy planetary missions including correction for nodal precession

Equations are derived by using branched trajectory optimization techniques and the maximum principle to maximize the payload capability of a reusable tug/expendable kickstage vehicle configuration for planetary missions. The two stages and the payload are launched into a low earth orbit by a single space shuttle. The analysis includes correction for precession of the orbit. This correction is done by the tug. The tug propels the payload and the kickstage to an energy beyond earth escape and returns within a specified time to the precessed orbit. After separating from the tug, the kickstage accelerates the payload to the required injection conditions. Planetary injection conditions are specified by the mission energy and a fixed declination and right ascension of the outgoing asymptote. The multipoint boundary value problem resulting from the analysis is solved by a Newton-Raphson iteration technique. Partial derivatives of the boundary conditions are obtained by perturbing the initial conditions one at a time, integrating the trajectory and adjoint equations, and observing the changes in boundary conditions. Maximum payload capability is derived for two typical mission energies. In addition, the variations of several mission and stage parameters are also examined.

Borsody, J.↗

The fully nonlinear development of Goertler vortices in growing boundary layers

The fully nonlinear development of small wavelength Goertler vortices in a growing boundary layer is investigated using a combination of asymptotic and numerical methods. The starting point for the analysis is the weakly nonlinear theory of Hall (1982b) who discussed the initial development of small amplitude vortices in a neighborhood of the location where they first become linearly unstable. That development is unusual in the context of nonlinear stability theory in that it is not described by the Stuart-Watson approach. In fact the development is governed by a pair of coupled nonlinear partial differential evolution equations for the vortex flow and the mean flow correction. Here the further development of this interaction is considered for vortices so large that the mean flow correction driven by them is as large as the basic state. Surprisingly it is found that such a nonlinear interaction can still be described by asymptotic means. It is shown that the vortices spread out across the boundary layer and effectively drive the boundary layer. In fact the system obtained by writing down the equations for the fundamental component of the vortex generate a differential equation for the basic state. Thus the mean flow adjusts so as to make these large amplitude vortices locally neutral. Moreover in the region where the vortices exist the mean flow has a square-root profile and the vortex velocity field can be written down in closed form. The upper and lower boundaries of the region of vortex activity are determined by a free-boundary problem involving the boundary layer equations. In general it is found that this region ultimately includes almost all of the original boundary layer and much of the free-stream. In this situation the mean flow has essentially no relationship to the flow which exists in the absence of the vortices.

Hall, Philip↗

Passive Roll Stabilization of the Near Earth Asteroid Scout Solar Sail Mission

The Near Earth Asteroid (NEA) Scout is a small satellite, solar sail mission set to launch on Artemis I. Analysis of the NEA Scout solar radiation pressure model determined that, for certain solar incidence angles, there exists at least one locally stable equilibrium point about the “roll” axis, or the axis normal to the solar sail plane. This analysis is extended to three other geometrically similar solar sail models, and a Lyapunov stability analysis is conducted demonstrating that the stable equilibria are locally stable for the majority of roll angles in a wide range of solar incidence angles. Under certain assumptions, it is possible that only two axis control is required for non-spinning solar sails. If the roll axis is controlled, an active rate damping controller can be designed to asymptotically stabilize the system. Simulation results are presented that follow the stability analysis, where in the homogenous case, the system enters a limit cycle and, in the actively damped case, the system asymptotically converges to the equilibrium point.

Ivan Rodrigues Bertaska↗

Method of Green’s Potentials for Elliptic PDEs in Domains with Random Apertures

Problems with topological uncertainties appear in many fields ranging from nano-device engineering to the design of bridges. In many of such problems, a part of the domains boundaries is subjected to random perturbations making inefficient conventional schemes that rely on discretization of the whole domain. Here, we study elliptic PDEs in domains with boundaries comprised of a deterministic part and random apertures, and apply the method of modified potentials with Green’s kernels defined on the deterministic part of the domain. This approach allows to reduce the dimension of the original differential problem by reformulating it as a boundary integral equation posed on the random apertures only. The multilevel Monte Carlo method is then applied to this modified integral equation and its optimal ϵ -2 asymptotical complexity is shown. Finally, we provide the qualitative analysis of the proposed technique and support it with numerical results.

97 MATHEMATICS AND COMPUTING↗