Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

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 505 records · Page 28

Characterization of DBD Plasma Actuators Performance without External Flow : Thrust-Voltage Quadratic Relationship in Logarithmic Space for Sinusoidal Excitation - Part I

We present results of thrust measurements of Dielectric Barrier Discharge (DBD) plasma actuators. We have used a test setup, measurement, and data processing methodology that we developed in prior work. The tests were conducted with High Density Polyethylene (HDPE) actuators of three thicknesses. The applied voltage driving the actuators was a pure sinusoidal waveform. The test setup was suspended actuators with a partial liquid interface. The tests were conducted at low ambient humidity. The thrust was measured with an analytical balance and the results were corrected for anti-thrust to isolate the plasma generated thrust. Applying this approach resulted in smooth and repeatable data. It also enabled curve fitting that yielded quadratic relations between the plasma thrust and voltage in log-log space at constant frequencies. The results contrast power law relationships developed in literature that appear to be a rough approximation over a limited voltage range.

Dielectric Barrier Discharge↗

SiAl Alloy Feedhorn Arrays: Material Properties, Feedhorn Design, and Astrophysical Applications

We present here a study of the use of the SiAl alloy CE7 for the packaging of silicon devices at cryogenic temperatures. We report on the development of baseplates and feedhorn arrays for millimeter wave bolometric detectors for astrophysics. Existing interfaces to such detectors are typically made either of metals, which are easy to machine but mismatched to the thermal contraction profile of Si devices, or of silicon, which avoids the mismatch but is difficult to directly machine. CE7 exhibits properties of both Si and Al, which makes it uniquely well suited for this application. We measure CE7 to a) superconduct below a critical transition temperature, T (sub c), 1.2 K, b) have a thermal contraction profile much closer to Si than metals, which enables simple mating, and c) have a low thermal conductivity which can be improved by Au-plating. Our investigations also demonstrate that CE7 can be machined well enough to fabricate small structures, such as #0-80 threaded holes, to tight tolerances (approximately 25 microns) in contrast with pure silicon and similar substrates. We have fabricated CE7 baseplates being deployed in the 93 gigahertz polarimetric focal planes used in the Cosmology Large Angular Scale Surveyor (CLASS). We also report on the development of smooth-walled feedhorn arrays made of CE7 that will be used in a focal plane of dichroic 150/220 gigahertz detectors for the CLASS High-Frequency camera.

Ali, Aamir M.↗

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

BCS-BEC crossover in a (t 2g ) 4 excitonic magnet

The condensation of spin-orbit-induced excitons in t 2 g 4 electronic systems is attracting considerable attention. At large Hubbard U , antiferromagnetism was proposed to emerge from the Bose-Einstein Condensation (BEC) of triplons ( J eff = 1 ). Furthermore, we show that even at intermediate U regimes, the spin-orbit exciton condensation is possible leading also to staggered magnetic order. The canonical electron-hole excitations (excitons) transform into local triplon excitations at large U , and this BEC strong coupling regime is smoothly connected to the intermediate U excitonic insulator region. We solved the degenerate three-orbital Hubbard model with spin-orbit coupling ( λ ) in one dimension using the density matrix renormalization group, while in two dimensions we use the Hartree-Fock approximation (HFA). Employing these techniques, we provide the full λ versus U phase diagrams for both one- and two-dimensional lattices. Our main result is that at intermediate Hubbard U , increasing λ at fixed U the system transitions from an incommensurate spin-density-wave metal to a Bardeen-Cooper-Schrieffer (BCS) excitonic insulator, with coherence length r coh of O ( a ) and O ( 10 a ) in 1 d and 2 d , respectively, with a being the lattice spacing. Further increasing λ , the system eventually crosses over to the BEC limit (with r coh ≪ a ).

1-dimensional systems↗

Tuning the Interpolation Basis in a Multigrid Decomposition for Local Error Control

In the compression of scientific data, error-controlled compressors enable to considerably decrease the size of the dataset while maintaining adequate levels of accuracy. In this paper, we note that multi-level refactoring scheme such as MGARD i) rely on an approximation of the data based on the interpolation of coefficients, ii) estimate the resulting error with global metrics on the dataset. To improve on these two aspects, we propose a method that aims to divide the original dataset into blocks based on their smoothness and refactors each block separately with the most relevant interpolation order. We show the relevance of such a method on tailored datasets and the benefits and challenges when applying it to large scientific data.

Vidal, Nicolas [ORNL]↗

X-ray observations of the southern cluster CA 0340-538 and the Horologium supercluster

X-ray observations of the Horologium supercluster are examined. Results show that the optically Coma-like cluster CA 0340-538 has a smooth X-ray structure which also resembles that of Coma. The density structure of its X-ray emitting gas, determined in a model-independent way, is found to be remarkably featureless, being consistent with a nearly isothermal polytrope out to approximately 2 Mpc from the cluster center. The total mass of the gas is determined to exceed 10% of the cluster virial mass and may contribute a substantial fraction of the cluster binding mass. Although CA 0329-527 is found to be an X-ray double cluster, its structure appears to be confused. CA 0325-539 is found to show no substantial X-ray emission. An upper limit of 1.3 x 10 to the 44th erg/s is determined for any diffuse 2-10 keV emission from the Horologium supercluster.

Ku, W. H.-M.↗

Tidal interactions between spherical galaxies

The impulse approximation and full N-body experiments are used to investigate the changes in mass and binding energy which result from gravitational encounters between spherical galaxies. Models with King (1966) and de Vaucouleurs (1948) surface density profiles, and radial, tangential, and isotropic velocity distributions. Results based on the impulse approximation depend on two collision parameters (the impact parameter and the ratio of the perturbing mass to the relative velocity at infinity); they agree well with N-body results. The asymptotic limit, valid for distant encounters, is found to approximate the correct results only for impact parameters larger than 10 effective radii. It is therefore of little use in most situations of interest. No simple 'linear' or 'quadratic' dependence of energy and mass loss on collision strength is found, but rather, a smooth intermediate behavior. Cross sections depend very little on the internal velocity structure of the galaxy, but quite strongly on the surface density profile in its outer parts. The reasons why these results disagree with the conclusions of previous workers are discussed.

Aguilar, L. A.↗

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid interface shape and convection patterns during solidification of pure Al and Al-4 wt. % Cu alloy. Steady state calculations were performed for different gravity levels and orientations. Furthermore, transient modeling was undertaken to investigate effect of the solidification velocity. The furnace configuration used in this analysis is the proposed International Space Station Furnace. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. The Solidification model was adopted from previous work and was based on the finite element code FIDAP. Thermocouple data and quenched interface shape from a pure Al sample flown on the Life Sciences and Microgravity Spacelab (LMS) mission, July 1996 was used for model validation. Good agreement was obtained between the predicted interface shape and that measured from the quenched LMS sample. It was found that the imposed temperature boundary condition must be sufficiently smooth for the heat fluxes in the model to be self-consistent. The model predicted that alloy sample Al-Cu is more sensitive to variations in the gravity level compared to pure Al. For Al-Cu alloy, solute diffusivity is approximately 10(exp 4) times smaller than thermal diffusivity of pure Al. Hence, in a microgravity environment the weak convection has no measurable effect on the heat fluxes, but is still strong enough to affect concentration distribution. Since the alloy melting temperature is determined by concentration the interface shape depends on the level of convection. This was found not to be the case for pure metal.

Bune, Andris V.↗

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗

Machine learning unifies flexibility and efficiency of spinodal structure generation for stochastic biomaterial design

Abstract Porous biomaterials design for bone repair is still largely limited to regular structures (e.g. rod-based lattices), due to their easy parameterization and high controllability. The capability of designing stochastic structure can redefine the boundary of our explorable structure–property space for synthesizing next-generation biomaterials. We hereby propose a convolutional neural network (CNN) approach for efficient generation and design of spinodal structure—an intriguing structure with stochastic yet interconnected, smooth, and constant pore channel conducive to bio-transport. Our CNN-based approach simultaneously possesses the tremendous flexibility of physics-based model in generating various spinodal structures (e.g. periodic, anisotropic, gradient, and arbitrarily large ones) and comparable computational efficiency to mathematical approximation model. We thus successfully design spinodal bone structures with target anisotropic elasticity via high-throughput screening, and directly generate large spinodal orthopedic implants with desired gradient porosity. This work significantly advances stochastic biomaterials development by offering an optimal solution to spinodal structure generation and design.

59 BASIC BIOLOGICAL SCIENCES↗

Stellar Abundance Maps of the Milky Way Disk

To understand the formation of the Milky Way's prominent bar it is important to know whether stars in the bar differ in the chemical element composition of their birth material as compared to disk stars. This requires stellar abundance measurements for large samples across the Milky Way's body. Such samples, e.g., luminous red giant stars observed by the Sloan Digital Sky Survey's APOGEE survey, will inevitably span a range of stellar parameters; as a consequence, both modeling imperfections and stellar evolution may preclude consistent and precise estimates of their chemical composition at a level of purported bar signatures, which has left current analyses of a chemically distinct bar inconclusive. Here, we develop a new self-calibration approach to eliminate both modeling and astrophysical abundance systematics among red giant branch (RGB) stars of different luminosities (and hence surface gravity $\mathrm{log} g$). We apply our method to 48,853 luminous APOGEE Data Release 16 RGB stars to construct spatial abundance maps of 20 chemical elements near the Milky Way's mid-plane, covering galactocentric radii of 0 kpc < R GC < 20 kpc. Our results indicate that there are no abundance variations whose geometry matches that of the bar, and that the mean abundance gradients vary smoothly and monotonically with galactocentric radius. We confirm that the high-α disk is chemically homogeneous, without spatial gradients. Furthermore, we present the most precise [Fe/H] versus R GC gradient to date with a slope of – 0.057 ±0.001 dex kpc –1 out to approximately 15 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗

Production of polydisperse sprays.

A device which produces a spray with a controllable drop size distribution is described. This system is especially well suited for generating nonsmooth (e.g., bimodal) distributions, although smooth distributions can also be produced. Droplets of uniform size result from the rapid growth of an oscillatory disturbance in a free liquid jet. Simultaneous generation of droplets of various diameters by this method produces the polydisperse spray. The range of drop diameters is approximately 300-3000 microns.

Pierce, T. H.↗

Two-dimensional magnetohydrodynamic turbulence - Cylindrical, non-dissipative model

Incompressible magnetohydrodynamic turbulence is treated in the presence of cylindrical boundaries which are perfectly conducting and rigidly smooth. The model treated is non-dissipative and two-dimensional, the variation of all quantities in the axial direction being ignored. Equilibrium Gibbs ensemble predictions are explored assuming the constraint of constant axial current (appropriate to tokamak operation). No small-amplitude approximations are made. The expectation value of the turbulent kinetic energy is found to approach zero for the state of maximum mean-square vector potential to energy ratio. These are the only states for which large velocity fluctuations are not expected.

Montgomery, D.↗

Determination of wind from NIMBUS 6 satellite sounding data

Objective methods of computing upper level and surface wind fields from NIMBUS 6 satellite sounding data are developed. These methods are evaluated by comparing satellite derived and rawinsonde wind fields on gridded constant pressure charts in four geographical regions. Satellite-derived and hourly observed surface wind fields are compared. Results indicate that the best satellite-derived wind on constant pressure charts is a geostrophic wind derived from highly smoothed fields of geopotential height. Satellite-derived winds computed in this manner and rawinsonde winds show similar circulation patterns except in areas of small height gradients. Magnitudes of the standard deviation of the differences between satellite derived and rawinsonde wind speeds range from approximately 3 to 12 m/sec on constant pressure charts and peak at the jet stream level. Fields of satellite-derived surface wind computed with the logarithmic wind law agree well with fields of observed surface wind in most regions. Magnitudes of the standard deviation of the differences in surface wind speed range from approximately 2 to 4 m/sec, and satellite derived surface winds are able to depict flow across a cold front and around a low pressure center.

Carle, W. E.↗

Theory for computing the field scattered from a smooth inflected surface

A theory is described for computing the reflected or scattered field from a smooth body with inflection points. These inflections occur in certain directions at each surface point for which the total (Gaussian) curvature is zero or negative. For surface illumination in one of these critical directions, the usual reflection formulas obtained by the high-frequency approximation are inapplicable, and a shadow zone exists in the reflected field. Scattering into the shadow zone is treated, as well as specular reflection. This theory should have a variety of applications such as for certain optics problems, computer graphics modeling of three-dimensional shapes, and the design and analysis of specialized microwave reflector antennas.

Barger, R. L.↗

Comparison of Hirs' equation of Moody's equation for determining rotordynamic coefficients of annular pressure seals

The rotordynamic coefficients of an incompressible-flow annular pressure seal were determined using a bulk-flow model in conjunction with two different friction factor relationships. The first, Hirs' equation, assumes the friction factor is a function of Reynolds number only. The second, Moody's equation, approximates Moody's diagram and assumes the friction factor is a function of both Reynolds number and relative roughness. For each value of relative roughness, Hirs' constants were determined so that both equations gave the same magnitude and slope of the friction factor. For smooth seals, both relationships give the same results. For rough seals (e/2 H sub 0 = 0.05) Moody's equation predicts 44% greater direct stiffness, 35% greater cross-coupled stiffness, 19% smaller cross-coupled damping, 59% smaller cross-coupled inertia, and nominally the same direct damping and direct inertia.

Nelson, Clayton C.↗

Comparison of Hirs' equation with Moody's equation for determining rotordynamic coefficients of annular pressure seals

The rotordynamic coefficients of an incompressible-flow annular pressure seal were determined using a bulk-flow model in conjunction with two different friction factor relationships. The first, Hirs' equation, assumes the friction factor is a function of Reynolds number only. The second, Moody's equation, approximates Moody's diagram and assumes the friction factor is a function of both Reynolds number and relative roughness. For each value of relative roughness, Hirs' constants were determined so that both equations gave the same magnitude and slope of the friction factor. For smooth seals, both relationships give the same results. For rough seals, Moody's equation predicts 44 percent greater direct stiffness, 35 percent greater cross-coupled stiffness, 19 percent smaller cross-coupled damping, 59 percent smaller cross-coupled inertia, and nominally the same direct damping and direct inertia.

Nelson, C. C.↗

Three-dimensional multigrid algorithms for the flux-split Euler equations

The Full Approximation Scheme (FAS) multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each of the splitting algorithms uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computation required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined.

Anderson, W. Kyle↗