Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearized distribution”

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,333 records · Page 74

Unified, Geometric Framework for Nonequilibrium Protocol Optimization

Controlling thermodynamic cycles to minimize the dissipated heat is a long-standing goal in thermodynamics, and more recently, a central challenge in stochastic thermodynamics for nanoscale systems. Here, we introduce a theoretical and computational framework for optimizing nonequilibrium control protocols that can transform a system between two distributions in a minimally dissipative fashion. These protocols optimally transport a system along paths through the space of probability distributions that minimize the dissipative cost of a transformation. Furthermore, we show that the thermodynamic metric—determined via a linear response approach—can be directly derived from the same objective function that is optimized in the optimal transport problem, thus providing a unified perspective on thermodynamic geometries. As a result, we investigate this unified geometric framework in two model systems and observe that our procedure for optimizing control protocols is robust beyond linear response.

36 MATERIALS SCIENCE↗

Reinforcement Learning via Gaussian Processes with Neural Network Dual Kernels

While deep neural networks (DNNs) and Gaussian Processes (GPs) are both popularly utilized to solve problems in reinforcement learning, both approaches feature undesirable drawbacks for challenging problems. DNNs learn complex non-linear embeddings, but do not naturally quantify uncertainty and are often data-inefficient to train. GPs infer posterior distributions over functions, but popular kernels exhibit limited expressivity on complex and high-dimensional data. Fortunately, recently discovered conjugate and neural tangent kernel functions encode the behavior of overparameterized neural networks in the kernel domain. We demonstrate that these kernels can be efficiently applied to regression and reinforcement learning problems by analyzing a baseline case study.We apply GPs with neural network dual kernels to solve reinforcement learning tasks for the first time. We demonstrate, using the well understood mountain-car problem, that GPs empowered with dual kernels perform at least as well as those using the conventional radial basis function kernel. Finally, we conjecture that by inheriting the probabilistic rigor of GPs and the powerful embedding properties of DNNs, GPs using NN dual kernels will empower future reinforcement learning models on difficult domains.

97 MATHEMATICS AND COMPUTING↗

Minimizing Timing Jitter’s Impact on Ground-Penetrating Radar Array Coupling Signals

This article presents a novel investigation into the effect of receiver timing jitter on the quality of imaging of ground-penetrating radar systems. We explicitly show the impact of timing jitter throughout an entire example processing pipeline. The process of coupling removal in systems whose receivers have a random Gaussian-distributed jitter with a high standard deviation leaves randomly distributed residue in the radar image. This residue, called jitter-induced coupling noise, extends even into the expected target domain which may cause false alarms. This residue approximately follows a sum of Gaussian and Gamma distributions based on the higher order derivatives of the coupling signal. Once migrated to create an image, the coupling residue is best described by a Weibull distribution. We show that simple filtering does not remove enough residue and propose two preventative design specifications called the linear slope and migrated probability specifications. Antenna and system designers should use these specifications to assess the true severity of coupling signals in a given jitter environment.

47 OTHER INSTRUMENTATION↗

Impact of beam–beam effects on absolute luminosity calibrations at the CERN Large Hadron Collider

At the Large Hadron Collider (LHC), absolute luminosity calibrations obtained by the van der Meer (vdM) method are affected by the mutual electromagnetic interaction of the two beams. The colliding bunches experience relative orbit shifts, as well as optical distortions akin to the dynamic-β effect, that both depend on the transverse beam separation and must therefore be corrected for when deriving the absolute luminosity scale. In the vdM regime, the beam–beam parameter is small enough that the orbit shift can be calculated analytically. The dynamic-β corrections to the luminometer calibrations, however, had until the end of Run 2 been estimated in the linear approximation only. In this report, the influence of beam–beam effects on the vdM-based luminosity scale is quantified, together with the associated systematic uncertainties, by means of simulations that fully take into account the non-linearity of the beam–beam force, as well as the resulting non-Gaussian distortions of the transverse beam distributions. Two independent multiparticle simulations, one limited to the weak-strong approximation and one that models strong-strong effects in a self-consistent manner, are found in excellent agreement; both predict a percent-level shift of the absolute pp-luminosity values with respect to those assumed until recently in the physics publications of the LHC experiments. These results also provide guidance regarding further studies aimed at reducing the beam–beam-related systematic uncertainty on beam–beam corrections to absolute luminosity calibrations by the van der Meer method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Review—Development of Highly Active and Durable Hybrid Compressive Platinum Lattice Catalysts for Polymer Electrolyte Membrane Fuel Cells: Mathematical Modeling and Experimental Work

This review provides a comprehensive overview on the development of highly active and durable platinum catalysts with ultra-low Pt loadings for polymer electrolyte membrane fuel cells (PEMFCs) through a combined mathematical modeling and experimental work. First, simulation techniques were applied to evaluate the validity of the Tafel approximation for the calculation of the mass activity (MA) and specific activity (SA). A one-dimensional agglomeration model was developed and solved to understand the effects of exchange current density, porosity, agglomerate size, Nafion ® film thickness, and Pt loading on the MA and SA. High porosity (> 60%) and agglomerations at high Pt loadings cause the loss of the Tafel approximation and consequently the decrease in MA and SA. A new structure parameter was introduced to estimate the real porous structure using the fractal theory. The volumetric catalyst density was corrected by the fractal dimension (measured by Hg porosimetry), which gave a good agreement with the experimental values. The loading-dependent Tafel equation was then derived, which contains both the utilization and the non-linear scaling factor. Second, activated carbon composite support (ACCS) with optimized surface area, porosity, pore size, and pore size distribution was developed. The hydrophilic/hydrophobic ratio, structural properties (amorphous/crystalline ratio), and the number of active sites were optimized through metal-catalyzed pyrolysis. Stability of ACCS and Pt/ACCS were evaluated using an accelerated stress test (AST). The results indicated that Pt/ACCS showed no significant loss of MA and power density after 5,000 cycles at 1.0–1.5 V, while the commercial Pt/C catalysts showed drastic losses of MA and power density. Finally, monolayers of compressed Pt (core–shell-type Pt 3 Co 1 ) catalysts were structured by diffusing Co atoms (previously embedded in ACCS) into Pt. Compressive Pt lattice (Pt * ) catalysts were synthesized through an annealing procedure developed at the University of South Carolina (USC). The Pt * /ACCS catalyst showed high initial power density (rated) of 0.174 g Pt kW −1 and high stability (24 mV loss) at 0.8 A cm −2 after 30,000 cycles (0.6–1.0 V). The outstanding performance of Pt * /ACCS is due to the synergistic effect of ACCS and compressive Pt * lattice.

25 ENERGY STORAGE↗

MoorPy (Quasi-Static Mooring Analysis in Python)

MoorPy is a quasi-static mooring model and a suite of associated functions for mooring system analysis. The core model supports quasi-static analysis of moored floating systems including any arrangement of mooring lines and floating platforms. It solves the distributed position and tension of each mooring line segment using standard catenary equations. Floating platforms can be represented with linear hydrostatic characteristics. MoorPy automatically computes a floating system's equilibrium state and can be queried to identify a mooring system's nonlinear force-displacement relationships. Linearized stiffness matrices are efficiently computed using semi-analytic Jacobians. MoorPy also includes plotting functions and a library of mooring component property and cost coefficients. MoorPy can be used directly from Python scripts to perform mooring design and analysis tasks, or it can be coupled with other tools to compute quasi-static mooring reactions as part of a larger simulation.

Hall, Mathew↗

Global weak solutions for a nonlocal multispecies Fokker–Planck–Landau system

The global-in-time existence of weak solutions to a spatially homogeneous multispecies Fokker–Planck–Landau system for plasmas in the three-dimensional whole space is shown. The Fokker–Planck–Landau system is a simplification of the Landau equations assuming a linearized, velocity-independent, and isotropic kernel. The resulting equations depend nonlocally and nonlinearly on the moments of the distribution functions via the multispecies local Maxwellians. Furthermore, the existence proof is based on a three-level approximation scheme, energy and entropy estimates, as well as compactness results, and it holds for both soft and hard potentials.

97 MATHEMATICS AND COMPUTING↗

Thermal elastoplastic structural analysis of non-metallic thermal protection systems

An incremental theory and numerical procedure to analyze a three-dimensional thermoelastoplastic structure subjected to high temperature, surface heat flux, and volume heat supply as well as mechanical loadings are presented. Heat conduction equations and equilibrium equations are derived by assuming a specific form of incremental free energy, entropy, stresses and heat flux together with the first and second laws of thermodynamics, von Mises yield criteria and Prandtl-Reuss flow rule. The finite element discretization using the linear isotropic three-dimensional element for the space domain and a difference operator corresponding to a linear variation of temperature within a small time increment for the time domain lead to systematic solutions of temperature distribution and displacement and stress fields. Various boundary conditions such as insulated surfaces and convection through uninsulated surface can be easily treated. To demonstrate effectiveness of the present formulation a number of example problems are presented.

Chung, T. J.↗

Advanced control concepts

The problems of excess control devices and insufficient trim control capability on shuttle ascent vehicles were investigated. The trim problem is solved at all time points of interest using Lagrangian multipliers and a Simplex based iterative algorithm developed as a result of the study. This algorithm has the capability to solve any bounded linear problem with physically realizable constraints, and to minimize any piecewise differentiable cost function. Both solution methods also automatically distribute the command torques to the control devices. It is shown that trim requirements are unrealizable if only the orbiter engines and the aerodynamic surfaces are used.

Sharp, J. B.↗

Energetic particles and ionization in the nighttime middle and low latitude ionosphere

Seven Nike Apache rockets, each equipped with an energetic particle spectrometer (12 E 80 keV) and electron-density experiments, were launched from Wallops Island, Virginia and Chilca, Peru, under varying geomagnetic conditions near midnight. At Wallops Island the energetic particle flux (E 40 keV) is found to be strongly dependent on Kp. The pitch-angle distribution is asymmetrical about a peak at 90 D signifying a predominately quasi-trapped flux and explaining the linear increase of count rate with altitute in the altitude region 120 to 200 km. The height-averaged ionization rates derived from the electron-density profiles are consistent with the rates calculated from the observed total particle flux for magnetic index Kp 3. In the region 90 to 110 km it is found that the nighttime ionization is primarily a result of Ly-beta radiation from the geocorona and interplanetary hydrogen for even very disturbed conditions. Below 90 km during rather disturbed conditions energetic electrons can be a significant ionization source. Two energetic particle precipitation zones have been identified at midlatitudes.

Voss, H. D.↗

The origin of the grooves on Phobos

Various theories for the long, linear depressions on the surface of Phobos are reviewed. Imagery from Viking Orbiters is used to map the surface distribution of the grooves, study their morphology, and date them by means of the density of superimposed impact craters. Data is presented which tends to support the hypothesis that the deep-seated fracturing was caused by a large, nearly catastrophic cratering event. It is suggested that the grooves were produced during the creation of the Stickney crater, rather than as the result of tidal stresses induced by Mars or by drag forces during the hypothetical capture of the satellite by Mars.

Thomas, P. C.↗

Subsonic 3-D surface panel method for rapid analysis of multiple geometry perturbations

For any baseline aircraft configuration, a matrix of partial derivatives of surface velocity potential with respect to geometry coordinates is calculated. Linear extrapolation is then used to analyze the subcritical potential flow corresponding to a series of arbitrary small geometry perturbations. Each perturbation analysis is more than an order of magnitude more efficient than a conventional panel method solution because no influence coefficients are calculated and no large system of linear algebraic equations is solved. Wing and wing-fuselage examples are presented to demonstrate that the predicted pressure distributions are nearly exact for large changes to wing camber, thickness, and leading edge radius.

Bristow, D. R.↗

Tortuosity of lightning return stroke channels

Data obtained from photographs of lightning are presented on the tortuosity of return stroke channels. The data were obtained by making piecewise linear fits to the channels, and recording the cartesian coordinates of the ends of each linear segment. The mean change between ends of the segments was nearly zero in the horizontal direction and was about eight meters in the vertical direction. Histograms of these changes are presented. These data were used to create model lightning channels and to predict the electric fields radiated during return strokes. This was done using a computer generated random walk in which linear segments were placed end-to-end to form a piecewise linear representation of the channel. The computer selected random numbers for the ends of the segments assuming a normal distribution with the measured statistics. Once the channels were simulated, the electric fields radiated during a return stroke were predicted using a transmission line model on each segment. It was found that realistic channels are obtained with this procedure, but only if the model includes two scales of tortuosity: fine scale irregularities corresponding to the local channel tortuosity which are superimposed on large scale horizontal drifts. The two scales of tortuosity are also necessary to obtain agreement between the electric fields computed mathematically from the simulated channels and the electric fields radiated from real return strokes. Without large scale drifts, the computed electric fields do not have the undulations characteristics of the data.

Levine, D. M.↗

VLBI observations of the RS Canum Venaticorum binary systems UX Arietis and HR 1099 at 1.65 GHz

VLBI observations of the RS CVn binaries UX Arietis and HR 1099 have been made at 1.65 GHz using a three-element array with a minimum fringe spacing of 11.5 milli-arcsec. Both sources were found to be unresolved within measurement uncertainties. In both cases, the derived upper limit to the source size was comparable to the overall size of each binary system. The lower limits to the brightness temperature were 1.4 x 10 to the 10th K for UX Arietis and 2.9 x 10 to the 10th K for HR 1099. Simultaneous polarization measurements at the VLA showed 4-8 percent circular polarization and less than 2 percent linear polarization. It is found that the data are consistent with gyrosynchrotron emission from a power-law energy distribution of electrons in a magnetic field B less than or approximately equal to 6 gauss.

Mutel, R. L.↗

Two-dimensional wind-tunnel interference from measurements on two contours

This paper describes how wall-induced velocities near a model in a two-dimensional wind tunnel can be estimated from upwash distributions measured along two contours surrounding a model. The method is applicable to flows that can be represented by linear theory. It was derived by applying the Schwarz Integral Formula separately to the two contours and by exploiting the free-air relationship between upwashes along the contours. Advantages of the method are that only one flow quantity need by measured and no representation of the model is required. A weakness of the method is that it assumes streamwise interference velocity vanishes far upstream of the model. This method was applied to a simple theoretical model of flow in a solid-wall wind tunnel. The theoretical interference velocities and the velocities computed using the method were in excellent agreement. The method was then used to analyze experimental data acquired during adaptive-wall experiments at Ames Research Center. This analysis confirmed that the wall adjustments reduced wall-induced velocities near the model.

Schairer, E. T.↗

Broadening of the Interplanetary Helium Cone Structure Due to Elastic Collisions of LISM Helium Atoms with Solar Wind Ions

Neutral interstellar particles penetrating into the heliosphere, besides being subject there to specific loss processes, suffer elastic collisions with KeV-solar wind ions. The momentum transfer to the neutrals connected with these collisions leads to a loss of angular momentum with respect to the Sun and to a fractional compensation of the effective solar gravity. The dynamical particle trajectories hence are changed into non-Keplerians leading to density and temperature distributions differing from those calculated in the past. This is found from a solution of the Boltzmann equation that linearizes the effect of this additional force. It is shown that the HeI-584A resonance glow of the heliospheric helium cone lead to substantially lower interstellar helium temperatures if re-interpreted on the basis of this revised theory. These temperatures now seem to be in accordance with the derived temperatures for interstellar hydrogen.

Fahr, H. J.↗

Morphological stability in the presence of fluid flow in the melt

Recent experiments have shown that the presence of a vertical buoyancy-driven flow adjacent to an initially cylindrical crystal-melt interface may produce a time-dependent helical deformation of the interface, with a rotation period ranging from several minutes to many hours, depending upon the width of the melt. The temperature distribution is such that the interface is expected to be morphologically stable in the absence of fluid flow. A linear stability analysis reveals that the instability is due to a coupling between a basic hydrodynamic instability in the buoyant flow and the deformable boundary separating the two phases. The crystal-melt interface lowers the critical Grashof number of an analogous rigid-walled system by an order of magnitude for succinonitrile with a Prandtl number P = 22.8; furthermore, the hydrodynamic mode that is actually destabilized by the interface is not the least stable mode in the rigid-walled system for P = 22.8. The results show that the instability may be regarded either as a rather large alteration of a basic hydrodynamic instability by the crystal-melt interface, or as a significant modification of the morphological stability of the interface by the presence of the buoyant flow.

Mcfadden, G. B.↗

Unified control/structure design and modeling research

To demonstrate the applicability of the control theory for distributed systems to large flexible space structures, research was focused on a model of a space antenna which consists of a rigid hub, flexible ribs, and a mesh reflecting surface. The space antenna model used is discussed along with the finite element approximation of the distributed model. The basic control problem is to design an optimal or near-optimal compensator to suppress the linear vibrations and rigid-body displacements of the structure. The application of an infinite dimensional Linear Quadratic Gaussian (LQG) control theory to flexible structure is discussed. Two basic approaches for robustness enhancement were investigated: loop transfer recovery and sensitivity optimization. A third approach synthesized from elements of these two basic approaches is currently under development. The control driven finite element approximation of flexible structures is discussed. Three sets of finite element basic vectors for computing functional control gains are compared. The possibility of constructing a finite element scheme to approximate the infinite dimensional Hamiltonian system directly, instead of indirectly is discussed.

Mingori, D. L.↗