Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearization”

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 307 records · Page 17

SU-E-T-384: Experimental Verification of a Monte Carlo Linear Accelerator Model Using a Radiochromic Film Stack Dosimeter

Purpose: To experimentally verify a Monte Carlo (MC) linear accelerator model for the simulation of intensity-modulated radiation therapy (IMRT) treatments of moving targets. Methods: A Varian Clinac™ 21EX linear accelerator was modeled using the EGSnrc user code BEAMnrc. The mean energy, radial-intensity distribution, and divergence of the electron beam incident on the bremsstrahlung target were adjusted to achieve agreement between simulated and measured percentage-depth-dose and transverse field profiles for a 6 MV beam. A seven-field step-and-shoot IMRT lung procedure was prepared using Varian Eclipse™ treatment planning software. The plan was delivered using a Clinac™ 21EX linear accelerator and measured with a Gafchromic™ EBT2 film stack dosimeter (FSD) in two separate static geometries: within a cylindrical water-equivalent-plastic phantom and within an anthropomorphic chest phantom. Two measurements were completed in each setup. The dose distribution for each geometry was simulated using the EGSnrc user code DOSXYZnrc. MC geometries of the treatment couch, cylindrical phantom, and chest phantom were developed by thresholding CT data sets using MATLAB™. The FSD was modeled as water. The measured and simulated dose distributions were normalized to the median dose within the FSD. Results: Using an electron beam with a mean energy of 6.05 MeV, a Gaussian radial-intensity distribution with a full width at half maximum of 1.5 mm, and a divergence of 0°, the measured and simulated dose profiles agree within 1.75% and 1 mm. Measured and simulated dose distributions within both the cylindrical and chest phantoms agree within 3% over 94% of the FSD volume. The overall uncertainty in the FSD measurements is 3.1% (k=1). Conclusion: MC simulations agree with FSD measurements within measurement uncertainty, thereby verifying the accuracy of the linear accelerator model for the simulation of IMRT treatments of static geometries. The experimental verification will be extended to treatments of targets undergoing three-dimensional motion.

60 APPLIED LIFE SCIENCES↗

SU-F-E-06: Dosimetric Characterization of Small Photons Beams of a Novel Linear Accelerator

Purpose: The aim of the present contribution was to measure the main dosimetric quantities of small fields produced by UNIQUE and evaluate its matching with the corresponding dosimetric data of one 21EX conventional linear accelerator (Varian) in operation at the same center. The second step was to evaluate comparative performance of the EDGE diode detector and the PinPoint micro-ionization chamber for dosimetry of small fields. Methods: UNIQUE is configured with MLC (120 leaves with 0.5 cm leaf width) and a single low photon energy of 6 MV. Beam data were measured with scanning EDGE diode detector (volume of 0.019 mm{sup 3}), a PinPoint micro-ionization chamber (PTW) and for larger fields (≥ 4×4cm{sup 2}) a PTW Semi flex chamber (0.125 cm{sup 3}) was used. The scanning system used was the 3D cylindrical tank manufactured by Sun Nuclear, Inc. The measurement of PDD and profiles were done at 100 cm SSD and 1.5 depth; the relative output factors were measured at 10 cm depth. Results: PDD and the profile data showed less than 1% variation between the two linear accelerators for fields size between 2×2 cm{sup 2} and 5×5cm{sup 2}. Output factor differences was less than 1% for field sizes between 3×3 cm{sup 2} and 10×10 cm{sup 2} and less of 1.5 % for fields of 1.5×1.5 cm{sup 2} and 2×2 cm{sup 2} respectively. The dmax value of the EDGE diode detector, measured from the PDD, was 8.347 mm for 0.5×0,5cm{sup 2} for UNIQUE. The performance of EDGE diode detector was comparable for all measurements in small fields. Conclusion: UNIQUE linear accelerator show similar dosimetrics characteristics as conventional 21EX Varian linear accelerator for small, medium and large field sizes.EDGE detector show good performance by measuring dosimetrics quantities in small fields typically used in IMRT and radiosurgery treatments.

60 APPLIED LIFE SCIENCES↗

Optimal Power Flow Derived Sparse Linear Solver Benchmarks

Due to the changing nature of the power grid, it is increasingly important to be able to solve a high-fidelity optimal power-flow models on large power networks. This high-fidelity problem, called AC Optimal Power Flow (ACOPF), is a nonlinear, nonconvex optimization problem. One of the few reliable ways of solving such a problem is interior point methods. These methods result in sparse linear systems where the coefficient matrix is symmetric, indefinite and nearly always ill-conditioned. As such, they are particularly challenging for sparse linear solvers and represent a considerable computational bottleneck in solving the ACOPF problem. In this paper, we introduce a repository of linear systems captured from ACOPF problems when solved by the open-source optimizer IPOPT. These matrices are meant to be used as a test suite for sparse linear solver development.

97 MATHEMATICS AND COMPUTING↗

Efficient Space–Time Reduced Order Model for Linear Dynamical Systems in Python Using Less than 120 Lines of Code

A classical reduced order model (ROM) for dynamical problems typically involves only the spatial reduction of a given problem. Recently, a novel space–time ROM for linear dynamical problems has been developed [Choi et al., Space–tume reduced order model for large-scale linear dynamical systems with application to Boltzmann transport problems, Journal of Computational Physics, 2020], which further reduces the problem size by introducing a temporal reduction in addition to a spatial reduction without much loss in accuracy. The authors show an order of a thousand speed-up with a relative error of less than 10−5 for a large-scale Boltzmann transport problem. In this work, we present for the first time the derivation of the space–time least-squares Petrov–Galerkin (LSPG) projection for linear dynamical systems and its corresponding block structures. Utilizing these block structures, we demonstrate the ease of construction of the space–time ROM method with two model problems: 2D diffusion and 2D convection diffusion, with and without a linear source term. For each problem, we demonstrate the entire process of generating the full order model (FOM) data, constructing the space–time ROM, and predicting the reduced-order solutions, all in less than 120 lines of Python code. We compare our LSPG method with the traditional Galerkin method and show that the space–time ROMs can achieve O(10−3) to O(10−4) relative errors for these problems. Depending on parameter–separability, online speed-ups may or may not be achieved. For the FOMs with parameter–separability, the space–time ROMs can achieve O(10) online speed-ups. Finally, we present an error analysis for the space–time LSPG projection and derive an error bound, which shows an improvement compared to traditional spatial Galerkin ROM methods.

97 MATHEMATICS AND COMPUTING↗

Sensitivity Study for Forecasting Variables of WRF-Solar Using a Tangent Linear Approach

Integrating solar generation in recent years has highlighted the need for improved accuracy in predicting solar power. Confidence in solar power forecasting can be achieved by designing an ensemble that provides reliable probabilistic information for solar radiation with reduced uncertainty and error. Ideally, ensemble members are created through the optimized perturbation of the initial conditions in numerical weather prediction (NWP) models. Tangent linear models are capable of efficiently investigating the sensitivity of solar radiation to model input parameters because they do not require individual perturbation of each variable. This sensitivity study using tangent linear models provide us the capability to identify the right variables to perturb in an ensemble prediction system. In this study, we developed tangent linear models for WRF-Solar modules that directly impact the computation of solar radiation and the simulation of cloud formation and dissipation including the Fast All-sky Model for Solar Applications (FARMS), the Noah land surface model (LSM), the Thompson microphysics, the Mello-Yamada-Nakanishi-Niino (MYNN) boundary layer parameterization, and the Deng scheme for a shallow-convection parameterization. A sensitivity analysis was conducted under various scenarios based on satellite observations and model simulations from the National Solar Radiation Data Base (NSRDB) and WRF-Solar, respectively. Critical forecasting variables that are highly sensitive to the forecasting of global horizontal irradiance (GHI), direct normal irradiance (DNI), cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes were determined using the relevant WRF-Solar module. This study will be used as a guidance on future research leading to high-quality probabilistic solar forecasting. In this presentation, we discuss the validation of tangent linear approach for WRF-Solar modules and illustrate how the sensitivity results are valuable in the improvement of probabilistic solar prediction.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Investigating the Linear Dependence of Direct and Indirect Radiative Forcing on Emission of Carbonaceous Aerosols in a Global Climate Model

Herein, the relationship between forcing and emission is investigated for black carbon (BC) and primary organic carbon (OC) emitted from North America and Asia. Direct and indirect radiative forcing (DRF and IRF) of BC and OC are simulated with CAM5.1. Two diagnostics are introduced to aid policy-relevant discussions: linearity and emission-normalized forcing. DRF is linearly related to emission for both BC and OC from two regions, and the linear relationship is similar, within 15%. IRF is linear in emissions when emissions are lower and regions are far from sources (North American BC and OC). Indirect radiative forcing is sublinear for strong sources and near-source regions (Asian OC). Emission-normalized IRF in North America is two to four times higher than that in Asia. The difference among regions and species is primarily caused by particle density as high density of BC results in fewer emitted particles and by the processes for accumulation mode particles to become cloud condensation nuclei and then to activate into cloud droplet. Lower emission-normalized IRF in more polluted regions means that reductions of OC in these regions would be relatively climate-neutral rather than causing significant warming via IRF reduction. An optimal aggregation area (30° × 30°) is identified for analysis of the forcing-to-emission relationship. For IRF, only 15–40% of the Earth’s surface is significantly affected by an emission region, but forcing in these regions comprises most of the global impact. Emission-normalized forcing can be used to estimate forcing changes due to emission reductions, as long as causes of nonlinearity are identified and considered.

54 ENVIRONMENTAL SCIENCES↗

The linear rheological responses of wedge‐type polymers

ABSTRACT The linear rheological responses of a series of specially designed wedge‐type polymers synthesized by the polymerization of large molecular weight monomers have been measured. These wedge polymers contained large side groups which contained three flexible branch chains per polymer chain unit. The master curves for these polymers were obtained by time temperature superposition of dynamic data at different temperatures from the terminal flow regime to well below the glass transition temperature, T g . While these polymers maintained a behavior similar to that of linear polymers, the influence of the large side group structure lead to low entanglement densities and extremely low rubbery plateau modulus values, being near to 13 kPa. The viscosity molecular weight dependence was also somewhat higher than that normally observed for linear polymers, tending toward a power law near to 4.2 rather than the typical 3.4 found in entangled linear chains. The glassy modulus of these branched polymers is also found to be extremely low, being less than 100 MPa at T g −60 °C. © 2015 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53 , 899–906

Hu, Miao↗

Linear two-pool models are insufficient to infer soil organic matter decomposition temperature sensitivity from incubations

Abstract Terrestrial carbon (C)-climate feedbacks depend strongly on how soil organic matter (SOM) decomposition responds to temperature. This dependency is often represented in land models by the parameter Q 10 , which quantifies the relative increase of microbial soil respiration per 10 °C temperature increase. Many studies have conducted paired laboratory soil incubations and inferred “active” and “slow” pool Q 10 values by fitting linear two-pool models to measured respiration time series. Using a recently published incubation study (Qin et al. in Sci Adv 5(7):eaau1218, 2019) as an example, here we first show that the very high parametric equifinality of the linear two-pool models may render such incubation-based Q 10 estimates unreliable. In particular, we show that, accompanied by the uncertain initial active pool size, the slow pool Q 10 can span a very wide range, including values as high as 100, although all parameter combinations are producing almost equally good model fit with respect to the observations. This result is robust whether or not interactions between the active and slow pools are considered (typically these interactions are not considered when interpreting incubation data, but are part of the predictive soil carbon models). This very large parametric equifinality in the context of interpreting incubation data is consistent with the poor temporal extrapolation capability of linear multi-pool models identified in recent studies. Next, using a microbe-explicit SOM model (RESOM), we show that the inferred two pools and their associated parameters (e.g., Q 10 ) could be artificial constructs and are therefore unreliable concepts for integration into predictive models. We finally discuss uncertainties in applying linear two-pool (or more generally multiple-pool) models to estimate SOM decomposition parameters such as temperature sensitivities from laboratory incubations. We also propose new observations and model structures that could enable better process understanding and more robust predictive capabilities of soil carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗

A homogeneous model for monotone mixed horizontal linear complementarity problems

In this paper, we propose a homogeneous model for the class of mixed horizontal linear complementarity problems. The proposed homogeneous model is always solvable and provides the solution of the original problem if it exists, or a certificate of infeasibility otherwise. Our formulation preserves the sparsity of the original formulation and does not reduce to the homogeneous model of the equivalent standard linear complementarity problem. We study the properties of the model and show that interior-point methods can be used efficiently for the numerical solutions of the homogeneous problem. Finally, numerical experiments show convincingly that it is more efficient to use the proposed homogeneous model for the mixed horizontal linear complementarity problem than to use known homogeneous models for the equivalent standard linear complementarity problem.

97 MATHEMATICS AND COMPUTING↗

Simple Expressions for the Design of Linear Tapers in Overmoded Corrugated Waveguides

In this paper, simple analytical formulae are presented for the design of linear tapers with very low mode conversion loss in overmoded corrugated waveguides. For tapers from waveguide radius a2 to a1, with a1 1 a 2 /λ. Here, λ is the wavelength of radiation. The fractional loss of the HE 11 mode in an optimized taper is 0.0293(a 2 -a 1 ) 4 /a$2\atop{1}$1a$2\atop{2}$. These formulae are accurate when a2≲2a 1 . Slightly more complex formulae, accurate for a 2 ≤4a 1 , are also presented in this paper. The loss in an overmoded corrugated linear taper is less than 1 % when a 2 ≤2.12a 1 and less than 0.1 % when a 2 ≤1.53a 1 . The present analytic results have been benchmarked against a rigorous mode matching code and have been found to be very accurate. The results for linear tapers are compared with the analogous expressions for parabolic tapers. Finally, parabolic tapers may provide lower loss, but linear tapers with moderate values of a 2 /a 1 may be attractive because of their simplicity of fabrication.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the Validity of Linear Response Theory in High-Dimensional Deterministic Dynamical Systems

This theoretical work considers the following conundrum: linear response theory is successfully used by scientists in numerous fields, but mathematicians have shown that typical low-dimensional dynamical systems violate the theory’s assumptions. Here we provide a proof of concept for the validity of linear response theory in high-dimensional deterministic systems for large-scale observables. We introduce an exemplary model in which observables of resolved degrees of freedom are weakly coupled to a large, inhomogeneous collection of unresolved chaotic degrees of freedom. By employing statistical limit laws we give conditions under which such systems obey linear response theory even if all the degrees of freedom individually violate linear response. We corroborate our result with numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spectral reflectance reconstruction with the locally weighted linear model

Spectral reflectance reconstruction by the local linear model has many advantages such as simple model structure, fast calculation speed and less chance to be over-fitting. However, local linear model often suffers from the outliers and prone to be under-fitting. In this paper, we propose a locally weighted linear model for spectral reflectance reconstruction. It improves the contribution of the local neighbors with higher similarity to the test point, which can reduce the influence of noisy points and redundant points. Experimental results show that the locally weighted linear model can effectively reduce the reconstruction error.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Precursor and Neutral Loss Scans in an RF Scanning Linear Quadrupole Ion Trap

Methodology for performing precursor and neutral loss scans in an RF scanning linear quadrupole ion trap is described and compared to the unconventional ac frequency scan technique. In the RF scanning variant, precursor ions are mass selectively excited by a fixed frequency resonance excitation signal at low Mathieu q while the RF amplitude is ramped linearly to pass ions through the point of excitation such that the excited ion’s m/z varies linearly with time. Ironically, a nonlinear ac frequency scan is still required for ejection of the product ions since their frequencies vary nonlinearly with the linearly varying RF amplitude. In the case of the precursor scan, the ejection frequency must be scanned so that it is fixed on a product ion m/z throughout the RF scan, whereas in the neutral loss scan, it must be scanned to maintain a constant mass offset from the excited precursor ions. Both simultaneous and sequential permutation scans are possible; only the former are demonstrated here. The scans described are performed on a variety of samples using different ionization sources: protonated amphetamine ions generated by nanoelectrospray ionization (nESI), explosives ionized by low-temperature plasma (LTP), and chemical warfare agent simulants sampled from a surface and analyzed with swab touch spray (TS). We lastly conclude that the ac frequency scan variant of these MS/MS scans is preferred due to electronic simplicity. In an accompanying manuscript, we thus describe the implementation of orthogonal double resonance precursor and neutral loss scans on the Mini 12 using constant RF voltage. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Generalized trigonometric Fourier-series method with automatic time step control for non-linear point kinetics equations

The trigonometric Fourier-series method (TFS) is generalized to provide approximate solutions for non-linear point kinetics equations with feedback using varying step sizes. This method can provide a very stable solution against the size of the discrete time step allowing much larger step sizes to be used. Systems of the point kinetics equations are solved using Fourier-series expansion over a partition of the total time interval. The approximate solution requires determining the series coefficients over each time step in that partition. These coefficients are determined using the high-order derivatives of the solution vector at the beginning of the time step introducing a system of linear algebraic equations to be solved at each step. This system is similar to the Vandermonde system. Two successive orders of the partial sums are used to estimate the local truncation error. This error and some other constrains are used to produce the largest step size allowable at each step while keeping the error within a specific tolerance. The process of calculating suitable step sizes should be automatic and inexpensive. Convergence and stability of the proposed method are discussed and a new formula is introduced to maintain stability. The proposed method solves the general linear and non-linear kinetics problems. The method has been applied to five different types of reactivities including step/ramp insertions with temperature feedback. The method is seemed to be valid for larger time intervals than those used in the conventional numerical integration, and is thus useful in reducing computing time. Computational results are found to be consistent with the analysis, they demonstrate that the convergence of the iteration scheme can be accelerated and the resulting computing times can be greatly reduced while maintain computational accuracy.

97 MATHEMATICS AND COMPUTING↗

Solving a system of linear equations: From centralized to distributed methods

For a wide range of control engineering applications, the problem of solving a system of linear equations is often encountered and has been well studied. Traditionally, this problem has been mainly solved in a centralized manner. However, for applications related to large-scale complex networked systems, centralized algorithms are often subjected to some practical issues due to limited computational power and communication bandwidth. As a promising and viable alternative, distributed algorithms can effectively address the issues associated with centralized algorithms by solving the problem efficiently in a multi-agent setting that accords with the distributed nature of networked systems. Distributed algorithms decompose the entire problem into many sub-problems that are solved by individual agents in a cooperative manner. In this survey paper, we provide a detailed overview of the state of the art relevant to distributed algorithms for solving a system of linear equations. We will first review basic distributed algorithms including both discrete-time and continuous-time algorithms. Then we will discuss the extended algorithms to achieve communication efficiency. Furthermore, we will also introduce distributed algorithms to obtain the minimum-norm solution for a system of linear equations with multiple solutions, as well as the least-squares solutions when there is no solution. Finally, the relationship of distributed algorithms for solving a system of linear equations to the existing distributed optimization algorithms is discussed.

Wang, Peng↗

High-order, linearly stable, partitioned solvers for general multiphysics problems based on implicit–explicit Runge–Kutta schemes

This work introduces a general framework for constructing high-order, linearly stable, partitioned solvers for multiphysics problems from a monolithic implicit–explicit Runge–Kutta (IMEX-RK) discretization of the semi-discrete equations. The generic multiphysics problem is modeled as a system of n systems of partial differential equations where the i th subsystem is coupled to the other subsystems through a coupling term that can depend on the state of all the other subsystems. This coupled system of partial differential equations reduces to a coupled system of ordinary differential equations via the method of lines where an appropriate spatial discretization is applied to each subsystem. The coupled system of ordinary differential equations is taken as a monolithic system and discretized using an IMEX-RK discretization with a specific implicit–explicit decomposition that introduces the concept of a predictor for the coupling term. We propose four coupling predictors that enable the monolithic system to be solved in a partitioned manner, i.e., subsystem-by-subsystem, and preserve the IMEX-RK structure and therefore the design order of accuracy of the monolithic scheme. The four partitioned solvers that result from these predictors are high-order accurate, allow for maximum re-use of existing single-physics software, and two of the four solvers allow the subsystems to be solved in parallel at a given stage and time step. We also analyze the stability of a coupled, linear model problem with a specific coupling structure and show that one of the partitioned solvers achieves unconditional linear stability for this problem, while the others are unconditionally stable only for certain values of the coupling strength. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems including a simple linear system of ODEs, advection–diffusion–reaction systems, fluid–structure interaction problems, and particle-laden flows, where we verify the design order of the IMEX schemes and study various stability properties.

97 MATHEMATICS AND COMPUTING↗

Lessons to be learned from a contentious challenge to mainstream radiobiological science (the linear no-threshold theory of genetic mutations)

There are both statistically valid and invalid reasons why scientists with differing default hypotheses can disagree in high-profile situations. Examples can be found in recent correspondence in this journal, which may offer lessons for resolving challenges to mainstream science, particularly when adherents of a minority view attempt to elevate the status of outlier studies and/or claim that self-interest explains the acceptance of the dominant theory. Edward J. Calabrese and I have been debating the historical origins of the linear no-threshold theory (LNT) of carcinogenesis and its use in the regulation of ionizing radiation. Professor Calabrese, a supporter of hormesis, has charged a committee of scientists with misconduct in their preparation of a 1956 report on the genetic effects of atomic radiation. Specifically he argues that the report mischaracterized the LNT research record and suppressed calculations of some committee members. After reviewing the available scientific literature, I found that the contemporaneous evidence overwhelmingly favored a (genetics) LNT and that no calculations were suppressed. Calabrese's claims about the scientific record do not hold up primarily because of lack of attention to statistical analysis. Ironically, outlier studies were more likely to favor supra-linearity, not sub-linearity. Finally, the claim of investigator bias, which underlies Calabrese's accusations about key studies, is based on misreading of text. Attention to ethics charges, early on, may help seed a counter narrative explaining the community's adoption of a default hypothesis and may help focus attention on valid evidence and any real weaknesses in the dominant paradigm. - Highlights: • Edward J Calabrese has made a contentious challenge to mainstream radiobiological science. • Such challenges should not be neglected, lest they enter the political arena without review. • Key genetic studies from the 1940s, challenged by Calabrese, were found consistent and unbiased. • A 1956 genetics report did not hide estimates and does not need investigation for misconduct. • The scientific record was strong for a no-threshold, linear genetic response to radiation.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Kernel reconstruction methods for Doppler broadening — Temperature interpolation by linear combination of reference cross sections at optimally chosen temperatures

This paper establishes a new family of methods to perform temperature interpolation of nuclear interactions cross sections, reaction rates, or cross sections times the energy. One of these quantities at temperature T is approximated as a linear combination of quantities at reference temperatures (T j ). The problem is formalized in a cross section independent fashion by considering the kernels of the different operators that convert cross section related quantities from a temperature T 0 to a higher temperature T — namely the Doppler broadening operation. Doppler broadening interpolation of nuclear cross sections is thus here performed by reconstructing the kernel of the operation at a given temperature T by means of linear combination of kernels at reference temperatures (T j ). The choice of the L 2 metric yields optimal linear interpolation coefficients in the form of the solutions of a linear algebraic system inversion. The optimization of the choice of reference temperatures (T j ) is then undertaken so as to best reconstruct, in the L∞ sense, the kernels over a given temperature range [T min ,T max ]. The performance of these kernel reconstruction methods is then assessed in light of previous temperature interpolation methods by testing them upon isotope 238 U. Temperature-optimized free Doppler kernel reconstruction significantly outperforms all previous interpolation-based methods, achieving 0.1% relative error on temperature interpolation of 238 U total cross section over the temperature range [300 K,3000 K] with only 9 reference temperatures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗