Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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 1,099 records · Page 61

Estimating and Forecasting Time-Varying Groundwater Recharge in Fractured Rock: A State-Space Formulation with Preferential and Diffuse Flow to the Water Table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

groundwater, Recharge, infiltration, recursive est↗

Multi-qubit nanoscale sensing with entanglement as a resource

Nitrogen vacancy (NV) centres in diamond are widely deployed as local magnetic sensors, using single-qubit control to measure both time-averaged fields and noise with nanoscale spatial resolution. Moving beyond single qubits to multi-qubit control enables new sensing modalities such as measuring nonlocal spatiotemporal correlators or using entangled states to enhance measurement sensitivity. Here, in this work, we describe protocols for using optically unresolved NV centre pairs and nuclear spins as multi-qubit sensors for measuring correlated noise at nanometre length scales. For noninteracting NV centres, we implement a phase-cycling protocol that disambiguates magnetic correlations from variance fluctuations, leveraging the presence of a third qubit, a 13 C nucleus, to effect coherent single-NV spin flips and enable phase cycling even for co-aligned NV centres that are spectrally unresolved. For length scales around 10 nm, we create maximally entangled Bell states through dipole–dipole coupling between two NV centres and use these entangled states to directly read out the magnetic field correlation, rather than reconstructing it from independent measurements of unentangled NV centres. Importantly, this changes the scaling of sensitivity with readout noise from quadratic to linear. For conventional off-resonant readout of the NV centre spin state (for which the readout noise is roughly 30 times the quantum projection limit), this results in more than an order of magnitude improvement in sensitivity. Finally, we demonstrate methods for detecting high spatial- and temporal-resolution correlators with pairs of strongly interacting NV centres.

condensed-matter physics↗

The growth of density perturbations in the last ~10 billion years from tomographic large-scale structure data

In order to investigate the origin of the ongoing tension between the amplitude of matter fluctuations measured by weak lensing experiments at low redshifts and the value inferred from the cosmic microwave background anisotropies, we reconstruct the evolution of this amplitude from z ~ 2 using existing large-scale structure data. To do so, we decouple the linear growth of density inhomogeneities from the background expansion, and constrain its redshift dependence making use of a combination of 6 different data sets, including cosmic shear, galaxy clustering and CMB lensing. We analyze these data under a consistent harmonic-space angular power spectrum-based pipeline. We show that current data constrain the amplitude of fluctuations mostly in the range 0.2 < z < 0.7, where it is lower than predicted by . This difference is mostly driven by current cosmic shear data, although the growth histories reconstructed from different data combinations are consistent with each other, and we find no evidence of systematic deviations in any particular experiment. In spite of the tension with , the data are well-described by the model, albeit with a lower value of S 8 ≡ σ 8 (Ω m /0.3) 0.5 . As part of our analysis, we find constraints on this parameter of S 8 = 0.7781 ± 0.0094 (68% confidence level), reaching almost percent-level errors comparable with CMB measurements, and 3.4σ away from the value found by .

79 ASTRONOMY AND ASTROPHYSICS↗

Galactic echoes

Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.

Galaxy: evolution↗

Microstructural pattern formation during liquid metal dealloying: Phase-field simulations and theoretical analyses

In recent years, liquid metal dealloying has emerged as a promising material processing method to generate micro- and nanoscale bicontinuous or porous structures. Most previous studies focused on the experimental characterization of the dealloying process and on the properties of the dealloyed materials, leaving the theoretical study incomplete to fully understand the fundamental mechanisms of the liquid metal dealloying process. In this paper, we use theoretical models and phase-field simulations to clarify the kinetics and pattern formation during liquid metal dealloying. Our investigation starts from a theoretical analysis of the 1D dissolution of a binary precursor alloy, which reveals that the 1D dissolution process involves two regimes. In the first regime, due to the low solubility of one of the elements in the melt, it accumulates at the solid-liquid interface, which reduces the dissolution kinetics. In the second regime, the interface kinetics reaches a stationary regime where both elements of the precursor alloy dissolve into the melt. Previous works revealed that in the early dealloying stage, the dealloying front is destabilized by an interfacial spinodal decomposition, which triggers the formation of interconnected ligaments. Here we extend this line of work by proposing a linear stability analysis able to predict the initial length scale of the ligaments formed in the initial stage of the dealloying. Combining this analysis with the 1D dissolution model proposed here enables us to better understand the initial conditions (composition of the precursor alloy and the melt) leading to a planar dissolution without interface destabilization. Finally, we report a strong influence of solid-state diffusion on dealloying. Although the solid-state diffusivity is four to five orders of magnitude smaller than in the liquid phase, it is found to affect both dissolution kinetics and ligament morphologies.

36 MATERIALS SCIENCE↗

QuantumNAT: quantum noise-aware training with noise injection, quantization and normalization

Parameterized Quantum Circuits (PQC) are promising towards quantum advantage on near-term quantum hardware. However, due to the large quantum noises (errors), the performance of PQC models has a severe degradation on real quantum devices. Take Quantum Neural Network (QNN) as an example, the accuracy gap between noise-free simulation and noisy results on IBMQ-Yorktown for MNIST-4 classification is over 60%. Existing noise mitigation methods are general ones without leveraging unique characteristics of PQC; on the other hand, existing PQC work does not consider noise effect. To this end, we present QuantumNAT, a PQC-specific framework to perform noise-aware optimizations in both training and inference stages to improve robustness. We experimentally observe that the effect of quantum noise to PQC measurement outcome is a linear map from noise-free outcome with a scaling and a shift factor. Motivated by that, we propose post-measurement normalization to mitigate the feature distribution differences between noise-free and noisy scenarios. Furthermore, to improve the robustness against noise, we propose noise injection to the training process by inserting quantum error gates to PQC according to realistic noise models of quantum hardware. Lastly, post-measurement quantization is introduced to quantize the measurement outcomes to discrete values, achieving the denoising effect. Extensive experiments on 8 classification tasks using 6 quantum devices demonstrate that QuantumNAT improves accuracy by up to 43%, and achieves over 94% 2-class, 80% 4-class, and 34% 10-class classification accuracy measured on real quantum computers. The code for construction and noise-aware training of PQC is available in the TorchQuantum library.

Wang, Hanrui↗

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗

Cosmic background radiation anisotropies in universes dominated by nonbaryonic dark matter

Detailed calculations of the temperature fluctuations in the cosmic background radiation for universes dominated by massive collisionless relics of the big bang are presented. An initially adiabatic constant curvature perturbation spectrum is assumed. In models with cold dark matter, the simplest hypothesis - that galaxies follow the mass distribution leads to small-scale anisotropies which exceed current observational limits if omega is less than 0.2 h to the -4/3. Since low values of omega are indicated by dynamical studies of galaxy clustering, cold particle models in which light traces mass are probably incorrect. Reheating of the pregalactic medium is unlikely to modify this conclusion. In cold particle or neutrino-dominated universes with omega = 1, presented predictions for small-scale and quadrupole anisotropies are below current limits. In all cases, the small-scale fluctuations are predicted to be about 10 percent linearly polarized.

Bond, J. R.↗

Simultaneous rocket and MST radar observation of an internal gravity wave breaking in the mesosphere

In June, 1983, the Structure and Atmospheric Turbulence Environment (STATE) rocket and Poker Flat Mesophere-Stratosphere-Troposphere radar campaign was conducted to measure the interaction between turbulence, electron density and electron density gradient that has produced unusually strong MST radar echoes from the summer mesosphere over Poker Flat, Alaska. Analysis or radar wind measurements and a concurrent wind and temperature profile obtained from a rocket probe carrying a three-axis accelerometer are given. The two data sets provide a fairly complete (and in some cases, redundant) picture of the breaking (or more correctly, the saturation) of a large-amplitude, low-frequency, long-wavelength internal gravity wave. The data show that small-scale turbulence and small-scale wave intensity is greatest at those altitudes where the large-scale wave-induced temperature lapse rate is most negative or most nearly unstable, but the wind shear due to the large-scale wave is a minimum. A brief review of linear gravity-wave theory is presented as an aid to the identification of the gravity-wave signature in the radar and rocket data. Analysis of the time and height cross sections of wind speed and turbulence intensity observed by the Poker Flat MST radar follows. Then, the vertical profile of temperature and winds measured by a rocket probe examined. Finally, the use of the independent data sets provided by the rocket and the radar are discussed and implications for theories of wave saturation are presented.

Smith, S. A.↗

Steady surficial core motions - An alternate method

A new method for deriving steady surficial core motions from geomagnetic field models is presented and applied. The original method determines the steady velocity field at the top of a frozen-flux core which best fits, in both the spatial and temporal linear least squares sense, a model of the large scale geomagnetic secular variation (SV) at the base of a source-free mantle (Voorhies, 1986). New method is based upon the same physical assumptions, but fits SV at earth's surface instead of the core-mantle boundary (CMB). The newly derived flow differs somewhat from prior solutions, and shows but slight westward drift; however, other key global properties are virtually unchanged. Over a 15 year interval the new method provides a better fit to the time-varying geomagnetic field at earth's surface than does the old.

Voorhies, Coerte V.↗

The current-voltage relationship in auroral current sheets

The current-voltage relation within narrow auroral current sheets is examined through the use of high-resolution data from the high altitude Dynamics Explorer 1 satellite. The north-south perpendicular electric field and the east-west magnetic field are shown for three cases in which there are large amplitude, oppositely directed paired electric fields and narrow current sheets. These data are shown to indicate that there is a linear Ohm's law relationship between the current density and the parallel potential drop within the narrow current sheets. This linear relationship had previously been verified for large-scale auroral formations greater than 20 km wide at the ionosphere. The evidence shown here extends our knowledge down to the scale size of discrete auroral arcs.

Weimer, D. R.↗

An analytically soluble problem in fully nonlinear statistical gravitational lensing

The amplification probability distribution p(I)dI for a point source behind a random star field which acts as the deflector exhibits a I exp-3 behavior for large amplification, as can be shown from the universality of the lens equation near critical lines. In this paper it is shown that the amplitude of the I exp-3 tail can be derived exactly for arbitrary mass distribution of the stars, surface mass density of stars and smoothly distributed matter, and large-scale shear. This is then compared with the corresponding linear result.

Schneider, P.↗

How many upper Eocene microspherule layers: More than we thought

The scientific controversy over the origin of upper Eocene tektites, microtektites and other microspherules cannot be logically resolved until it is determined just how many events are involved. The microspherule-bearing beds in marine sediments have been dated using standard biozonal techniques. Although a powerful stratigraphic tool, zonal biostratigraph has its limitations. One is that if an event, such as a microspherule occurrence, is observed to occur in a zone at one locality and then a similar event observed in the same zone at another locality, it still may be unwarranted to conclude that these events exactly correlate. To be in a zone a sample only need be between the fossil events that define the zone boundaries. It is often very difficult to accurately determine where within a zone one might be. Further, the zone defining events do not everywhere occur at the same points in time. That is, the ranges of the defining taxa are not always filled. Thus, the length of time represented by a zone (but not, of course, its chronozone) can vary from place to place. These problems can be offset by use of chronostratigraphic modelling techniques such as Graphic Correlation. This technique was used to build a Cretaceous and Cenozoic model containing fossil, magnetopolarity, and other events. The scale of the model can be demonstrated to be linear with time. This model was used to determine the chronostratigraphic position of upper Eocene microspherule layers.

Hazel, Joseph E.↗

The resonance of twin supersonic jets

This paper presents an analytical study of the resonant interaction between twin supersonic jets. An instability wave model is used to describe the large scale coherent structures in the jet mixing layers. A linearized shock cell model is also given for the jets when operating off design. The problem's geometry admits four types of normal modes associated with each azimuthal mode number in the single jet. The stability of these modes is examined for both a vortex sheet model of the jet and a jet flow represented by realistic profiles. The growth rates of each mode number and type are found to vary with jet separation and mixing layer thickness and Strouhal number. Contours of equal pressure level are obtained for each mode. The region close to the symmetry axis is found to have the greatest pressure fluctuation amplitude.

Morris, Philip J.↗

Modal techniques for analyzing airplane dynamics

A series of techniques are presented for analyzing airplane control and output characteristics. State-space matrix equations describing the linear perturbation dynamics are transformed from physical coordinates into scaled, modal coordinates. Techniques for analyzing the impacts of system inputs on the fundamental modes of motion and the appearance of these modes in the system outputs are explained. Scaled modal matrices are used to compute steady-state control inputs which optimize the steady-state response of selected system outputs. Graphics which promote quick understanding of the analysis are presented to display the resulting vectors and matrices. The defined analysis techniques are applied to an example airplane model, illustrating the insight which can be acquired using the described modal techniques.

Arbuckle, P. Douglas↗

Linear tearing modes of a forced current-sheet equilibrium

Linear tearing modes are studied in a nearly singular forced current-sheet equilibrium, such as could result from global magnetic forces in the solar corona. Growth rates for the tearing modes, determined by solving the linearized reduced MHD (Strauss) equations numerically, were found to scale as (gamma)tau(d) = about S exp 4/5 k(y) exp 4/5, where S is the Lundquist number, k(y) is the wavenumber, and tau(d) is the classical resistive diffusion time. This scaling is in agreement with predictions from analytical theory. Because of the faster S scaling of these modes compared to the tearing modes of a diffuse current-sheet equilibrium, the modes have much higher growth rates (by a factor of about 10,000) for coronal values of S (about 10 to the 12th). For coronal parameters, the growth times of these new modes are estimated to be on the order of several hours to days, as compared to growth times of months to years for the tearing modes in a diffuse current-shear equilibrium. The growth times are comparable to reconnection times scales required in models of coronal heating by magnetic field dissipation.

Liewer, Paulett C.↗

A controls engineering approach for analyzing airplane input-output characteristics

An engineering approach for analyzing airplane control and output characteristics is presented. State-space matrix equations describing the linear perturbation dynamics are transformed from physical coordinates into scaled coordinates. The scaling is accomplished by applying various transformations to the system to employ prior engineering knowledge of the airplane physics. Two different analysis techniques are then explained. Modal analysis techniques calculate the influence of each system input on each fundamental mode of motion and the distribution of each mode among the system outputs. The optimal steady state response technique computes the blending of steady state control inputs that optimize the steady state response of selected system outputs. Analysis of an example airplane model is presented to demonstrate the described engineering approach.

Arbuckle, P. Douglas↗

Computation of wind tunnel wall effects for complex models using a low-order panel method

A technique for determining wind tunnel wall effects for complex models using the low-order, three dimensional panel method PMARC (Panel Method Ames Research Center) has been developed. Initial validation of the technique was performed using lift-coefficient data in the linear lift range from tests of a large-scale STOVL fighter model in the National Full-Scale Aerodynamics Complex (NFAC) facility. The data from these tests served as an ideal database for validating the technique because the same model was tested in two wind tunnel test sections with widely different dimensions. The lift-coefficient data obtained for the same model configuration in the two test sections were different, indicating a significant influence of the presence of the tunnel walls and mounting hardware on the lift coefficient in at least one of the two test sections. The wind tunnel wall effects were computed using PMARC and then subtracted from the measured data to yield corrected lift-coefficient versus angle-of-attack curves. The corrected lift-coefficient curves from the two wind tunnel test sections matched very well. Detailed pressure distributions computed by PMARC on the wing lower surface helped identify the source of large strut interference effects in one of the wind tunnel test sections. Extension of the technique to analysis of wind tunnel wall effects on the lift coefficient in the nonlinear lift range and on drag coefficient will require the addition of boundary-layer and separated-flow models to PMARC.

Ashby, Dale L.↗