Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linearity preservation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Six-dimensional matching of intense beam with linear accelerating structure

Beam matching is a common technique that is routinely employed in accelerator design with the aim of minimizing beam losses and preservation of beam brightness. Despite being widely used, a full theoretical understanding of beam matching in 6D remains elusive. In this work, we present an analytical treatment of 6D beam matching of a high-intensity beam onto an RF structure. We begin our analysis within the framework of a linear model, and apply the averaging method to a set of 3D beam envelope equations. Accordingly, we obtain a matched solution that is comprised of smoothed envelopes and periodic terms, describing envelope oscillations with the period of the focusing structure. We then consider the nonlinear regime, where the beam size is comparable with the separatrix size. Stating with a Hamiltonian analysis in 6D phase space, we attain a self-consistent beam profile and show that it is significantly different from the commonly used ellipsoidal shape. Subsequently, we analyze the special case of an equilibrium with equal space charge depression between all degrees of freedom. Comparison of beam dynamics for equipartitioned, equal space charge depression, and equal emittances beams is given. Finally, we present experimental results on beam matching in the LANSCE linac.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Is Mini Beam Ready for Human Trials? Results of Randomized Study of Treating De-Novo Brain Tumors in Canines Using Linear Accelerator Generated Mini Beams

The main challenge in treating malignant brain neoplasms lies in eradicating the tumor while minimizing treatmentrelated damage. Conventional radiation treatments are associated with considerable side effects. Synchrotron generated micro-beam radiation (SMBRT) has shown to preserve brain architecture while killing tumor cells, however physical characteristics and limited facility access restrict its use. We have created a new clinical device which produces mini beams on a linear accelerator, to provide a new type of treatment called mini-beam radiation therapy (MBRT). The objective of this study is to compare the treatment outcomes of linear accelerator based MBRT versus standard radiation treatment (SRT), to evaluate the tumor response and the treatment-related changes in the normal brain with respect to each treatment type. Pet dogs with de-novo brain tumors were accrued for treatment. Dogs were randomized between standard fractionated stereotactic (9 Gy in 3 fractions) radiation treatment vs. a single fraction of MBRT (26 Gy mean dose). Dogs were monitored after treatment for clinical assessment and imaging. When the dogs were euthanized, a veterinary pathologist assessed the radiation changes and tumor response. We accrued 16 dogs, 8 dogs in each treatment arm. In the MBRT arm, 71% dogs achieved complete pathological remission. The radiation-related changes were all confined to the target region. Structural damage was not observed in the beam path outside of the target region. In contrast, none of the dogs in control group achieved remission and the treatment related damage was more extensive. Therapeutic superiority was observed with MBRT, including both tumor control and the normal structural preservation. The MBRT findings are suggestive of an immune related mechanism which is absent in standard treatment. These findings together with the widespread availability of clinical linear accelerators make MBRT a promising research topic to explore further treatment and clinical trial opportunities.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1991

Forward flight trim calculation and frequency response validation of a high-order helicopter simulation model

This paper describes a new trim procedure, that includes the calculation of the steady-state response of the rotor blades, and that is applicable to straight flight and steady coordinated turns. This paper also describes the results of a validation study for a high order linearized model of helicopter flight dynamics, that includes rotor, inflow, and actuator dynamics. The model is obtained by numerical perturbations of a nonlinear, blade element type mathematical model. Predicted responses are compared with flight test data for two values of flight speed. The comparison is carried out in the frequency domain. Numerical simulations show that the trim algorithm is very accurate, and preserves the periodicity of the aircraft states. The results also indicate that the predictions of the linearized model are in good agreement with flight test data, especially at medium and high frequencies.

Kim, Frederick D.↗
DOE OSTI2022

Structure of Shock Waves and Inelasticity in Shock-Compressed Cemented Tungsten Carbides

In this study, shock wave experiments are conducted on General Carbide cemented tungsten carbide with 3.7wt.% cobalt binder to determine its shock-induced compression behavior up to 100 GPa. The measured wave profiles indicate the cemented tungsten carbide to undergo elastic-plastic deformation during shock compression. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot Elastic Limit (HEL), an elastic-plastic ramp indicating substantial post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, the shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relation, and the longitudinal stress (σ x ) vs. specific volume (V) curve for the samples. The HEL of the material was determined to lie between 4.41 GPa and 4.58 GPa. The U s - u p relation was determined to be U s = 4.97 + 1. 457u p for particle velocities greater than 0.75 km/s. The measured plastic shock velocities for particle velocities less than 0.7 km/s were found to be larger than those predicted using the linear U s - u p Hugoniot relationship, indicating the cemented WC samples to preserve substantial shear strength in the post-yield deformation region. No phase transformation was observed up to 100 GPa.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Celestial amplitudes as AdS-Witten diagrams

Both celestial and momentum space amplitudes in four dimensions are beset by divergences resulting from spacetime translation and sometimes scale invariance. In this paper we consider a (linearized) marginal deformation of the celestial CFT for Yang-Mills theory which preserves 2D conformal invariance but breaks both spacetime translation and scale invariance and involves a chirally coupled massive scalar. The resulting MHV celestial amplitudes are completely finite (apart from the usual soft and collinear divergences and isolated poles in the sum of the weights) and take the canonical CFT form. Moreover, we show they can be simply rewritten in terms of AdS 3 -Witten contact diagrams which evaluate to the well-known D-functions, thereby forging a direct connection between flat and AdS holography.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Circularly Polarized Attosecond Pulses Enabled by an Azimuthal Phase and Polarization Grating

High-harmonic generation (HHG) is an extreme nonlinear optical process that can map the properties of an infrared driving laser beam onto short wavelength attosecond pulse trains. However, current techniques for generating circularly polarized high harmonics for probing magnetic materials and chiral systems have limitations: two-color collinear counter-rotating driving lasers result in a low cutoff photon energy, while single-color non-collinear counter-rotating schemes suffer from low conversion efficiency. In this work, we generate circularly polarized attosecond pulse trains by using a structured laser driver which has a rotating polarization and phase grating along the azimuthal coordinate. Furthermore, our experimental and numerical results demonstrate the production of left and right circularly polarized harmonics, which naturally separate upon propagation. Our approach uses a single laser color in a collinear geometry, that can be scaled for high efficiency. Simulations show this scheme can extend into the soft x-ray region when driven by mid-infrared driving lasers, while preserving the same high phase-matching cutoff photon energy as for linearly-polarized high harmonics.

attosecond science↗
DOE OSTI2023

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot elastic limit (HEL), an elastic-plastic ramp indicating post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Polarization aberration analysis of the Advanced X-ray Astrophysics Facility Telescope assembly

The Advanced X-ray Astrophysics Facility (AXAF) Telescope consists of six concentric paraboloid-hyperboloid pairs of mirrors that operate near grazing incidence. Because of the substantial polarization effects at large angles of incidence there has been concern regarding the feasibility of doing polarimetry near the telescope focal plane. The primary mirror is here shown to act as a tangentially directed half-wave linear retarder and almost completely to depolarize the linearly polarized component of the light. The secondary mirror introduces and additional half-wave of linear retardance. The two-mirror telescope assembly acts as a tangentially directed one-wave linear retarder. Each mirror depolarizes alone but together the two-mirror assembly preserves the polarization state. The net instrumental polarization effects are small and polarimetry is feasible with AXAF.

Chipman, Russell A.↗
NASA NTRS1995

Investigation of dispersion-relation-preserving scheme and spectral analysis methods for acoustic waves

Important characteristics of the aeroacoustic wave propagation are mostly encoded in their dispersion relations. Hence, a computational aeroacoustic (CAA) algorithm, which reasonably preserves these relations, was investigated. It was derived using an optimization procedure to ensure, that the numerical derivatives preserved the wave number and angular frequency of the differential terms in the linearized, 2-D Euler equations. Then, simulations were performed to validate the scheme and a compatible set of discretized boundary conditions. The computational results were found to agree favorably with the exact solutions. The boundary conditions were transparent to the outgoing waves, except when the disturbance source was close to a boundary. The time-domain data generated by such CAA solutions were often intractable until their spectra was analyzed. Therefore, the relative merits of three different methods were included in the study. For simple, periodic waves, the periodogram method produced better estimates of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window was more effective when used with the weighted-overlapped-segment-averaging and Blackman-Tukey methods gave better results than the periodogram method. Finally, it was demonstrated that the representation of time domain-data was significantly dependent on the particular spectral analysis method employed.

Vanel, Florence O.↗
DOE OSTI2023

Optimization of transformer ratio and beam loading in a plasma wakefield accelerator with a structure-exploiting algorithm

Plasma-based acceleration has emerged as a promising candidate as an accelerator technology for a future linear collider or a next-generation light source. We consider the plasma wakefield accelerator (PWFA) concept where a plasma wave wake is excited by a particle beam and a trailing beam surfs on the wake. For a linear collider, the energy transfer efficiency from the drive beam to the wake and from the wake to the trailing beam must be large, while the emittance and energy spread of the trailing bunch must be preserved. One way to simultaneously achieve this when accelerating electrons is to use longitudinally shaped bunches and nonlinear wakes. In the linear regime, there is an analytical formalism to obtain the optimal shapes. In the nonlinear regime, however, the optimal shape of the driver to maximize the energy transfer efficiency cannot be precisely obtained because currently no theory describes the wake structure and excitation process for all degrees of nonlinearity. In addition, the ion channel radius is not well defined at the front of the wake where the plasma electrons are not fully blown out by the drive beam. We present results using a novel optimization method to effectively determine a current profile for the drive and trailing beam in PWFA that provides low energy spread, low emittance, and high acceleration efficiency. We parameterize the longitudinal beam current profile as a piecewise-linear function and define optimization objectives. For the trailing beam, the algorithm converges quickly to a nearly inverse trapezoidal trailing beam current profile similar to that predicted by the ultrarelativistic limit of the nonlinear wakefield theory. For the drive beam, the beam profile found by the optimization in the nonlinear regime that maximizes the transformer ratio also resembles that predicted by linear theory. Furthermore, the current profiles found from the optimization method provide higher transformer ratios compared with the linear ramp predicted by the relativistic limit of the nonlinear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Geometry of Non-Diffusive Tracer Transport in Gridded Atmospheric Models

A first-order linear and numerically non-diffusive Eulerian transport algorithm, minVAR, was recently developed for preservation of correlations between interrelated tracers during advective transport. The present study extends this work by: (1) providing further investigation of several interesting geometric constructions found in contours of constant minVAR, short for minimum variance, through extension to three coordinate dimensions. These contours capture point-by-point representations of thousands of individual atmospheric aerosol and/or cloud particles as they evolve and are rendered on Eulerian grids at a level of sub-grid resolution limited only by numerical precision; and (2) exploration of geometric similarities between the Arakawa C-grid, used to obtain interpolated values of the wind field at grid scale and minVAR. In particular, we consider interpolation of the u and v horizontal components of wind velocity from grid to sub-grid scales. The last results are motivated by recent applications of the Weather Research and Forecasting (WRF) model applied in the coastal Houston region, where the recent TRacking Aerosol Convection Interactions ExpeRiment (TRACER) field campaign was organized. A unique and fully consistent mapping is obtained between particles moving along meteorological wind trajectories and the non-diffusive, non-dispersive representation of such trajectories on an Eulerian grid.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2024

Parton distributions from boosted fields in the Coulomb gauge

We propose a new method to calculate parton distribution functions (PDFs) from lattice correlations of boosted quarks and gluons in the Coulomb gauge. Compared to the widely used gauge-invariant Wilson-line operators, these correlations greatly simplify the renormalization thanks to the absence of linear power divergence. Additionally, they enable access to larger off-axis momenta under preserved 3D rotational symmetry, as well as enhanced long-range precision that facilitates the Fourier transform. We verify the factorization formula that relates this new observable to the quark PDF at one-loop order in perturbation theory. Moreover, through a lattice calculation of the pion valence quark PDF, we demonstrate the aforementioned advantage and features of the Coulomb gauge correlation and show that it yields consistent results with the gauge-invariant method. This opens the door to a more efficient way to calculate parton physics on the lattice.

Baryons↗
NASA NTRS2013

High-Order Space-Time Methods for Conservation Laws

Current high-order methods such as discontinuous Galerkin and/or flux reconstruction can provide effective discretization for the spatial derivatives. Together with a time discretization, such methods result in either too small a time step size in the case of an explicit scheme or a very large system in the case of an implicit one. To tackle these problems, two new high-order space-time schemes for conservation laws are introduced: the first is explicit and the second, implicit. The explicit method here, also called the moment scheme, achieves a Courant-Friedrichs-Lewy (CFL) condition of 1 for the case of one-spatial dimension regardless of the degree of the polynomial approximation. (For standard explicit methods, if the spatial approximation is of degree p, then the time step sizes are typically proportional to 1/p(exp 2)). Fourier analyses for the one and two-dimensional cases are carried out. The property of super accuracy (or super convergence) is discussed. The implicit method is a simplified but optimal version of the discontinuous Galerkin scheme applied to time. It reduces to a collocation implicit Runge-Kutta (RK) method for ordinary differential equations (ODE) called Radau IIA. The explicit and implicit schemes are closely related since they employ the same intermediate time levels, and the former can serve as a key building block in an iterative procedure for the latter. A limiting technique for the piecewise linear scheme is also discussed. The technique can suppress oscillations near a discontinuity while preserving accuracy near extrema. Preliminary numerical results are shown

Huynh, H. T.↗
DOE OSTI2024

A structure-preserving machine learning framework for accurate prediction of structural dynamics for systems with isolated nonlinearities

The nonlinearities present in structural systems are often found in isolated regions within the structure, such as those containing joints or interfaces. However, despite the localized nature of these nonlinearities their presence serves to couple together the modes of the underlying linear system and significantly complicate the development of appropriate reduced-order models; the localized nonlinearities have a global effect on the dynamics of the system. Further, in the presence of evolving structural health the nonlinearities can arise from accumulating damage, with dynamics distinct from those observed in the healthy state. The present work develops a data-driven formulation to identify and include the contributions of the isolated nonlinearities on the dynamics of the underlying linear structure. A novel coordinate separation is developed that decomposes those nonlinearities restricted to the isolated subdomain from the known linear system defined over the entire domain, and the influence of the isolated nonlinearities is reintroduced as an appropriately identified traction at the boundary of the isolated subdomain, referred to as the deviatoric force. In the region exterior to the nonlinear subdomain the response of the ideal linear system recovers that of the original nonlinear system. In this work, the deviatoric force component is predicted using a structure-preserving multilayer perceptron, based only on measured responses at the boundary of the isolated subdomain. Therefore introduction of the perceptron is able to bypass the direct numerical simulation of the nonlinearities within the isolated subdomain. This approach is illustrated through a simple structural system in which an interior region contains cubic nonlinearities and hysteretic damping. Once trained, the machine learning system is able to accurately predict the deviatoric force so that the ideal system recovers the response of the original system in the region outside the isolated nonlinear subdomain. Moreover, the data-driven approach is able to accurately predict the response when the system is subject to differing initial conditions and external excitation without the need for retraining, so that the proposed approach provides a robust description of the structural dynamics of the overall system.

Machine learning↗
NASA NTRS2015

Comparison of the Tangent Linear Properties of Tracer Transport Schemes Applied to Geophysical Problems.

A number of geophysical applications require the use of the linearized version of the full model. One such example is in numerical weather prediction, where the tangent linear and adjoint versions of the atmospheric model are required for the 4DVAR inverse problem. The part of the model that represents the resolved scale processes of the atmosphere is known as the dynamical core. Advection, or transport, is performed by the dynamical core. It is a central process in many geophysical applications and is a process that often has a quasi-linear underlying behavior. However, over the decades since the advent of numerical modelling, significant effort has gone into developing many flavors of high-order, shape preserving, nonoscillatory, positive definite advection schemes. These schemes are excellent in terms of transporting the quantities of interest in the dynamical core, but they introduce nonlinearity through the use of nonlinear limiters. The linearity of the transport schemes used in Goddard Earth Observing System version 5 (GEOS-5), as well as a number of other schemes, is analyzed using a simple 1D setup. The linearized version of GEOS-5 is then tested using a linear third order scheme in the tangent linear version.

4DVAR↗
DOE OSTI2026

Magnetic brightening and nanoscale imaging of spin-polarized helical edge modes in ZrTe 5

Efficient charge transport remains a fundamental challenge for nanoelectronic devices, as their performance is constrained by high dissipation and the impedance mismatch between high-frequency signal sources and nanoscale circuit interfaces. Although topologically protected helical edge modes offer a dissipationless and backscattering-resilient alternative, achieving nanoscale control over these modes requires direct visualization of their spatial electrodynamics, especially under high-frequency operation. Here, by utilizing cryogenic magneto-infrared scattering-type scanning near-field optical microscopy (cm-IR-sSNOM), we image spin-polarized helical edge channels in ZrTe 5 at the nanoscale, revealing magnetic-field-induced brightening as a high-frequency electrodynamic signature of robust topological edge modes. Operating at 1.8 K and under a magnetic field of up to 5 T, we observe the emergence of edge-state polarizability at infrared frequencies that is notably resilient to the magnetic gaps that typically quench d.c. and microwave edge transport. Our results reveal a topological ‘two-lane’ spatial reorganization in which an external magnetic field induces a spin-population imbalance between counterpropagating edge modes. Favoured helical branches are confined against physical boundaries to activate a net infrared near-field contrast. This electrodynamic response scales linearly with the number of atomic layers, which confirms that individual layers in ZrTe 5 preserve their discrete quantum spin Hall identities. These findings may be useful for developing topological spintronic devices, as the magnetic infrared tunability of helical edge channels provides a pathway for low-loss nanoscale interconnects and high-speed information processing.

36 MATERIALS SCIENCE↗