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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Limitations of Fault-Tolerant Quantum Linear System Solvers for Quantum Power Flow

Quantum computers hold promise for solving problems intractable for classical computers, especially those with high time or space complexity. Practical quantum advantage can be said to exist for such problems when the end-to-end time for solving such a problem using a classical algorithm exceeds that required by a quantum algorithm. Reducing the power flow (PF) problem into a linear system of equations allows for the formulation of quantum PF (QPF) algorithms, which are based on solving methods for quantum linear systems such as the Harrow-Hassidim-Lloyd (HHL) algorithm. Speedup from using QPF algorithms is often claimed to be exponential when compared to classical PF solved by state-of-the-art algorithms. Here, we investigate the potential for practical quantum advantage in solving QPF compared to classical methods on gate-based quantum computers. Notably, this paper does not present a new QPF solving algorithm but scrutinizes the end-to-end complexity of the QPF approach, providing a nuanced evaluation of the purported quantum speedup in this problem. Our analysis establishes a best-case bound for the HHL-based quantum power flow complexity, conclusively demonstrating that the HHL-based method has higher runtime complexity compared to the classical algorithm for solving the direct current power flow (DCPF) and fast decoupled load flow (FDLF) problem. Notably, our analysis and conclusions can be extended to any quantum linear system solver with rigorous performance guarantees, based on the known complexity lower bounds for this problem. Additionally, we establish that for potential practical quantum advantage (PQA) to exist it is necessary to consider DCPF-type problems with a very narrow range of condition number values and readout requirements.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2014

Simulations and Analysis of the Reshocked Inclined Interface Richtmyer–Meshkov Instability for Linear and Nonlinear Interface Perturbations

A computational study of the Richtmyer–Meshkov instability (RMI) is presented for an inclined interface perturbation in support of experiments being performed at the Texas A&M shock tube facility. The study is comprised of 2D, viscous, diffusive, compressible simulations performed using the arbitrary Lagrange Eulerian code, ARES, developed at Lawrence Livermore National Laboratory. These simulations were performed to late times after reshock with two initial interface perturbations, in the linear and nonlinear regimes each, prescribed by the interface inclination angle. The interaction of the interface with the reshock wave produced a complex 2D set of compressible wave interactions including expansion waves, which also interacted with the interface. Distinct differences in the interface growth rates prior to reshock were found in previous work. The current work provides in-depth analysis of the vorticity and enstrophy fields to elucidate the physics of reshock for the inclined interface RMI. After reshock, the two cases exhibit some similarities in integral measurements despite their disparate initial conditions but also show different vorticity decay trends, power law decay for the nonlinear and linear decay for the linear perturbation case.

42 ENGINEERING↗
DOE OSTI2021

Comparison of Rotary and Linear Cutting Methodology in Determining Specific Cutting Energy of Granite

Abstract Single polycrystalline diamond compact (PDC) cutting is a practical technique to understand the rock–tool interactions in drag-bit-type geothermal drilling operations. This paper introduces a rotary cutting method to determine specific cutting energy (SCE) and compares it with the conventional linear cutting method. In this work, granite is selected to represent hard rock formations in geothermal drilling. Cutting tests are conducted on a Computer Numerical Control (CNC) machine with a realistic cutting speed of 12.7 m/min and several chip loads ranging from 0.08 to 0.25 mm. The cutting force is measured using a dynamometer and then converted to specific cutting energy (SCE). The results show that the rotary method produces an inverse relationship between SCE and chip load, whereas the linear method shows the opposite. As a result, the produced SCE by the rotary method tends to be lower than that of the linear method at a higher chip load at and over 0.16 mm. The difference may be attributed to the cutting configuration and associated force components.

Engineering↗
DOE OSTI2020

High-resolution 3D microwave imaging using an optically tracked linear array

Active three-dimensional microwave or millimeter-wave imaging is typically performed using a linear array that is mechanically scanned over a uniform planar or cylindrical aperture. However, an image can also be formed with a linear array that is scanned over a non-uniform irregular aperture. This requires that the precise position of the array be known for each radar sample. This paper explores high-resolution 3D microwave imaging of stationary targets using an optically tracked manually scanned linear array. The image reconstruction techniques necessary to account for a non-uniform irregular aperture are also developed.

Clark, Richard T.↗
DOE OSTI2014

SU-E-T-526: On the Linearity, Stability and Beam Energy Dependence of CdSe Quantum Dots as Scintillating Probes

Purpose: Structure and energy transfer mechanisms confer colloidal quantum dots (cQDs) interesting properties, among them their potential as scintillators. CdSe multi-shell cQDs in powder were investigated under photons irradiation. The purpose of this work is to characterize signal to dose linearity, stability with time and to quantify the dependence of their light output with beam energy. Methods: The cQDs are placed at the extremity of a non-scintillating plastic collecting fiber, with the other extremity connected to an Apogee U2000C CCD camera. The CCD camera collects the fluorescence light from irradiated cQDs from which the delivered dose is extracted. This signal is corrected for Cerenkov contamination at MV energies using the chromatic technique. The detector was irradiated with two devices: Xstrahl 200 orthovoltage unit for 120, 180 and 220 kVp and a Varian Clinac iX for 6 and 23 MV. Results: Linear output response with varying dose is observed for all beam energies with R2 factors > 0,999. Reproducibility measurements were performed at 120 kVp: the same set-up was irradiated at different time intervals (one week and three months). The results showed only a small relative decrease of light output of 3,2 % after a combine deposited dose of approximately 95 Gy. CdSe nanocrystals response has been studied as a function of beam energy. The output increases with decreasing energy from 120 kVp to 6 MV and increase again for 23 MV. This behavior could be explained in part by the cQDs high-Z composition. Conclusion: The fluorescence light output of CdSe cQDs was found to be linear as a function of dose. The results suggest stability of the scintillation output of cQDs over time. The specific composition of cQDs is the main cause of the observed energy dependence. We will further look into particle beam dependence of the cQDs. Bourse d'excellence aux etudes graduees du CRC (Centre de Recherche sur le Cancer, Universite Laval) Bourse d'excellence aux etudes graduees du CRC (Centre de Recherche sur le Cancer, Universite Laval)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Poster — Thur Eve — 44: Linearization of Compartmental Models for More Robust Estimates of Regional Hemodynamic, Metabolic and Functional Parameters using DCE-CT/PET Imaging

Quantitative analysis of dynamic positron emission tomography (PET) data usually involves minimizing a cost function with nonlinear regression, wherein the choice of starting parameter values and the presence of local minima affect the bias and variability of the estimated kinetic parameters. These nonlinear methods can also require lengthy computation time, making them unsuitable for use in clinical settings. Kinetic modeling of PET aims to estimate the rate parameter k{sub 3}, which is the binding affinity of the tracer to a biological process of interest and is highly susceptible to noise inherent in PET image acquisition. We have developed linearized kinetic models for kinetic analysis of dynamic contrast enhanced computed tomography (DCE-CT)/PET imaging, including a 2-compartment model for DCE-CT and a 3-compartment model for PET. Use of kinetic parameters estimated from DCE-CT can stabilize the kinetic analysis of dynamic PET data, allowing for more robust estimation of k{sub 3}. Furthermore, these linearized models are solved with a non-negative least squares algorithm and together they provide other advantages including: 1) only one possible solution and they do not require a choice of starting parameter values, 2) parameter estimates are comparable in accuracy to those from nonlinear models, 3) significantly reduced computational time. Our simulated data show that when blood volume and permeability are estimated with DCE-CT, the bias of k{sub 3} estimation with our linearized model is 1.97 ± 38.5% for 1,000 runs with a signal-to-noise ratio of 10. In summary, we have developed a computationally efficient technique for accurate estimation of k{sub 3} from noisy dynamic PET data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-F-BRD-16: Relative Biological Effectiveness of Double-Strand Break Induction for Modeling Cell Survival in Pristine Proton Beams of Different Dose-Averaged Linear Energy Transfers

Purpose: High throughput in vitro experiments assessing cell survival following proton radiation indicate that both the alpha and the beta parameters of the linear quadratic model increase with increasing proton linear energy transfer (LET). We investigated the relative biological effectiveness (RBE) of double-strand break (DSB) induction as a means of explaining the experimental results. Methods: Experiments were performed with two lung cancer cell lines and a range of proton LET values (0.94 – 19.4 keV/µm) using an experimental apparatus designed to irradiate cells in a 96 well plate such that each column encounters protons of different dose-averaged LET (LETd). Traditional linear quadratic survival curve fitting was performed, and alpha, beta, and RBE values obtained. Survival curves were also fit with a model incorporating RBE of DSB induction as the sole fit parameter. Fitted values of the RBE of DSB induction were then compared to values obtained using Monte Carlo Damage Simulation (MCDS) software and energy spectra calculated with Geant4. Other parameters including alpha, beta, and number of DSBs were compared to those obtained from traditional fitting. Results: Survival curve fitting with RBE of DSB induction yielded alpha and beta parameters that increase with proton LETd, which follows from the standard method of fitting; however, relying on a single fit parameter provided more consistent trends. The fitted values of RBE of DSB induction increased beyond what is predicted from MCDS data above proton LETd of approximately 10 keV/µm. Conclusion: In order to accurately model in vitro proton irradiation experiments performed with high throughput methods, the RBE of DSB induction must increase more rapidly than predicted by MCDS above LETd of 10 keV/µm. This can be explained by considering the increased complexity of DSBs or the nature of intra-track pairwise DSB interactions in this range of LETd values. NIH Grant 2U19CA021239-35.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A linear, separable two-parameter model for dual energy CT imaging of proton stopping power computation

Purpose: To evaluate the accuracy and robustness of a simple, linear, separable, two-parameter model (basis vector model, BVM) in mapping proton stopping powers via dual energy computed tomography (DECT) imaging. Methods: The BVM assumes that photon cross sections (attenuation coefficients) of unknown materials are linear combinations of the corresponding radiological quantities of dissimilar basis substances (i.e., polystyrene, CaCl{sub 2} aqueous solution, and water). The authors have extended this approach to the estimation of electron density and mean excitation energy, which are required parameters for computing proton stopping powers via the Bethe–Bloch equation. The authors compared the stopping power estimation accuracy of the BVM with that of a nonlinear, nonseparable photon cross section Torikoshi parametric fit model (VCU tPFM) as implemented by the authors and by Yang et al. [“Theoretical variance analysis of single- and dual-energy computed tomography methods for calculating proton stopping power ratios of biological tissues,” Phys. Med. Biol. 55, 1343–1362 (2010)]. Using an idealized monoenergetic DECT imaging model, proton ranges estimated by the BVM, VCU tPFM, and Yang tPFM were compared to International Commission on Radiation Units and Measurements (ICRU) published reference values. The robustness of the stopping power prediction accuracy of tissue composition variations was assessed for both of the BVM and VCU tPFM. The sensitivity of accuracy to CT image uncertainty was also evaluated. Results: Based on the authors’ idealized, error-free DECT imaging model, the root-mean-square error of BVM proton stopping power estimation for 175 MeV protons relative to ICRU reference values for 34 ICRU standard tissues is 0.20%, compared to 0.23% and 0.68% for the Yang and VCU tPFM models, respectively. The range estimation errors were less than 1 mm for the BVM and Yang tPFM models, respectively. The BVM estimation accuracy is not dependent on tissue type and proton energy range. The BVM is slightly more vulnerable to CT image intensity uncertainties than the tPFM models. Both the BVM and tPFM prediction accuracies were robust to uncertainties of tissue composition and independent of the choice of reference values. This reported accuracy does not include the impacts of I-value uncertainties and imaging artifacts and may not be achievable on current clinical CT scanners. Conclusions: The proton stopping power estimation accuracy of the proposed linear, separable BVM model is comparable to or better than that of the nonseparable tPFM models proposed by other groups. In contrast to the tPFM, the BVM does not require an iterative solving for effective atomic number and electron density at every voxel; this improves the computational efficiency of DECT imaging when iterative, model-based image reconstruction algorithms are used to minimize noise and systematic imaging artifacts of CT images.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-F-J-41: Experimental Validation of a Cascaded Linear System Model for MVCBCT with a Multi-Layer EPID

Purpose: The purpose of this study was to validate the use of a cascaded linear system model for MV cone-beam CT (CBCT) using a multi-layer (MLI) electronic portal imaging device (EPID) and provide experimental insight into image formation. A validated 3D model provides insight into salient factors affecting reconstructed image quality, allowing potential for optimizing detector design for CBCT applications. Methods: A cascaded linear system model was developed to investigate the potential improvement in reconstructed image quality for MV CBCT using an MLI EPID. Inputs to the three-dimensional (3D) model include projection space MTF and NPS. Experimental validation was performed on a prototype MLI detector installed on the portal imaging arm of a Varian TrueBeam radiotherapy system. CBCT scans of up to 898 projections over 360 degrees were acquired at exposures of 16 and 64 MU. Image volumes were reconstructed using a Feldkamp-type (FDK) filtered backprojection (FBP) algorithm. Flat field images and scans of a Catphan model 604 phantom were acquired. The effect of 2×2 and 4×4 detector binning was also examined. Results: Using projection flat fields as an input, examination of the modeled and measured NPS in the axial plane exhibits good agreement. Binning projection images was shown to improve axial slice SDNR by a factor of approximately 1.4. This improvement is largely driven by a decrease in image noise of roughly 20%. However, this effect is accompanied by a subsequent loss in image resolution. Conclusion: The measured axial NPS shows good agreement with the theoretical calculation using a linear system model. Binning of projection images improves SNR of large objects on the Catphan phantom by decreasing noise. Specific imaging tasks will dictate the implementation image binning to two-dimensional projection images. The project was partially supported by a grant from Varian Medical Systems, Inc. and grant No. R01CA188446-01 from the National Cancer Institute.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-71: Improving CT Quality for Radiation Therapy Planning and Delivery Guidance Using a Non-Linear Contrast Enhancement Technique

Purpose: Advanced image post-processing techniques which enhance soft-tissue contrast in CT have not been widely employed for RT planning or delivery guidance. The purpose of this work is to assess the soft-tissue contrast enhancement from non-linear contrast enhancing filters and its impact in RT. The contrast enhancement reduces patient alignment uncertainties. Methods: Non-linear contrast enhancing methods, such as Best Contrast (Siemens), amplify small differences in X-ray attenuation between two adjacent structure without significantly increasing noise. Best Contrast (BC) separates a CT into two frequency bands. The low frequency band is modified by a non-linear scaling function before recombination with the high frequency band. CT data collected using a CT-on-rails (Definition AS Open, Siemens) during daily CT-guided RT for 6 prostate cancer patients and an image quality phantom (The Phantom Laboratory) were analyzed. Images acquired with a standard protocol (120 kVp, 0.6 pitch, 18 mGy CTDIvol) were processed before comparison to the unaltered images. Contrast and noise were measured in the the phantom. Inter-observer variation was assessed by placing prostate contours on the 12 CT study sets, 6 enhanced and 6 unaltered, in a blinded study involving 8 observers. Results: The phantom data demonstrate that BC increased the contrast between the 1.0% supra-slice element and the background substrate by 46.5 HU while noise increased by only 2.3 HU. Thus the contrast to noise ratio increased from 1.28 to 6.71. Furthermore, the variation in centroid position of the prostate contours was decreased from 1.3±0.4 mm to 0.8±0.3 mm. Thus the CTV-to-PTV margin was reduced by 1.1 mm. The uncertainty in delineation of the prostate/rectum edge decreased by 0.5 mm. Conclusion: As demonstrated in phantom and patient scans the BC filter accentuates soft-tissue contrast. This enhancement leads to reduced inter-observer variation, which should improve RT planning and delivery. Supported by Siemens.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-138: An Extension of PCA-Based Respiratory Deformation Modeling Via Multi-Linear Decomposition

Purpose: To address and lift the limited degree of freedom (DoF) of globally bilinear motion components such as those based on principal components analysis (PCA), for encoding and modeling volumetric deformation motion. Methods: We provide a systematic approach to obtaining a multi-linear decomposition (MLD) and associated motion model from deformation vector field (DVF) data. We had previously introduced MLD for capturing multi-way relationships between DVF variables, without being restricted by the bilinear component format of PCA-based models. PCA-based modeling is commonly used for encoding patient-specific deformation as per planning 4D-CT images, and aiding on-board motion estimation during radiotherapy. However, the bilinear space-time decomposition inherently limits the DoF of such models by the small number of respiratory phases. While this limit is not reached in model studies using analytical or digital phantoms with low-rank motion, it compromises modeling power in the presence of relative motion, asymmetries and hysteresis, etc, which are often observed in patient data. Specifically, a low-DoF model will spuriously couple incoherent motion components, compromising its adaptability to on-board deformation changes. By the multi-linear format of extracted motion components, MLD-based models can encode higher-DoF deformation structure. Results: We conduct mathematical and experimental comparisons between PCA- and MLD-based models. A set of temporally-sampled analytical trajectories provides a synthetic, high-rank DVF; trajectories correspond to respiratory and cardiac motion factors, including different relative frequencies and spatial variations. Additionally, a digital XCAT phantom is used to simulate a lung lesion deforming incoherently with respect to the body, which adheres to a simple respiratory trend. In both cases, coupling of incoherent motion components due to a low model DoF is clearly demonstrated. Conclusion: Multi-linear decomposition can enable decoupling of distinct motion factors in high-rank DVF measurements. This may improve motion model expressiveness and adaptability to on-board deformation, aiding model-based image reconstruction for target verification. NIH Grant No. R01-184173.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-44: A Comparison of the Pre-Plan, Intra-Operative Plan, and Post-Implant Dosimetry for a Prostate Implant Case Using Prefabricated Linear Polymer-Encapsulated Pd-103

Purpose: To investigate the reproducibility and limitations of Pd-103 prostate brachytherapy using fixed length linear sources (CivaString). Methods: An LDR prostate brachytherapy case which was preplanned on MR images with prefabricated linear polymer-encapsulated Pd-103 sources (CivaString) was studied and compared with ultrasound based intra-operative planning and CT based post-implant dosimetry. We evaluated the following parameters among the three studies: prostate geometry (volume and cross sectional area), needle position and alignment deviations, and dosimetry parameters (D90). Results: The prostate volumes and axial cross sectional areas at center of prostate were measured as 41.8, 39.3 and 36.8 cc, and 14.9, 14.3, and 11.3 respectively on pre-plan MR, inter-op US, and post-implant CT studies. The deviation of prostate volumes and axial cross sectional areas measured on pre-planning MR and intra-operative US were within 5%. 17 out of 19 pre-planned needles were positioned within 5mm (the template grid size). One needle location was adjusted intra-operatively and another needle was removed due to proximity to urethra. The needle pathways were not always parallel to the trans-rectal probe due to the flexibility of CivaString. The angle of deviation was up to 10 degrees. Two pairs of needles were exchanged to better fit the length of prostate at the time of implant. This resulted in a prostate D90 of 153.8 Gy (124%) and 131.4 Gy (106.7%) for intra-op and PID respectively. Conclusion: Preplanning is a necessary part of implants performed with prefabricated linear polymer sources. However, as is often the case, there were real-time deviations from the pre-plan. Intra-op planning provides the ability conform to anatomy at the time of implant. Therefore, we propose to develop a systematic way to order extra strings of different length to provide the flexibility to perform intra-operative planning with fixed length strands.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-130: [18F]-FDG Uptake Dose Response in Lung Correlates Linearly with Proton Therapy Dose

Purpose: Analysis of clinical outcomes in lung cancer patients treated with protons using 18F-FDG uptake in lung as a measure of dose response. Methods: A test case lung cancer patient was selected in an unbiased way. The test patient’s treatment planning and post treatment positron emission tomography (PET) were collected from picture archiving and communication system at the UT M.D. Anderson Cancer Center. Average computerized tomography scan was registered with post PET/CT through both rigid and deformable registrations for selected region of interest (ROI) via VelocityAI imaging informatics software. For the voxels in the ROI, a system that extracts the Standard Uptake Value (SUV) from PET was developed, and the corresponding relative biological effectiveness (RBE) weighted (both variable and constant) dose was computed using the Monte Carlo (MC) methods. The treatment planning system (TPS) dose was also obtained. Using histogram analysis, the voxel average normalized SUV vs. 3 different doses was obtained and linear regression fit was performed. Results: From the registration process, there were some regions that showed significant artifacts near the diaphragm and heart region, which yielded poor r-squared values when the linear regression fit was performed on normalized SUV vs. dose. Excluding these values, TPS fit yielded mean r-squared value of 0.79 (range 0.61–0.95), constant RBE fit yielded 0.79 (range 0.52–0.94), and variable RBE fit yielded 0.80 (range 0.52–0.94). Conclusion: A system that extracts SUV from PET to correlate between normalized SUV and various dose calculations was developed. A linear relation between normalized SUV and all three different doses was found.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-BRA-08: Diaphragm Motion Tracking Based On KV CBCT Projections with a Constrained Linear Regression Optimization

Purpose: To develop a novel strategy to extract the respiratory motion of the thoracic diaphragm from kilovoltage cone beam computed tomography (CBCT) projections by a constrained linear regression optimization technique. Methods: A parabolic function was identified as the geometric model and was employed to fit the shape of the diaphragm on the CBCT projections. The search was initialized by five manually placed seeds on a pre-selected projection image. Temporal redundancies, the enabling phenomenology in video compression and encoding techniques, inherent in the dynamic properties of the diaphragm motion together with the geometrical shape of the diaphragm boundary and the associated algebraic constraint that significantly reduced the searching space of viable parabolic parameters was integrated, which can be effectively optimized by a constrained linear regression approach on the subsequent projections. The innovative algebraic constraints stipulating the kinetic range of the motion and the spatial constraint preventing any unphysical deviations was able to obtain the optimal contour of the diaphragm with minimal initialization. The algorithm was assessed by a fluoroscopic movie acquired at anteriorposterior fixed direction and kilovoltage CBCT projection image sets from four lung and two liver patients. The automatic tracing by the proposed algorithm and manual tracking by a human operator were compared in both space and frequency domains. Results: The error between the estimated and manual detections for the fluoroscopic movie was 0.54mm with standard deviation (SD) of 0.45mm, while the average error for the CBCT projections was 0.79mm with SD of 0.64mm for all enrolled patients. The submillimeter accuracy outcome exhibits the promise of the proposed constrained linear regression approach to track the diaphragm motion on rotational projection images. Conclusion: The new algorithm will provide a potential solution to rendering diaphragm motion and ultimately improving tumor motion management for radiation therapy of cancer patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Cyclone–anticyclone vortex asymmetry mechanism and linear Ekman friction

Allowance for the linear Ekman friction has been found to ensure a threshold (in rotation frequency) realization of the linear dissipative–centrifugal instability and the related chiral symmetry breaking in the dynamics of Lagrangian particles, which leads to the cyclone–anticyclone vortex asymmetry. An excess of the fluid rotation rate ω{sub 0} over some threshold value determined by the fluid eigenfrequency ω (i.e., ω{sub 0} > ω) is shown to be a condition for the realization of such an instability. A new generalization of the solution of the Karman problem to determine the steady-state velocity field in a viscous incompressible fluid above a rotating solid disk of large radius, in which the linear Ekman friction was additionally taken into account, has been obtained. A correspondence of this solution and the conditions for the realization of the dissipative–centrifugal instability of a chiral-symmetric vortex state and the corresponding cyclone–anticyclone vortex asymmetry has been shown. A generalization of the well-known spiral velocity distribution in an “Ekman layer” near a solid surface has been established for the case where the fluid rotation frequency far from the disk ω differs from the disk rotation frequency ω{sub 0}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Analysis of the processes occurring in a submicrosecond discharge with a linear current density of up to 3 MA/cm through a thick-wall stainless-steel electrode

The state of conductors carrying a megampere current from the generator to the load is studied experimentally. It is found that the plasma produced from cylindrical stainless-steel tubes during the passage of a submicrosecond current pulse with a linear density of 3 MA/cm expands with a velocity of 5.5 km/s. Numerical results on the diffusion of the magnetic field induced by a current with a linear density of 1–3MA/cm into metal electrodes agree with the experimental data on the penetration time of the magnetic field. For a linear current density of 3.1 MA/cm, the experimentally determined electric field strength on the inner surface of the tube is 4 kV/cm. The calculated electric field strength on the inner surface of the tube turns out to be two times higher, which can be explained by plasma production on the outer and inner surfaces of the electrode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Quasi-linear gyrokinetic predictions of the Coriolis momentum pinch in NSTX

This paper presents quasi-linear gyrokinetic predictions of the Coriolis momentum pinch for low aspect-ratio NSTX H-modes where previous experimental measurements were focused. Local, linear calculations predict that in the region of interest (just outside the mid-radius) of these relatively high-beta plasmas, profiles are most unstable to microtearing modes that are only effective in transporting electron energy. However, sub-dominant electromagnetic and electrostatic ballooning modes are also unstable, which are effective at transporting energy, particles and momentum. The quasi-linear prediction of transport from these weaker ballooning modes, assuming they contribute transport in addition to that from microtearing modes in a nonlinear turbulent state, leads to a very small or outward convection of momentum, inconsistent with the experimentally measured inward pinch, and opposite to predictions in conventional aspect ratio tokamaks. Additional predictions of a low beta L-mode plasma, unstable to more traditional electrostatic ion temperature gradient-trapped electron mode instability, show that the Coriolis pinch is inward but remains relatively weak and insensitive to many parameter variations. The weak or outward pinch predicted in NSTX plasmas appears to be at least partially correlated to changes in the parallel mode structure that occur at finite beta and low aspect ratio, as discussed in previous theories. The only conditions identified where a stronger inward pinch is predicted occur either in the purely electrostatic limit or if the aspect ratio is increased. As the Coriolis pinch cannot explain the measured momentum pinch, additional theoretical momentum transport mechanisms are discussed that may be potentially important.

Tokamak Turbulence Transport Momentum pinch↗
DOE OSTI2019

Electron inertial effects on linearly polarized electromagnetic ion cyclotron waves at Earth's magnetosphere

We discuss a role of the electron inertial effect on linearly polarized electromagnetic ioncyclotron (EMIC) waves at Earth. The linearly polarized EMIC waves have been previouslysuggested to be generated via mode conversion from the fast compressional wave at the ion-ionhybrid (IIH) resonance. When the electron inertial effects are neglected, the wave normal angleof the mode-converted IIH waves is 90 degrees because the wavevector perpendicular to themagnetic field becomes infinite at the IIH resonance. When the electron inertial effect isconsidered, the mode-converted IIH waves can propagate across the magnetic field lines and thewavelength perpendicular to the magnetic field approaches the electron inertial length scale nearthe Buchsbaum resonance. These waves are referred to as electron inertial waves. Due to theelectron inertial effect, the perpendicular wavenumber to the ambient magnetic field near theIIH resonance remains finite and the wave normal angle is less than 90 degrees. The wave normalangle where the maximum absorption occurs in a dipole magnetic field is 30-80 degrees, whichis consistent with the observed values near the magnetic equator. Therefore, the numericalresults suggest that the linearly polarized EMIC wave generated via mode conversion near theIIH resonance can be detected in between the Buchsbaum and the IIH resonance frequencies,and these waves can have normal angle less than 90 degrees.

Kim, Eun-Hwa↗