Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shape 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 73 records · Page 4

Adversarial sampling of unknown and high-dimensional conditional distributions

Many engineering problems require the prediction of realization-to-realization variability or a refined description of modeled quantities. In that case, it is necessary to sample elements from unknown high-dimensional spaces with possibly millions of degrees of freedom. While there exist methods able to sample elements from probability density functions (PDF) with known shapes, several approximations need to be made when the distribution is unknown. In this paper the sampling method, as well as the inference of the underlying distribution, are both handled with a data-driven method known as generative adversarial networks (GAN), which trains two competing neural networks to produce a network that can effectively generate samples from the training set distribution. In practice, it is often necessary to draw samples from conditional distributions. When the conditional variables are continuous, only one (if any) data point corresponding to a particular value of a conditioning variable may be available, which is not sufficient to estimate the conditional distribution. This work handles this problem using an a priori estimation of the conditional moments of a PDF. Herein, two approaches, stochastic estimation, and an external neural network are compared for computing these moments; however, any preferred method can be used. The algorithm is demonstrated in the case of the deconvolution of a filtered turbulent flow field. It is shown that all the versions of the proposed algorithm effectively sample the target conditional distribution with minimal impact on the quality of the samples compared to state-of-the-art methods. Additionally, the procedure can be used as a metric for the diversity of samples generated by a conditional GAN (cGAN) conditioned with continuous variables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗

Exploring transport mechanisms in atomic precision advanced manufacturing enabled pn junctions

We investigate the different transport mechanisms that can occur in pn junction devices made using atomic precision advanced manufacturing at temperatures ranging from cryogenic to room temperature. We first elucidate the potential cause of the anomalous behavior observed in the forward-bias response of these devices in recent cryogenic temperature measurements, which deviates from the theoretical response of a silicon Esaki diode. These anomalous behaviors include current suppression at low voltages in the forward-bias response and a much lower valley voltage at cryogenic temperatures than theoretically expected for a silicon diode. To investigate the potential causes of these anomalies, we studied the effects of a few possible transport mechanisms, including band-to-band tunneling, bandgap narrowing, potential impact of non-Ohmic contacts, band quantization, impact of leakage, and inelastic trap-assisted tunneling, through semi-classical simulations. We find that a combination of two sets of band-to-band tunneling (BTBT) parameters can qualitatively approximate the shape of the tunneling current at low bias. This can arise from band quantization and realignment due to the strong potential confinement in δ-layers. We also find that the lower-than-theoretically-expected valley voltage can be attributed to modifications in the electronic band structure within the δ-layer regions, leading to a significant bandgap narrowing induced by the high density of dopants. Finally, we extend our analyses to room temperature operation and predict that trap-assisted tunneling (TAT) facilitated by phonon interactions may become significant, leading to a complex superposition of BTBT and TAT transport mechanisms in the electrical measurements.

band gap↗

Progress towards understanding ultranonlocality through the wave-vector and frequency dependence of approximate exchange-correlation kernels

In the framework of time-dependent density functional theory (TDDFT), the exact exchange-correlation (xc) kernel f xc (n, q, ω) determines the ground-state energy, excited-state energies, lifetimes, and the time-dependent linear density response of any many-electron system. The recently developed MCP07 xc kernel f xc (n, q, ω) of A. Ruzsinszky et al. [Phys. Rev. B 101, 245135 (2020)] yields excellent uniform electron gas (UEG) ground-state energies and plausible plasmon lifetimes. As MCP07 is constructed to describe f xc of the UEG, it cannot capture optical properties of real materials. To verify this claim, we follow Nazarov et al. [Phys. Rev. Lett. 102, 113001 (2009)] to construct the long-range, dynamic xc kernel, lim q→0 f xc (n, q, ω) = -α(ω)e 2 /q 2 , of a weakly inhomogeneous electron gas, using MCP07 and other common xc kernels. The strong wavevector and frequency dependence of the “ultranonlocality” coefficient α(ω) is demonstrated for a variety of simple metals and semiconductors. We examine how imposing exact constraints on an approximate kernel shapes α(ω). Comparisons to kernels derived from correlated-wavefunction calculations are drawn.

36 MATERIALS SCIENCE↗

Robustness of the Galactic Center excess morphology against masking

The Galactic Center excess (GCE) remains an enduring mystery, with leading explanations being annihilating dark matter or an unresolved population of millisecond pulsars. Analyzing the morphology of the GCE provides critical clues to identify its exact origin. We investigate the robustness of the inferred GCE morphology against the effects of masking, an important step in the analysis where the gamma-ray emission from point sources and the galactic disk are excluded. Using different masks constructed from Fermi point source catalogs and a wavelet method, we find that the GCE morphology, particularly its ellipticity and cuspiness, is relatively independent of the choice of mask for energies above 2–3 GeV. The GCE morphology systematically favors an approximately spherical shape, as expected for dark matter annihilation. Compared to various stellar bulge profiles, a spherical dark matter annihilation profile better fits the data across different masks and galactic diffuse emission backgrounds, except for the stellar bulge profile which provides a similar fit to the data. Modeling the GCE with two components, one from dark matter annihilation and one tracing the Coleman bulge, we find this two-component model outperforms any single component or combinations of dark matter annihilation and other stellar bulge profiles. Uncertainty remains about the exact fraction contributed by each component across different background models and masks. Furthermore, when the Coleman bulge dominates, its corresponding spectrum lacks characteristics typically associated with millisecond pulsars, suggesting that it mostly models the emission from other sources instead of the GCE that is still present and spherically symmetric.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Surrogate Distributed Radiological Sources—Part II: Aerial Measurement Campaign

In this second part of a multipaper series, we present results from outdoor aerial measurements of surrogate distributed gamma-ray sources. Here we detail the design, manufacture, and testing of 300 individual ~7 mCi Cu-64 sealed sources at the Washington State University (WSU) research reactor and their deployment in various source patterns (each comprising up to 100 point sources) during the aerial measurement campaign. We show the results of two such measurements, in which approximate source shapes and qualitative source intensities can be seen from the count rate versus position plots, even without performing reconstructions. We also detail our efforts in ground-truthing the deployed sources and comparing measured gamma-ray data to model predictions. In particular, we compare measured versus expected count data using the Poisson deviance formalism of Part I to evaluate whether the fielded surrogate point-source arrays “look like” their truly continuous distributed source analogs. More generally, we find that the point-source array technique provides high source placement accuracy, relative ease of quantifying the true source configuration, scalability to source dimensions of ≲100 m, ease of reconfiguration and removal, and relatively low dose to personnel. Finally, we consider potential improvements and generalizations of the point-source array technique for future measurement campaigns.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Multi-Factor-Coupled, Ahead-of-Time Aggregation of Power Flexibility Under Forecast Uncertainty

The increasing penetration of distributed energy resources (DERs) is significantly reshaping the role of distribution systems under active energy management. To aggregate the active-reactive power flexibility of DERs dispersed at the feeder and provide capacity support to the transmission system, it is essential to efficiently identify feasible substation power injection trajectories. This paper introduces a novel ahead-of-time flexibility characterization method to address it. First, a polyhedral non-feeder-level power flexibility region (PFR) is constructed, accounting for various time-dependent, power-coupled, and forecast error uncertainties. Then, a polyhedral feeder-level PFR is analytically derived through a coordinate transformation, which can reveal the uncertainty propagation path, i.e., how uncertainty applies to the feeder-level PFR. To facilitate the high-level application, a tractable chance-constrained Chebyshev centering optimization model is further developed to find a ball-shaped inner approximation of the feeder-level PFR. Finally, the proposed method is validated on a modified IEEE 123-bus test system. Here, both theoretical and experimental results show that, with appropriate robustness parameter settings, the proposed method can make the approximated PFR less conservative with abundant robustness against forecast error uncertainty.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Implementation of the Glued Sphere Discrete Element Model for Non-Spherical Particles in MFiX Software

To enhance solver capabilities, simulation flexibility and model validation within the MFiX software, the U.S. Department of Energy (DOE) is funding efforts to develop and integrate the glued-sphere discrete element method into the latest version of MFiX as a dedicated computational module. The glued-sphere discrete element method is a numerical technique to depict the behavior of non-spherical particles in granular flows or particulate systems by representing them as a collection of component spheres. These spheres are bonded together to approximate the shape and mechanical/chemical properties of a more complex particle. The method effectively reuses the existing sphere-sphere collision algorithm, interphase momentum and heat transfer calculations utilized in the traditional discrete element method, extending these capabilities to non-spherical particles. Additionally, this method explicitly resolves intra-particle temperature and species distributions. The MFiX glued-sphere computational module includes tools for generating glued sphere configurations, a dedicated solver, and visualization capabilities in post-processing. More specifically within the computational module, collision detection and calculations were first performed on component spheres and then mapped onto non-spherical particles. The linear spring-dashpot model was utilized to simulate the sphere-sphere interactions.

Ke, Renjie↗

Three-dimensional elastic stress and displacement analysis of finite circular geometry solids containing cracks

A seminumerical method is presented for solving a set of coupled partial differential equations subject to mixed and coupled boundary conditions. The use of this method is illustrated by obtaining solutions for two circular geometry and mixed boundary value problems in three-dimensional elasticity. Stress and displacement distributions are calculated in an axisymmetric, circular bar of finite dimensions containing a penny-shaped crack. Approximate results for an annular plate containing internal surface cracks are also presented.

Gyekenyesi, J. P.↗

Design of FIR digital filters for pulse shaping and channel equalization using time-domain optimization

Three algorithms are developed for designing finite impulse response digital filters to be used for pulse shaping and channel equalization. The first is the Minimax algorithm which uses linear programming to design a frequency-sampling filter with a pulse shape that approximates the specification in a minimax sense. Design examples are included which accurately approximate a specified impulse response with a maximum error of 0.03 using only six resonators. The second algorithm is an extension of the Minimax algorithm to design preset equalizers for channels with known impulse responses. Both transversal and frequency-sampling equalizer structures are designed to produce a minimax approximation of a specified channel output waveform. Examples of these designs are compared as to the accuracy of the approximation, the resultant intersymbol interference (ISI), and the required transmitted energy. While the transversal designs are slightly more accurate, the frequency-sampling designs using six resonators have smaller ISI and energy values.

Houts, R. C.↗

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗

Transonic airfoil design code

Program aids in design of shockless airfoils, assists development of fuel-conserving, supercritical wings. Algorithm calculates approximate airfoil shape given prescribed pressure distribution. This allows design of families of transonic airfoils for use in aircraft wings or turbine and compressor blades. Program is written in FORTRAN IV for batch execution on CDC-6000.

Bauer, F.↗

Complete supersonic flowfields over blunt bodies in a generalized orthogonal coordinate system

A general orthogonal coordinate system is used to describe various axisymmetric and two-dimensional shapes. Close approximations to planetary probe configurations are possible. The full Navier-Stokes equations are discretized in this coordinate system in a manner based on Allen and Cheng's numerical procedure. The blow shock is treated as a discontinuity which floats between grid points. Completely coupled flows over the forebody, base, and near wake were calculated over a cylinder, sphere, and an approximation to the Viking Aeroshell. Some problem areas in determining the base flow for increasing Reynolds number are discussed. In particular, it is found that the mean free path of the fluid near the wall immediately below the corner of the Viking Aeroshell, which experiences a severe expansion, can become greater than the local mesh size required to resolve the boundary layer in the forebody.

Gnoffo, P. A.↗

A measurement of forward-flight effects on the noise from a JT15D-1 turbofan engine in the NASA-Ames 40- by 80-Foot Wind Tunnel

A Pratt and Whitney JT15D-1 turbofan engine was tested in two facilities at Ames Research Center: the outdoor Static Test Facility and the 40- by 80-Foot Wind Tunnel. The primary purposes of the test were to determine the effects of forward velocity on the turbofan spectra in the forward quadrant for the cruise inlet and to compare these wind-tunnel spectra with outdoor spectra to determine the possibility of simulating forward-velocity effects from purely outdoor measurements. The wind-tunnel data show a reduction in the blade-passage frequency tones of the order of 10 dB with increasing forward velocity at subsonic fan-tip speeds. No forward-velocity variation was observed at supersonic tip speeds. Comparison of in-duct spectra for the cruise inlet at forward velocity, with spectra from outdoor tests with a distortion-control inlet shows excellent agreement for the in-duct data when allowance is made for different in-duct volumes. This is also reflected in good agreement for the far-field spectra at small forward angles. The comparisons of wind-tunnel and outdoor data also indicate that at least for the JT15D-1, it may be possible to approximate the shape of the far-field spectra at large directivity angles from an outdoor measurement with the cruise inlet, providing an effective inflow control device is used.

Ahtye, W. F.↗

Variation of Local Liquid-Water Concentration About and Ellipsoid of Fineness Ratio 5 Moving in a Droplet Field

Trajectories of water droplets about an ellipsoid of revolution with a fineness ratio of 5 (which often approximates the shape of an aircraft fuselage or missile) were computed with the aid of a differential analyzer. Analyses of these trajectories indicate that the local concentration of liquid water at various points about an ellipsoid in flight through a droplet field varies considerably and under some conditions may be several times the free-stream concentration. Curves of the local concentration factor as a function of spatial position were obtained and are presented in terms of dimensionless parameters Re(sub 0) (free-stream Reynolds number) and K (inertia), which contain flight and atmospheric conditions. These curves show that the local concentration factor at any point is very sensitive to change in the dimensionless parameters Re(sub 0) and K. These data indicate that the expected local concentration factors should be considered when choosing the location of, or when determining antiicing heat requirements for, water- or ice-sensitive devices that protrude into the stream from an aircraft fuselage or missile. Similarly, the concentration factor should be considered when choosing the location on an aircraft of instruments that measure liquid-water content or droplet-size distribution in the atmosphere.

Dorsch, Robert G.↗

Theory, computation, and application of exponential splines

A generalization of the semiclassical cubic spline known in the literature as the exponential spline is discussed. In actuality, the exponential spline represents a continuum of interpolants ranging from the cubic spline to the linear spline. A particular member of this family is uniquely specified by the choice of certain tension parameters. The theoretical underpinnings of the exponential spline are outlined. This development roughly parallels the existing theory for cubic splines. The primary extension lies in the ability of the exponential spline to preserve convexity and monotonicity present in the data. Next, the numerical computation of the exponential spline is discussed. A variety of numerical devices are employed to produce a stable and robust algorithm. An algorithm for the selection of tension parameters that will produce a shape preserving approximant is developed. A sequence of selected curve-fitting examples are presented which clearly demonstrate the advantages of exponential splines over cubic splines.

Mccartin, B. J.↗

A quantitative assessment of resampling errors

Applications associated with digital geographic imagery are subject to great diversity in required cell size, cartographic projection, etc. The need for resampling remote sensing scaner data is evident in all but the most undemanding cases. It is shown that proper resampling of such data is dependent in important ways on the detailed knowledge of the original scanner's effective point-spread function and to the desired point-spread function of resampled data. When both of these are known, it is relatively straightforward to compute the resampling coefficients which do the best job of approximating the shape and position of the synthesized point-spread function. The resulting synthesized psf are compared with an ideal psf located at various interpixel positions and any differences observed as errors.

Dye, R. H.↗

IUE observations of reflection nebulae

Low-resolution IUE spectra in the SWP and LWR ranges were obtained of several reflection nebulae in the vicinity of their respective illuminating stars. Data for NGC 7023 (HD200775), NGC 1435 (Merope, 23 Tau), NGC 1432 (Maia, 20 Tau) and the Electra (17 Tau) nebula were found to display significant differences in the shape of their normalized nebular spectra. If the dust in these nebulae is similar from one object to the next, these differences can be most readily explained by differences in the line-of-sight dust distribution, coupled with a wavelength-dependent phase function, changing from a strongly forward throwing form (g = 0.6-0.7) in the visible to a more nearly isotropic shape (g approximately 0.25) at 1400A.

Witt, A. N.↗