Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Continuous derivative”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dynamic probabilistic risk assessment and game theory for cyber security risk analysis in nuclear power plants

Nuclear Power Plants and energy systems have become more prone to cyber-attacks with their digitalization and the increased use of smart equipment. Hence, it is important to quantify the risk associated with cyber-attacks in such systems. Dynamic Probabilistic Risk Assessment which involves studying the evolution of a system due to random events and operator and attacker actions during a cyber-attack by employing a physics-based model of the system is a suitable framework to quantify cybersecurity risk in nuclear power plants. In addition to the plant dynamics, it is also important to model the strategies of the attackers and plant operators for an effective cybersecurity risk assessment. Game theory provides a set of necessary tools to model such strategic interactions. In this research, a framework that integrates dynamic probabilistic risk assessment with game theory for cybersecurity risk analysis in nuclear power plants is presented. The mathematical formulation is derived based on the theory of continuous event trees. We propose a game theory based action model, that utilizes physics-based rewards to define the strategies of attackers and operators at every decision epoch. As a case study, the risk associated with cyber-attacks on the digital components in the secondary side of a pressurized water reactor is studied using a reduced order model. A set of attacker actions and a set of operator actions are defined for the system. The operator and attacker interactions were modelled using simultaneous game, their action policies were computed using the concept of mixed strategy Nash equilibrium and the evolution of the system was studied.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Laser Ablation-Capillary Absorption Spectroscopy: A novel approach for high throughput and increased spatial resolution measurements of δ 13 C in plant-soil systems

Spatial and temporal heterogeneity of nutrient exchange within the rhizosphere is a topic of increasing interest, although challenging to study due to limits in existing analytical capabilities. Here, we developed and demonstrated a new approach applying laser ablation sample introduction to capillary absorption spectroscopy (LA-CAS) to characterize carbon isotopic distribution within plant tissues, rhizosphere, and soil. We exposed switchgrass plants to 13 CO 2 to allow tracing of 13 C-labelled photosynthates within plant biomass and into the associated soil. The LA-CAS methods we describe leverage continuous measurements of a sample stream derived from laser ablation line scans (10-25 µm in width) over a sample surface which enables the user to produce an isotope map with a) higher data density or b) over larger spatial areas versus previously existing techniques. This versatility of LA-CAS is assessed through testing of a range of laser parameters (spot size, scan rates) on various materials (soil, plant biomass/tissues, and rhizosphere). Here, we demonstrate the ability of LA-CAS to provide near instantaneous $δ$ 13 C measurements over isotopically distinct surfaces to enable high spatially resolved mapping of 13 C-labelled material within the rhizosphere. Applying LA-CAS analysis to plant biomass, we observed higher $δ$ 13 C values concentrated within phloem structures, consistent with localized photosynthate transport. When mapping across the rhizosphere, 13 C-enriched soil was typically present within 5-10µm of root boundaries with a steep spatial increase in $δ$ 13 C when the scan approached the middle of the root. As with all LA approaches, care is required to ensure accurate results as phenomena linked to ablation, combustion, and isotopic measurement can impart artifacts if not carefully controlled. Still, taken as a whole, our demonstrations highlight the increased sample throughput, improved data density, and enhanced $δ$ 13 C capability of LA-CAS versus other LA techniques and emphasize the role this method can play in future plant and rhizosphere related studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2011

Current drive in recombining plasma

The Langevin equations describing the average collisional dynamics of suprathermal particles in nonstationary plasma remarkably admit an exact analytical solution in the case of recombining plasma. The current density produced by arbitrary particle fluxes is derived including the influence of charge recombination. Since recombination has the effect of lowering the charge density of the plasma, thus reducing the charged particle collisional frequencies, the evolution of the current density can be modified substantially compared to plasma with fixed charge density. The current drive efficiency is derived and optimized for discrete and continuous pulses of current, leading to the discovery of a nonzero ''residual'' current density that persists indefinitely under certain conditions, a feature not present in stationary plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential applications in specialized areas such as segmented inverse beta decay neutrino detectors, astronomy, machine learning, and more. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Physics↗
DOE OSTI2025

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part IV: PFNS and $\bar{ν}$ Uncertainty Propagation

An analytic benchmark with continuous-energy cross sections was previously derived to validate criticality calculations. Here, to extend the utility of the analytic benchmark to verify the implementation of $\bar{ν}$ and prompt fission neutron spectrum (PFNS) uncertainty propagation methods, new simplified forms that are dependent on the incident (fission-causing) neutron energy, as well as the outgoing neutron energy for the PFNS, are introduced in this work. The analytical forms for the flux and adjoint flux are derived for the extended benchmark and used to determine the 𝑘-eigenvalue sensitivity to $\bar{ν}$ and PFNS. The 𝑘-eigenvalue uncertainty due to $\bar{ν}$ and PFNS is calculated for the analytic benchmark using simplified$\bar{ν}$ and PFNS representations based on the ENDF-B/VIII.0 239 Pu evaluation. Because of the low sensitivity of the analytic benchmark to the physical PFNS, a nonphysical high-sensitivity PFNS is also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

STELLAR PULSATIONS AND PERIOD CHANGES IN THE SX PHOENICIS STAR XX CYGNI

Time-series photometric observations were made for the SX Phoenicis star XX Cyg between 2007 and 2011 at the Xinglong Station of National Astronomical Observatories of China. With the light curves derived from the new observations, we do not detect any secondary maximum in the descending portion of the light curves of XX Cyg, as reported in some previous work. Frequency analysis of the light curves confirms a fundamental frequency f{sub 0} = 7.4148 cycles day{sup -1} and up to 19 harmonics, 11 of which are newly detected. However, no secondary mode of pulsation is detected from the light curves. The O-C diagram, produced from 46 newly determined times of maximum light combined with those derived from the literature, reveals a continuous period increase with the rate of (1/P)(dP/dt) = 1.19(13) Multiplication-Sign 10{sup -8} yr{sup -1}. Theoretical rates of period change due to the stellar evolution were calculated with a modeling code. The result shows that the observed rate of period change is fully consistent with period change caused by evolutionary behavior predicted by standard theoretical models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Most general neutron decay correlations: Standard model recoil and radiative corrections

To continue from our previous work, we derive the full Standard Model prediction of the most general free neutron differential decay rate with all massive particles (neutron, proton, and electron) polarized, including the $\mathcal{O}$⁢(1/$\mathcal{m}$ N ) recoil corrections and $\mathcal{O}$⁢(α/π) radiative corrections. For the latter we adopt the newly developed pseudoneutrino formalism which is compatible to realistic experimental setups, in which neutrinos and photons are not detected. We also provide readily executable Mathematica notebooks to evaluate these corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Current drive in recombining plasma

The Langevin equations describing the average collisional dynamics of suprathermal particles in nonstationary plasma remarkably admit an exact analytical solution in the case of recombining plasma. The current density produced by arbitrary particle fluxes is derived including the effect of charge recombination. Since recombination has the effect of lowering the charge density of the plasma, thus reducing the charged particle collisional frequencies, the evolution of the current density can be modified substantially compared to plasma with fixed charge density. The current drive efficiency is derived and optimized for discrete and continuous pulses of current, leading to the discovery of a nonzero “residual” current density that persists indefinitely under certain conditions, a feature not present in stationary plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Data Analysis, Statistics and Probability (physics↗
DOE OSTI2026

Towards retrieving cloud top entrainment velocities from MISR cloud motion vectors

Although important, direct retrievals of entrainment rates in cloud-topped planetary boundary layer (PBL) remain elusive. Here we present a novel technique for retrieving cloud-top entrainment velocities using only Multi-angle Imaging Spectro-Radiometer (MISR) stereoscopic retrievals of cloud-motion vectors (CMVs) and cloud-top heights (CTHs). Mesoscale vertical air velocity at CTH is diagnosed from the continuity equation and then used to derive entrainment velocities from the PBL mass-budget equation. The algorithm is demonstrated through a case of marine stratocumulus deck off the California coast, with comparisons made against data from the European Centre for Medium-range Weather Forecasts (ECMWF) reanalysis (ERA5) and the data from other satellites. MISR low-cloud CTH for this case were lower than the ERA5 reported PBL depth by 189 ± 87 m. These differences in cloud top heights partly modulate the differences in the ERA5 and MISR horizontal winds, with larger differences in meridional over zonal wind components. Average difference between ERA5 and MISR derived mesoscale vertical air motion at cloud top was 0.14 ± 0.73 cm s −1 , while the same for entrainment rate was −0.09 ± 0.46 cm s −1 . The uncertainties in the utilized CTHs and CMVs are propagated to derive systematic and random retrieval uncertainties. Fractional uncertainty is lower than 25 % when the retrieved mesoscale vertical air motion is stronger than ±0.04 cm s −1 and entrainment velocities are stronger than ±0.03 cm s −1 . These results showcase the ability to derive mesoscale vertical air motion and entrainment rates from MISR observations and motivate its extension to generate a global climatology leveraging its full 23-year record (2000–2022). Nonetheless comprehensive validation of the retrievals is warranted through comparisons with estimates from an independent dataset across diverse weather conditions.

Mitra, Arka [Argonne National Laboratory (ANL), Ar↗
NASA NTRS1972

Analog simulation of a gas pressure regulator.

The behavior characteristics of direct-acting gas pressure regulators are discussed in terms of spring force, flow force, and friction force. The dynamics of the metering valve and bellows assembly are described by equating all forces on the valve to zero. Continuity equations for the regulator are derived. Simulations are carried out to determine effects of sensing orifice diameter, poppet stem friction, output volume, and supply pressure on regulator behavior. Several regulator configurations are studied.

Dustin, M. O.↗
NASA NTRS1973

Coordinated adaptive washout for motion simulators.

This paper introduces a new method of providing motion cues to a moving base six-degree-of-freedom flight simulator utilizing nonlinear filters. Coordinated adaptive filters, used to coordinate translational and rotational motion, are derived based on the method of continuous steepest descent, and the basic concept of the digital controllers used for the uncoordinated heave and yaw cues is also presented. The coordinated adaptive washout method is illustrated by an application in a six-degree-of-freedom fixed-base environment.

Parrish, R. V.↗
NASA NTRS1975

Self-regulating proportionally controlled heating apparatus and technique

A self-regulating proportionally controlled heating apparatus and technique is provided wherein a single electrical resistance heating element having a temperature coefficient of resistance serves simultaneously as a heater and temperature sensor. The heating element is current-driven and the voltage drop across the heating element is monitored and a component extracted which is attributable to a change in actual temperature of the heating element from a desired reference temperature, so as to produce a resulting error signal. The error signal is utilized to control the level of the heater drive current and the actual heater temperature in a direction to reduce the noted temperature difference. The continuous nature of the process for deriving the error signal feedback information results in true proportional control of the heating element without the necessity for current-switching which may interfere with nearby sensitive circuits, and with no cyclical variation in the controlled temperature.

Strange, M. G.↗
NASA NTRS1980

Sporadic-E layers and unstable wind shears

Electron density profiles of sporadic-E layers have been observed with good height resolution using rocket-borne probes. These generally show a simple shape consistent with the effect of a linear wind shear acting on metallic ions. Occasionally more complex shapes have been recorded, including double peaks and, on one occasion, a nearly rectangular profile. A direct method of obtaining the wind profile from the concentration profile of metallic ions has been developed. The metallic ion concentration profile itself is obtained from the electron density profile. Both procedures derive from the steady-state continuity equation. For linear wind shears it is found that the maximum value of the shear is about 50 m/s/km which corresponds to a Richardson number of 1/4. Layers of complex shape are associated with non-linear wind shears in which the maximum shear considerably exceeds this value. It is concluded that the complex profiles of sporadic-E layers can be interpreted as an effect of unstable wind shears.

Smith, L. G.↗
NASA NTRS1982

Implementation of structural response sensitivity calculations in a large-scale finite-element analysis system

The implementation includes a generalized method for specifying element cross-sectional dimensions as design variables that can be used in analytically calculating derivatives of output quantities from static stress, vibration, and buckling analyses for both membrane and bending elements. Limited sample results for static displacements and stresses are presented to indicate the advantages of analytically calclating response derivatives compared to finite difference methods. Continuing developments to implement these procedures into an enhanced version of the system are also discussed.

Giles, G. L.↗