Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear projection”

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 145 records · Page 8

Potential Non-Linearities in the High Latitude Circulation and Ozone Response to Stratospheric Aerosol Injection

The impacts of Stratospheric Aerosol Injection (SAI) on the atmosphere and surface climate depend on when and where the sulfate aerosol precursors are injected, as well as on how much surface cooling is to be achieved. We use a set of CESM2(WACCM6) SAI simulations achieving three different levels of global mean surface cooling and demonstrate that unlike some direct surface climate impacts driven by the reflection of solar radiation by sulfate aerosols, the SAI-induced changes in the high latitude circulation and ozone are more complex and could be non-linear. This manifests in our simulations by disproportionally larger Antarctic springtime ozone loss, significantly larger intra-ensemble spread of the Arctic stratospheric jet and ozone responses, and non-linear impacts on the extratropical modes of surface climate variability under the strongest-cooling SAI scenario compared to the weakest one. These potential non-linearities may add to uncertainties in projections of regional surface impacts under SAI.

54 ENVIRONMENTAL SCIENCES↗

Optimal Energy Dispatch of Distributed PVs for the Next Generation of Distribution Management Systems

Advanced Distribution Management Systems (ADMS) are being widely adopted by electric utilities for managing and optimizing the operations of their distribution systems. Distributed photovoltaic (DPV) systems with smart inverters can be controlled to adjust active power and reactive power outputs, and they are envisioned to become a part of (centrally or distributed) controllable assets managed by the ADMS for optimizing grid operations. This paper proposes an optimal energy dispatch strategy controlling DPV systems for regulating distribution voltages and achieving conservation voltage reduction. A convex optimization model is proposed with the use of linearized power flow, and the gradient projection algorithm is used to solve the optimal active power and reactive power outputs of smart inverters. The proposed optimal energy dispatch is implemented using an open-source ADMS platform, and simulation results have demonstrated the effectiveness of the proposed approach on improving distribution grid operations.

42 ENGINEERING↗

Preliminary Design of a LAMP DTL

A preliminary Drift Tube Linac (DTL) layout was designed to create an algorithm for developing the conceptual design of the proposed LANSCE Accelerator Modernization Project (LAMP) final section: the proton linear accelerator from 3 MeV to 100 MeV. Preceding reports describe the proposed layouts of the LEBT, the RFQ, and the MEBT subsections of the linac. Initial estimates of the needed RF power for the DTL are also presented in the report. Present report describes the initial layout of the DTL and the first longitudinal beam dynamics results of the simulations n simplified models. The follow-up reports will include the transverse focusing scheme and beam dynamics in details. Present report in the last section will only describe the planned focusing scheme, based on the existing LANCE DTL.

43 PARTICLE ACCELERATORS↗

Preliminary Design of a LAMP DTL (Update)

Present report is a continuation of the previously published report LA-UR-23-26600 from 19-Jun-2023 and incorporates the mentioned report as its part. All the additional information is contained in the part III of this report, the DTL Transverse Focusing Scheme. A preliminary Drift Tube Linac (DTL) layout was designed to create an algorithm for developing the conceptual design of the proposed LANSCE Accelerator Modernization Project (LAMP) final section: the proton linear accelerator from 3 MeV to 100 MeV. Preceding reports describe the proposed layouts of the LEBT, the RFQ, and the MEBT subsections of the linac. Initial estimates of the needed RF power for the DTL are also presented in the report. Present report describes the initial layout of the DTL and the first longitudinal beam dynamics (BD) results of the simulations in simplified models. The part III of this report includes the transverse focusing scheme and transverse BD in details.

43 PARTICLE ACCELERATORS↗

On long-term periodicities in the sunspot record

Sunspot records are systematically maintained, with the knowledge that an 11 year average period exists since about 1850. Thus, the sunspot record of highest quality and considered to be the most reliable is that of cycle eight through the present. On the basis of cycles 8 through 20, various combinations of sine curves were used to approximate the observed R sub MAX values (where R sub MAX is the smoothed sunspot number at cycle maximum). It is found that a three component sinusoidal function, having an 11 cycle and a 2 cycle variation on a 90 cycle periodicity, yields computed R sub MAX values which fit, reasonably well, observed R sub MAX values for the modern sunspot cycles. Extrapolation of the empirical functions forward in time allows for the projection of values of R sub MAX for cycles 21 and 22. For cycle 21, the function projects a value of 157.3, very close to the actually observed value of 164.5. For cycle 22, the function projects a value of about 107. Linear regressions applied to cycle 22 indicate a long-period cycle (cycle duration 132 months). An extensive bibliography on techniques used to estimate the time dependent behavior of sunspot cycles is provided.

Wilson, R. M.↗

Tensor-GMRES method for large sparse systems of nonlinear equations

This paper introduces a tensor-Krylov method, the tensor-GMRES method, for large sparse systems of nonlinear equations. This method is a coupling of tensor model formation and solution techniques for nonlinear equations with Krylov subspace projection techniques for unsymmetric systems of linear equations. Traditional tensor methods for nonlinear equations are based on a quadratic model of the nonlinear function, a standard linear model augmented by a simple second order term. These methods are shown to be significantly more efficient than standard methods both on nonsingular problems and on problems where the Jacobian matrix at the solution is singular. A major disadvantage of the traditional tensor methods is that the solution of the tensor model requires the factorization of the Jacobian matrix, which may not be suitable for problems where the Jacobian matrix is large and has a 'bad' sparsity structure for an efficient factorization. We overcome this difficulty by forming and solving the tensor model using an extension of a Newton-GMRES scheme. Like traditional tensor methods, we show that the new tensor method has significant computational advantages over the analogous Newton counterpart. Consistent with Krylov subspace based methods, the new tensor method does not depend on the factorization of the Jacobian matrix. As a matter of fact, the Jacobian matrix is never needed explicitly.

Feng, Dan↗

Wavelet Approximation in Data Assimilation

Estimation of the state of the atmosphere with the Kalman filter remains a distant goal because of high computational cost of evolving the error covariance for both linear and nonlinear systems. Wavelet approximation is presented here as a possible solution that efficiently compresses both global and local covariance information. We demonstrate the compression characteristics on the the error correlation field from a global two-dimensional chemical constituent assimilation, and implement an adaptive wavelet approximation scheme on the assimilation of the one-dimensional Burger's equation. In the former problem, we show that 99%, of the error correlation can be represented by just 3% of the wavelet coefficients, with good representation of localized features. In the Burger's equation assimilation, the discrete linearized equations (tangent linear model) and analysis covariance are projected onto a wavelet basis and truncated to just 6%, of the coefficients. A nearly optimal forecast is achieved and we show that errors due to truncation of the dynamics are no greater than the errors due to covariance truncation.

Tangborn, Andrew↗

An Update on the Effectiveness of Postmission Disposal in LEO

The commonly‐adopted orbital debris mitigation measures were developed to reduce the growth of the future debris population. A major component in debris mitigation is post-mission disposal (PMD). The key PMD element for LEO satellites is the 25‐year rule. It is intended to limit the long‐term presence of rocket bodies (R/Bs) and spacecraft (S/C), as well as mission-related debris, in the environment. The effectiveness of PMD has been demonstrated and documented since the development of mitigation measures began in the 1990s. This paper summarizes an updated study, based on the current environment, using the NASA LEGEND model. The study focused on the > or = 10 cm population in LEO. The historical simulation covered 1957 through 2011 and followed the recorded launches and known breakup events. The future projection was carried out for 200 years. An eight‐year launch traffic, 2004 - 2011, was repeated during the projection period. An eight‐year mission lifetime was assumed for future S/C. No stationkeeping and no collision avoidance maneuver were implemented. Only objects 10 cm and larger were included in collision consideration. No explosion was allowed for R/Bs and S/C launched after 2011. The 25‐year PMD rule success rates were set at 0%, 10%, 50%, 75%, and 95%, respectively, for the 5 study scenarios. Results of the simulations were analyzed to quantify the differences among the different compliance rates. As expected, the 0% PMD projection followed a rapid and non‐linear increase in the next 200 years. The LEO population, on average, more than tripled at the end of the simulations. With a 50% compliance of the 25‐year rule, the population growth was reduced approximately by half. However, even with a 95% compliance of the 25‐year rule, the LEO debris population would still increase by an average of more than 50% in 200 years. These simulation results provide an updated assessment of the effectiveness of the 25‐year rule. It is the first and the most cost‐effective defense against future population growth. In addition, the results also confirm the instability of the LEO population and lay the foundation for the need to consider environment remediation in the future.

Liou, J.-C.↗

Solar Signals in CMIP-5 Simulations: The Stratospheric Pathway

The 11 year solar-cycle component of climate variability is assessed in historical simulations of models taken from the Coupled Model Intercomparison Project, phase 5 (CMIP-5). Multiple linear regression is applied to estimate the zonal temperature, wind and annular mode responses to a typical solar cycle, with a focus on both the stratosphere and the stratospheric influence on the surface over the period approximately 1850-2005. The analysis is performed on all CMIP-5 models but focuses on the 13 CMIP-5 models that resolve the stratosphere (high-top models) and compares the simulated solar cycle signature with reanalysis data. The 11 year solar cycle component of climate variability is found to be weaker in terms of magnitude and latitudinal gradient around the stratopause in the models than in the reanalysis. The peak in temperature in the lower equatorial stratosphere (approximately 70 hPa) reported in some studies is found in the models to depend on the length of the analysis period, with the last 30 years yielding the strongest response. A modification of the Polar Jet Oscillation (PJO) in response to the 11 year solar cycle is not robust across all models, but is more apparent in models with high spectral resolution in the short-wave region. The PJO evolution is slower in these models, leading to a stronger response during February, whereas observations indicate it to be weaker. In early winter, the magnitude of the modeled response is more consistent with observations when only data from 1979-2005 are considered. The observed North Pacific high-pressure surface response during the solar maximum is only simulated in some models, for which there are no distinguishing model characteristics. The lagged North Atlantic surface response is reproduced in both high- and low-top models, but is more prevalent in the former. In both cases, the magnitude of the response is generally lower than in observations.

stratosphere↗

STUDY OF THE η'π 0 SYSTEM IN GLUEX

The η'π is an important decay channel for the potentially exotic hybrid meson π 1 (1600) that has been detected in the COMPASS, VES and E852 experiments. These previous experiments used a charged pion beam as a probe. In contrast, the GlueX experiment is designed to search for the exotic hybrid mesons using a linearly polarized photon beam. This thesis project analyzes the η'π 0 system produced in GlueX by the reaction γp → η'π 0 p. After optimizing the event selection to minimize background contributions, a mass independent partial wave analysis is performed in search of an exotic signature. Within the currently available statistical precision, no signal was found.

Dotel, Rupesh↗

Radiation Shielding Analysis for the PIP-II Linac at Fermilab

The Proton Improvement Plan-II (PIP-II) [1] has been developed at Fermilab to provide powerful proton beams to the laboratory’s experiments. An 800-MeV superconducting linear accelerator—a centerpiece of the project—is currently under construction in Batavia, Illinois (USA). After completion, the superconducting linac will be the starting point for the 1.2 MW (Phase 1) and 2.4-MW (Phase 2) proton beam that is needed for the Long-Baseline Neutrino Facility (LBNF) at Fermilab [2]. Due to unavoidable loss of a fraction of the beam in the accelerator components, a certain level of radiation will be generated in the accelerator tunnel both during normal operation and at accidents. This work deals with radiation shielding design for the accelerator facility.

Rakhno, Igor↗

A life cycle cost economics model for automation projects with uniformly varying operating costs

The described mathematical model calculates life-cycle costs for projects with operating costs increasing or decreasing linearly with time. The cost factors involved in the life-cycle cost are considered, and the errors resulting from the assumption of constant rather than uniformly varying operating costs are examined. Parameters in the study range from 2 to 30 years, for project life; 0 to 15% per year, for interest rate; and 5 to 90% of the initial operating cost, for the operating cost gradient. A numerical example is presented.

Remer, D. S.↗

Identification of motion parameters of a rigid body from its orthogonal and perspective projections

An estimate is made of the motion parameters, namely, linear and angular velocities, of a rigid body rotating and translating in three-space. The authors assume that the velocities are constant and that the motion is not completely observable. They consider two separate cases of partial observations corresponding to the orthogonal and the perspective projections, respectively. If (x, y, z) is the Cartesian coordinate of the three-space, the authors assume in the first case that the projection of the motion on the x-y plane is observed. If (r, theta, phi) is the polar coordinates of the three-space, they assume in the second case that the parameter vector (theta, phi) is observed. The use of both of these cases to estimate the motion parameters is discussed.

Ganguly, S.↗

Otolith-Canal Convergence In Vestibular Nuclei Neurons

The current final report covers the period from June 1, 1999 to May 31, 2002. The primary objective of the investigation was to determine how information regarding head movements and head position relative to gravity is received and processed by central vestibular nuclei neurons in the brainstem. Specialized receptors in the vestibular labyrinths of the inner ear function to detect angular and linear accelerations of the head, with receptors located in the semicircular canals transducing rotational head movements and receptors located in the otolith organs transducing changes in head position relative to gravity or linear accelerations of the head. The information from these different receptors is then transmitted to central vestibular nuclei neurons which process the input signals, then project the appropriate output information to the eye, head, and body musculature motor neurons to control compensatory reflexes. Although a number of studies have reported on the responsiveness of vestibular nuclei neurons, it has not yet been possible to determine precisely how these cells combine the information from the different angular and linear acceleration receptors into a correct neural output signal. In the present project, rotational and linear motion stimuli were separately delivered while recording responses from vestibular nuclei neurons that were characterized according to direct input from the labyrinth and eye movement sensitivity. Responses from neurons receiving convergent input from the semicircular canals and otolith organs were quantified and compared to non-convergent neurons.

Dickman, J. David↗

Metal colloids and semiconductor quantum dots: Linear and nonlinear optical properties

One aspect of this project involves a collaborative effort with the Solid State Division of ORNL. The thrust behind this research is to develop ion implantion for synthesizing novel materials (quantum dots wires and wells, and metal colloids) for applications in all optical switching devices, up conversion, and the synthesis of novel refractory materials. In general the host material is typically a glass such as optical grade silica. The ions of interest are Au, Ag, Cd, Se, In, P, Sb, Ga and As. An emphasis is placed on host guest interactions between the matrix and the implanted ion and how the matrix effects and implantation parameters can be used to obtain designer level optical devices tailored for specific applications. The specific materials of interest are: CdSe, CdTe, InAs, GaAs, InP, GaP, InSb, GaSb and InGaAs. A second aspect of this research program involves using porous glass (25-200 A) for fabricating materials of finite size. In this part of the program, we are particularly interested in characterizing the thermodynamic and optical properties of these non-composite materials. We also address how phase diagram of the confined material is altered by the interfacial properties between the confined material and the pore wall.

Henderson, D. O.↗

Bayesian Vision for Shape Recovery

We present a new Bayesian vision technique that aims at recovering a shape from two or more noisy observations taken under similar lighting conditions. The shape is parametrized by a piecewise linear height field, textured by a piecewise linear irradiance field, and we assume Gaussian Markovian priors for both shape vertices and irradiance variables. The observation process. also known as rendering, is modeled by a non-affine projection (e.g. perspective projection) followed by a convolution with a piecewise linear point spread function. and contamination by additive Gaussian noise. We assume that the observation parameters are calibrated beforehand. The major novelty of the proposed method consists of marginalizing out the irradiances considered as nuisance parameters, which is achieved by Laplace approximations. This reduces the inference to minimizing an energy that only depends on the shape vertices, and therefore allows an efficient Iterated Conditional Mode (ICM) optimization scheme to be implemented. A Gaussian approximation of the posterior shape density is computed, thus providing estimates both the geometry and its uncertainty. We illustrate the effectiveness of the new method by shape reconstruction results in a 2D case. A 3D version is currently under development and aims at recovering a surface from multiple images, reconstructing the topography by marginalizing out both albedo and shading.

Jalobeanu, Andre↗

Projection techniques for iterative solution of A(bar x) equals (bar b) with successive right-hand sides

Two projection techniques for computing approximate solutions to linear systems of the form A(bar x)(sup n) = (bar b)(sup n), for a sequence n = 1, 2, ..., e.g., such as arises from time discretization of a partial differential equation, are presented. The inexpensive approximate solutions can be used as initial guesses for iterative solution of the system, resulting in significantly reduced computational expense. Examples of two- and three-dimensional incompressible Navier-Stokes calculations are presented in which x represents the pressure, and A is a discrete Poisson operator. In flows containing significant dynamic activity, these projection techniques lead to as much as a two-fold reduction in solution time.

Fischer, Paul F.↗