Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ASYMPTOTIC FUNCTION”

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 145 records · Page 8

Modes of interconnected lattice trusses using continuum models, part 1

This represents a continuing systematic attempt to explore the use of continuum models--in contrast to the Finite Element Models currently universally in use--to develop feedback control laws for stability enhancement of structures, particularly large structures, for deployment in space. We shall show that for the control objective, continuum models do offer unique advantages. It must be admitted of course that developing continuum models for arbitrary structures is no easy task. In this paper we take advantage of the special nature of current Large Space Structures--typified by the NASA-LaRC Evolutionary Model which will be our main concern--which consists of interconnected orthogonal lattice trusses each with identical bays. Using an equivalent one-dimensional Timoshenko beam model, we develop an almost complete continuum model for the evolutionary structure. We do this in stages, beginning only with the main bus as flexible and then going on to make all the appendages also flexible-except for the antenna structure. Based on these models we proceed to develop formulas for mode frequencies and shapes. These are shown to be the roots of the determinant of a matrix of small dimension compared with mode calculations using Finite Element Models, even though the matrix involves transcendental functions. The formulas allow us to study asymptotic properties of the modes and how they evolve as we increase the number of bodies which are treated as flexible. The asymptotics, in fact, become simpler.

Balakrishnan, A. V.↗

Scaling properties of induction times in heterogeneous nucleation

The heterogeneous-to-homogeneous induction time ratio is obtained as a function of the contact angle in the asymptotic limit of a high nucleation barrier. Model-dependent corrections to t(ind) are investigated, particularly in cases of the Turnbull-Fisher model used in numerical simulations by Greer et al. (1990).

Shneidman, Vitaly A.↗

Dynamics and Control of Articulated Anisotropic Timoshenko Beams

The paper illustrates the use of continuum models in control design for stabilizing flexible structures. A 6-DOF anisotropic Timoshenko beam with discrete nodes where lumped masses or actuators are located provides a sufficiently rich model to be of interest for mathematical theory as well as practical application. We develop concepts and tools to help answer engineering questions without having to resort to ad hoc heuristic ("physical") arguments or faith. In this sense the paper is more mathematically oriented than engineering papers and vice versa at the same time. For instance we make precise time-domain solutions using the theory of semigroups of operators rather than formal "inverse Laplace transforms." We show that the modes arise as eigenvalues of the generator of the semigroup, which are then related to the eigenvalues of the stiffness operator. With the feedback control, the modes are no longer orthogonal and the question naturally arises as to whether there is still a modal expansion. Here we prove that the eigenfunctions yield a biorthogonal Riesz basis and indicate the corresponding expansion. We prove mathematically that the number of eigenvalues is nonfinite, based on the theory of zeros of entire functions. We make precise the notion of asymptotic modes and indicate how to calculate them. Although limited by space, we do consider the root locus problem and show for instance that the damping at first increases as the control gain increases but starts to decrease at a critical value, and goes to zero as the gain increases without bound. The undamped oscillatory modes remain oscillatory and the rigid-body modes go over into deadbeat modes. The Timoshenko model dynamics are translated into a canonical wave equation in a Hilbert space. The solution is shown to require the use of an "energy" norm which is no more than the total energy: potential plus kinetic. We show that, under an appropriate extension of the notion of controllability, rate feedback with a collocated sensor can stabilize the structure in the sense that all modes are damped and the energy decays to zero. An example, non-numeric, is worked out in some detail illustrating the concepts and theory developed.

Balakrishnan, A. V.↗

Understanding how molecular hydrogen impacts the total dose and dose rate response of linear bipolar circuits

Recent Enhanced Low Dose Rate Sensitivity (ELDRS) investigations carried out by RLP Research, Arizona State University (ASU) and Jet Propulsion Laboratory (JPL) have shown significant differences in the degradation of bipolar micro-circuits with total dose in the presence of molecular hydrogen (H2) in packages. This has a significant impact on radiation hardness assurance and opens up opportunities to improve device performance. The general objectives of this program are to: 1. Determine the extent to which hydrogen contamination affects the total dose and dose rate response of linear bipolar circuits; 2. Develop a model that will enable the prediction of high dose rate (HDR) and low dose rate (LDR) response asymptotes and transition dose rates as a function of total dose, temperature, pressurized hydrogen, defect precursors, and other process dependent variables; 3. Explore the possibility of an accelerated hardness assurance method and possible hardening approaches; and 4. Extend the work to other technologies that have total dose response affected by hydrogen contamination.In this document, we xperimentally demonstrate with test transistors and circuits that hydrogen is correlated with ELDRS in bipolar linear circuits. We show that the amount of hydrogen determines: 1) the total dose response versus dose rate and 2) the transition dose rates between the high and low dose rate responses. The experimental results are supported with a steady state drift/diffusion analytical model as well as modeling calculations using COMSOL Multiphysics.

Pease, Ron↗

Asymptotic Expansions and Estimates for the Capillary Problem

The asymptotic properties for the small Bond number B of the equilibrium capillary interface interior to a circular cylindrical tube vertically dipped in an infinite reservoir of liquid are discussed. (The Bond number B is a dimensionless parameter which is the ratio of gravitational to capillary forces.) The formal expansion in powers of B of the solution to the differential equation describing the equilibrium surface (as can be obtained by standard perturbation methods) is proved to be truly asymptotic -- to all orders and uniformly in the variable and parameter gamma, the contact angle. Sequences of general estimates, in closed form, from above and from below, are also given for the solution and related functions. The M-th term in these sequences are asymptotically exact to order m. An idiosyncrasy of the problem, crucial in obtaining these estimates, is the absolute monotonicity of the structural function of the system in integral form.

Brulois, F. P.↗

Theory and modeling of atmospheric turbulence, part 2

Two dimensional geostrophic turbulence driven by a random force is investigated. Based on the Liouville equation, which simulates the primitive hydrodynamical equations, a group-kinetic theory of turbulence is developed and the kinetic equation of the scaled singlet distribution is derived. The kinetic equation is transformed into an equation of spectral balance in the equilibrium and non-equilibrium states. Comparison is made between the propagators and the Green's functions in the case of the non-asymptotic quasi-linear equation to prove the equivalence of both kinds of approximations used to describe perturbed trajectories of plasma turbulence. The microdynamical state of fluid turbulence is described by a hydrodynamical system and transformed into a master equation analogous to the Vlasov equation for plasma turbulence. The spectral balance for the velocity fluctuations of individual components shows that the scaled pressure strain correlation and the cascade transfer are two transport functions that play the most important roles.

Chen, C. M.↗

A uniform geometrical optics and an extended uniform geometrical theory of diffraction for evaluating high frequency EM fields near smooth caustics and composite shadow boundaries

A uniform geometrical optics (UGO) and an extended uniform geometrical theory of diffraction (EUTD) are developed for evaluating high frequency electromagnetic (EM) fields within transition regions associated with a two and three dimensional smooth caustic of reflected rays and a composite shadow boundary formed by the caustic termination or the confluence of the caustic with the reflection shadow boundary (RSB). The UGO is a uniform version of the classic geometrical optics (GO). It retains the simple ray optical expressions of classic GO and employs a new set of uniform reflection coefficients. The UGO also includes a uniform version of the complex GO ray field that exists on the dark side of the smooth caustic. The EUTD is an extension of the classic uniform geometrical theory of diffraction (UTD) and accounts for the non-ray optical behavior of the UGO reflected field near caustics by using a two-variable transition function in the expressions for the edge diffraction coefficients. It also uniformly recovers the classic UTD behavior of the edge diffracted field outside the composite shadow boundary transition region. The approach employed for constructing the UGO/EUTD solution is based on a spatial domain physical optics (PO) radiation integral representation for the fields which is then reduced using uniform asymptotic procedures. The UGO/EUTD analysis is also employed to investigate the far-zone RCS problem of plane wave scattering from two and three dimensional polynomial defined surfaces, and uniform reflection, zero-curvature, and edge diffraction coefficients are derived. Numerical results for the scattering and diffraction from cubic and fourth order polynomial strips are also shown and the UGO/EUTD solution is validated by comparison to an independent moment method (MM) solution. The UGO/EUTD solution is also compared with the classic GO/UTD solution. The failure of the classic techniques near caustics and composite shadow boundaries is clearly demonstrated and it is shown that the UGO/EUTD results remain valid and uniformly reduce to the classic results away from the transition regions. Mathematical details on the asymptotic properties and efficient numerical evaluation of the canonical functions involved in the UGO/EUTD expressions are also provided.

Constantinides, E. D.↗

Generative Vulnerability Assessment for Cyber-Physical Systems

Cyber-physical systems (CPS) are highly susceptible to malicious attacks due to their complex dynamics and interconnectivity. A comprehensive understanding of their vulnerabilities is essential for designing effective resilience measures. This paper presents a data-driven attack generative system for evaluating the vulnerability of CPS. The proposed approach formulates the vulnerability assessment problem as determining the feasibility of a specific attack set based on two boundary functions that represent the effectiveness and stealthiness of attacks. The attack generative model is trained using a custom loss function, with two universal approximators designed to learn the effectiveness and stealthiness functions simultaneously. Theoretical results for successful generation and asymptotic convergence of the resulting training algorithm are given. As a result, the proposed approach is evaluated via numerical simulation of an IEEE 14-bus system and gas pipeline systems, demonstrating its viability in learning how to attack nonlinear CPS and identify potential vulnerabilities.

Computer systems organization↗

Stability of large-scale systems with stable and unstable subsystems.

The purpose of this paper is to develop new methods for constructing vector Liapunov functions and broaden the application of Liapunov's theory to stability analysis of large-scale dynamic systems. The application, so far limited by the assumption that the large-scale systems are composed of exponentially stable subsystems, is extended via the general concept of comparison functions to systems which can be decomposed into asymptotically stable subsystems. Asymptotic stability of the composite system is tested by a simple algebraic criterion. With minor technical adjustments, the same criterion can be used to determine connective asymptotic stability of large-scale systems subject to structural perturbations. By redefining the constraints imposed on the interconnections among the subsystems, the considered class of systems is broadened in an essential way to include composite systems with unstable subsystems. In this way, the theory is brought substantially closer to reality since stability of all subsystems is no longer a necessary assumption in establishing stability of the overall composite system.

Grujic, Lj. T.↗

A direct method for the solution of unsteady two-dimensional incompressible Navier-Stokes equations

The unsteady incompressible Navier-Stokes equations are formulated in terms of vorticity and stream function in generalized curvilinear orthogonal coordinates to facilitiate analysis of flow configurations with general geometries. The numerical method developed solves the conservative form of the transport equation using the alternating-direction implicit method, whereas the stream-function equation is solved by direct block Gaussian elimination. The method is applied to a model problem of flow over a back-step in a doubly infinite channel, using clustered conformal coordinates. One-dimensional stretching functions, dependent on the Reynolds number and the asymptotic behavior of the flow, are used to provide suitable grid distribution in the separation and reattachment regions, as well as in the inflow and outflow regions. The optimum grid distribution selected attempts to honor the multiple length scales of the separated-flow model problem. The asymptotic behavior of the finite-differenced transport equation near infinity is examined and the numerical method is carefully developed so as to lead to spatially second-order accurate wiggle-free solutions, i.e., with minimum dispersive error. Results have been obtained in the entire laminar range for the backstep channel and are in good agreement with the available experimental data for this flow problem.

Ghia, K. N.↗

Theory of Ostwald ripening in a two-component system

When a two-component system is cooled below the minimum temperature for its stability, it separates into two or more immiscible phases. The initial nucleation produces grains (if solid) or droplets (if liquid) of one of the phases dispersed in the other. The dynamics by which these nuclei proceed toward equilibrium is called Ostwald ripening. The dynamics of growth of the droplets depends upon the following factors: (1) The solubility of the droplet depends upon its radius and the interfacial energy between it and the surrounding (continuous) phase. There is a critical radius determined by the supersaturation in the continuous phase. Droplets with radii smaller than critical dissolve, while droplets with radii larger grow. (2) The droplets concentrate one component and reject the other. The rate at which this occurs is assumed to be determined by the interdiffusion of the two components in the continuous phase. (3) The Ostwald ripening is constrained by conservation of mass; e.g., the amount of materials in the droplet phase plus the remaining supersaturation in the continuous phase must equal the supersaturation available at the start. (4) There is a distribution of droplet sizes associated with a mean droplet radius, which grows continuously with time. This distribution function satisfies a continuity equation, which is solved asymptotically by a similarity transformation method.

Baird, J. K.↗

On the efficient calculation of the incomplete Airy function with application to edge diffraction

When the geometrical theory of diffraction (GTD) and its extensions are used to calculate fields across transition regions or in focal areas, various special functions are necessary to create uniform representations. One such transition function is the incomplete Airy function, which is applicable when two reflection points are arbitrarily close to each other and the endpoint of the domain considered. An efficient calculation of the incomplete Airy function is presented and shown to equal known asymptotic expansions when its argument is away from the critical point. By using these expansions, uniform field expressions are found which are consistent with the GTD and maintain its computational efficiency.

Cwik, Tom↗

Undamped plasma waves

This paper describes small-amplitude nonlinear plasma wave solutions to the one-dimensional Vlasov-Maxwell equations. A sufficient condition for waves of a given phase velocity to exist arbitrarily close to a given spatially uniform Vlasov equilibrium is developed, and sufficient analytical information for the construction of approximate expressions for the electric potential and distribution functions is derived, with exact knowledge of the asymptotic behavior of the error terms. These results have a very surprising physical implication: the Landau damping of small-amplitude waves is not inevitable. Instead, there exist plasma waves that trap particles even at arbitrarily small amplitude and do not damp.

Holloway, James P.↗

Supersonic flow calculation using a Reynolds-stress and an eddy thermal diffusivity turbulence model

A second-order model for the velocity field and a two-equation model for the temperature field are used to calculate supersonic boundary layers assuming negligible real gas effects. The modeled equations are formulated on the basis of an incompressible assumption and then extended to supersonic flows by invoking Morkovin's hypothesis, which proposes that compressibility effects are completely accounted for by mean density variations alone. In order to calculate the near-wall flow accurately, correction functions are proposed to render the modeled equations asymptotically consistent with the behavior of the exact equations near a wall and, at the same time, display the proper dependence on the molecular Prandtl number. Thus formulated, the near-wall second order turbulence model for heat transfer is applicable to supersonic flows with different Prandtl numbers. The model is validated against flows with different Prandtl numbers and supersonic flows with free-stream Mach numbers as high as 10 and wall temperature ratios as low as 0.3. Among the flow cases considered, the momentum thickness Reynolds number varies from approximately 4,000 to approximately 21,000. Good correlation with measurements of mean velocity, temperature, and its variance is obtained. Discernible improvements in the law-of-the-wall are observed, especially in the range where the big-law applies.

Sommer, T. P.↗

Two Improved Algorithms for Envelope and Wavefront Reduction

Two algorithms for reordering sparse, symmetric matrices or undirected graphs to reduce envelope and wavefront are considered. The first is a combinatorial algorithm introduced by Sloan and further developed by Duff, Reid, and Scott; we describe enhancements to the Sloan algorithm that improve its quality and reduce its run time. Our test problems fall into two classes with differing asymptotic behavior of their envelope parameters as a function of the weights in the Sloan algorithm. We describe an efficient 0(nlogn + m) time implementation of the Sloan algorithm, where n is the number of rows (vertices), and m is the number of nonzeros (edges). On a collection of test problems, the improved Sloan algorithm required, on the average, only twice the time required by the simpler Reverse Cuthill-Mckee algorithm while improving the mean square wavefront by a factor of three. The second algorithm is a hybrid that combines a spectral algorithm for envelope and wavefront reduction with a refinement step that uses a modified Sloan algorithm. The hybrid algorithm reduces the envelope size and mean square wavefront obtained from the Sloan algorithm at the cost of greater running times. We illustrate how these reductions translate into tangible benefits for frontal Cholesky factorization and incomplete factorization preconditioning.

Kumfert, Gary↗

Spectral characteristics of digital phase-modulated signals

The power spectral density of two-level digital phase modulated signals is determined using various polar pulse shaping functions, the only restriction being that the pulses are independent and do not overlap. Rectangular pulses and pulses having finite rise and decay times are considered. It is shown that a substantial improvement in terms of minimizing spectral occupancy on a power basis over a given frequency band can be achieved by using linear and raised cosine pulse shaping functions. Further, the rate of decrease of the power spectral density is given for the asymptotic limit as f goes to infinity. The use of these and possibly other well chosen pulse shaping functions can then provide an additional aid to overall efficient spectrum utilization and management.

Mcgregor, D. N.↗

Embedded function methods for supersonic turbulent boundary layers

The development of embedded functions to represent the mean velocity and total enthalpy distributions in the wall layer of a supersonic turbulent boundary layer is considered. The asymptotic scaling laws (in the limit of large Reynolds number) for high speed compressible flows are obtained to facilitate eventual implementation of the embedded functions in a general prediction method. A self-consistent asymptotic structure is derived, as well as a compressible law of the wall in which the velocity and total enthalpy are logarithmic within the overlap zone, but in the Howarth-Dorodnitsyn variable. Simple outer region turbulence models are proposed (some of which are modifications of existing incompressible models) to reflect the effects of compressibility. As a test of the methodology and the new turbulence models, a set of self-similar outer region profiles is obtained for constant pressure flow; these are then coupled with embedded functions in the wall layer. The composite profiles thus obtained are compared directly with experimental data and good agreement is obtained for flows with Mach numbers up to 10.

He, J.↗