Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent discretization”

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 253 records · Page 14

Non-Markovian relaxation spectroscopy of fluxonium qubits

Recent studies have shown that parasitic two-level systems (TLS) in superconducting qubits, which are a leading source of decoherence, can have relaxation times longer than the qubits themselves. However, the standard techniques used to characterize qubit relaxation is only valid for measuring T 1 under the Born-Markov approximation and could mask environmental memory effects in practice. Here, we introduce two-timescale relaxometry, a technique to probe the qubit and environment relaxation simultaneously and efficiently. We apply it to high-coherence fluxonium qubits over a frequency range of 0.1-0.4 GHz, and reveal a discrete spectrum of TLS with millisecond lifetimes. Our analysis of the spectrum is consistent with a random distribution of TLS in the aluminum oxide tunnel barrier of the Josephson junction chain of the fluxonium, with a spectral and volumetric density and average electric dipole similar to previous TLS studies at much higher frequencies. Our study suggests that investigating and mitigating TLS in the junction chain is crucial to the development of various types of noise-protected qubits in circuit QED.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Catalytic enhancement in the performance of the microscopic two-stroke heat engine

Here we consider a model of heat engine operating in the microscopic regime: the two-stroke engine. It produces work and exchanges heat in two discrete strokes that are separated in time. The working body of the engine consists of two d-level systems initialized in thermal states at two distinct temperatures. Additionally, an auxiliary nonequilibrium system called catalyst may be incorporated with the working body of the engine, provided the state of the catalyst remains unchanged after the completion of a thermodynamic cycle. This ensures that the work produced by the engine arises solely from the temperature difference. Upon establishing the rigorous thermodynamic framework, we characterize twofold improvement stemming from the inclusion of a catalyst. Firstly, we prove that in the noncatalytic scenario, the optimal efficiency of the two-stroke heat engine with a working body composed of two-level systems is given by the Otto efficiency, which can be surpassed by incorporating a catalyst with the working body. Secondly, we show that incorporating a catalyst allows the engine to operate in frequency and temperature regimes that are not accessible for noncatalytic two-stroke engines. We conclude with a general conjecture about the advantage brought by a catalyst: including the catalyst with the working body always allows to improve efficiency over the noncatalytic scenario for any microscopic two-stroke heat engines. We prove this conjecture for two-stroke engines where the working body is composed of two d-level systems initialized in thermal states at two distinct temperatures, as long as the final joint state leading to optimal efficiency in the noncatalytic scenario is not a product state, or at least one of the d-level system is not thermal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Far-ultraviolet imagery of the Barnard Loop Nebula

An electrographic Schmidt camera carried on a sounding rocket has yielded far-ultraviolet (1050-2000 A and 1230-2000 A) images of the Barnard Loop Nebula and of the general background in the Orion region due to scattering of ultraviolet starlight by interstellar dust particles. The total intensity in the Barnard Loop region agrees well with OAO-2 measurements, but the discrete Loop structure contributes only some 15% of the total. The measurements are consistent with a relatively high albedo for the dust grains in the far-ultraviolet.

Carruthers, G. R.↗

Radiative cooling efficiencies and predicted spectra of species of the Io plasma torus

Calculations of the physical condition of the Io plasma torus have been made based on the recent Voyager EUV observations. The calculations represent an assumed thin plasma collisional ionization equilibrium among the states within each species. The observations of the torus are all consistent with this condition. The major energy loss mechanism is radiative cooling in discrete transitions. Calculations of radiative cooling efficiencies of the identified species leads to an estimated energy loss rate of at least 1.5 x 10 to the 12th watts. The mean electron temperature and density of the plasma are estimated to be 100,000 K and 2100/cu cm. The estimated number densities of S III, S IV, and O III are roughly 95, 80, and 190-740/cu cm. Upper limits have been placed on a number of other species based on the first published Voyager EUV spectrum of the torus. The assumption that energy is supplied to the torus through injection of neutral particles from Io leads to the conclusion that ion loss rates are controlled by diffusion, and relative species abundances consequently are not controlled by collisional ionization equilibrium.

Shemansky, D. E.↗

Unconditional instability of inflow dependent boundary conditions in difference approximations to hyperbolic systems

The stability of finite difference approximations to initial boundary hyperbolic systems is studied. As is well known, a proper specification of boundary conditions for such systems is essential for their solutions to be well defined. A discrete analogue of the above is proved - if the numerical boundary conditions are consistent with an inflow part of the problem, they render the overall computation unstable. An example of the inviscid gasdynamics equations is considered.

Tadmor, E.↗

The unconditional instability of inflow-dependent boundary conditions in difference approximations to hyperbolic systems

The stability of finite difference approximations to initial boundary hyperbolic systems is studied. As is well known, a proper specification of boundary conditions for such systems is essential for their solutions to be well defined. A discrete analogue of the above is proved - if the numerical boundary conditions are consistent with an inflow part of the problem, they render the overall computation unstable. An example of the inviscid gasdynamics equations is considered. Previously announced in STAR as N81-33874

Tadmor, E.↗

Calibration of the 7- to 14-micron brightness spectra of Uranus and Neptune

An independent measure of the absolute brightness of the disk-averaged spectrum measured by Orton et al. (1987) is attempted via observations of Uranus and Neptune at discrete wavelengths between 8.17 and 13.0 microns. These observations are found to be consistent with the suggested upward recalibrations of the Uranus and Neptune spectra by factors of 1.6 and 1.12, respectively. Both graphic and tabular results are presented.

Orton, Glenn S.↗

Application of the boundary element method to transient heat conduction

An advanced boundary element method (BEM) is presented for the transient heat conduction analysis of engineering components. The numerical implementation necessarily includes higher-order conforming elements, self-adaptive integration and a multiregion capability. Planar, three-dimensional and axisymmetric analyses are all addressed with a consistent time-domain convolution approach, which completely eliminates the need for volume discretization for most practical analyses. The resulting general purpose algorithm establishes BEM as an attractive alternative to the more familiar finite difference and finite element methods for this class of problems. Several detailed numerical examples are included to emphasize the accuracy, stability and generality of the present BEM. Furthermore, a new efficient treatment is introduced for bodies with embedded holes. This development provides a powerful analytical tool for transient solutions of components, such as casting moulds and turbine blades, which are cumbersome to model when employing the conventional domain-based methods.

Dargush, G. F.↗

Sensitivity analysis for the dynamic response of viscoplastic shells of revolution

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric response of viscoplastic shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. Spatial discretization is performed using finite elements, with discontinuous stress resultants across element interfaces. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) with an implicit constitutive update. The sensitivity coefficients are evaluated using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading, and a circular plate subjected to impulsive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to thickness, mass density, Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients are presented, along with spatial distributions of these quantities at selected times.

Kulkarni, Makarand↗

On spurious behavior of CFD simulations

Spurious behavior in underresolved grids and/or semi-implicit temporal discretizations for four computational fluid dynamics (CFD) simulations are studied. The numerical simulations consist of (a) a 1-D chemically relaxed nonequilibrium model, (b) the direct numerical simulation (DNS) of 2-D incompressible flow over a backward facing step, (c) a loosely-coupled approach for a 2-D fluid-structure interaction, and (d) a 3-D compressible unsteady flow simulation of vortex breakdown in delta wings. Using knowledge from dynamical systems theory, various types of spurious behaviors that are numerical artifacts were systematically identified. These studies revealed the various possible dangers of misinterpreting numerical simulation of realistic complex flows that are constrained by the available computing power. In large scale computations underresolved grids, semi-implicit procedures, loosely-coupled implicit procedures, and insufficiently long time integration in DNS are most often unavoidable. Consequently, care must be taken in both computation and in interpretation of the numerical data. The results presented confirm the important role that dynamical systems theory can play in the understanding of the nonlinear behavior of numerical algorithms and in aiding the identification of the sources of numerical uncertainties in CFD.

Numerical Solution↗

Titan's Atmospheric Dynamics and Meteorology

Titan, after Venus, is the second example of an atmosphere with a global cyclostrophic circulation in the solar system, but a circulation that has a strong seasonal modulation in the middle atmosphere. Direct measurement of Titan's winds, particularly observations tracking the Huygens probe at 10degS, indicate that the zonal winds are generally in the sense of the satellites rotation. They become cyclostrophic approx. 35 km above the surface and generally increase with altitude, with the exception of a sharp minimum centered near 75 km, where the wind velocity decreases to nearly zero. Zonal winds derived from the temperature field retrieved from Cassini measurements, using the thermal wind equation, indicate a strong winter circumpolar vortex, with maximum winds at mid northern latitudes of 190 ms-' near 300 km. Above this level, the vortex decays. Curiously, the zonal winds and temperatures are symmetric about a pole that is offset from the surface pole by approx.4 degrees. The cause of this is not well understood, but it may reflect the response of a cyclostrophic circulation to the offset between the equator, where the distance to the rotation axis is greatest, and the solar equator. The mean meridional circulation can be inferred from the temperature field and the meridional distribution of organic molecules and condensates and hazes. Both the warm temperatures in the north polar region near 400 km and the enhanced concentration of several organic molecules suggests subsidence there during winter and early spring. Stratospheric condensates are localized at high northern latitudes, with a sharp cut-off near 50degN. Titan's winter polar vortex appears to share many of the same characteristics of winter vortices on Earth-the ozone holes. Global mapping of temperatures, winds, and composition in he troposphere, by contrast, is incomplete. The few suitable discrete clouds that have bee found for tracking indicate smaller velocities than aloft, consistent with the Huygens measurements. At low latitudes the zonal winds near the surface appear not to be westward as on Earth, but eastward. Because the net zonal-mean time-averaged torq exerted by the surface on the atmosphere should vanish, this implies westward flow o part of the surface; the question is where. The latitude contrast in tropospheric temperatures, deduced from radio occultations at low, mid, and high latitudes, is small approx.5 K at the tropopause and approx.3 K at the surface.

Flasar, F. M.↗

State-Chart Autocoder

A computer program translates Unified Modeling Language (UML) representations of state charts into source code in the C, C++, and Python computing languages. ( State charts signifies graphical descriptions of states and state transitions of a spacecraft or other complex system.) The UML representations constituting the input to this program are generated by using a UML-compliant graphical design program to draw the state charts. The generated source code is consistent with the "quantum programming" approach, which is so named because it involves discrete states and state transitions that have features in common with states and state transitions in quantum mechanics. Quantum programming enables efficient implementation of state charts, suitable for real-time embedded flight software. In addition to source code, the autocoder program generates a graphical-user-interface (GUI) program that, in turn, generates a display of state transitions in response to events triggered by the user. The GUI program is wrapped around, and can be used to exercise the state-chart behavior of, the generated source code. Once the expected state-chart behavior is confirmed, the generated source code can be augmented with a software interface to the rest of the software with which the source code is required to interact.

Clark, Kenneth↗

Navier-Stokes Dynamics by a Discrete Boltzmann Model

This work investigates the possibility of particle-based algorithms for the Navier-Stokes equations and higher order continuum approximations of the Boltzmann equation; such algorithms would generalize the well-known Pullin scheme for the Euler equations. One such method is proposed in the context of a discrete velocity model of the Boltzmann equation. Preliminary results on shock structure are consistent with the expectation that the shock should be much broader than the near discontinuity predicted by the Pullin scheme, yet narrower than the prediction of the Boltzmann equation. We discuss the extension of this essentially deterministic method to a stochastic particle method that, like DSMC, samples the distribution function rather than resolving it completely.

Rubinstein, Robet↗

Chapter 13. Atmospheric Dynamics and Meteorology

Titan, after Venus, is the second example in the solar system of an atmosphere with a global cyclostrophic circulation, but in this case a circulation that has a strong seasonal modulation in the middle atmosphere. Direct measurement of Titan's winds, particularly observations tracking the Huygens probe at 10 deg S, indicate that the zonal winds are mostly in the sense of the satellite's rotation. They generally increase with altitude and become cyclostrophic near 35 km above the surface. An exception to this is a sharp minimum centered near 75 km, where the wind velocity decreases to nearly zero. Zonal winds derived from temperatures retrieved from Cassini orbiter measurements, using the thermal wind equation, indicate a strong winter circumpolar vortex, with maximum winds of 190 m/s at mid northern latitudes near 300 km. Above this level, the vortex decays. Curiously, the stratospheric zonal winds and temperatures in both hemispheres are symmetric about a pole that is offset from the surface pole by about 4 deg. The cause of this is not well understood, but it may reflect the response of a cyclostrophic circulation to the onset between the equator, where the distance to the rotation axis is greatest, and the seasonally varying subsolar latitude. The mean meridional circulation can be inferred from the temperature field and the meridional distribution of organic molecules and condensates and hazes. Both the warm temperatures near 400 km and the enhanced concentration of several organic molecules suggest subsidence in the north polar region during winter and early spring. Stratospheric condensates are localized at high northern latitudes, with a sharp cut-off near 50 deg N. Titan's winter polar vortex appears to share many of the same characteristics of isolating high and low-latitude air masses as do the winter polar vortices on Earth that envelop the ozone holes. Global mapping of temperatures, winds, and composition in the troposphere, by contrast, is incomplete. The few suitable discrete clouds that have been found for tracking indicate smaller velocities than aloft, consistent: with the Huygens measurements, Along the descent trajectory, the Huygens measurements indicate eastward zonal winds down to 7 km, where they shift westward, and then eastward again below 1 km dawn to the surface. The low-latitude dune fields seen in Cassini RADAR images have been interpreted as longitudinal dunes occurring in a mean eastward zonal wind. This is not like Earth, where the low-latitude winds are westward above the surface. Because the net zonal-mean time-averaged torque exerted by the surface on the atmosphere should vanish, there must be westward flow over part of the surface; the question is where and when. The meridional contrast in tropospheric temperatures deduced from radio occultations at low, mid, and high latitudes. is small, approximately 5 K at the tropopause and approximately 3 K at the surface. This implies efficient heat transport, probably by axisymmetric meridional circulations. The effect of the methane "hydrological" cycle on the atmospheric circulation is not well constrained by existing measurements, Understanding the mature of the surface-atmosphere coupling will be critical to elucidating the atmospheric transports of momentum, heat, and volatiles.

Flasar, F. M.↗

Helicopter Pilot Performance for Discrete-maneuver Flight Tasks

This paper describes a current study of several basic helicopter flight maneuvers. The data base consists of in-flight measurements from instrumented helicopters using experienced pilots. The analysis technique is simple enough to apply without automatic data processing, and the results can be used to build quantitative matah models of the flight task and some aspects of the pilot control strategy. In addition to describing the performance measurement technqiue, some results are presented which define the aggressiveness and amplitude of maneuvering for several lateral maneuvers including turns and sidesteps.

Heffley, R. K.↗

Optimal application of Morrison's iterative noise removal for deconvolution. Appendices

Morrison's iterative method of noise removal, or Morrison's smoothing, is applied in a simulation to noise-added data sets of various noise levels to determine its optimum use. Morrison's smoothing is applied for noise removal alone, and for noise removal prior to deconvolution. For the latter, an accurate method is analyzed to provide confidence in the optimization. The method consists of convolving the data with an inverse filter calculated by taking the inverse discrete Fourier transform of the reciprocal of the transform of the response of the system. Various length filters are calculated for the narrow and wide Gaussian response functions used. Deconvolution of non-noisy data is performed, and the error in each deconvolution calculated. Plots are produced of error versus filter length; and from these plots the most accurate length filters determined. The statistical methodologies employed in the optimizations of Morrison's method are similar. A typical peak-type input is selected and convolved with the two response functions to produce the data sets to be analyzed. Both constant and ordinate-dependent Gaussian distributed noise is added to the data, where the noise levels of the data are characterized by their signal-to-noise ratios. The error measures employed in the optimizations are the L1 and L2 norms. Results of the optimizations for both Gaussians, both noise types, and both norms include figures of optimum iteration number and error improvement versus signal-to-noise ratio, and tables of results. The statistical variation of all quantities considered is also given.

Ioup, George E.↗

Computation of rotor aerodynamic loads with a constant vorticity contour free wake model

An analytical method is presented which facilitates the study of isolated rotors with an improved approach to wake simulation. Vortex filaments are simulated along contours of constant sheet strength for the sheet of vorticity resulting from each rotor blade. Curved vortex elements comprise the filaments which can be distorted by the local velocity field. Called the Constant Vorticity Contour wake model, the approach permits the simulation of the blades' wakes corresponding to the full span of the rotor blade. The discretization of the wake of the rotor blade produces spacing and structure that are consistent with the spatial and temporal variations in the loading. A vortex-lattice aerodynamic model of the blade is also included which introduces a finite-element structural model of the blade and consideration of the force and moment trim analysis. Results of the present version of the simulation, called RotorCRAFT, are found to correlate well with H-34 flight-test data.

Quackenbush, Todd R.↗

Effect of mesh distortion on the accuracy of transverse shear stresses and their sensitivity coefficients in multilayered composites

A study is made of the effect of mesh distortion on the accuracy of transverse shear stresses and their first-order and second-order sensitivity coefficients in multilayered composite panels subjected to mechanical and thermal loads. The panels are discretized by using a two-field degenerate solid element, with the fundamental unknowns consisting of both displacement and strain components, and the displacement components having a linear variation throughout the thickness of the laminate. A two-step computational procedure is used for evaluating the transverse shear stresses. In the first step, the in-plane stresses in the different layers are calculated at the numerical quadrature points for each element. In the second step, the transverse shear stresses are evaluated by using piecewise integration, in the thickness direction, of the three-dimensional equilibrium equations. The same procedure is used for evaluating the sensitivity coefficients of transverse shear stresses. Numerical results are presented showing no noticeable degradation in the accuracy of the in-plane stresses and their sensitivity coefficients with mesh distortion. However, such degradation is observed for the transverse shear stresses and their sensitivity coefficients. The standard of comparison is taken to be the exact solution of the three-dimensional thermoelasticity equations of the panel.

Noor, Ahmed K.↗