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 existence of time-periodic stochastic motions of an incompressible fluid is obtained. Here the fluid is subject to a time-periodic body force and an additional time-periodic stochastic force that is produced by a rigid body moves periodically stochastically with the same period in the fluid.
We consider a quantum mechanical system represented in phase space (referred to hereafter as “Wigner space”), coupled to a harmonic oscillator bath. We derive quantum hierarchal Fokker-Planck (QHFP) equations not only in real time but also in imaginary time, which represents an inverse temperature. This is an extension of a previous work, in which we studied a spin-boson system, to a Brownian system. It is shown that the QHFP in real time obtained from a correlated thermal equilibrium state of the total system possesses the same form as those obtained from a factorized initial state. A modified terminator for the hierarchal equations of motion is introduced to treat the non-Markovian case more efficiently. Using the imaginary-time QHFP, numerous thermodynamic quantities, including the free energy, entropy, internal energy, heat capacity, and susceptibility, can be evaluated for any potential. These equations allow us to treat non-Markovian, non-perturbative system-bath interactions at finite temperature. Through numerical integration of the real-time QHFP for a harmonic system, we obtain the equilibrium distributions, the auto-correlation function, and the first- and second-order response functions. These results are compared with analytically exact results for the same quantities. This provides a critical test of the formalism for a non-factorized thermal state and elucidates the roles of fluctuation, dissipation, non-Markovian effects, and system-bath coherence. Employing numerical solutions of the imaginary-time QHFP, we demonstrate the capability of this method to obtain thermodynamic quantities for any potential surface. It is shown that both types of QHFP equations can produce numerical results of any desired accuracy. The FORTRAN source codes that we developed, which allow for the treatment of Wigner space dynamics with any potential form (TanimuranFP15 and ImTanimuranFP15), are provided as the supplementary material.
The correct folding of proteins is of paramount importance for their function, and protein misfolding is believed to be the primary cause of a wide range of diseases. Protein folding has been investigated with time-averaged methods and time-resolved spectroscopy, but observing the structural dynamics of the unfolding process in real-time is challenging. Here, we demonstrate an approach to directly reveal the structural changes in the unfolding reaction. We use nano- to millisecond time-resolved x-ray solution scattering to probe the unfolding of apomyoglobin. The unfolding reaction was triggered using a temperature jump, which was induced by a nanosecond laser pulse. We demonstrate a new strategy to interpret time-resolved x-ray solution scattering data, which evaluates ensembles of structures obtained from molecular dynamics simulations. We find that apomyoglobin passes three states when unfolding, which we characterize as native, molten globule, and unfolded. The molten globule dominates the population under the conditions investigated herein, whereas native and unfolded structures primarily contribute before the laser jump and 30 μs after it, respectively. The molten globule retains much of the native structure but shows a dynamic pattern of inter-residue contacts. Our study demonstrates a new strategy to directly observe structural changes over the cause of the unfolding reaction, providing time- and spatially resolved atomic details of the folding mechanism of globular proteins.
A new time-resolved opacity spectrometer (OpSpecTR) is currently under development for the National Ignition Facility (NIF) opacity campaign. Here, the spectrometer utilizes Icarus version 2 (IV2) hybridized complementary metal–oxide–semiconductor sensors to collect gated data at the time of the opacity transmission signal, unlocking the ability to collect higher-temperature measurements on NIF. Experimental conditions to achieve higher temperatures are feasible; however, backgrounds will dominate the data collected by the current time-integrating opacity spectrometer. The shortest available OpSpecTR integration time of ~2 ns is predicted to reduce self-emission and other late-time backgrounds by up to 80%. Initially, three Icarus sensors will be used to collect data in the self-emission, backlighter, and absorption regions of the transmission spectrum, with plans to upgrade to five Daedalus sensors in future implementations with integration times of ~1.3 ns. We present the details of the diagnostic design along with recent characterization results of the IV2 sensors.
Time-resolved measurements of materials provide a wealth of information on quasiparticle dynamics, and have been the focus of optical studies for decades. In this paper, we develop a theory for explicitly evaluating time-resolved resonant inelastic x-ray scattering (tr-RIXS). We apply the theory to a noninteracting electronic system and reveal the particle-hole spectrum and its evolution during the pump pulse. With a high-frequency pump, the frequency and amplitude dependence analysis of the spectra agrees well with the steady-state assumptions and Floquet excitations. When the pump frequency is low, the spectrum extracts real-time dynamics of the particle-hole continuum in momentum space. These results verify the correctness of our theory and demonstrate the breadth of physical problems that tr-RIXS could shed light on.
Data driven science is becoming increasingly more common, complex, and is placing tremendous stresses on visualization and analysis frameworks. Data sources producing 10GB per second (and more) are becoming increasingly commonplace in both simulation, sensor and experimental sciences. These data sources, which are often distributed around the world, must be analyzed by teams of scientists that are also distributed. Enabling scientists to view, query and interact with such large volumes of data in near-real-time requires a rich fusion of visualization and analysis techniques, middleware and workflow systems. Here, this paper discusses initial research into visualization and analysis of distributed data workflows that enables scientists to make near-real-time decisions of large volumes of time varying data.
The modern power grid is a cyber-physical system. While the grid is becoming more intelligent with emerging sensing and communication techniques, new vulnerabilities are introduced and cyber security becomes a major concern. One type of cyber-attacks - False Data Injection Attacks (FDIAs) - exploits the limitations in traditional power system state estimation, and modifies system states without being detected. In this paper, we propose a physics-guided deep learning (PGDL) approach to defend against FDIAs. The PGDL takes real-time measurements as inputs to neural networks, outputs the estimated states, and reconstructs measurements considering power system physics. A deep recurrent neural network - Long Short-Term Memory (LSTM) - is employed to learn the temporal correlations among states. This hybrid learning model leads to a time-series state estimation method to defend against FDIAs. The simulation results using IEEE 14-bus test system demonstrate the accuracy and robustness of the proposed time-series state estimation under FDIAs.
The conditions and dynamics of neon gas puff Z-pinch plasmas at pinch time are studied on a pulsed power generator with a current rise time of approximately 200-ns and 0.9-mA peak current. Radial tailoring of the gas puff massdensity profile using a triple-nozzle coaxial valve (two annular gas liners and a central jet) allows production of both more stable and less stable (with regard to the magneto-Rayleigh– Taylor instability) Z-pinch implosions. A 526.5-nm, 10-J Thomson scattering diagnostic laser enables probing of the flow dynamics and plasma conditions of these implosions with both spatial and temporal resolutions. The 2.2-ns laser pulse scatters from the plasma electrons and is carried by one optical fiber to a visible light streak camera, and by a bundle of optical fibers to a time-gated camera, both after spectral dispersion by 750-mm spectrometers. The streak camera, with a 10-ns full streak time, provides subnanosecond resolution of the evolution of the pinch plasma parameters through stagnation. The timegated camera provides spatially resolved spectra (across a field of view of 6.1 mm) at the same time as the streak. Scattering spectra suggest that temperatures are high at stagnation, with the ion temperature as much as three times higher than the electron temperature. Furthermore, we consider the possibility of nonthermal explanations for the broad scattering spectra and high effective ion temperature, including collisionality, implosion velocity distributions (velocity gradients), and small-scale hydrodynamic motion.
Purpose: To develop a real-time applicator position monitoring system (RAPS) for intracavitary brachytherapy using an infrared camera and reflective markers. Methods: 3D image-guided brachytherapy requires high accuracy of applicator localization; however, applicator displacement can happen during patient transfer for imaging and treatment delivery. No continuous applicator position monitoring system is currently available. The RAPS system was developed for real-time applicator position monitoring without additional radiation dose to patients. It includes an infrared camera, reflective markers, an infrared illuminator, and image processing software. After reflective markers are firmly attached to the applicator and the patient body, applicator displacement can be measured by computing the relative change in distance between the markers. The reflective markers are magnetic resonance imaging (MRI) compatible, which is suitable for MRI-guided HDR brachytherapy paradigm. In our prototype, a Microsoft Kinect sensor with a resolution of 640 by 480 pixels is used as an infrared camera. A phantom study was carried out to compare RAPS' measurements with known displacements ranging from −15 to +15 mm. A reproducibility test was also conducted. Results: The RAPS can achieve 4 frames/s using a laptop with Intel{sup ®} Core™2 Duo processor. When the pixel size is 0.95 mm, the difference between RAPS' measurements and known shift values varied from 0 to 0.8 mm with the mean value of 0.1 mm and a standard deviation of 0.44 mm. The system reproducibility was within 0.6 mm after ten reposition trials. Conclusions: This work demonstrates the feasibility of a real-time infrared camera based gynecologic intracavitary brachytherapy applicator monitoring system. Less than 1 mm accuracy is achieved when using an off-the-shelf infrared camera.
In this work, we study rotating stratified turbulence (RST) making use of numerical data stemming from a large parametric study varying the Reynolds, Froude and Rossby numbers, Re, Fr and Ro in a broad range of values. The computations are performed using periodic boundary conditions on grids of 10243 points, with no modeling of the small scales, no forcing and with large-scale random initial conditions for the velocity field only, and there are altogether 65 runs analyzed in this paper. The buoyancy Reynolds number defined as R B = ReFr 2 varies from negligible values to ≈ 105, approaching atmospheric or oceanic regimes. This preliminary analysis deals with the variation of characteristic time scales of RST with dimensionless parameters, focusing on the role played by the partition of energy between the kinetic and potential modes, as a key ingredient for modeling the dynamics of such flows. We find that neither rotation nor the ratio of the Brunt-Väisälä frequency to the inertial frequency seem to play a major role in the absence of forcing in the global dynamics of the small-scale kinetic and potential modes. Specifically, in these computations, mostly in regimes of wave turbulence, characteristic times based on the ratio of energy to dissipation of the velocity and temperature fluctuations, T V and T P , vary substantially with parameters. Their ratio γ=T V /T P follows roughly a bell-shaped curve in terms of Richardson number Ri. It reaches a plateau --on which time scales become comparable, γ≈0.6 -- when the turbulence has significantly strengthened, leading to numerous destabilization events together with a tendency towards an isotropization of the flow.
In recent years, streaming-based data processing has been gaining substantial traction for dealing with overwhelming data generated by real-time applications, from both enterprise sources and scientific computing. In this work, however, we look at an emerging class of scientific data with Near Real-Time (NRT) requirement, in which data is typically generated in a bursty fashion with the near real-time constraints being applied primarily between bursts, rather than within a stream. A key challenge for this types of data sources is that the processing time per data element is not uniform, and not always feasible to predict. Given the observations on the increasing unpredictability of compute load and system dynamics, this work looks to adapt streaming-based approach to the context of this new class of large experiments and simulations that have complex run-time control and analysis issues. In particular, we deploy a novel two-tier scheme for handling the increasing unpredictability of runtime behaviors: Instead of relying on determining what and where to run the scientific workflows beforehand or partial dynamically, the decision will also be adaptively enhanced online according to system runtime status. This is enabled by embedding workflow along with data streams. Specifically, we break data outputs generated from experiments or simulations into multiple self-describing "chunks", which we call active data objects. As such, if there is a transient hotspot observed, a data object with unfinished workflow pipeline can break its previous schedule and search for a least loaded location to continue the execution. Our preliminary experiment results based on synthetic workloads demonstrate the proposed active workflow system as a very promising solution by outperforming the state-of-the-art semi-dynamic workflow schedulers with an improved workflow completion time, as well as a good scalability.
The Fermilab PIP-II (proton improvement plan - II) project is being constructed at Fermilab to deliver $800\,MeV$ protons of $>1\,MW$ beam power to replace the present LINAC and provide protons to the remainder of the existing accelerator complex. The new LINAC consists of a warm front end, 23 superconducting RF cryomodules, and a beam transfer line to the existing complex. The cryomodules (CMs) are to be tested at Fermilab's CryoModule Test Facility (CMTF).An important measurement in cryogenic testing is the heat load of each CM. Traditionally, at Fermilab, these measurements were made collecting archived data offline and analyzing it. The new control system for PIP-II is being developed with the EPICS (Experimental Physics and Industrial Control System) framework, which allows us to compute the heat load in real time using the HePak library.We are exploring other $He$ properties, such as flow, where flow meters are not available, which can also be calculated in real time and fed back to the cryogenics engineers.This paper details the real time heat load calculation and $He$ flow software developed for CM testing at CMTF, as well as the first results from the prototype HB650 CM. Future plans for 2-phase $LHe$ flow will also be outlined.
This study examines the dynamics and control of an engine operated with late intake valve closure (LIVC) timings in a dual-fuel combustion mode. The engine features a fuel delivery system in which diesel is direct-injected and natural gas is port-injected. Despite the benefits of LIVC and dual-fuel strategy, combining these two techniques resulted in efficiency losses due to the variability of the combustion process across cylinders. The difference in power production across cylinders ranges from 9% at an IVC of 570°ATDC* to 38% at an IVC of 620 °ATDC and indicates an increasingly uneven fuel distribution as the intake valve remains open longer in the compression stroke. This paper describes an approach for controlling the amount of fuel injected into each cylinders’ port of an inline six- cylinder heavy-duty dual-fuel engine to minimize the variations in fuel distribution across cylinder. Since measuring the actual fuel reaching each cylinder is not a practical option on stock engines, this study aimed to produce a feedforward control technique. Using a model developed on an engine simulation software, the optimal fuel injection amounts at each cylinder port across 125 operating conditions (consisting of sweeps of five engine speeds, five engine loads, and five IVC timings) were determined by using PI controllers that adjust the amount of fuel injected at each port based on the amount of fuel reaching the cylinders. The results from the PI controller were expressed as a function of the operating condition (engine speed, engine load, and IVC timing); and consequently, could be reproduced without relying on a measurement of the amount of fuel reaching each cylinder. The feedforward control strategy, unlike the PI controller, is easily implementable on a stock engine and significantly reduces the variations in fuel distribution across cylinders.
A magnetic particle time-of-flight (MagPTOF) diagnostic has been designed to measure shock- and compression-bang time using DHe-fusion protons and DD-fusion neutrons, respectively, at the National Ignition Facility (NIF). This capability, in combination with shock-burn weighted areal density measurements, will significantly constrain the modeling of the implosion dynamics. This design is an upgrade to the existing particle time-of-flight (pTOF) diagnostic, which records bang times using DD or DT neutrons with an accuracy better than ±70 ps [H. G. Rinderknecht et al., Rev. Sci. Instrum. 83, 10D902 (2012)]. The inclusion of a deflecting magnet will increase DHe-proton signal-to-background by a factor of 1000, allowing for the first time simultaneous measurements of shock and compression-bang times in DHe-filled surrogate implosions at the NIF.
A method is presented for factoring a two-impulse orbital transfer into a three- or four-impulse transfer which solves the rendezvous problem and satisfies an intermediate timing constraint. Both the time of rendezvous and the intermediate time of a alinement are formulated as any element of a finite sequence of times. These times are integer multiples of a constant plus an additive constant. The rendezvous condition is an equality constraint, whereas the intermediate alinement is an inequality constraint. The two timing constraints are satisfied by factoring the impulses into collinear parts that vectorially sum to the original impulse and by varying the resultant period differences and the number of revolutions in each orbit. Five different types of solutions arise by considering factoring either or both of the two impulses into two or three parts with a limit for four total impulses. The impulse-factoring technique may be applied to any two-impulse transfer which has distinct orbital periods.
The proceedings of a conference on Precise Time and Time Interval Planning are presented. The subjects discussed include the following: (1) satellite timing techniques, precision frequency sources, and very long baseline interferometry, (2) frequency stabilities and communications, and (3) very low frequency and ultrahigh frequency propagation and use. Emphasis is placed on the accuracy of time discrimination obtained with time measuring equipment and specific applications of time measurement to military operations and civilian research projects.
Subjects in a time-free environment for 14 days estimated the hour and day several times a day. Half of the subjects were under a heavy exercise regime. During the waking hours, the no-exercise group showed no difference between estimated and real time, whereas the exercise group showed significantly shorter estimated than real time. Neither group showed a difference after the sleeping periods. However, the mean accumulated error for the two groups was 48.73 hours and was strongly related to the displacements of sleep/waking behavior. It is concluded that behavioral cues are the primary determinants of time estimates in time-free environments.
The development of a high precision time interval digitizer is described. The time digitizer is a 10 psec resolution stop watch covering a range of up to 340 msec. The measured time interval is determined as a separation between leading edges of a pair of pulses applied externally to the start input and the stop input of the digitizer. Employing an interpolation techniques and a 50 MHz high precision master oscillator, the equivalent of a 100 GHz clock frequency standard is achieved. Absolute accuracy and stability of the digitizer are determined by the external 50 MHz master oscillator, which serves as a standard time marker. The start and stop pulses are fast 1 nsec rise time signals, according to the Nuclear Instrument means of tunnel diode discriminators. Firing level of the discriminator define start and stop points between which the time interval is digitized.