Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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,153 records · Page 64

Series FACTS Devices for Increasing Resiliency in Severe Weather Conditions

Severe weather conditions are low-probability, high-impact events that affect grid operations. The majority of power outages are caused by severe weather conditions. Grid resiliency to weather events can be enhanced by decreasing the reliance on its affected sections. One way to do this is to reduce the power flow through lines vulnerable to severe weather. If a line is disconnected, its initial power flow is distributed through the neighbor lines, which may cause congestion in the grid. FACTS devices can be used to control the power flow of lines that have a higher chance of power outages. Most previous works do not consider weather events in power flow control. In this work, a linearized optimal power flow (OPF)–based algorithm is developed to minimize the real power flow of vulnerable lines considering the thermal limits of lines to prevent infeasible solutions; the simulation is fast, making it suitable for large-scale systems. The proposed optimization problem is presented as a mixed-integer linear program (MILP), making it capable of using short-term load forecasting due to its high solution speed. The proposed optimization problem considers multiple lines with different outage probabilities and the uncertainties of the weather forecast. Moreover, it estimates the power reduction in vulnerable lines due to changes in the series FACTS devices. The performance of the proposed optimization problem is tested on IEEE 14-, 30-, and 118-bus systems for several scenarios. The results are validated with the AC power flow results from MATPOWER.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Aeroacoustic computation of gust-blade interaction

To better understand and address the challenges faced in computing the acoustics of flow fields, test problems must be considered. In the present study, the sound radiated by the interaction of a flat plate with an oncoming gust containing a two component, mean velocity is computed. The gust has a uniform mean flow in x with Mach number M(infinity) equal to 0.5. The gust's mean velocity in y is of smaller amplitude and is given by: v = 0.1 sin(pi/8(x/M(sub infinity) - t)). This problem has been posed for an upcoming ICASE/LaRC workshop on benchmark problems in computational aeroacoustics. A plate with a length of 30 units in x is used. The plate is assumed to be infinitesimally thin and is centered at the origin. All variables are made dimensionless using the scales specified. Acoustic quantities are obtained by numerically integrating the linearized Euler equations. Integration is performed on the computational domain -100.0 less than or equal to x less than or equal to 100.0, -100.0 less than or equal to y less than or equal to 100.0, using unit length grid spacing in x and in y. An integration scheme is sought which will provide accurate solution to the small quantities of interest at a minimal computational expense. Results indicate that with the given discretization a scheme of minimal fourth order accuracy might be adequate to approximate the waves within the given flow. Thus, a variation of the MacCormack scheme with fourth order accuracy in space and second order accuracy in time was chosen. A scheme with sixth order accuracy in space has also been implemented and results compared with those of the fourth order accurate scheme. To ensure no mass flux, zero normal velocity is assigned at the plate. This condition will induce a discontinuity in the pressure across the plate location. Values for the perturbation pressure p' along the surface of the plate are obtained using a one-sided, third order Taylor expansion, such that p'(sub y) = O. In accordance with the Kutta condition, perturbation pressure at the trailing edge is assigned to zero. In the far field, radiation boundary conditions have been implemented. The effectiveness of the far field conditions are validated by computing in a larger computational domain and comparing the results. Early in time, sound waves begin to radiate from the plate. A Doppler effect is observed. After the initial transients disappear, the strongest waves leave the trailing edge at an approximate 45 degree angle. The intensity pattern of pressure fluctuations shows five lobes (of increasing magnitude with increasing downstream direction) emerging from the plate. Undesirable short wave contaminants are observed in the computed pressure distribution along the plate surface. For a more accurate solution at small scales, a more refined discretization will be required.

Martin, James E.↗

Equilibrium Separation of Siloxanes in Metal–Organic Frameworks

We present an in silico assessment of metal–organic frameworks (MOFs) for the equilibrium separation of linear and cyclic siloxanes. Using a combination of configurational bias/continuous fractional component Monte Carlo (CB/CFCMC) simulations and Ideal Adsorbed Solution Theory (IAST), we investigated the adsorption of both equimolar and nonequimolar mixtures of linear and cyclic siloxanes in a selection of synthesizable MOFs with medium to large pore volumes. We showed that configurational entropy effects drive the preferential adsorption of linear siloxanes over cyclic siloxanes. Based on synthesizability metrics we identified ZIF-70 as a promising adsorbent for the separation of linear and cyclic siloxanes. Self-diffusivities of linear and cyclic siloxanes in ZIF-70 calculated from molecular dynamics simulations show that equilibrium can be reached on reasonable time scales. We explored vacuum-temperature swing adsorption (VTSA) as a potential method for the recovery of adsorbed linear siloxanes from ZIF-70, achieving a 75% reduction in adsorbed phase concentration. Additionally, we demonstrated that supercritical CO 2 offers an alternative desorption strategy by displacing adsorbed linear siloxanes in ZIF-70 at pressures above 200 bar driven by entropy effects.

Adsorption↗

An improved shooting method for solving minimum-time maneuver problems

The minimum-time attitude maneuvers of a rigid spacecraft are considered. Generally, the associated nonsingular, nonlinear two-point boundary-value problem, derived by using Pontryagin's Maximum Principle, can be solved through shooting methods to find the switching times for the bang-bang control. However, a good initial guess for the missing initial costates is important because the convergence range of the initial costates is very narrow, especially when a multi-bang-bang control strategy is needed. This paper proposes a numerical approach to generate the initial costates so that they are very close to the converged values. The initial costates are first generated by using a recently developed quasi-linearization technique for the near-minimum-time problem and then reduced in magnitude, through some properly designed scales, before the shooting method is used. Numerical examples indicate the closeness of the results from these two techniques (quasi-linearization and shooting).

Li, Feiyue↗

Development of the Transport Class Model (TCM) Aircraft Simulation From a Sub-Scale Generic Transport Model (GTM) Simulation

A six degree-of-freedom, flat-earth dynamics, non-linear, and non-proprietary aircraft simulation was developed that is representative of a generic mid-sized twin-jet transport aircraft. The simulation was developed from a non-proprietary, publicly available, subscale twin-jet transport aircraft simulation using scaling relationships and a modified aerodynamic database. The simulation has an extended aerodynamics database with aero data outside the normal transport-operating envelope (large angle-of-attack and sideslip values). The simulation has representative transport aircraft surface actuator models with variable rate-limits and generally fixed position limits. The simulation contains a generic 40,000 lb sea level thrust engine model. The engine model is a first order dynamic model with a variable time constant that changes according to simulation conditions. The simulation provides a means for interfacing a flight control system to use the simulation sensor variables and to command the surface actuators and throttle position of the engine model.

Hueschen, Richard M.↗

Dispersive interactions between standard and Dirac materials and the role of dimensionality

The van der Waals (vdW) interaction plays a prominent role between neutral objects at separations where short ranged chemical forces are negligible. This type of dispersive coupling is determined by the interplay between geometry and response properties of the materials making up the objects. Here, we investigate the vdW interaction between 1D, 2D, and 3D standard and Dirac materials within the Random Phase Approximation, which takes into account collective excitations originating from the electronic Coulomb potential. A comprehensive understanding of characteristic functionalities and scaling laws are obtained for systems with parabolic energy dispersion (standard materials) and crossing linear bands (Dirac materials). By comparing the quantum mechanical and thermal limits the onset of thermal fluctuations in the vdW interaction is discussed showing that thermal effects are significantly pronounced at smaller scales in reduced dimensions.

36 MATERIALS SCIENCE↗

The saturation of the electron-cyclotron maser instability and the interpretation of solar millisecond spikes

A self-consistent numeric two-dimensional code of the kinetic wave-particle equations developed to investigate the maser dynamics in the solar context is applied to solar millisecond-spike observations in order to improve the diagnostic capabilities of the theory of the electron-cyclotron maser instablitity. Attention is given to the inhomogeneity of the magnetic field selecting magneto-ionic modes with relatively short saturation lengths and suppressing mechanisms such as collisional deflection, free-free absorption, and gyroresonance absorption. The time scales of maser saturation in respect to time scales of global particle changes in a magnetic loop are covered, relevant observations of solar millisecond spikes are described, and the interpretation in terms of physical parameters deduced from the quasi-linear maser simulations are presented. It is demonstrated that the quasi-linear simulations make it possible to constrain the physical parameters from the observed time scale and frequency.

Aschwanden, M. J.↗

Optimally Scaled H(sub infinity) Full Information Control Synthesis with Real Uncertainty

This paper presents an algorithm to synthesize optimal controllers for the scaled H(sub infinity). full information problem with real and complex uncertainty. The control problem is reduced to a linear matrix inequality which can be solved via a finite dimensional convex optimization. This technique is compared with the optimal scaled H(sub infinity). full information with only complex uncertainty and D - K iteration control design to synthesize controllers for a missile autopilot. Directly including real parametric uncertainty into the control design results in improved robust performance of the missile autopilot. The controller synthesized via D - K iteration achieves results similar to the optimal designs.

Balas, Gary J.↗

A two-time-scale autopilot for high-performance aircraft

A two-time-scale autopilot is proposed for the Aircraft Controls Design Challenge problem. This control law uses a nonlinear aircraft model constructed from the given vehicle simulation. The vehicle model is partitioned into slow translational dynamics and fast rotational dynamics. Feedback linearization is then employed to synthesize control laws for these two-time scales. Due to the nature of the synthesis, the control law is suitable for automatic trajectory following, and also for pilot control.

Menon, P. K. A.↗

Data-Driven Smoothers for Extreme-Scale Computing

Patch-based relaxation refers to a family of methods for solving linear systems which partitions the matrix into smaller pieces often corresponding to groups of adjacent degrees of freedom residing within patches of the computational domain. The two most common families of patch-based methods are block-Jacobi and Schwarz methods, where the former typically corresponds to non-overlapping domains and the later implies some overlap. We focus on cases where each patch consists of the degrees of freedom on a finite element method mesh cell. Patch methods often capture complex local physics much more effectively than simpler point-smoothers such as Jacobi; however, forming, inverting, and applying each patch can be prohibitively expensive in terms of both storage and computation time. To this end, we propose several approaches for performing analysis on these patches and constructing a reduced representation. The compression techniques rely on either matrix norm comparisons or unsupervised learning via a clustering approach. We illustrate how it is frequently possible to retain/factor less than 5% of all patches and still develop a method that converges only a little slower than when all patches are stored/factored.

97 MATHEMATICS AND COMPUTING↗

Scaling, elasticity, and CLPT

The first few viewgraphs describe the general solution properties of linear elasticity theory which are given by the following two statements: (1) for stress B.C. on S(sub sigma) and zero displacement B.C. on S(sub u) the altered displacements u(sub i)(*) and the actual stresses tau(sub ij) are elastically dependent on Poisson's ratio nu alone: thus the actual displacements are given by u(sub i) = mu(exp -1)u(sub i)(*); and (2) for zero stress B.C. on S(sub sigma) and displacement B.C. on S(sub u) the actual displacements u(sub i) and the altered stresses tau(sub ij)(*) are elastically dependent on Poisson's ratio nu alone: thus the actual stresses are given by tau(sub ij) = E tau(sub ij)(*). The remaining viewgraphs describe the minimum parameter formulation of the general classical laminate theory plate problem as follows: The general CLT plate problem is expressed as a 3 x 3 system of differential equations in the displacements u, v, and w. The eighteen (six each) A(sub ij), B(sub ij), and D(sub ij) system coefficients are ply-weighted sums of the transformed reduced stiffnesses (bar-Q(sub ij))(sub k); the (bar-Q(sub ij))(sub k) in turn depend on six reduced stiffnesses (Q(sub ij))(sub k) and the material and geometry properties of the k(sup th) layer. This paper develops a method for redefining the system coefficients, the displacement components (u,v,w), and the position components (x,y) such that a minimum parameter formulation is possible. The pivotal steps in this method are (1) the reduction of (bar-Q(sub ij))(sub k) dependencies to just two constants Q(*) = (Q(12) + 2Q(66))/(Q(11)Q(22))(exp 1/2) and F(*) - (Q(22)/Q(11))(exp 1/2) in terms of ply-independent reference values Q(sub ij); (2) the reduction of the remaining portions of the A, B, and D coefficients to nondimensional ply-weighted sums (with 0 to 1 ranges) that are independent of Q(*) and F(*); and (3) the introduction of simple coordinate stretchings for u, v, w and x,y such that the process is neatly completed.

Brunelle, Eugene J.↗

30 years of semiconductor nanowire research: A Personal Journey

This project is about flux line pinning in high T c superconductors using columnar defects as the pinning centers. The basic idea is that linear structural defects can pin the flux lines within HTSCs so that their critical current density (the maximal current density that these HTSCs can carry with zero resistance) can be enhanced. The question is how to create such linear structural defects within crystalline cuprate thin films. Scientists have already successfully demonstrated the concept by irradiating HTSC crystals with high-energy ions at synchrotron beamline facilities. However, such heavy ion irradiation is practically not viable for large-scale processes. Finally, we reasoned that maybe we could introduce some sort of linear nanoscopic objects into HTSC to form a composite, and then these linear nanoscopic objects would serve as the magnetic flux line pinning centers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Porous Flow Modeling of Axial Gas Redistribution in Fragmented LWR Fuel Rods using MOOSE

Understanding how gas axially redistributes within fragmented fuel pellets is crucial for predicting the behavior of Light Water Reactor (LWR) fuel rods, particularly during transient and accidental scenarios. The time scale of this phenomenon plays a fundamental role in determining the progression and hazard of a Loss Of Coolant Accident (LOCA), especially when high burn-up fuel in a severe state of fragmentation is involved. Here, this study presents a Computational Fluid Dynamics (CFD) model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) to predict the time-scale of plenum depressurization in Light-Water Reactor (LWR) fuel rods driven by axial gas transport through fragmented pellets. The model examines the effects of incorporating non-linearities in the friction term by comparing the results with experimental data. The experiment employed surrogate fuel rods containing pellets subjected to mechanical and/or thermal loadings to simulate various severity of cracking, and aimed at studying the influence of fuel conditions on axial gas redistribution. The results of this analysis indicate that under certain flow regime conditions - determined by the value of an equivalent Reynolds number - accounting for the non-linear friction term in Navier-Stokes equations guarantees better predictions for the time-scale of plenum depressurization. Also, the model enabled the simulation of the pressure decay by assigning distinct permeability values to each pellet instead of a single uniform value. Multiple simulations were run across all possible pellet position combinations, having each pellet assigned with values of permeability extracted from the experimental data. This allows to quantify the impact of the considering various non-uniform distributions of permeability on the dynamics of axial gas redistribution. The present work findings enhance the understanding of axial gas transport, and provide valuable insights for the integration of a model for predicting the axial gas redistribution during a LOCA scenario into the BISON fuel performance code.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon

Event Horizon Telescope (EHT) observations at 230 GHz have now imaged polarized emission around the supermassive black hole in M87 on event-horizon scales. This polarized synchrotron radiation probes the structure of magnetic fields and the plasma properties near the black hole. Here we compare the resolved polarization structure observed by the EHT, along with simultaneous unresolved observations with the Atacama Large Millimeter/submillimeter Array, to expectations from theoretical models. The low fractional linear polarization in the resolved image suggests that the polarization is scrambled on scales smaller than the EHT beam, which we attribute to Faraday rotation internal to the emission region. We estimate the average density n_e of order 10^4-7 cm-3, magnetic field strength B of order 1-30 G, and electron temperature Te of order (1-12) x 10^10 K of the radiating plasma in a simple one-zone emission model. We show that the net azimuthal linear polarization pattern may result from organized, poloidal magnetic fields in the emission region. In a quantitative comparison with a large library of simulated polarimetric images from general relativistic magnetohydrodynamic (GRMHD) simulations, we identify a subset of physical models that can explain critical features of the polarimetric EHT observations while producing a relativistic jet of sufficient power. The consistent GRMHD models are all of magnetically arrested accretion disks, where near-horizon magnetic fields are dynamically important. We use the models to infer a mass accretion rate onto the black hole in M87 of (3-20) x 10^-4 Msun yr-1.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic-ballooning-limited pedestals in spherical tokamak plasmas

Abstract A theoretical model is presented that for the first time matches experimental measurements of the pedestal width-height Diallo scaling in the low-aspect-ratio high- β tokamak NSTX. Combining linear gyrokinetics with self-consistent pedestal equilibrium variation, kinetic-ballooning, rather than ideal-ballooning plasma instability, is shown to limit achievable confinement in spherical tokamak pedestals. Simulations are used to find the novel Gyrokinetic Critical Pedestal constraint, which determines the steepest pressure profile a pedestal can sustain subject to gyrokinetic instability. Gyrokinetic width-height scaling expressions for NSTX pedestals with varying density and temperature profiles are obtained. These scalings for STs depart significantly from that of conventional aspect ratio tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Theoretical Basis for the Scaling Law of Broadband Shock Noise Intensity in Supersonic Jets

A theoretical basis for the scaling of broadband shock noise intensity In supersonic jets was formulated considering linear shock-shear wave interaction. Modeling of broadband shock noise with the aid of shock-turbulence interaction with special reference to linear theories is briefly reviewed. An hypothesis has been postulated that the peak angle of incidence (closer to the critical angle) for the shear wave primarily governs the generation of sound in the interaction process with the noise generation contribution from off-peak incident angles being relatively unimportant. The proposed hypothesis satisfactorily explains the well-known scaling law for the broadband shock-associated noise in supersonic jets.

Kandula, Max↗

Improved importance sampling technique for efficient simulation of digital communication systems

A new, improved importance sampling (IIS) approach to simulation is considered. Some basic concepts of IS are introduced, and detailed evolutions of simulation estimation variances for Monte Carlo (MC) and IS simulations are given. The general results obtained from these evolutions are applied to the specific previously known conventional importance sampling (CIS) technique and the new IIS technique. The derivation for a linear system with no signal random memory is considered in some detail. For the CIS technique, the optimum input scaling parameter is found, while for the IIS technique, the optimum translation parameter is found. The results are generalized to a linear system with memory and signals. Specific numerical and simulation results are given which show the advantages of CIS over MC and IIS over CIS for simulations of digital communications systems.

Lu, Dingqing↗

Atmospheric Parameter Climatologies from AIRS: Monitoring Short-, and Longer-Term Climate Variabilities and 'Trends'

The AIRS instrument is currently the best space-based tool to simultaneously monitor the vertical distribution of key climatically important atmospheric parameters as well as surface properties, and has provided high quality data for more than 5 years. AIRS analysis results produced at the GODDARD/DAAC, based on Versions 4 & 5 of the AIRS retrieval algorithm, are currently available for public use. Here, first we present an assessment of interrelationships of anomalies (proxies of climate variability based on 5 full years, since Sept. 2002) of various climate parameters at different spatial scales. We also present AIRS-retrievals-based global, regional and 1x1 degree grid-scale "trend"-analyses of important atmospheric parameters for this 5-year period. Note that here "trend" simply means the linear fit to the anomaly (relative the mean seasonal cycle) time series of various parameters at the above-mentioned spatial scales, and we present these to illustrate the usefulness of continuing AIRS-based climate observations. Preliminary validation efforts, in terms of intercomparisons of interannual variabilities with other available satellite data analysis results, will also be addressed. For example, we show that the outgoing longwave radiation (OLR) interannual spatial variabilities from the available state-of-the-art CERES measurements and from the AIRS computations are in remarkably good agreement. Version 6 of the AIRS retrieval scheme (currently under development) promises to further improve bias agreements for the absolute values by implementing a more accurate radiative transfer model for the OLR computations and by improving surface emissivity retrievals.

Molnar, Gyula↗