Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear response”

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 325 records · Page 18

Electron ptychography achieves atomic-resolution limits set by lattice vibrations

Transmission electron microscopes use electrons with wavelengths of a few picometers, potentially capable of imaging individual atoms in solids at a resolution ultimately set by the intrinsic size of an atom. However, owing to lens aberrations and multiple scattering of electrons in the sample, the image resolution is reduced by a factor of 3 to 10. By inversely solving the multiple scattering problem and overcoming the electron-probe aberrations using electron ptychography, we demonstrate an instrumental blurring of less than 20 picometers and a linear phase response in thick samples. The measured widths of atomic columns are limited by thermal fluctuations of the atoms. As a result, our method is also capable of locating embedded atomic dopant atoms in all three dimensions with subnanometer precision from only a single projection measurement.

42 ENGINEERING↗

Simple strategy for the simulation of axially symmetric large-area metasurfaces

Metalenses are composed of nanostructures for focusing light and have been widely explored in many exciting applications. However, their expanding dimensions pose simulation challenges. We propose a method to simulate metalenses in a timely manner using vectorial wave and ray tracing models. We sample the metalens’s radial phase gradient and locally approximate the phase profile by a linear phase response. Each sampling point is modeled as a binary blazed grating, employing the chosen nanostructure, to build a transfer function set. The metalens transmission or reflection is then obtained by applying the corresponding transfer function to the incoming field on the regions surrounding each sampling point. Fourier optics is used to calculate the scattered fields under arbitrary illumination for the vectorial wave method, and a Monte Carlo algorithm is used in the ray tracing formalism. We validated our method against finite-difference time domain simulations at 632 nm, and we were able to simulate metalenses larger than 3000 wavelengths in diameter on a personal computer.

Martins, Augusto (ORCID:0000000295546481)↗

Evaluation of Low Cost Infrasound Sensor Packages

Sandia National Laboratories has tested and evaluated the performance of the following five models of low-cost infrasound sensors and sensor packages: Camas microphone, Gem Infrasound Logger, InfraBSU sensor, Raspberry Boom, and the Samsung S10 smartphone utilizing the Redvox app. The purpose of this infrasound sensor evaluation is to measure the performance characteristics in such areas as power consumption, sensitivity, self-noise, dynamic range, response, passband, linearity, sensitivity variation due to changes in static pressure and temperature, and sensitivity to vertical acceleration. The infrasound monitoring community has leveraged such sensors and integrated packages in novel ways; better understanding the performance of these units serves the geophysical monitoring community.

42 ENGINEERING↗

The Stochastic Gause Predator-Prey model: Noise-induced extinctions and invariance

We consider the Gause predator-prey with general bounded or sub‑linear functional responses, – which includes those of Holling types Ⅰ–Ⅳ. – and multiplicative Gaussian noise. In contrast to previous studies, the prey in our model follows logistic dynamics while the predator's population is solely regulated by consumption of the prey. To ensure well-posedeness, we derive explicit Lyapunov‐type criteria ensuring global positivity and moment boundedness of solutions. We find conditions for noise‑induced extinctions, proving that stochasticity can drive either population to collapse even when the deterministic analogue predicts stable coexistence. In the case when the predator becomes extinct, we establish a limiting distribution for the predator's population. Last, for functional responses of Holling type Ⅰ, we provide sufficient conditions on the intensity of the noise for the existence and uniqueness of a stationary distribution.

Gause model↗

Trace Chemical Detection Using Intercalated MXenes as a Signal Enhancing Substrate in Optical Probes

MXenes are 2D materials composed of layered transition metal nitrides or carbides. These materials are synthesized by HF exfoliation from MAX phases (Ti{sub 3}AlC{sub 2}). The 2D nanomaterial was synthesized by the removal of the 'A' element, resulting in a Mxene product (Ti{sub 3}C{sub 2}). MXenes have the general formula M{sub n+1}X{sub n}T{sub x}, where M is an early transition metal, X is Carbon and/or Nitrogen, and T accounts for surface terminated functional groups such as Fluoride, hydroxyl, and oxygen. These materials have very unique properties, similar to graphene, that allows them to be applied in a variety of trace detection techniques including surface-enhanced Raman spectroscopy (SERS). MXenes have also been demonstrated to selectively uptake uranyl ion, UO{sub 2}{sup 2+}. If this property can be combined with SERS or fluorescence detection, it may be possible to use MXenes as the basis for an alternative method to kinetic phosphorescence analysis (KPA) for trace uranyl measurements. Objectives: To confirm that MXene Nano materials are suitable substrates for SERS and sensor development by enhancing Raman signaling. To determine if certain MXene preparation methods yield materials that are more suitable for trace sensing methods. To determine uranyl uptake properties of these MXene materials and test them for analytical signals. Sample Preparation: Preparation of Ti{sub 3}C{sub 2}MXene (at FSU). MXenes were prepared by etching Al from Ti{sub 3}AlC{sub 2} (MAX phase)material. Two etching techniques yield different MXene products: LiF/HCl: Milder reaction, larger MXene flakes. HF: Harsher reaction, smaller flakes, larger layer separation. Products washed to remove etchant, vacuum filtered, and dried. Dried MXene flakes are air-stable. Film preparation for sensor testing (at SRNL): Suspend powder in diH{sub 2}O, purge with Ar, sonicate for 30 min. Centrifuge and collect supernate with suspended particles. Observed LiF-etched Mxene yielded a higher density of particles and darker collected solution. Drop-cast (4 ml) supernate onto slides and dried with Ar. For Rhodamine B (RhB) testing, drop-cast 4 ml drops onto Mxene spots and dried with Ar. Scanning Electron Microscopy conditions: 10 kV Beam energy, high vacuum; Working distance of 8 mm; beam penetration depth appx. 4 microns, beam spot size appx. 2 nanometers. Results: Detection of aluminum correlates with bright spots on image. Presence of aluminum shows that LiF/HCl etching was less thorough than HF etching. Trace Cl detection in LiF images suggests incomplete rinsing. HF has smaller feature size, more layer structure, and increased homogeneity, consistent with expectations. Macroscopic Raman spectroscopy measurements: 532 nm excitation, ∼50 mW with a ∼100 micron spot size (InPhotonics RPB probe). Kaiser Optical Holospec f/1.8 spectrometer with cooled (-60 deg.C) Andor iDus OE420 CCD. LiF 1x supernate showed good signal for trace measurements of Rhodamine B. HF and 1/4x LiF supernates showed little Mxene or Rhodamine B signal. Low deposition densities led to excess background signal from glass slides. For LiF film, response is linear with Rhodamine B concentration over range tested. Will retest with Raman microscope (∼1 micron spot size) to characterize SERS of more dilute LiF and HF etched Mxenes. Conclusions: The LiF etched material was more suitable for macroscopic SERS measurements because it was more concentrated, resulting in a thicker film than the HF etched Mxene and diluted LiF sample. However, the other materials may give greater SERS enhancements, which we hope to determine from measurements with the Raman microscope. From characterization with SEM we concluded that the HF etched Mxene is more uniform/homogenous and has smaller particle size than the LiF etched Mxene. There is still aluminum present in both samples indicating that etching wasn't complete, but the removal of the aluminum was more efficient in the HF method. Path Forward: Observe SERS with Raman microscopy, to obtain better signals for the more diluted samples and be able to compare enhancement effects for the different MXenes. Characterize uranyl sorption into MXene films and test Raman and fluorescence signals. Revisit the etching conditions to improve removal of aluminum. FSU and SRNL will continue to collaborate to create and characterize different Mxene materials and test their usefulness for sensor applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Amplitude- and frequency-dependent activation of layer II/III neurons by intracortical microstimulation

Intracortical microstimulation (ICMS) has been used for the development of brain machine interfaces. However, further understanding about the spatiotemporal responses of neurons to different electrical stimulation parameters is necessary to inform the design of optimal therapies. In this study, we employed in vivo electrophysiological recording, two-photon calcium imaging, and electric field simulation to evaluate the acute effect of ICMS on layer II/III neurons. Our results show that stimulation frequency non-linearly modulates neuronal responses, whereas the magnitude of responses is linearly correlated to the electric field strength and stimulation amplitude before reaching a steady state. Temporal dynamics of neurons’ responses depends more on stimulation frequency and their distance to the stimulation electrode. In addition, amplitude-dependent post-stimulation suppression was observed within ~500 μm of the stimulation electrode, as evidenced by both calcium imaging and local field potentials. These findings provide insights for selecting stimulation parameters to achieve desirable spatiotemporal specificity of ICMS.

59 BASIC BIOLOGICAL SCIENCES↗

Quasi-linear toroidal simulations of resonant magnetic perturbations in eight ITER H-mode scenarios

Abstract Both linear and quasi-linear aspects of the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for various H-mode scenarios in ITER, covering the pre-fusion power operation and the fusion power operation phases. Linear response computations for eight ITER scenarios, with varying plasma current and toroidal magnetic field, reveal that the best coil current phasing for controlling the type-I edge localized modes (ELMs) scales roughly linearly with the edge safety factor. The coil phasing is defined as the relative toroidal phase of the coil currents between different rows, for a given toroidal harmonic. Quasi-linear initial value simulation, which is the focus of the present study, shows that application of the n = 3 ( n is the toroidal mode number) RMP field has a minimum side effect on the plasma core momentum confinement but potentially a large effect on the global particle transport. Generally, the RMP field with the best (worst) coil phasing for ELM control produces the strongest (weakest) effect on the plasma edge flow and the overall density. This robustly holds for all eight ITER scenarios. Consequently, in order to minimize the RMP induced side effects while achieving ELM control (suppression) in ITER, a compromise is necessary in choosing the coil current configuration.

Physics↗

Coastal Wetland Carbon Sequestration in a Warmer Climate (Final Report)

Coastal wetlands are global hotspots of carbon storage and locations where carbon and nitrogen cycles have a disproportionately large impact on land, water, and air in comparison to the area they occupy. The extremely high rates of carbon sequestration in these systems are the result of complex feedbacks between vegetation and the physical environment. The plant and microbial interactions that drive these feedbacks are absent in ecosystem- and global-scale Earth System models. The goal of this proposal was to address the uncertain response of the coastal carbon sink to climate change, a critical knowledge gap in coastal carbon research and a significant barrier to incorporating coastal wetlands into Earth System Models. The Salt Marsh Accretion Response to Temperature eXperiment (SMARTX) was designed to understand the ecosystem-scale consequences of warming and elevated CO₂ (eCO₂) in tidal wetlands at the coastal terrestrial-aquatic interface. We successfully designed and built a novel whole-ecosystem warming and eCO₂ experiment in a coastal wetland. The gradient design of the warming treatments (+0, +1.7, +3.4, +5.1 ⁰C) allowed us to discover unexpected non-linear and non-additive responses arising from plant-microbe interactions. We concluded that these non-linear responses to warming were the result of plants and microbes responding to temperature at different thresholds, with important consequences for forecasting terrestrial ecosystem responses to climate change and understanding global trends.

54 ENVIRONMENTAL SCIENCES↗

Nonlinear photoconductivities and quantum geometry of chiral multifold fermions

Here, chiral multifold fermions are quasiparticles that appear only in chiral crystals such as transition metal silicides in the cubic B20 structure (i.e., the CoSi family), and they may show exotic physical properties. Here we study the injection and shift photoconductivities and also the related geometrical quantities for several types of chiral multifold fermions, including spin-1/2 as well as pseudospin-1 and -3/2 fermions, dubbed as Kramers Weyl, triple-point, and Rarita-Schwinger-Weyl (RSW) fermions, respectively. We utilize the minimal symmorphic model to describe the triple-point fermions (TPF). We also consider the more realistic model Hamiltonian for the CoSi family including both linear and quadratic terms. We find that injection currents due to circularly polarized light are quantized as a result of the Chern numbers carried by the multifold fermions within the linear models. Surprisingly, we discover that in the TPF model, the linear shift conductivities, responsible for the shift current generation by linearly polarized light, are proportional to the pseudo spin-orbit coupling and independent of photon frequency. In contrast, for the RSW and Kramer Weyl fermions, the linear shift conductivity is linearly proportional to photon frequency. The numerical results agree with the power-counting analysis for quadratic Hamiltonians. The frequency independence of the linear shift conductivity could be attributed to the strong resonant symplectic Christoffel symbols of the flat bands. Moreover, the calculated symplectic Christoffel symbols show significant peaks at the nodes, suggesting that the shift currents are due to the strong geometrical response near the topological nodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Gauge invariance of light-matter interactions in first-principle tight-binding models

In this work, we study the different ways of introducing light-matter interaction in first-principle tight-binding (TB) models. The standard way of describing optical properties is the velocity gauge, defined by linear coupling to the vector potential. In finite systems a transformation to represent the electromagnetic radiation by the electric field instead is possible, albeit subtleties arise in periodic systems. The resulting dipole gauge is a multi-orbital generalization of the Peierls substitution. In this work we investigate the accuracy for both pathways, with particular emphasis on gauge invariance, for TB models constructed from maximally localized Wannier functions. This approach accurately captures the light-matter interaction close to the Fermi level. Focusing on paradigmatic two-dimensional materials, we construct first-principle models and calculate the response to electromagnetic fields in linear response and for strong excitations. Benchmarks against fully converged first-principle calculations allow for ascertaining the accuracy of the TB models. We find that the dipole gauge provides a more accurate description than the velocity gauge in all cases. The main deficiency of the velocity gauge is an imperfect cancellation of paramagnetic and diamagnetic current. Formulating a corresponding sum rule however provides a way to explicitly enforce this cancellation. This procedure corrects the TB models in the velocity gauge, yielding excellent agreement with dipole gauge and thus restoring gauge invariance.

2-dimensional systems↗

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the Resistive DCON code [A.H. Glasser, Z.R. Wang and J.-K. Park, Physics of Plasmas, \textbf{23}, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Since the solution of ideal outer region is critical to the asymptotic matching and is challenging to be solved in the toroidal geometry due to the singular power series solution at the resonant surfaces, systematic verification of the outer region $\Delta^\prime$ matrix is made by reproducing the well known analytical $\Delta^{\prime}$ result in [H.P. Furth, P.H. Rutherford and H. Selberg, The Physics of Fluids, \textbf{16}, 1054-1063 (1073)] as well as by making a quantitative benchmark with the PEST3 code [A. Pletzer and R.L. Dewar, J. Plasma Physics, \textbf{45}, 427-451 (1991)]. Finally, the reconstructed numerical solution of resistive plasma response from the toroidal matching matrix is presented. Comparing with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

MHD simulation↗

Ignition Delay Measurements of Four Component Model Gasolines Exploring the Impacts of Biofuels and Aromatics

This study explores the impacts of combinations of biofuel (ethanol, isobutanol and 2-methyl furan) and aromatic (toluene) compounds in a four component fuel blend, at fixed research octane number (RON) on ignition delay measured in an advanced fuel ignition delay analyzer (AFIDA 2805). Ignition delay measurements were performed over a range of temperatures from 400 to 725 °C (673 to 998 K) and two chamber pressures of 10 and 20 bar. The four component mixtures are compared to primary reference fuels at RON values of 90 and 100. The ignition delay measurements show that as the aromatic and biofuel concentrations increased, two stage ignition behavior was suppressed, at both initial chamber pressures. But both RON 100 (isooctane) and RON 90 reference fuels showed two stage ignition behavior, as did fuel mixtures with low biofuel and aromatic content. RON 90 fuels showed stronger two stage ignition behavior than RON 100 fuels, as expected. Depending on the type of biofuel in the mixture, the ignition delay at low chamber temperatures could be far greater than for the reference fuels. In particular, for the RON 100 mixtures at either 10 or 20 bar initial chamber pressure, the ignition delay at 400 °C (673 K) for the high level blend of 2-methyl furan and toluene (30 vol% of each) exhibited an ignition delay that was 10 times longer than for neat isooctane. The results show the strong non-linear octane blending response of these three biofuel compounds, especially in concert with the kinetic antagonism that toluene is known to display in mixtures with isooctane. These results have implications for the formulation of biofuel mixtures for spark ignition and advanced compression ignition engines, where this non-linear octane blending response could be exploited to improve knock resistance, or modulate the autoignition process.

2-methyl furan↗

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving the linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the resistive DCON code [A. H. Glasser, Z. R. Wang, and J.-K. Park, Phys. Plasmas 23, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Additionally, the solution of the ideal outer region is critical to the asymptotic matching and is challenging in toroidal geometry due to the singular power series solution at the resonant surfaces. Thus, a systematic verification of the outer region Δ' matrix is made by reproducing the well-known analytical Δ' result in Furth et al. [Phys. Fluids 16, 1054–1063 (1073)] and by making a quantitative benchmark with the PEST3 code [A. Pletzer and R. L. Dewar, J. Plasma Phys. 45, 427–451 (1991)]. Finally, the reconstructed numerical solution of the resistive plasma response from the toroidal matching matrix is presented. Compared with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cyber physical attack detection

A cyber-security threat detection system and method stores physical data measurements from a cyber-physical system and extracts synchronized measurement vectors synchronized to one or more timing pulses. The system and method synthesizes data integrity attacks in response to the physical data measurements and applies alternating parametrized linear and non-linear operations in response to the synthesized data integrity attacks. The synthesis renders optimized model parameters used to detect multiple cyber-attacks.

97 MATHEMATICS AND COMPUTING↗

Cyber physical attack detection

A cyber-security threat detection system and method stores physical data measurements from a cyber-physical system and extracts synchronized measurement vectors synchronized to one or more timing pulses. The system and method synthesize data integrity attacks in response to the physical data measurements and applies alternating parameterized linear and non-linear operations in response to the synthesized data integrity attacks. The synthesis renders optimized model parameters used to detect multiple cyber-attacks.

Ferragut, Erik M.↗

A data-driven linear formulation of the optimal demand response scheduling problem for an industrial air separation unit

Demand response (DR) has become a key element in balancing the power grid as the contribution of time-varying renewable power generation increases. Chemical plants are appealing candidates for DR programs as they offer large, concentrated and flexible loads. DR participation calls for frequent production rate changes over time scales that overlap with the dominant dynamics of the plant. Production scheduling should therefore consider the process dynamics explicitly. Here we present a data-driven approach for modelling the scheduling-relevant dynamics based on historical closed-loop operating data using autoregressive with extra inputs (ARX) models. We introduce a new, linear scheduling problem formulation based on the ARX representation, and demonstrate its implementation on an industrial air separation unit.

42 ENGINEERING↗

Performance Demonstration of Self-Powered Neutron Detectors for Steady-State Reactor Operations

The irradiation testing of sensors in reactors is a crucial step towards calibrating and qualifying sensors prior to their deployment in experiments. This report details the process toward qualifying and calibrating custom-designed rhodium-based self-power neutron detectors (Rh-SPNDs) for steady-state reactor irradiations. This process serves to both demonstrate the performance capabilities of Rh-SPNDs as well as to provide experimental data for development of a sensor sensitivity model. Two designs of Rh-SPNDs were tested in various reactors to demonstrate: detection resolution in a low neutron flux environment, a delayed-response compensation technique, output linearity in a large range of neutron flux, and measurement accuracy verified with dosimetry. The detection resolution and compensation technique was demonstrated in the AGN-201m reactor at Idaho State University. The irradiation confirmed the sensors’ capability to perform steady-state operations in a low neutron flux of ~2E8 n/cm 2 -sec. Sensor output linearity coupled with the delayed-response compensation was investigated at the neutron radiography reactor at Idaho National Laboratory. A Rh-SPND was irradiated to neutron fluxes ranging from 2E8 to 2E13 n/cm2-sec range. The measured data demonstrated a wide and linear range of operation with a measured linear sensitivity of 1.0129 ×10 -13 A/W with a correlation-squared value of r 2 =0.9927. The measurement accuracy was investigated at the Advanced Test Reactor Critical reactor. The SPNDs were inserted into a test vehicle with collocated flux wires. Two irradiations with different flux levels were performed, and the SPNDs relative measurement between the two irradiations was calculated to be 1.2613 ± 0.0153 for the small SPND design and 1.1809 ± 0.0108 for the large SPND design. Both SPND measurements fell between the co-axial dosimetry result, which reported 1.218 ± 0.047. Additionally, the preliminary MCNP model for calculating SPND sensitivity was developed in parallel to this work. Modeled neutron spectrum with measured magnitude was used for inputs to determine the simulated SPND output. The results showed an overestimation of signal strength by a factor of 5, which was expected because of model simplification. This leads to future modeling work to account for signal losses from additional physical properties, including high temperature environments for FY-21.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗