Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation analysis”

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,477 records · Page 82

Efficient algorithms for use in probabilistic finite element analysis

This paper investigates the use of Fast Probability Integration (FPI) algorithms in a Finite Element environment. A method allowing the representation of correlated fields in terms of a vector of uncorrelated transformed variables, based on the spectral decomposition of the variance-covariance matrix is developed. The response of the deterministic model corresponding to selected perturbations of these uncorrelated variables is then obtained via a Newton-type iterative scheme. The results of the perturbed problems are used to construct a local representation of the model's behavior in the neighborhood of the deterministic state, which the FPI algorithm will use to estimate the reliability of the system. Although the proposed strategy has thus far only been applied to linear elastostatics, the extension of the method to a broader class of problems appears to be feasible.

Dias, J. B.↗

Measured Effects of Turbulence on the Loudness and Waveforms of Conventional and Shaped Minimized Sonic Booms

Turbulence has two distinctive effects on sonic booms: there is distortion in the form of random perturbations that appear behind the shock waves, and shock rise times are increased randomly. A first scattering theory by S.C. Crow in the late 1960s quantified the random distortions, and Crow's theory was shown to agree with available flight test data. A variety of theories for the shock thickness have been presented, all supporting the role of turbulence in increasing rise time above that of a basic molecular-relaxation structure. The net effect of these phenomena on the loudness of shaped minimized booms is of significant interest. Initial analysis suggests that there would be no change to average loudness, but this had not been experimentally investigated. The January 2004 flight test of the Shaped Sonic Boom Demonstrator (SSBD), together with a reference unmodified F-5E, included a 12500- foot linear ground sensor array with 28 digitally recorded sensor sites. This data set provides an opportunity to re-test Crow's theory for the post-shock perturbations, and to examine the net effect of turbulence on the loudness of shaped sonic booms.

Plotkin, Kenneth J.↗

Flavor fragmentation function factorization

A definition of partonic jet flavor that is both theoretically well-defined and experimentally robust would have profound implications for measurements and predictions especially for heavy flavor applications. Recently, a definition of jet flavor was introduced as the net flavor flowing along the direction of the Winner-Take-All axis of a jet which is soft safe to all orders, but not collinear safe. Here, we exploit the lack of collinear safety and propose a factorization theorem of perturbative flavor fragmentation functions that resum collinear divergences and describe the evolution of flavor from the short distance of jet production to the long distance at which hadronization occurs. Collinear flavor evolution is governed by a small modification of the DGLAP equations. We present a detailed all-orders analysis and identify exact relations that must hold amongst the various anomalous dimensions by probability conservation and the existence of fixed points of the renormalization group flow. We explicitly validate the factorization theorem at one-loop order, and demonstrate its consistency at two loops in particular flavor channels. Starting at two-loops, constraints on phase space imposed by flavor measurements potentially allow for non-trivial soft contributions, but we demonstrate that they are scaleless and so explicitly vanish, ensuring that soft particles are summed inclusively and all divergences are exclusively collinear in nature. This factorization theorem opens the door to precision calculations with identified flavor in the infrared.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Error field measurements with rotating RMP fields for DIII-D H-mode

3D magnetic sensors are employed to identify the amplitude and toroidal phase of error fields (EF) by analyzing the torque balance for magnetic islands entrained by rotating resonant magnetic perturbations (RMPs) in DIII-D H-mode plasmas. This technique of torque balance allows for efficient error field identification, offering a valuable tool for scenario-specific and optimized error field compensation (EFC) and requires only magnetic diagnostics. The torque balance used in this work includes the contribution from electromagnetic torque due to error fields, wall response, RMP fields, and viscous. Results show that viscous torque plays a crucial role, particularly during locked modes and H-mode plasmas, ensuring accurate data fits with lower residuals. The torque balance technique reveals that the L- and H-mode plasmas have distinct EF configurations, and consistent with a model-based EF analysis including MHD response in IPEC and the SURFMN EF simulation. This technique shows great robustness in measuring the intrinsic EF amplitude regardless of its amplitude or toroidal phase. Repeated discharges with EFC disparities exhibit consistent results of intrinsic error field within a reasonable range near the “standard” error field compensation. Additionally, the use of a rotating n = 1 resonant magnetic perturbation offers the advantage of reducing disruption risks by entraining saturated magnetic islands. These findings are instrumental for optimizing EF correction in fusion devices, thereby enhancing tearing mode suppression and overall plasma stability.

3D magnetic sensors↗

Surrogate-driven Variance-based Sensitivity Analysis of Thermal Storage Tanks in Integrated Energy Systems

Sensitivity analysis and uncertainty quantification are essential steps for enhancing the accuracy of computational models by identifying and mitigating uncertainties. This study focuses on these steps for the Thermal Energy Delivery System at Idaho National Laboratory, specifically targeting the thermocline tank. Using a Modelica/Dymola simulation model, the study perturbed various design parameters and boundary conditions, including shape factor, porosity, outlet temperature, inlet mass flow rate, and system pressure, to predict and quantify uncertainty in the tank’s ax- ial temperature. A dataset of over 1,000 simulations was generated, and surrogate models were developed using the pyMAISE (Michigan Artificial Intelligence Standard Environment) library, which is an Automatic Machine Learning library for nuclear engineering applications. The optimal model, a feedforward neural network with two hidden layers, achieved an R2 score above 0.99 and a mean absolute error below 1 Kelvin. Sensitivity analyses using Sobol indices and Fourier amplitude sensitivity testing methods on this surrogate model revealed that the inlet mass flow rate at initial timestamps and porosity significantly impacts predicted temperatures across all sensors and time steps.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Linear statistical analysis of maneuver optimization techniques

It is shown how a method for maneuver design during the Jovian satellite tour phase of the Galileo mission requires a full-scale redesign of the tour at each maneuver. Software was developed to provide insight into the statistical properties of this process; it is based on the use of linear perturbations to provide trajectory propagation and a least squares representation of the optimization. The result is a Monte Carlo program, making it possible to evaluate the proposed strategy.

Maize, Earl H.↗

Ameloblastin binding to biomimetic models of cell membranes – A continuum of intrinsic disorder

A 37-residue amino acid sequence corresponding to the segment encoded by exon-5 of murine ameloblastin (Ambn), AB2 (Y67-Q103), has been implicated with membrane association, ameloblastin self-assembly, and amelogenin-binding. Here, our aim was to characterize, at the residue level, the structural behavior of AB2 bound to chemical mimics of biological membranes using NMR spectroscopy. To better define the structure of AB2 using NMR-based methods, recombinant 13 C- and 15 N-labelled AB2 (*AB2) was prepared and data collected free in solution and with deuterated dodecylphosphocholine (dPC) micelles, deuterated bicelles, and both small and large unilamellar vesicles. Amide chemical shift and intensity perturbations observed in 1 H- 15 N HSQC spectra of *AB2 in the presence of bicelles and dPC micelles suggest that a region of *AB2, S6-E36 (murine Ambn S68 – E98), associates with the membrane biomimetics. A CSI-3 analysis of the NMR chemical shift assignments for *AB2 free in solution and bound to dPC micelles indicated the peptide remains disordered except for the adoption of a short, 12-residue α-helix, F10-G21 (murine Ambn F72-G83). In dPC micelles, the NOE NMR data was void of patterns characteristic of long-lived helical structure indicating this helix was transient in nature. A continuum of intrinsic disorder in the membrane-bound state may be responsible for ameloblastin’s ability to dynamically interact with multiple partners at the same site during amelogenesis.

59 BASIC BIOLOGICAL SCIENCES↗

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Erosion and impurity transport for the edge localized mode suppression window in KSTAR

A 3D analysis of plasma wall interactions and global impurity transport for the edge localized mode suppression window in KSTAR during H-Mode discharges has shown that carbon erosion at the divertor plates is a strong function of the resonant magnetic perturbation (n = 1) coil current and relative phasing. The Generalized Perturbation Equilibrium Code was used to determine a realistic initial perturbed plasma equilibrium, and EMC3-EIRENE was used to calculate the resulting scrape-off layer plasma used in this study as a fixed background for the ERO2.0 plasma–material interaction model. The resulting transport leads to deposition of impurities along the targets positioned at the high-field side of the device. An attempt at calculating the resulting effective charge state has demonstrated a similar dependence on the perturbation coil current and has been able to determine a window for the experimentally observed values of Zeff by including contributions of all ionized carbon charge states and deuterium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exponential stability of large-scale discrete systems

The concept of vector Liapunov functions is used to obtain conditions for the exponential stability of large-scale discrete systems which can be decomposed into a number of interconnected subsystems with the same stability property. Both the structurally invariant composite systems and the large-scale systems under structural perturbations are considered. Connective absolute stability of a large-scale system composed of the interconnected Lur'e-type subsystems is defined and resolved in this context, resulting in a computationally and conceptually attractive alternative to a straightforward stability analysis of the system by frequency-domain criteria.

Grujic, L. T.↗

Revised Perturbation Statistics for the Global Scale Atmospheric Model

Magnitudes and scales of atmospheric perturbations about the monthly mean for the thermodynamic variables and wind components are presented by month at various latitudes. These perturbation statistics are a revision of the random perturbation data required for the global scale atmospheric model program and are from meteorological rocket network statistical summaries in the 22 to 65 km height range and NASA grenade and pitot tube data summaries in the region up to 90 km. The observed perturbations in the thermodynamic variables were adjusted to make them consistent with constraints required by the perfect gas law and the hydrostatic equation. Vertical scales were evaluated by Buell's depth of pressure system equation and from vertical structure function analysis. Tables of magnitudes and vertical scales are presented for each month at latitude 10, 30, 50, 70, and 90 degrees.

Justus, C. G.↗

Geographical representation of radial orbit perturbations due to ocean tides: Implications for satellite altimetry

In analogy to the geographical representation of the zeroth-order radial orbit perturbations due to the static geopotential, similar relationships have been derived for radial orbit perturbations due to the ocean tides. At each location these perturbations are seen to be coherent with the tide height variations. The study of this singularity is of obvious importance to the estimation of ocean tides from satellite altimeter data. We derive analytical expressions for the sensitivity of altimeter derived ocean tide models to the ocean tide force model induced errors in the orbits of the altimeter satellite. In particular, we focus on characterizing and quantifying the nonresonant tidal orbit perturbations, which cannot be adjusted into the empirical accelerations or radial perturbation adjustments commonly used during orbit determination and in altimeter data processing. As an illustration of the utility of this technique, we study the differences between a TOPEX/POSEIDON-derived ocean tide model and the Cartwright and Ray 1991 Geosat model. This analysis shows that nearly 60% of the variance of this difference for M(sub 2) can be explained by the Geosat radial orbit eror due to the omission of coefficients from the GEM-T2 background ocean tide model. For O(sub 1), K(sub 1), S(sub 2), and K(sub 2) the orbital effects account for approximately 10 to 40% of the variances of these differences. The utility of this technique to assessment of the ocean tide induced errors in the TOPEX/POSEIDON-derived tide models is also discussed.

Bettadpur, Srinivas V.↗

Effects of Increasing Drag on Conjunction Assessment

Conjunction Assessment Risk Analysis relies heavily on the computation of the Probability of Collision (Pc) and the understanding of the sensitivity of this calculation to the position errors as defined by the covariance. In Low Earth Orbit (LEO), covariance is predominantly driven by perturbations due to atmospheric drag. This paper describes the effects of increasing atmospheric drag through Solar Cycle 24 on Pc calculations. The process of determining these effects is found through analyzing solar flux predictions on Energy Dissipation Rate (EDR), historical relationship between EDR and covariance, and the sensitivity of Pc to covariance. It is discovered that while all LEO satellites will be affected by the increase in solar activity, the relative effect is more significant in the LEO regime around 700 kilometers in altitude compared to 400 kilometers. Furthermore, it is shown that higher Pc values can be expected at larger close approach miss distances. Understanding these counter-intuitive results is important to setting Owner/Operator expectations concerning conjunctions as solar maximum approaches.

Frigm, Ryan Clayton↗

Land Processes Can Substantially Impact the Mean Climate State

Terrestrial processes influence the atmosphere by controlling land-to-atmosphere fluxes of energy, water, and carbon. Prior research has demonstrated that parameter uncertainty drives uncertainty in land surface fluxes. However, the influence of land process uncertainty on the climate system remains underexplored. Here, we quantify how assumptions about land processes impact climate using a perturbed parameter ensemble for 18 land parameters in the Community Earth System Model version 2 under preindustrial conditions. We find that an observationally-informed range of land parameters generate biogeophysical feedbacks that significantly influence the mean climate state, largely by modifying evapotranspiration. Global mean land surface temperature ranges by 2.2°C across our ensemble (σ = 0.5°C) and precipitation changes were significant and spatially variable. Our analysis demonstrates that the impacts of land parameter uncertainty on surface fluxes propagate to the entire Earth system, and provides insights into where and how land process uncertainty influences climate.

54 ENVIRONMENTAL SCIENCES↗

Correlative Surface Imaging Reveals Chemical Signatures for Bacterial Hotspots on Plant Roots

The rhizosphere is arguably the most complex microbial habitat on earth, comprising an integrated network of plant roots, soil and a highly diverse microbial community (the rhizosphere microbiome)1. Understanding, predicting and controlling plant-microbe interactions in the rhizosphere will allow us to harness the plant microbiome as a means to increase or restore plant ecosystem productivity, improve plant responses to a wide range of environmental perturbations, and mitigate effects of climate change by designing ecosystems for long-term soil carbon storage1. To this end, it is imperative to develop new molecular approaches with high spatial resolution to capture interactions at the plant-microbe, microbe-microbe, and plant-plant interfaces. In this work, we designed an imaging sample holder that allows integrated surface imaging tools to map the same locations of a plant root-microbe interface with submicron lateral resolutions, providing novel in vivo analysis of root-microbe interactions. Specifically, confocal fluorescence microscopy, time-of-flight secondary ion mass spectrometry (ToF-SIMS), X-ray photoelectron spectroscopy (XPS), and scanning electron microscope (SEM) were used for the first time for correlative imaging of the Brachypodium distachyon root and its interaction with Pseudomonas SW25, a typical plant growth-promoting soil bacterium. Imaging data suggest that the root surface is inhomogeneous and that the interaction between Pseudomonas and Brachypodium roots was confined to only few spots along the sampled root segments and that the bacterial attachment spots were enriched in Na- and S-related and high-mass organic species. We conclude that attachment of the Pseudomonas cells to the root surface is outcompeted by strong root-soil mineral interactions but facilitated by formation of extracellular polymeric substances (EPS).

Plant root, Plant growth-promoting bacteria, soil,↗

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗

Grassmannian Shape Representations for Aerodynamic Applications: Preprint

Airfoil shape design is a classical problem in engineering, science, and manufacturing. Our motivation is to combine principled physics-based considerations for the shape design problem with modern computational techniques informed by a data-driven approach. Traditional analyses of airfoil shapes emphasize a flow-based sensitivity to deformations which can be represented generally by affine transformations (rotation, scaling, shearing, shifting). We present a novel representation of shapes which decouples affine-style deformations from a rich set of data-driven deformations over a submanifold of the Grassmannian. The Grassmannian representation, informed by a database of physically relevant airfoils, offers (i) a rich set of novel 2D airfoil deformations not previously captured in the data, (ii) improved low-dimensional parameter domain for inferential statistics, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal shapes.

blade representation↗