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 541 records · Page 30

Linear Actuator Has Long Stroke and High Resolution

The term precision linear actuator, direct drive ( PLADD ) refers to a robust linear actuator designed to be capable of repeatedly performing, over a lifetime of the order of 5 to 10 years, positioning maneuvers that include, variously, submicron increments or slews of the order of a centimeter. The PLADD is capable of both long stroke (120 mm) and high resolution (repeatable increments of 20 nm). Unlike precise linear actuators of prior design, the PLADD contains no gears, levers, or hydraulic converters. The PLADD, now at the prototype stage of development, is intended for original use as a coarse-positioning actuator in a spaceborne interferometer. The PLADD could also be adapted to terrestrial applications in which there are requirements for long stroke and high resolution: potential applications include medical imaging and fabrication of semiconductor devices. The PLADD (see figure) includes a commercially available ball-screw actuator driven directly by a commercially available three-phase brushless DC motor. The ball-screw actuator comprises a spring-preloaded ball nut on a ball screw that is restrained against rotation as described below. The motor is coupled directly (that is, without an intervening gear train) to a drive link that, in turn, is coupled to the ball nut. By eliminating the gear train, the direct-drive design eliminates the complexity, backlash, and potential for misalignment associated with a gear train. To prevent inadvertent movement, there is a brake that includes flexured levers compressed against the drive link by preload springs. This is a power-off brake: There are also piezoelectric stacks that can be activated to oppose the springs and push the levers away from the drive link. Hence, power must be applied to the piezoelectric stacks to release the drive link from braking. To help ensure long operational life, all of the mechanical drive components are immersed in an oil bath within hermetically sealed bellows. The outer end of the bellows holds the outer end of the ball screw, thereby preventing rotation of the ball screw. Positioning is controlled by an electronic control system that includes digital and analog subsystems that interact with the motor and brake and with two sensor/encoder units: a Hall-effect-sensor rotation encoder and a linear glass-scale encoder. This system implements a proportional + integral + derivative control algorithm that results in variation of voltage commands to each of the three pairs of windings of the brushless DC motor. In one of two alternative control modes, the voltages are applied to the windings in a trapezoidal commutation scheme on the basis of timing signals obtained from the Hall-effect sensors; this scheme yields relatively coarse positioning - 24 steps per motor revolution. The second control mode involves a sinusoidal commutation scheme in which the output of the linear glass-scale encoder is transposed to rotational increments to yield much finer position feedback - more than 400,000 steps per revolution.

Cook, Brant T.↗

Inclusion of Linearized Moist Physics in Nasa's Goddard Earth Observing System Data Assimilation Tools

Inclusion of moist physics in the linearized version of a weather forecast model is beneficial in terms of variational data assimilation. Further, it improves the capability of important tools, such as adjoint-based observation impacts and sensitivity studies. A linearized version of the relaxed Arakawa-Schubert (RAS) convection scheme has been developed and tested in NASA's Goddard Earth Observing System data assimilation tools. A previous study of the RAS scheme showed it to exhibit reasonable linearity and stability. This motivates the development of a linearization of a near-exact version of the RAS scheme. Linearized large-scale condensation is included through simple conversion of supersaturation into precipitation. The linearization of moist physics is validated against the full nonlinear model for 6- and 24-h intervals, relevant to variational data assimilation and observation impacts, respectively. For a small number of profiles, sudden large growth in the perturbation trajectory is encountered. Efficient filtering of these profiles is achieved by diagnosis of steep gradients in a reduced version of the operator of the tangent linear model. With filtering turned on, the inclusion of linearized moist physics increases the correlation between the nonlinear perturbation trajectory and the linear approximation of the perturbation trajectory. A month-long observation impact experiment is performed and the effect of including moist physics on the impacts is discussed. Impacts from moist-sensitive instruments and channels are increased. The effect of including moist physics is examined for adjoint sensitivity studies. A case study examining an intensifying Northern Hemisphere Atlantic storm is presented. The results show a significant sensitivity with respect to moisture.

Convection↗

PURE: Scalable Phase Unwrapping with Spatial Redundant Arcs

Phase unwrapping is a key problem in many coherent imaging systems, such as syntheticapertureradar(SAR)interferometry. Ageneralformulationforredundant integration of finite differences for phase unwrapping (Costantini et al., 2010) was shown to produce a more reliable solution by exploiting redundant differential estimates. However, this technique requires a commercial linear programming solver for large-scale problems. For a linear cost function, we propose a method based on Dual Decomposition that breaks the given problem defined over a nonplanar graph into tractable sub-problems over planar subgraphs. We also propose a decomposition technique that exploits the underlying graph structure for solving thesub-problemsefficientlyandguaranteesasymptoticconvergencetotheglobally optimal solution. The experimental results demonstrate that the proposed approach is comparable to the existing state-of-the-art methods in terms of the estimate with a better runtime and memory footprint.

Lanka, Ravi↗

Large-scale harmonic balance simulations with Krylov subspace and preconditioner recycling

The multi-harmonic balance method combined with numerical continuation provides an efficient framework to compute a family of time-periodic solutions, or response curves, for large-scale, nonlinear mechanical systems. The predictor and corrector steps repeatedly solve a sequence of linear systems that scale by the model size and number of harmonics in the assumed Fourier series approximation. In this paper, a novel Newton–Krylov iterative method is embedded within the multi-harmonic balance and continuation algorithm to efficiently compute the approximate solutions from the sequence of linear systems that arise during the prediction and correction steps. Further, the method recycles, or reuses, both the preconditioner and the Krylov subspace generated by previous linear systems in the solution sequence. A delayed frequency preconditioner refactorizes the preconditioner only when the performance of the iterative solver deteriorates. The GCRO-DR iterative solver recycles a subset of harmonic Ritz vectors to initialize the solution subspace for the next linear system in the sequence. The performance of the iterative solver is demonstrated on two exemplars with contact-type nonlinearities and benchmarked against a direct solver with traditional Newton–Raphson iterations.

97 MATHEMATICS AND COMPUTING↗

Linear regulator design for stochastic systems by a multiple time scales method

A hierarchically-structured, suboptimal controller for a linear stochastic system composed of fast and slow subsystems is considered. The controller is optimal in the limit as the separation of time scales of the subsystems becomes infinite. The methodology is illustrated by design of a controller to suppress the phugoid and short period modes of the longitudinal dynamics of the F-8 aircraft.

Teneketzis, D.↗

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗

Time–Concentration Superposition for Linear Viscoelasticity of Polymer Solutions

The concentration dependence of linear viscoelastic properties of polymer solutions is a well-studied topic in polymer physics. Dynamic scaling theories allow qualitative predictions of polymer solution rheology, but quantitative predictions are still limited to model polymers. Meanwhile, the scaling properties of non-model polymer solutions must be determined experimentally. In present paper, the time–concentration superposition (TCS) of experimental data is shown to be a robust procedure for studying the concentration scaling properties of binary and ternary polymer solutions. TCS can not only identify whether power law scaling may exist or not, and over which concentration range, but also unambiguously estimate the concentration scaling exponents of linear viscoelastic properties for a range of non-model polymer solutions.

42 ENGINEERING↗

A Lagrangian perturbation theory in the presence of massive neutrinos

We develop a Lagrangian Perturbation Theory (LPT) framework to study the clustering of cold dark matter (CDM) in cosmologies with massive neutrinos. We follow the trajectories of CDM particles with Lagrangian displacements fields up to third order in perturbation theory. Once the neutrinos become non-relativistic, their density fluctuations are modeled as being proportional to the CDM density fluctuations, with a scale-dependent proportionality factor. This yields a gravitational back-reaction that introduces additional scales to the linear growth function, which is accounted for in the higher order LPT kernels. Through non-linear mappings from Eulerian to Lagrangian frames, we ensure that our theory has a well behaved large scale behavior free of unwanted UV divergences, which are common when neutrino and CDM densities are not treated on an equal footing, and in resummation schemes that manifestly break Galilean invariance. We use our theory to construct correlation functions for both the underlying matter field, as well as for biased tracers using ConvolutionLPT. Redshift-space distortions effects are modeled using the Gaussian Streaming Model. When comparing our analytical results to simulated data from the Quijote 1 simulation suite, we find good accuracy down to r = 20 Mpc h -1 at redshift z = 0.5, for the real space and redshift space monopole particle correlation functions with no free parameters. The same accuracy is reached for the redshift space quadrupole if we additionally consider an effective field theory parameter that shifts the pairwise velocity dispersion. For modeling the correlation functions of tracers we adopt a simple Lagrangian biasing scheme with only density and curvature operators, which we find sufficient to reach down to r = 20 Mpc h -1 when comparing to simulated halos.

79 ASTRONOMY AND ASTROPHYSICS↗

Towards Energy Scale Calibration and Drift Correction of TES Detectors for Athena X-IFU

The Athena X-Ray Integral Field Unit (X-IFU) comprises a 2376-pixel array of transition edge sensors (TES) read out with time-division multiplexing (TDM). X-IFU will provide spatially resolved, high-resolution spectroscopy (2.5 eV full-width-half-maximum up to 7 keV) over the energy range 0.2 to 12 keV, with an absolute energy scale accuracy of 0.4 eV. The energy scale function maps the optimally filtered pulse height, in arbitrary engineering units, to real calibrated energy. Uncertainties in the calibration can result from imperfect fitting of the energy scale between the known calibration points. Furthermore, temporal changes in the TES operating environment, such as heat-sink temperature, magnetic field and bias voltage, can cause significant variations in the detector gain function over time. If not properly corrected, this can result in degradation of the energy resolution, and systematic errors in the absolute energy scale. The non-linear nature of TES detectors, coupled with the possibility of multiple simultaneously occurring sources of drift, can make effective corrections over the full bandpass of the instrument extremely challenging. Athena X-IFU will employ an on-board calibration source that provides known reference x-ray lines. This provides real-time monitoring of the gain stability of the detector system and information that can be used to correct for gain drifts. For X-IFU the baseline approach is to measure a series of calibration curves under different environmental conditions, which bound the expected drifts the instrument is predicted to see over the course of the mission. Using the information from the in-flight calibration source, these energy scale functions can be interpolated to generate a new corrected energy scale as a function of time. In this paper we discuss progress towards demonstrating that the X-IFU energy scale requirements can be met. We present measurements on ~ 200 pixels in a prototype X-IFU array read out with 8-column x 32-row TDM. We use a rotating target source containing 12 fluorescent targets to generate x-ray lines covering the energy range 4 keV (Sc-Kα) to 12 keV (Br-Kα). We present measurements of the non-linear energy scale function and show how variations in heat-sink temperature, TES bias voltage and magnetic field affect the shape of TES energy scale differently and introduce different residual gain errors over the bandpass. We explore different drift correction algorithms that use either a single or multiple referential lines to track and correct the gain from these various sources of drift. In addition to the pulse-height, the DC ‘baseline’ level of the TES can contain information about its bias conditions. Thus, we test a multi-parameter gain correction algorithm that attempts to incorporate both the pulse height and the additional baseline information into the algorithm.

Stephen J Smith↗

Complexity-theoretic foundations of BosonSampling with a linear number of modes

BosonSampling is the leading candidate for demonstrating quantum computational advantage in photonic systems. While we have recently seen many impressive experimental demonstrations, there is still a formidable distance between the complexity-theoretic hardness arguments and current experiments. One of the largest gaps involves the ratio of {particles} to modes -- all current hardness evidence assumes a dilute regime in which the number of linear optical modes scales at least quadratically in the number of particles. By contrast, current experiments operate in a saturated regime with a linear number of modes. In this paper we bridge this gap, bringing the hardness evidence for experiments in the saturated regime to the same level as had been previously established for the dilute regime. This involves proving a new worst-to-average-case reduction for computing the Permanent which is robust to both large numbers of row repetitions and also to distributions over matrices with correlated entries. We also apply similar arguments to give evidence for hardness of Gaussian BosonSampling in the saturated regime.

Computational Complexity (cs.CC)↗

Computation of Large-Scale Structure Jet Noise Sources With Weak Nonlinear Effects Using Linear Euler

An approximate technique is presented for the prediction of the large-scale turbulent structure sound source in a supersonic jet. A linearized Euler equations code is used to solve for the flow disturbances within and near a jet with a given mean flow. Assuming a normal mode composition for the wave-like disturbances, the linear radial profiles are used in an integration of the Navier-Stokes equations. This results in a set of ordinary differential equations representing the weakly nonlinear self-interactions of the modes along with their interaction with the mean flow. Solutions are then used to correct the amplitude of the disturbances that represent the source of large-scale turbulent structure sound in the jet.

Dahl, Milo D.↗

Flame zone of a composite propellant expanded by a laser source

Technique scales flame structure linearly with gas kinetic mean free path, which increases two to three orders of magnitude as pressure decreases like amount. Kinetic and transport time scales expand in proportion so that regression rates for laser-induced flames are two to three orders of magnitude slower.

Friedman, R.↗

Linear regulator design for stochastic systems by a multiple time scales method

This paper develops a hierarchically-structured, suboptimal controller for a linear stochastic system composed of fast and slow subsystems. It is proved that the controller is optimal in the limit as the separation of time scales of the subsystems becomes infinite. The methodology is illustrated by design of a controller to suppress the phugoid and short period modes of the longitudinal dynamics of the F-8 aircraft.

Teneketzis, D.↗

Modeling Temperature Profiles in the Pedestal of NSTX with Reduced Models

This paper describes new modeling capabilities for predicting H-mode pedestal profiles in spherical tokamaks. Temperature profiles for NSTX discharges 132543 and 132588 are modeled by coupling the \textsc{astra} transport solver with neoclassical transport and gyrokinetic-based reduced models for electron temperature gradient (ETG) and kinetic ballooning mode (KBM) instabilities. A quasi-linear surrogate model for ion-scale transport is developed using linear \textsc{gene} simulations, requiring only a single free parameter calibrated to one discharge. Time-evolving the temperatures with fixed density yields good agreement with experiments for both discharges. Systematic analysis of the transport mechanisms reveals that neoclassical transport is huge across the entire pedestal region for the ion channel. ETG turbulence is large in the plasma edge and low density gradient region, contributing substantially to the electron channel. However, KBM/MHD-like modes also drive significant transport in both the ion and electron thermal channels, making them essential for accurate pedestal modeling. Further refinements, including explicit E×B shear suppression and scaled ETG transport, produce quantitative but not qualitative improvements. This work lays the foundation for predictive modeling of future devices. This paper is on arxiv and has been submitted to Nuclear Fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improved perfluoroalkylether fluid development

The objective of this program was to optimize and scale up the linear perfluoroalkylether stabilization process and to provide test data regarding the fluids' thermal oxidative stability in the presence of metal alloys. The stabilization of Fomblin Z-25 was scaled up to 300 g of fluid. The modified fluid was stable at 316 C in oxygen in the presence of M-50 alloy for more than 24 hrs but less than 40 hrs; the amount of volatiles produced after 24 hrs was 5.5 mg/g. In the presence of Ti(4Al,4Mn) alloy, under the above conditions, following an exposure of 24 hrs, the amount of volatiles formed was 6.2 mg/g; 56 hrs exposure yielded 13.9 mg/g. The commercial fluid at 288 C (in oxygen) in the presence of M-50 after 15 hrs of exposure decomposed extensively, 342 mg/g; in the presence of Ti(4Al,4Mn) alloy after only 8 hrs at 288 C, the amount of volatiles was 191 mg/g. Formulation of the commercial fluid with C2PN3 additive was not as effective as the stabilization processing. All the perfluoroalkylether fluids studied were stable in nitrogen at 343 C. The thermal oxidative stability in the absence of metal alloys varied, with Aflunox exhibiting the best behavior. All the fluids were degraded in oxygen at 316 C during 24 hrs exposure to Ti(4Al,4Mn) alloy with the exception of a perfluoroalkylether substituted triazine and the modified Z-25.

Paciorek, K. L.↗

The HyMeX Special Observation Period in Central Italy: Precipitation Measurements, Retrieval Techniques and Preliminary Results

The Mediterranean area concentrates the major natural risks related to the water cycle, including heavy precipitation and flash-flooding during the fall season. The capability to predict such high-impact events remains weak because of the contribution of very fine-scale processes and their non-linear interactions with the larger scale processes. These societal and science issues motivate the HyMeX (Hydrological cycle in the Mediterranean Experiment, http://www.hymex.orgl) experimental programme. HyMeX aims at a better quantification and understanding of the water cycle in the Mediterranean with emphasis on intense events. The observation strategy of HyMEX is organized in a long-term (4 years) Enhanced Observation Periods (EOP) and short-term (2 months) Special Observation Periods (SOP). HyMEX has identified 3 main Mediterranean target areas: North-West (NW), Adriatic (A) and South-East (SE). Within each target area several hydrometeorological sites for heavy rainfall and flash flooding have been set up. The hydrometeorological sire in Central Italy (CI) is interested by both western and eastern fronts coming from the Atlantic Ocean and Siberia, respectively. Orographic precipitations play an important role due to the central Apennine range, which reaches nearly 3000 m (Gran Sasso peak). Moreover, convective systems commonly develop in CI during late summer and beginning of autumn, often causing localized hailstorms with cluster organized cells. Western fronts may heavily hit the Tiber basin crossing large urban areas (Rome), whereas eastern fronts can cause flash floods along the Adriatic coastline. Two major basins are involved within Cl region: Tiber basin (1000 km long) and its tributary Aniene and the Aterno-Pescara basin (300 km long). The first HyMeX SOP1.1 was carried out from Sept. till Nov. 2012 in the NW target area The Italian SOP1.1 was coordinated by the Centre of Excellence CETEMPS, University of L'Aquila, a city located in the CI heart. The CI area was covered by a uniquely dense meteorological instrumentation thanks to a synergy between Italian institutions and NASA-GSFC. The following RADARs were operated: a Doppler single-polarization C-band radar located at Mt Midia; the Polar 55C Doppler dual-polarization C-band radar located in Rome; a Doppler C-hand polarimetric radar located at Il Monte (Abnazo); a polarimetric X-band mini-radar in L' Aquila; a polarimetric X-hand portable mini-radar in Rome; a single-polarization X-band mini-radar in Rome. DISDROMETERs were also deployed: 4 Parsivel optical disdrometers in Rome (at Sapienza, CNR-ISAC and CNR-INSEAN); 1 2D-video disdrometer in Rome; 3 Parsivels optical disdrometer respectively in L'Aquila (Abnazo), Avezzano (Abruzzo) and Pescara (Abnazo). Other INSTRUMENTS were available: 1 K-band vertically-pointing micro rain-radar (MRR), 2 Pludix X-band disdrometers, 1 VLF lightning sensor, 1 microwave radiometer at 23-31 GHz in Rome (at Sapienza); the raingauge network with more than 200 stations in Central Italy. Three overpasses in CI were also performed by the Falcon 20 aircraft equipped with the 950Hz cloud radar RASTA Analysis of the SOP1.1 main events in CI will be described by focusing on the raindrop size distribution statistics and its geographical variability. Intercomparison of rainfall estimates from disdrometers, raingauges and radars will be illustrated with the aim to provide a quality-controlled and physically consistent rainfall dataset for meteorological modeling validation and assimilation purposes.

Gatlin, Patrick↗