SEARCH · Engineering Papers
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.
The main goal of this is to implement a digital filtering method to generate turbulence upstream and a sub-grid model (Localized dynamic k-equation model (LDKM)) in the framework of LES which together enable reproduction of the effect of inflow turbulence on the transition seen in the VSPT experimental data at the cruise condition (incidence angle of 40o and Tu=0.5%, 5%, 10% and 15%.).
As the design complexity of modern accelerators grows, there is more interest in using advanced simulations that have fast execution time or yield additional insights like gradients. The FAST/IOTA facility has been working on implementing and experimentally validating an end-to-end digital twin that is both fast and gradient-aware, allowing for rapid prototyping of new software and experiments with minimal beam time costs. Our framework integrates physics and ML codes for linac and ring simulation through a set of generic interfaces between surrogate and physics-based sections. To reproduce device inputs and outputs, system state is exposed as a deterministic event loop in a specialized discrete event simulator architecture. Because Fermilab is undergoing control system transition, several APIs were implemented as final user interfaces - a fully asynchronous EPICS soft IOC, a gRPC-based Data Pool Manager (DPM), and legacy ACNET protocols. We discuss implementation details as well as challenges handling live data assimilation and future plans to extend modelling to main complex proton accelerators like PIPII and Booster.
To investigate the scope of ferroelectric behavior in La-substituted BiFeO 3 films, La x Bi 1−x FeO 3 epitaxial films were synthesized using off-axis co-sputtering on SrTiO 3 (001) and DyScO 3 (110) substrates with a SrRuO 3 bottom electrode layer. A digital-doping deposition method was used to enable precise control and continuous tuning of La concentration in high-quality LaxBi 1−x FeO 3 films across a wide range of x = 0.05–0.60, which was systematically investigated using piezoresponse force microscopy. Robust and reversible out-of-plane ferroelectric switching has been observed up to x = 0.35, while films with x ≥ 0.37 exhibit no measurable ferroelectric behavior, indicating a sharp ferroelectric-to-paraelectric phase transition between x = 0.35 and 0.37. This represents the highest reported La concentration in LaxBi 1−x FeO 3 films that retains ferroelectric ordering, highlighting opportunities to engineer ferroelectric and multiferroic properties in complex oxide heterostructures.
Widely recognized as one of the worst natural disaster in Puerto Rico's history, hurricane María made landfall on September 20, 2017 in southeast Puerto Rico as a high-end category 4 hurricane on the Saffir-Simpson scale causing widespread destruction, fatalities and forest disturbance. This study focused on hurricane María's effect on Puerto Rico's forests as well as the effect of landform and forest characteristics on observed disturbance patterns. In this work, we used Google Earth Engine (GEE) to assess the severity of forest disturbance using a disturbance metric based on Landsat 8 satellite data composites with pre and post-hurricane María. Forest structure, tree phenology characteristics, and landforms were obtained from satellite data products, including digital elevation model and global forest canopy height. Our analyses showed that forest structure, and characteristics such as forest age and forest type affected patterns of forest disturbance. Among forest types, highest disturbance values were found in sierra palm, transitional, and tall cloud forests; seasonal evergreen forests with coconut palm; and mangrove forests. For landforms, greatest disturbance metrics was found at high elevations, steeper slopes, and windward surfaces. As expected, high levels of disturbance were also found close to the hurricane track, with disturbance less severe as hurricane María moved inland. Results demonstrated that forest and landform characteristics accounted for 34% of the variation in spatial forest spectral disturbance patterns. This study demonstrated an informative regional approach, combining remote sensing with statistical analyses to investigate factors that result in variability in hurricane effects on forest ecosystems.
Digital mechanics explores information processing through binary, mechanical circuits. This work demonstrates a flexural, mechanical integrated circuit (m-IC) that achieves reversible, non-reciprocal signal propagation through integrated AND logic and memory. Our approach exploits sequential bistable transitions with symmetric energy wells, tunable stiffness, impedance matching, and AND gate non-linearity, to enable signal propagation, repeatability, and reversibility. We present a generalized model of logic kinematics and energetics, validated experimentally, to study energy flows, quantify energetic limits, and identify operating regimes for reversible logic. Macro-scale experiments confirm propagation dynamics, and new fabrication methods extend the architecture to micro-scale devices. By achieving controlled, reversible signal transmission across interconnected logic and memory, this work establishes a scalable platform for robust mechanical computing and adaptive sensing.
A simple mathematical jet engine description is presented where the measured transition functions of the engine thrust can be simulated in a quasi-stationary operation as well as in acceleration or deceleration schedules. Because of its simplicity and high fidelity, it is especially suited for representing jet engines in digital simulation programs.
Ability to lock onto signal increased with minimal equipment. Technique proposed to improve acquisition of spread-spectrum (pseudonoise) coded signal modulated by digital data signal. Divides correlation time into subintervals, integrates over subintervals, passes integrator outputs through square-law detectors, and adds detector outputs over all subintervals to reach decision regarding detection of code. Reduces effect of data-bit transitions at price of non-coherent-combining loss.
The PIP-II project will significantly enhance neutrino production for DUNE, Fermilab’s flagship long-baseline neutrino oscillation experiment, by doubling the beam power delivered by the accelerator complex. To achieve this, the total charge injected from the new PIP-II linac into the Booster rapid cycling synchrotron (RCS) will increase from $4.5\times 10^{12}$ to $6.5\times 10^{12}$ protons per pulse. Simultaneously, the Booster’s ramp rate will rise from 15 to 20 Hz, while the injection energy will go from 400 MeV to 800 MeV. The Booster accelerates the beam to 8 GeV and crosses transition at $\gamma_t=5.45$. In current operations, no dedicated $\gamma_t$ jump system is employed; instead, longitudinal emittance growth is controlled via an active quadrupole-mode feedback system. Downstream, the Recycler Ring uses slip-stacking to accumulate beam and increase bunch intensity. However, this process imposes a constraint on the longitudinal emittance at Booster extraction, limiting it to 0.1 eV-s (95\%) to avoid excessive particle loss. Because collective effects scale with intensity, additional measures to mitigate transition crossing may be required to stay below this limit. One potential approach is the so-called triple phase-jump technique. While the method has known limitations, it offers some advantages: it can be implemented using the existing digital low-level RF (LLRF) system, requires no additional magnets or pulsed power supplies, and remains compatible with quadrupole-mode feedback.
The PIP-II project will significantly enhance neutrino production for DUNE, Fermilab s flagship long-baseline neutrino oscillation experiment, by doubling the beam power delivered by the accelerator complex. To achieve this, the total charge injected from the new PIP-II linac into the Booster rapid cycling synchrotron (RCS) will increase from $4.5\times 10^{12}$ to $6.5\times 10^{12}$ protons per pulse. Simultaneously, the Booster s ramp rate will rise from 15 to 20 Hz, while the injection energy will go from 400 MeV to 800 MeV. The Booster accelerates the beam to 8 GeV and crosses transition at $\gamma_t=5.45$. In current operations, no dedicated $\gamma_t$ jump system is employed; instead, longitudinal emittance growth is controlled via an active quadrupole-mode feedback system. Downstream, the Recycler Ring uses slip-stacking to accumulate beam and increase bunch intensity. However, this process imposes a constraint on the longitudinal emittance at Booster extraction, limiting it to 0.1 eV-s (95\%) to avoid excessive particle loss. Because collective effects scale with intensity, additional measures to mitigate transition crossing may be required to stay below this limit. One potential approach is the so-called triple phase-jump technique. While the method has known limitations, it offers some advantages: it can be implemented using the existing digital low-level RF (LLRF) system, requires no additional magnets or pulsed power supplies, and remains compatible with quadrupole-mode feedback.
Abstract. Vegetation phenology plays a key role in controlling the seasonality of ecosystem processes that modulate carbon, water and energy fluxes between the biosphere and atmosphere. Accurate modelling of vegetation phenology in the interplay of Earth's surface and the atmosphere is thus crucial to understand how the coupled system will respond to and shape climatic changes. Phenology is controlled by meteorological conditions at different timescales: on the one hand, changes in key meteorological variables (temperature, water, radiation) can have immediate effects on the vegetation development; on the other hand, phenological changes can be driven by past environmental conditions, known as memory effects. However, the processes governing meteorological memory effects on phenology are not completely understood, resulting in their limited performance of vegetation phenology represented in land surface models. A deep learning model, specifically a long short-term memory network (LSTM), has the potential to capture and model the meteorological memory effects on vegetation phenology. Here, we apply the LSTM to model the vegetation phenology using meteorological drivers and high-temporal-resolution canopy greenness observations through digital repeat photography by the PhenoCam network. We compare a multiple linear regression model, a no-memory-effect LSTM model and a full-memory-effect LSTM model to predict the whole seasonal greenness trajectory and the corresponding phenological transition dates across 50 sites and 317 site years during 2009–2018, covering deciduous broadleaf forests, evergreen needleleaf forests and grasslands. Results show that the deep learning model outperforms the multiple linear regression model, and the full-memory-effect LSTM model performs better than the no-memory-effect model for all three plant function types (median R2 of 0.878, 0.957 and 0.955 for broadleaf forests, evergreen needleleaf forests and grasslands). We also find that the full-memory-effect LSTM model is capable of predicting the seasonal dynamic variations of canopy greenness and reproducing trends in shifting phenological transition dates. We also performed a sensitivity analysis of the full-memory-effect LSTM model to assess its plausibility, revealing its coherence with established knowledge of vegetation phenology sensitivity to meteorological conditions, particularly changes in temperature. Our study highlights that (1) multi-variate meteorological memory effects play a crucial role in vegetation phenology, and (2) deep learning opens up new avenues for improving the representation of vegetation phenological processes in land surface models via a hybrid modelling approach.
Multilevel amplitude-shift-keying (MASK) and quadrature amplitude-shift-keying (QASK) as signaling techniques for multilevel digital communications systems, and the problem of providing symbol synchronization in the receivers of such systems are discussed. A technique is presented for extracting symbol sync from an MASK or QASK signal. The scheme is a generalization of the data transition tracking loop used in PSK systems. The performance of the loop was analyzed in terms of its mean-squared jitter and its effects on the data detection process in MASK and QASK systems.
Sub-scale water tests were undertaken to assess the viability of utilizing high resolution, high frame-rate digital video recordings of a liquid rocket engine turbopump axial inducer to characterize cavitation instabilities. These high speed video (HSV) images of various cavitation phenomena, including higher order cavitation, rotating cavitation, alternating blade cavitation, and asymmetric cavitation, as well as non-cavitating flows for comparison, were recorded from various orientations through an acrylic tunnel using one and two cameras at digital recording rates ranging from 6,000 to 15,700 frames per second. The physical characteristics of these cavitation forms, including the mechanisms that define the cavitation frequency, were identified. Additionally, these images showed how the cavitation forms changed and transitioned from one type (tip vortex) to another (sheet cavitation) as the inducer boundary conditions (inlet pressures) were changed. Image processing techniques were developed which tracked the formation and collapse of cavitating fluid in a specified target area, both in the temporal and frequency domains, in order to characterize the cavitation instability frequency. The accuracy of the analysis techniques was found to be very dependent on target size for higher order cavitation, but much less so for the other phenomena. Tunnel-mounted piezoelectric, dynamic pressure transducers were present throughout these tests and were used as references in correlating the results obtained by image processing. Results showed good agreement between image processing and dynamic pressure spectral data. The test set-up, test program, and test results including H-Q and suction performance, dynamic environment and cavitation characterization, and image processing techniques and results will be discussed.
Work on the dynamics of cable-connected spinning spacecraft was completed by formulating the equations of motion by both the canonical equations and Lagrange's equations and programming them for numerical solution on a digital computer. These energy-based formulations will permit future addition of the effect of cable mass. Comparative runs indicate that the canonical formulation requires less computer time. Available literature on urban mass transportation was surveyed. Areas of the private rapid transit concept of urban transportation are also studied.
A submicrosecond time dissemination system was developed. The system is composed of an experimental transit improvement satellite (TRIAD), TRANET Satellite Tracking System, and special satellite signal receivers each tied to separate cesium clocks. The 400 MHz signal from the satellite is phase modulated plus or minus 45 deg by a pseudo-random noise (PRN) digital code with a bit rate of 5/3 MHz. Evaluation of this system has shown that satellite oscillator frequency instabilities are the major source of error. Experiments indicate that Global time transfers can be made with errors less than 100 nanoseconds.
The design of a data acquisition system often involves the integration of a Field Programmable Gate Array (FPGA) with analog front-end components to achieve precise timing and control. Reuse of these hardware systems can be difficult since they need to be tightly coupled to the communications interface and timing requirements of the specific ADC used. A hybrid design exploring the use of FPGA as a coprocessor to a traditional CPU in a dataflow architecture is presented. Reduction in the volume of data and gradual transitioning of data processing away from a hard real-time environment are both discussed. Chief design concerns, including data throughput and precise synchronization with external stimuli, are addressed. The discussion is illustrated by the implementation of a multi-channel digital integrator, a device based entirely on commercial off-the-shelf (COTS) equipment.
Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS 2 channel can be significantly altered. Further, we found that larger cations with weaker hydration energy (K + , Rb + , and Cs + ) exhibited a distinctive CN decrease under confinement, while smaller cations (Li + and Na + ) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K + , Rb + , and Cs + ) in MoS 2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS 2 channel.
The problem of the parameter identification of large scale dynamic systems involving a system matrix characterized by approximately 200 elements is addressed. By using phase variable transformations, a mathematical model of an aeroelastic airplane is described in a form that is amenable to partial or piecemeal acceptance of parameters estimated from flight data. A mathematical model of the U.S. Air Force Total In-Flight Simulator was computed using the FLEXSTAB digital computer program. As data became available during the progress of the flight test program, this data was processed and substituted in the mathematical model for parameters analytically obtained from the FLEXSTAB program. The results tend to show a progressive and orderly transition from an analytically defined mathematical model to one obtained from the flight tests of the actual aircraft.
Methods of adaptive control have been applied to suppress a potentially violent flutter condition of a half-span model of a lightweight figher aircraft. This marked the confluence of several technologies with active flutter suppression, digital control and adaptive control theory the primary contributors. The control algorithm was required to adapt both to slowly varying changes, corresponding to changes in the flight condition or fuel loading and to rapid changes, corresponding to a store release or the transition from a stable to an unstable flight condition. The development of the adaptive control methods was followed by a simulation and checkout of the complete system and a wind tunnel demonstration. As part of the test, a store was released from the model wing tip, transforming the model abruptly from a stable configuration to a violent flutter condition. The adaptive algorithm recognized the unstable nature of the resulting configuration and implemented a stabilizing control law in a fraction of a second. The algorithm was also shown to provide system stability over a range of wind tunnel Mach numbers and dynamic pressures.