Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential mode”

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 217 records · Page 12

Flutter of a uniform wing with an arbitrarily placed mass according to a differential-equation analysis and a comparison with experiment

A method is presented for the calculation of the flutter speed of a uniform wing carrying an arbitrarily placed concentrated mass. The method, an extension of recently published work by Goland and Luke, involves the solution of the differential equations of motion of the wing at flutter speed and therefore does not require the assumption of specific normal modes of vibration. The order of the flutter determinant to be solved by this method depends upon the order of the system of differential equations and not upon the number of modes of vibration involved. The differential equations are solved by operational methods, and a brief discussion of operational methods as applied to boundary-value problems is included in one of two appendixes. A comparison is made with experiment for a wing with a large eccentrically mounted weight and good agreement is obtained. Sample calculations are presented to illustrate the method; and curves of amplitudes of displacement, torque, and shear for a particular case are compared with corresponding curves computed from the first uncoupled normal modes.

Runyan, Harry L↗

Global and Local Helioseismic Studies of Solar Convection Zone Dynamics Using SOI-MDI on SOHO

Our joint collaborative analyses of global mode data to characterize the solar differential rotation (e.g. Thompson et al. 1996, Schou et al. 1998), and most recently to detect and analyze temporal variations in angular velocity Omega profiles both within the convection zone and in the deeper radiative interior (e.g. Howe et al 2000a,b; Toomre et al. 2000), have led to a series of fascinating discoveries. These should be pursued further as the solar cycle continues. The physical deductions being made from these studies have been greatly strengthened by utilizing both SOI-MDI and GONG data in order to have two independent observational realizations of Doppler images spanning a five-year interval, using two separate procedures to determine global mode splittings, and then analyzing those splitting data sets using both RLS and SOLA inversion procedures. There are considerable subtleties in the effects of instrumental response functions and calibrations, sensitivity of peak finding algorithms and their mode leakage estimates, and stochastic variations in mode amplitudes that can all contribute to apparent changes in the Omega profiles being inferred from sequences of helioseismic data. We have come to understand the implications of many of these calibration and analysis steps, greatly aided by frequent multi-week collaborative working sessions in our Helioseismic Analysis Facility (HAF) at JILA involving many members of the SO1 dynamics and inversion team, including most of our Co-Is during the summer months when we hold intensive working sessions. Considerable further focused attention is required in a collaborative setting on such global mode issues as we continue studying the changing sun.

Toomre, Juri↗

Nonlinear interaction between a pair of oblique modes in a supersonic mixing layer: Long-wave limit

The nonlinear interaction between a pair of symmetric, oblique, and spatial instability modes is studied in the long-wave limit using asymptotic methods. The base flow is taken to be a supersonic mixing layer whose Mach number is such that the corresponding vortex sheet is marginally stable according to Miles' criterion. It is shown that the amplitude of the mode obeys a nonlinear integro-differential equation. Numerical solutions of this equation show that, when the obliqueness angle is less than pi/4, the effect of the nonlinearity is to enhance the growth rate of the instability. The solution terminates in a singularity at a finite streamwise location. This result is reminiscent of that obtained in the vicinity of the neutral point by other authors in several different types of flows. On the other hand, when the obliqueness angle is more than pi/4, the streamwise development of the amplitude is characterized by a series of modulations. This arises from the fact that the nonlinear term in the amplitude equation may be either stabilizing or destabilizing, depending on the value of the streamwise coordinate. However, even in this case the amplitude of the disturbance increases, though not as rapidly as in the case for which the angle is less than pi/4. Quite generally then, the nonlinear interaction between two oblique modes in a supersonic mixing layer enhances the growth of the disturbance.

Balsa, Thomas F.↗

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗

Snekmer: a scalable pipeline for protein sequence fingerprinting based on amino acid recoding

Abstract Motivation The vast expansion of sequence data generated from single organisms and microbiomes has precipitated the need for faster and more sensitive methods to assess evolutionary and functional relationships between proteins. Representing proteins as sets of short peptide sequences (kmers) has been used for rapid, accurate classification of proteins into functional categories; however, this approach employs an exact-match methodology and thus may be limited in terms of sensitivity and coverage. We have previously used similarity groupings, based on the chemical properties of amino acids, to form reduced character sets and recode proteins. This amino acid recoding (AAR) approach simplifies the construction of protein representations in the form of kmer vectors, which can link sequences with distant sequence similarity and provide accurate classification of problematic protein families. Results Here, we describe Snekmer, a software tool for recoding proteins into AAR kmer vectors and performing either (i) construction of supervised classification models trained on input protein families or (ii) clustering for de novo determination of protein families. We provide examples of the operation of the tool against a set of nitrogen cycling families originally collected using both standard hidden Markov models and a larger set of proteins from Uniprot and demonstrate that our method accurately differentiates these sequences in both operation modes. Availability and implementation Snekmer is written in Python using Snakemake. Code and data used in this article, along with tutorial notebooks, are available at http://github.com/PNNL-CompBio/Snekmer under an open-source BSD-3 license. Supplementary information Supplementary data are available at Bioinformatics Advances online.

59 BASIC BIOLOGICAL SCIENCES↗

PLZT matrix-type block data composers

Description of the characteristics of PLZT-matrix type block data composers operated in the strain-biased, scattering, differential phase, and edge-effect modes. The operation peculiarities of each of these four modes are discussed, and their respective merits are compared.

Drake, M. D.↗

A factored implicit scheme for numerical weather prediction

An implicit method is proposed to factor the nonlinear partial differential equations governing fast and slow modes of dynamic motion in numerical weather prediction schemes. The method permits separate factorization of the slow and fast modes of the implicit operator. A simple two-dimensional version of the system of three-dimensional equations governing atmospheric dynamics over shallow water was analyzed to assess the accuracy of the proposed method. It is shown that the method has a small error which is comparable to other discretization errors in the overall scheme.

Augenbaum, J. M.↗

Progress of 2-micron Detectors for Application to Lidar Remote Sensing

AlGaAsSb/InGaAsSb heterojunction phototransistors were developed at Astropower, Inc under Laser Risk Reduction Program (LRRP) for operation in the 2-micron region. These phototransistors were optimized for 2-micron detection and have high quantum efficiency (>60%), high gain (>10(exp 3)) and low noise-equivalent- power (<5x10(exp -14) W/Hz), while operating at low bias voltage. One of these phototransistors was tested in lidar mode using the 2-micron CO2 Differential Absorption Lidar (DIAL) system currently under development under the Instrument Incubator Program (IIP) at NASA Langley. Lidar measurements included detecting atmospheric structures consisting of thin clouds in the mid-altitude and near-field boundary layer. These test results are very promising for the application of phototransistors for the two-micron lidar remote sensing. In addition, HgCdTe avalanche photodiodes (APD) acquired from Raytheon were used in atmospheric testing at 2-microns. A discussion of these measurements is also presented in this paper.

Abedin, M. N.↗

Analytic quasi-steady evolution of marginally unstable kinetic plasma instabilities in the presence of drag and scattering

The 1D bump-on-tail problem is studied in order to determine the influence of drag on quasi-steady solutions near marginal stability (1-γ d /γ L << 1) when effective collisions are much larger than the instability growth rate (ν>>γ). In this common tokamak regime, it is rigorously shown that the paradigmatic Berk–Breizman cubic equation for the nonlinear mode evolution reduces to a much simpler differential equation, dubbed the time-local cubic equation, which can be solved directly. It is found that in addition to increasing the saturation amplitude, drag introduces a shift in the apparent oscillation frequency by modulating the saturated wave envelope. Furthermore, excellent agreement is found between the analytic solution for the mode evolution and both the numerically integrated Berk–Breizman cubic equation and fully nonlinear 1D Vlasov simulations. Experimentally isolating the contribution of drag to the saturated mode amplitude for verification purposes is explored but complicated by the reality that the amount of drag cannot be varied independently of other key parameters in realistic scenarios. While the influence of drag is modest when the ratio of effective drag to effective scattering α/ν is very small, it can become substantial when α/ν ≳ 0.5, suggesting that drag should be accounted for in quantitative models of fast-ion-driven instabilities in fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Test of lepton flavor universality and search for lepton flavor violation in B → Kℓℓ decays

We present measurements of the branching fractions for the decays B → Kμ + μ - and B → Ke + e - , and their ratio (R K ), using a data sample of 711 fb -1 that contains 772 × 10 6 B$\bar{B}$ events. The data were collected at the Υ(4S) resonance with the Belle detector at the KEKB asymmetric-energy e + e - collider. The ratio R K is measured in five bins of dilepton invariant-mass-squared (q 2 ): q 2 $\epsilon$ (0.1, 4.0), (4.00, 8.12), (1.0, 6.0), (10.2, 12.8) and (> 14.18) GeV 2 /c 4 , along with the whole q 2 region. The R K value for q 2 $\epsilon$ (1.0, 6.0) GeV 2 /c 4 is 1.03$^{+0.28}_{-0.24}$ ± 0.01. The first and second uncertainties listed are statistical and systematic, respectively. All results for RK are consistent with Standard Model predictions. We also measure CP-averaged isospin asymmetries in the same q 2 bins. The results are consistent with a null asymmetry, with the largest difference of 2.6 standard deviations occurring for the q 2 $\epsilon$ (1.0, 6.0) GeV 2 /c 4 bin in the mode with muon final states. The measured differential branching fractions, dB/dq 2 , are consistent with theoretical predictions for charged B decays, while the corresponding values are below the expectations for neutral B decays. We have also searched for lepton-flavor-violating B → Kμ ± e ∓ decays and set 90% confidence-level upper limits on the branching fraction in the range of 10 -8 for B + → K + μ ± e ∓ , and B 0 → K 0 μ ± e ∓ modes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Kinetic theory for electrostatic waves due to transverse velocity shears

A kinetic theory in the form of an integral equation is provided to study the electrostatic oscillations in a collisionless plasma immersed in a uniform magnetic field and a nonuniform transverse electric field. In the low temperature limit the dispersion differential equation is recovered for the transverse Kelvin-Helmholtz modes for arbitrary values of K parallel, where K parallel is the component of the wave vector in the direction of the external magnetic field assumed in the z direction. For higher temperatures the ion-cyclotron-like modes described earlier in the literature by Ganguli, Lee and Plamadesso are recovered. In this article, the integral equation is reduced to a second-order differential equation and a study is made of the kinetic Kelvin-Helmholtz and ion-cyclotron-like modes that constitute the two branches of oscillation in a magnetized plasma including a transverse inhomogeneous dc electric field.

Ganguli, G.↗

Ultrafast method for scanning tunneling spectroscopy

A scanning tunneling microscope (STM) combines unique capabilities in imaging and spectroscopy with atomic precision, and it can obtain energy-resolved spectroscopic data with atomic resolution. In this paper, we utilize a recently proposed modification to the STM feedback control loop to acquire high quality d 2 I/dV 2 images. We have developed a constant differential conductance imaging method by closing the STM feedback loop with a high precision dI/dV measurement. In this mode, the tip’s vertical position is adjusted so as to keep the differential conductance constant during raster scanning of the surface. Based on this imaging mode, we propose a new technique to acquire fast and reliable scanning tunneling spectroscopy (STS) data simultaneously with the imaging.

47 OTHER INSTRUMENTATION↗

Static Feed Water Electrolysis Subsystem Testing and Component Development

A program was carried out to develop and test advanced electrochemical cells/modules and critical electromechanical components for a static feed (alkaline electrolyte) water electrolysis oxygen generation subsystem. The accomplishments were refurbishment of a previously developed subsystem and successful demonstration for a total of 2980 hours of normal operation; achievement of sustained one-person level oxygen generation performance with state-of-the-art cell voltages averaging 1.61 V at 191 ASF for an operating temperature of 128F (equivalent to 1.51V when normalized to 180F); endurance testing and demonstration of reliable performance of the three-fluid pressure controller for 8650 hours; design and development of a fluid control assembly for this subsystem and demonstration of its performance; development and demonstration at the single cell and module levels of a unitized core composite cell that provides expanded differential pressure tolerance capability; fabrication and evaluation of a feed water electrolyte elimination five-cell module; and successful demonstration of an electrolysis module pressurization technique that can be used in place of nitrogen gas during the standby mode of operation to maintain system pressure and differential pressures.

Koszenski, E. P.↗

Regularizing the r-mode Problem for Nonbarotropic Relativistic Stars

We present results for r-modes of relativistic nonbarotropic stars. We show that the main differential equation, which is formally singular at lowest order in the slow-rotation expansion, can be regularized if one considers the initial value problem rather than the normal mode problem. However, a more physically motivated way to regularize the problem is to include higher order terms. This allows us to develop a practical approach for solving the problem and we provide results that support earlier conclusions obtained for uniform density stars. In particular, we show that there will exist a single r-mode for each permissible combination of 1 and m. We discuss these results and provide some caveats regarding their usefulness for estimates of gravitational-radiation reaction timescales. The close connection between the seemingly singular relativistic r-mode problem and issues arising because of the presence of co-rotation points in differentially rotating stars is also clarified.

Lockitch, Keith H.↗

Axisymmetric perturbations of thin gaseous disks. I - Unstable convective modes and their consequences for the solar nebula

The axisymmetric perturbations of a thin, differentially rotating gas disk in which the vertical temperature stratification is superadiabatic are analyzed. The growth rate of adiabatic convective normal modes is calculated, departures from simple polytropic disk models are briefly discussed, and the detailed vertical structure of the eigenfunctions is analyzed. The stabilizing effect of radiative diffusion on convective modes is considered. It is found that rotation and compressibility tend to reduce the rate of growth of the disturbances, and that the growth rate increases without limit and is proportional to the square root of the radial wavenumber in polytropic equilibria. The maximum radial size of convective eddies scales like the square root of the degree of superadiabaticity times the size of the convective zone. There is significant convective penetration into radiatively stable layers only for the fundamental and low-order harmonic modes, whose vertical wavelength is comparable to the size of the convective layer.

Ruden, Steven P.↗

A spherical harmonic decomposition technique for analysing steady photospheric flows

Steady flows in the photosphere, including differential rotation, meridional circulation and convection, can be represented by a spectrum of modes. A technique is described in which the spherical harmonic functions are used to determine the characteristics of this spectrum and the nature of these flows. Some information about the spectrum is lost because only one hemisphere is seen, only the line of sight velocity is measured and the measurements contain noise. This produces an apparent mixing between spectral modes. By analyzing synthetic data it is found that, in spite of this mixing, differential rotation can be accurately measured, meridional circulations with small amplitudes can be extracted from the data and giant cells might be separated from supergranules.

Hathaway, David H.↗

A self-consistent linear-mode model of stellar convection

A normal-mode expansion of the linearized fluid equations in terms of small subset of spherical harmonics can provide a foundation for a physically motivated, self-consistent description of a solar-type convection zone. In the absence of dissipation, a second-order differential equation governs the radial dependence of the modes, so that interpretation of the effects on convection quantities of the normal-form 'potential well' is straightforward. The philosophy is quite different from the more recent work of Narasimha and Antia (1982): all envelopes presented here differ substantially from MLT envelopes, and therefore, from theirs, which are constructed to be consistent with MLT. The amplitude of all modes is set by a Kelvin-Helmholtz-('shear'-) instability argument unrelated to solar observations, with the result that the convection description may be considered to arise from 'first-hueristic-principles'. The thermodynamics modelled vaguely resemble the sun's, and more vigorously convective envelopes show some phenomena qualitatively like solar observations (e.g., atmospheric velocity spectra).

Macauslan, J.↗

Distinguishing fissile uranium isotopes using an active well neutron coincidence counter

Proposed thorium-based nuclear fuel cycles are likely to require quantification and verification of 233 U within nuclear material. Because of their similar fission cross sections, active neutron nondestructive assay (NDA) systems may respond similarly to 233 U and 235 U. Traditional safeguards equipment has been optimized for 235 U and 238 U quantification associated with conventional uranium/plutonium fuel cycles and may not be directly applicable to 233 U quantification when mixed with other actinides. This work used models of the large volume active well coincidence counter (LV-AWCC) at Oak Ridge National Laboratory to evaluate the performance of this neutron NDA system to differentiate fissile uranium isotopes. The models were developed to simulate NDA system performance in response to a number of triangular radiation signature training device sources within the central cavity or well. This work predicted that the LV-AWCC can effectively differentiate 233 U from 235 U in certain modes of operation. In active mode, the LV-AWCC with the cadmium liner results in different doubles count rates between the fissile isotopes for a given fissile uranium mass. Without the cadmium liner, the uranium isotopes provide a statistically indistinguishable doubles count rate response for the fissile masses considered in this work (up to approximately 150 g). The cadmium liner serves to harden the neutron interrogation spectrum, which better exploits the notable difference in the 233 U and 235 U fission cross sections at approximately 1 eV. In passive mode, the two fissile isotopes exhibit different doubles and singles count rates regardless of liner presence because the passive source strength of 233 U is approximately 2 orders of magnitude stronger than that of 235 U due to the shorter half-life and correspondingly higher (α, n) yield. We conclude that using neutron interrogation in the LV-AWCC, two measurements are needed to quantify 233 U content in mixed uranium items. The first measurement is used to determine the total fissile uranium mass using a mode that cannot distinguish fissile isotopes (i.e., where a similar response is observed for both fissile uranium isotopes such as active doubles without cadmium or using a thermal neutron interrogation source). In conclusion, the second measurement is used to determine the 233 U content by using a differentiating technique (e.g., passive doubles, passive doubles to singles ratio, active doubles with cadmium).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗