Search NASA⌕ Search

SEARCH · Search NASA

Results for “sensitivity gradient”

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 163 records · Page 9

Handling Iterative Solvers in an Algorithmic Differentiation Framework Using Implicit Methods

Differentiable programming is a powerful concept as it enables the seemly propagation of gradients through functions, algorithms, and/or whole physics simulations. These gradients are useful for a wide variety of applications, including sensitivity studies and machine learning, but one of particular interest is optimization. Gradient-based optimization, enabled through automatic/algorithmic differentiation (AD), can be used on predictive physical models to efficiently optimize a set of design variables. AD methods are a particularly promising approach to complex physics simulations because they can be shown to scale well with an increasing number of design variables; however, care must be taken when coupling between different models or different states of a single model.

algorithmic differentiation↗

Capturing the fractocohesive length scale in elastomers through a statistical mechanics-based gradient enhanced damage model

Here, this study aims to examine modeling flaw sensitivity in elastomers. The direct incorporation of polymer chain statistical mechanics considerations into a continuum stretch-based gradient-enhanced damage formulation, in turn, allows a representation of diffuse chain damage and fracture events that align with known micromechanical mechanisms. Through a series of numerical experiments, we simulate crack propagation and extract the fracture energy as an output of the model, while keeping track of the micromechanical signatures of diffuse chain damage that accommodate fracture propagation and eventually influence flaw sensitivity. Finally, by combining the fracture toughness and the work to rupture, we identify a fractocohesive length of the material, corresponding to the full width of the damage process zone. As the damage-to-fracture cascade in the proposed GED model is influenced by the introduction of a length scale associated with network imperfection and long-range load transfer, the emerging relationship of the two length scales is discussed, providing a potential link between microscopic damage mechanisms and the observed macroscopic fracture response.

36 MATERIALS SCIENCE↗

Gyrokinetic simulations of core turbulence and thermal transport in the high-β P discharge on EAST

The properties of core turbulence and thermal transport are investigated for EAST high-β P (β P ~ 3.1) plasmas with dominant electron heating (Te/Ti > 1) via gyrokinetic simulation with the NLT code. Linear simulations identify that the electrostatic η e -driven trapped electron mode (η e -TEM) dominates in the core region (ρ< 0.7) and the ion temperature gradient (ITG) mode dominates in the out region (ρ≥ 0.7), consistent with the linear threshold analysis of micro-instabilities. Sensitivity analysis shows that the normalized electron density gradient (R/L ne ) and ITG (R/L Ti ) are two effective parameters to stabilize TEM instability. Nonlinear simulations are also carried out and compared with the experimental results, which show that the electron thermal internal transport barrier (ITB) in EAST high-β P plasma is determined by the TEM-induced turbulence. A higher zonal flow shearing rate is observed in the ITB region (0.2 < ρ < 0.34), which can regulate energy transport induced by TEM turbulence and facilitate the formation of e-ITB. A plausible positive feedback mechanism to mitigate turbulence transport and improve energy confinement via enhanced ion heating is proposed for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating profile stiffness and critical gradients in shaped TCV discharges using local gyrokinetic simulations of turbulent transport

The experimental observation made on the TCV tokamak of a significant confinement improvement in plasmas with negative triangularity ( δ < 0) compared to those with standard positive triangularity has been interpreted in terms of different degrees of profile stiffness (Sauter et al. Phys. Plasmas 21 055906, 2014) and/or different critical gradients. Employing the Eulerian gyrokinetic code GENE (Jenko et al. Phys. Plasmas 7 1904, 2000), profile stiffness and critical gradients are studied under TCV relevant conditions. For the considered experimental discharges, trapped electron modes (TEMs) and electron temperature gradient (ETG) modes are the dominant microinstabilities, with the latter providing a significant contribution to the non-linear electron heat fluxes near the plasma edge. Two series of simulations with different levels of realism are performed, addressing the question of profile stiffness at various radial locations. Retaining finite collisionality, impurities and electromagnetic effects, as well as the physical electron-to-ion mass ratio are all necessary in order to approach the experimental flux measurements. However, flux-tube simulations are unable to fully reproduce the TCV results, pointing towards the need to carry out radially nonlocal (global) simulations, i.e. retaining finite machine size effects, in a future study. Some conclusions about the effect of triangularity can nevertheless be drawn based on the flux-tube results. In particular, the importance of considering the sensitivity to both temperature and density gradient is shown. The flux tube results show an increase of the critical gradients towards the edge, further enhanced when δ < 0, and they also appear to indicate a reduction of profile stiffness towards plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observational Signatures of Transverse Magnetohydrodynamic Waves and Associated Dynamic Instabilities in Coronal Flux Tubes

Magnetohydrodynamic (MHD) waves permeate the solar atmosphere and constitute potential coronal heating agents. Yet, the waves detected so far may be but a small subset of the true existing wave power. Detection is limited by instrumental constraints but also by wave processes that localize the wave power in undetectable spatial scales. In this study, we conduct 3D MHD simulations and forward modeling of standing transverse MHD waves in coronal loops with uniform and non-uniform temperature variation in the perpendicular cross-section. The observed signatures are largely dominated by the combination of the Kelvin–Helmholtz instability (KHI), resonant absorption, and phase mixing. In the presence of a cross-loop temperature gradient, we find that emission lines sensitive to the loop core catch different signatures compared to those that are more sensitive to the loop boundary and the surrounding corona, leading to an out-of-phase intensity and Doppler velocity modulation produced by KHI mixing. In all of the considered models, common signatures include an intensity and loop width modulation at half the kink period, a fine strand-like structure, a characteristic arrow-shaped structure in the Doppler maps, and overall line broadening in time but particularly at the loop edges. For our model, most of these features can be captured with a spatial resolution of 0.″33 and a spectral resolution of 25 km s{sup −1}, although we do obtain severe over-estimation of the line width. Resonant absorption leads to a significant decrease of the observed kinetic energy from Doppler motions over time, which is not recovered by a corresponding increase in the line width from phase mixing and KHI motions. We estimate this hidden wave energy to be a factor of 5–10 of the observed value.

79 ASTRONOMY AND ASTROPHYSICS↗

Integration of high confinement, high poloidal beta plasma with dual radiated power and detachment controls for divertor protection and ELM suppression

Divertor detachment without serious core confinement quality loss in DIII-D’s high poloidal β scenario has been combined with impurity-induced ELM mitigation without disruption. Use of Ne previously granted access to a detached, non-ELMing regime that retained high confinement quality due to stimulation of Internal Transport Barrier (ITB) growth, but suffered from on-going core fuel dilution and high disruptivity. Excess Ne accumulation in the core plasma has now been avoided by feeding back core radiated power ($P_{rad}$) measurements to control Ne seeding, rather than using attachment fraction ($A_{frac}$) control with Ne; this also reduces disruptivity. At the same time, N 2 seeding is used in a feedback loop with $A_{frac}$ measurements, which previously posed low disruption risk. In this way, the effect of Ne in the core is managed while avoiding excess seeding, and N 2 acts to correct for any excess heat exhaust that might interfere with detachment. The average Ne flow rate was 38% of what was used in pure Ne $A_{frac}$ control, plus average N 2 flow that was 43% of pure N 2 $A_{frac}$ control, even while meeting an even deeper detachment target. Meanwhile, a steep pressure gradient in the core plasma reduces sensitivity to impurity-induced degradation of the pedestal and reasonable confinement quality was maintained despite operational challenges that blocked formation of an ITB in these experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct Observation of the Local Reaction Environment during the Electrochemical Reduction of CO 2

The electrochemical reduction of carbon dioxide is sensitive to electrolyte polarization, which causes gradients in pH and the concentration of carbon dioxide to form near the cathode surface. It is desirable to measure the concentration of reaction-relevant species in the immediate vicinity of the cathode because the intrinsic kinetics of carbon dioxide reduction depend on the composition of the local reaction environment. Meeting this objective has proven difficult because conventional analytical methods only sample products from the bulk electrolyte. In this study, we describe the use of differential electrochemical mass spectrometry to measure the concentration of carbon dioxide and reaction products in the immediate vicinity of the cathode surface. This capability is achieved by coating the electrocatalyst directly onto the pervaporation membrane used to transfer volatile species into the mass spectrometer, thereby enabling species to be sampled directly from the electrode–electrolyte interface. This approach has been used to investigate hydrogen evolution and carbon dioxide reduction over Ag and Cu. We find that the measured CO 2 reduction activity of Ag agrees well with what is measured by gas chromatography of the effluent from an H-cell operated with the same catalyst and electrolyte. A distinct advantage of our approach is that it enables observation of the depletion of carbon dioxide near the cathode surface due to reaction with hydroxyl anions evolved at the cathode surface, something that cannot be done using conventional analytical techniques. We also demonstrate that the influence of this relatively slow chemical reaction can be minimized by evaluating electrocatalytic activity during a rapid potential sweep, thereby enabling measurement of the intrinsic kinetics. For CO 2 reduction over Cu, nine products can be observed simultaneously in real time. A notable finding is that the abundance of aldehydes relative to alcohols near the cathode surface is much higher than that observed in the bulk electrolyte. It is also observed that for increasingly cathodic potentials the relative abundance of ethanol increases at the expense of propionaldehyde. Finally, these findings suggest that acetaldehyde is a precursor to ethanol and propionaldehyde and that propionaldehyde is a precursor to n-propanol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Close coupling of the peeling ballooning and external kink modes in high poloidal-beta discharges with strong shape and large pedestal width

In plasmas with strong shape and large pedestal width, like DIII-D high poloidal-beta discharges, the eigenvalue spectrum of the peeling ballooning mode (PBM) is shown to shift to the very low-n regime near one, where n is the toroidal mode number. This peeling-type eigenvalue spectrum has a further shift to n = 1, thus being smoothly connected to the n = 1 external kink mode (EKM), as the normalized beta (β N ) increases. Once this connection occurs, the mode takes a mixture form of the PBM and EKM with its mode structure varying from the PBM-like to the EKM-like one as β N increases. The mode stability also becomes sensitive to both the local pedestal gradient and global β N , thus allowing an anti-correlation between the two driving forces. Further, these results appear to provide a qualitative explanation of the two unexpected features observed in the DIII-D high poloidal-beta discharges, that is, the dominance of the n = 1 mode in the edge-localized-modes and the negative correlation between the pedestal height and the internal transport barrier strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

$Z$-Pinch Interferometry Analysis With the Fourier-Based TNT Code

Here, we present the analysis of interferometry diagnostics with the user-friendly Talbot Numerical Tool (TNT), a Fourier-based postprocessing code that enables real-time assessment of plasma systems. TNT performance was explored with visible and infrared interferometry in pulsed-power-driven Z -pinch configurations to expand its capabilities beyond Talbot X-ray interferometry in the high-intensity laser environment. TNT enabled accurate electron density characterization of magnetically driven plasma flows and shocks through phase-retrieval methods that did not require data modification or masking. TNT demonstrated enhanced resolution, detecting below 4 % fringe shift, which corresponds to 8.7 × 10 15 cm –2 within 28 μ m, approaching the laser probing system limit. TNT was tested against a well-known interferometry analysis software, delivering an average resolving power nearly ten times better (~28 μ m versus ~ 210 μ m) when resolving plasma ablation features. TNT demonstrated higher sensitivity when probing sharp electron density gradients in supersonic shocks. A maximum electron areal density of 4.1 × 10 17 cm –2 was measured in the shocked plasma region, and a minimum electron density detection of ~ 1.0 × 10 15 cm –2 was achieved. When probing colliding plasma flows, the calculations of the effective adiabatic index and the associated errors were improved from γ* = 2.6 ± 1.6 – 1.4 ± 0.2 with TNT postprocessing, contributing valuable data for the interpretation of radiative transport. Additional applications of TNT in the characterization of pulsed-power plasmas and beyond are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Feasibility study of three-dimensional holographic interferometry for aerodynamics

Laser holographic interferometry was tried as a nonintrusive diagnostic tool for studying unsteady two dimensional flows. A NACA 0012 airfoil was tested, while undergoing dynamic stall, over a range of Mach numbers of 0.3 to 0.5, Reynolds number of 0.5 x 10 to the 6th power to 2 x 10 to the 6th power, and at reduced frequencies of 0.015 to 0.15. It was found that both quantitative and qualitative data could be obtained by the technique. Surface pressures on the airfoil can be measured to within 1% of those measured with orifices and pressure transducers when the flow is attached. Velocity profiles were measured near the wake region, and they compared very well with laser velocimeter data for attached flows. For separated flows with large scale vortices, densities can be measured, but pressures and velocities cannot be deduced with the assumption of constant pressure gradient in the normal direction. The sensitivity of the interferograms was good at a Mach number of 0.4 and a Reynolds number of 4 x 10 to the 6th power/ft; the sensitivity worsened at smaller Mach numbers and Reynolds numbers, and improved at larger ones.

Craig, J. E.↗

Investigation of the influence of wind shear on the aerodynamic characteristics of aircraft using a vortex-lattice method

The objective was to investigate and characterize the aerodynamic effect of shear flow through a series of sensitivity studies of the wind velocity gradients and wing planform geometry parameters. The wind shear effect was computed using a modified vortex-lattice computer program and characterized through the formulation of wind shear aerodynamic coefficients. The magnitude of the aerodynamic effect was demonstrated by computing the resultant change in the aerodynamics of a conventional wing and tail combination on a fixed flight path through a simulated microburst. The results of the study indicate that a significant amount of the control authority of an airplane may be required to counteract the wind shear induced forces and moments in the microburst environment.

Vicroy, Dan D.↗

Experience with duplex bearings in narrow angle oscillating applications

Duplex ball bearings are matched pairs on which the abutting faces of the rings have been accurately ground so that when the rings are clamped together, a controlled amount of interference (preload) exists across the balls. These bearings are vulnerable to radial temperature gradients, blocking in oscillation and increased sensitivity to contamination. These conditions decrease the service life of these bearings. It was decided that an accelerated thermal vacuum life test should be conducted. The test apparatus and results are described and the rationale is presented for reducing a multiyear life test on oil lubricated bearings to less than a year.

Phinney, D. D.↗

Influence of wind shear on the aerodynamic characteristics of airplanes

The aerodynamic effect of shear flow through a series of sensitivity studies of the wind velocity gradients and wing planform geometry parameters is investigated and characterized. The wind shear effect is computed using a modified vortex-lattice computer program and characterized through the formulation of wind shear aerodynamic coefficients. The magnitudes if the aerodynamic effects are demonstrated by computation of the resultant change in the aerodynamics of a conventional wing and horizontal stability configuration on a fixed flight path through a simulated microburst. The results indicate that as much as 20 percent of the control authority of the airplane may be required to counteract the wind-shear-induced forces and moments in the microburst environment.

Vicroy, Dan D.↗

Active control law synthesis for flexible aircraft

An application of an active control law synthesis procedure is presented, to meet multiple design requirements for a flexible aircraft modeled by a large order state space system. In this synthesis procedure, a linear quadratic Gaussian type cost function is minimized by updating the free parameters of the control law, while satisfying a set of constraints on the design loads, responses and stability margins. Analytical expressions for gradients of the cost function and the constraints, with respect to the control law design variables are used to facilitate rapid numerical convergence. These gradients can also be used for sensitivity study. A stable classical control law as well as an estimator-based full or reduced order control law can be modified, in order to meet individual root-mean-square response limitations as well as minimum singular value restrictions. Both analog and digital control laws can be optimized. Low order, robust control laws were synthesized for flutter suppression of a flexible aircraft.

Mukhopadhyay, V.↗

The chromospheres of late-type stars. I - Eridani as a test case of multiline modelling

A new model of the lower chromosphere of the dwarf K2 star Epsilon Eridani is derived by matching flux profiles of the Ca IR triplet lines 8498 and 8542 A H-alpha and H-beta lines and the Na D lines (all observed simultaneously at the AAT), and the Ca II K line. The coupled non-LTE equations of statistical equilibrium and radiative transfer are solved under the constraint of hydrostatic equilibrium using the Carlsson (1986) code. Within the framework of the model, the Na D lines are an important photospheric diagnostic, and the Ca IR triplet lines can be used to locate the temperature minimum. The computed H-alpha and H-beta depths are highly sensitive constraints on the transition zone gradients and base pressures allowing us to derive a pressure at the base of the transition zone of 0.9 dyn/cm.

Thatcher, John D.↗

On the use of forbidden lines as density diagnostics in stratified media

The influence of density and ionization gradients on the observed density-sensitive forbidden-line ratios from a spatially unresolved source is investigated by means of a simple model. In order to assess the role played by the source geometry, both spherical and cylindrical systems where all spatial dependences are given by power laws are considered. It is shown that the present ratios can be used to estimate the spatial dependence of density and ionization as well as the source geometry.

Safier, Pedro N.↗

Parametric Deformation of Discrete Geometry for Aerodynamic Shape Design

We present a versatile discrete geometry manipulation platform for aerospace vehicle shape optimization. The platform is based on the geometry kernel of an open-source modeling tool called Blender and offers access to four parametric deformation techniques: lattice, cage-based, skeletal, and direct manipulation. Custom deformation methods are implemented as plugins, and the kernel is controlled through a scripting interface. Surface sensitivities are provided to support gradient-based optimization. The platform architecture allows the use of geometry pipelines, where multiple modelers are used in sequence, enabling manipulation difficult or impossible to achieve with a constructive modeler or deformer alone. We implement an intuitive custom deformation method in which a set of surface points serve as the design variables and user-specified constraints are intrinsically satisfied. We test our geometry platform on several design examples using an aerodynamic design framework based on Cartesian grids. We examine inverse airfoil design and shape matching and perform lift-constrained drag minimization on an airfoil with thickness constraints. A transport wing-fuselage integration problem demonstrates the approach in 3D. In a final example, our platform is pipelined with a constructive modeler to parabolically sweep a wingtip while applying a 1-G loading deformation across the wingspan. This work is an important first step towards the larger goal of leveraging the investment of the graphics industry to improve the state-of-the-art in aerospace geometry tools.

Anderson, George R.↗

TPSAS-NF1676L-31914-DND

Analysis of aircraft systems involves many disciplines: aerodynamics, structural dynamics, aeroacoustics, etc. - rotorcraft aeromechanics is a typical example. Gradient-based, multidisciplinary optimization (MDO) requires sensitivity analysis for all disciplines. Adjoint-based sensitivity analysis capabilities for: - Coupled aero/structure FUN3D/DYMORE system; - Coupled aero/acoustics FUN3D/ANOPP2 system.

Li Wang↗