Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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 1,153 records · Page 64

Analysis of the track- and dose-averaged LET and LET spectra in proton therapy using the GEANT4 Monte Carlo code

Purpose: The motivation of this study was to find and eliminate the cause of errors in dose-averaged linear energy transfer (LET) calculations from therapeutic protons in small targets, such as biological cell layers, calculated using the GEANT 4 Monte Carlo code. Furthermore, the purpose was also to provide a recommendation to select an appropriate LET quantity from GEANT 4 simulations to correlate with biological effectiveness of therapeutic protons. Methods: The authors developed a particle tracking step based strategy to calculate the average LET quantities (track-averaged LET, LET{sub t} and dose-averaged LET, LET{sub d}) using GEANT 4 for different tracking step size limits. A step size limit refers to the maximally allowable tracking step length. The authors investigated how the tracking step size limit influenced the calculated LET{sub t} and LET{sub d} of protons with six different step limits ranging from 1 to 500 μm in a water phantom irradiated by a 79.7-MeV clinical proton beam. In addition, the authors analyzed the detailed stochastic energy deposition information including fluence spectra and dose spectra of the energy-deposition-per-step of protons. As a reference, the authors also calculated the averaged LET and analyzed the LET spectra combining the Monte Carlo method and the deterministic method. Relative biological effectiveness (RBE) calculations were performed to illustrate the impact of different LET calculation methods on the RBE-weighted dose. Results: Simulation results showed that the step limit effect was small for LET{sub t} but significant for LET{sub d}. This resulted from differences in the energy-deposition-per-step between the fluence spectra and dose spectra at different depths in the phantom. Using the Monte Carlo particle tracking method in GEANT 4 can result in incorrect LET{sub d} calculation results in the dose plateau region for small step limits. The erroneous LET{sub d} results can be attributed to the algorithm to determine fluctuations in energy deposition along the tracking step in GEANT 4. The incorrect LET{sub d} values lead to substantial differences in the calculated RBE. Conclusions: When the GEANT 4 particle tracking method is used to calculate the average LET values within targets with a small step limit, such as smaller than 500 μm, the authors recommend the use of LET{sub t} in the dose plateau region and LET{sub d} around the Bragg peak. For a large step limit, i.e., 500 μm, LET{sub d} is recommended along the whole Bragg curve. The transition point depends on beam parameters and can be found by determining the location where the gradient of the ratio of LET{sub d} and LET{sub t} becomes positive.

60 APPLIED LIFE SCIENCES↗

Heat flow control and segregation in directional solidification

Optimization of the vertical Bridgman technique for growth of electronic materials in single crystal form was investigated. The limitations of the crystal growth configuration were experimentally determined and heat transfer related deficiencies identified. Design of an alternate system was based on the use of heat pipes separated by a gradient region. Heat transfer analyses based on one and two dimensional models indicated the necessity of a flexible gradient zone configuration. Directional melting of binary systems as encountered during seeding in melt growth was analysed for concurrent compositional changes at the crystal-metal interface, and the theoretical treatment numerically applied to HgCdTe and Ga doped germanium. A theoretical and experimental study of the thermal effects associated with current flow was conducted. It was found that experimental measurements of dc induced growth during crystal pulling can be used for the precise determination of the Peltier coefficient.

Witt, A. F.↗

Active control of bearing preload using piezoelectric translators

In many spacecraft applications, mechanisms are required to perform precision pointing operations or to sometimes dither about or track a moving object. These mechanisms perform in a predictable and repeatable manner in benign temperature environments. Severe thermal gradients experienced in actual space applications however, cause assemblies to expand and contract around their bearings. This results in unpredictable changes in bearing preload, and hence bearing friction. This becomes a limitation for servos controlling pointing accuracy. Likewise, uncontrollable vibrations may couple into fixed preload (hence, fixed stiffness) mechanisms and limit pointing accuracy. Consequently, a complex problem faced today is how to design mechanisms that remain insensitive to changing thermal and vibrational spacecraft environments. Research presented involves the simplified modeling and test results of an actuator module that used piezoelectrically preload controlled bearings. The feasibility of actively controlling bearing preload was demonstrated. Because bearing friction is related to preload, a thermally active system designed with aluminum components and a 440 C bearing, was friction tested at temperatures ranging from 0 to 70 C (32 to 158 F). Effectiveness of the translators were demonstrated by mapping a controllable friction range throughout tested temperatures. It was learned that constant preload for this system could be maintained over an approximate 44 C (79 F) temperature span. From testing, it was also discovered that at the more deviate temperatures, expansions were so large that radial clearances were taken up and the duplex bearing became radially preloaded. Thus, active control of bearing preload is feasible but may be limited by inherent geometry constraints and materials used in the system.

Nye, Ted W.↗

Hypersonic flow computations around re-entry vehicles

The development of an algorithm for the solution of the compressible Euler equations at high Mach numbers on unstructured tetrahedral meshes is described. The basic algorithm is constructed in the form of a central difference scheme plus an explicit added artificial viscosity based upon fourth order differences of the solution. The stability of the solution in the vicinity of strong gradients is preserved by the incorporation of an additional artificial viscosity based upon a second order difference. Higher order accuracy is regained by using the ideas of flux corrected transport to limit the amount of added viscosity. The solution is advanced to steady state by means of an explicit multi-stage time-stepping method. The computational efficiency of the complete process is improved by incorporating an unstructured multigrid acceleration procedure. A number of flows of practical interest are analyzed to demonstrate the numerical performance of the proposed approach.

Peraire, J.↗

Conditions for Aeronomic Applicability of the Classical Electron Heat Conduction Formula

Conditions for the applicability of the classical formula for heat conduction in the electrons in ionized gas are investigated. In a fully ionised gas ( V(sub en) much greater than V(sub ei)), when the mean free path for electron-electron (or electron-ion) collisions is much larger than the characteristic thermal scale length of the observed system, the conditions for applicability break down. In the case of the Venus ionosphere this breakdown is indicated for a large fraction of the electron temperature data from altitudes greater than 180 km, for electron densities less than 10(exp 4)/cc cm. In a partially ionised gas such that V(sub en) much greater than V(sub ei) there is breakdown of the formula not only when the mean free path of electrons greatly exceeds the thermal scale length, but also when the gradient of neutral particle density exceeds the electron thermal gradient. It is shown that electron heat conduction may be neglected in estimating the temperature of joule heated electrons by observed strong 100 Hz electric fields when the conduction flux is limited by the saturation flux. The results of this paper support our earlier aeronomical arguments against the hypothesis of planetary scale whistlers for the 100 Hz electric field signal. In turn this means that data from the 100 Hz signal may not be used to support the case for lightning on Venus.

Cole, K. D.↗

Gravitational Influences on Flame Propagation through Non-Uniform, Premixed Gas Systems

Flame propagation through non-uniformly premixed (or layered) gases has importance both in useful combustion systems and in unintentional fires. As summarized previously, non-uniform premixed gas combustion receives scant attention compared to the more usual limiting cases of diffusion or uniformly premixed flames, especially regarding the role gravity plays. This paper summarizes our progress on furthering the knowledge of layered combustion, in which a fuel concentration gradient exists normal to the direction of flame spread. We present experimental and numerical results for flame spread through propanol-air layers formed near the flash point temperature (25 C) or near the stoichiometric temperature (33 C). Both the model and experimental results show that the removal of gravity results in a faster spreading flame, by as much as 80% depending on conditions. This is exactly the opposite effect as that predicted by an earlier model reported. We also found that having a gallery lid results in faster flame spread, an effect more pronounced at normal gravity, demonstrating the importance of enclosure geometry. Also reported here is the beginning of our spectroscopic measurements of fuel vapor.

Miller, Fletcher J.↗

Estimation of diastolic intraventricular pressure gradients by Doppler M-mode echocardiography

Previous studies have shown that small intraventricular pressure gradients (IVPG) are important for efficient filling of the left ventricle (LV) and as a sensitive marker for ischemia. Unfortunately, there has previously been no way of measuring these noninvasively, severely limiting their research and clinical utility. Color Doppler M-mode (CMM) echocardiography provides a spatiotemporal velocity distribution along the inflow tract throughout diastole, which we hypothesized would allow direct estimation of IVPG by using the Euler equation. Digital CMM images, obtained simultaneously with intracardiac pressure waveforms in six dogs, were processed by numerical differentiation for the Euler equation, then integrated to estimate IVPG and the total (left atrial to left ventricular apex) pressure drop. CMM-derived estimates agreed well with invasive measurements (IVPG: y = 0.87x + 0.22, r = 0.96, P < 0.001, standard error of the estimate = 0.35 mmHg). Quantitative processing of CMM data allows accurate estimation of IVPG and tracking of changes induced by beta-adrenergic stimulation. This novel approach provides unique information on LV filling dynamics in an entirely noninvasive way that has previously not been available for assessment of diastolic filling and function.

Non-NASA Center↗

Transition from ITG to MTM linear instabilities near pedestals of high density plasmas

Here, investigation of linear gyrokinetic ion-scale modes ( k θ ρ s = 0.3) finds that a transition from ion temperature gradient to microtearing mode (MTM) dominance occurs as the density is increased near the pedestal region of a parameterized DIII-D sized tokamak. H-modes profile densities, temperatures, and equilibria are parameterized utilizing the OMFIT PRO_create module. With these profiles, linear gyrokinetic ion-scale instabilities are predicted with CGYRO. This transition (nMTM) has a weak dependence on radial location in the region near the top of the pedestal ( ρ = 0.7 – 0.9), which allows simulating single radii to examine the approximate scaling of nMTM with global parameters. The critical nMTM is found to scale with plasma current. Additionally, increasing the minor radius by decreasing the aspect ratio and increasing the major radius are found to reduce nMTM. However, any relationship between nMTM and density limit physics remains unclear as nMTM increases relative to the Greenwald density with larger minor radius and with larger magnetic field, suggesting that the transport due to MTM may be less important for a reactor. Additionally, nMTM is sensitive to the pedestal temperature, the local electron and ion gradients, the ratio of ion to electron temperature T i / T e, and the current profile. MTMs are predicted to be the dominant instability in the core at similar Greenwald fractions for DIII-D, NSTX, and NSTX-U H-mode experiments, supporting the results of the parameterized study. Additionally, MTMs continue to be the dominant linear instability in a DIII-D L-mode after an H–L transition as the plasma approaches a density limit disruption despite the large change in plasma profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-E-T-262: Planning for Proton Pencil Beam Scanning (PBS): Applications of Gradient Optimization for Field Matching

Purpose: To propose the gradient optimization(GO) approach in planning for matching proton PBS fields and present two commonly used applications in our institution. Methods: GO is employed for PBS field matching in the scenarios that when the size of the target is beyond the field size limit of the beam delivery system or matching is required for beams from different angles to either improve the sparing of important organs or to pass through a short and simple beam path. Overlap is designed between adjacent fields and in the overlapped junction, the dose was optimized such that it gradually decreases in one field and the decrease is compensated by increase from another field. Clinical applications of this approach on craniospinal irradiation(CSI) and whole pelvis treatment were presented. Mathematical model was developed to study the relationships between dose errors, setup errors and junction lengths. Results: Uniform and conformal dose coverage to the entire target volumes was achieved for both applications using GO approach. For CSI, the gradient matching (6.7cm junction) between fields overcame the complexity of planning associated with feathering match lines. A slow dose gradient in the junction area significantly reduced the sensitivity of the treatment to setup errors. For whole pelvis, gradient matching (4cm junction) between posterior fields for superior target and bilateral fields for inferior target provided dose sparing to organs such as bowel, bladder and rectum. For a setup error of 3 mm in longitudinal direction from one field, mathematical model predicted dose errors of 10%, 6% and 4.3% for junction length of 3, 5 and 7cm. Conclusion: This GO approach improves the quality of the PBS treatment plan with matching fields while maintaining the safety of treatment delivery relative to potential misalignments.

60 APPLIED LIFE SCIENCES↗

Overview of density pedestal structure: Role fueling versus transport

How much of the density pedestal structure profile is determined by fueling versus transport is still a valid question. In most current tokamaks, the Scrape-Off Layer is not opaque enough to screen the neutrals and limit their penetration past the separatrix. As such recycled neutrals as well as additional gas puffing can directly affect the pedestal structure through ionization. Yet, in ITER and other future burning plasma magnetic confinement devices, this ionization source will be pushed further out into the Scrape-Off Layer. This poses the question, how much of the pedestal density structure is governed by this edge ionization source versus plasma transport effects. Theoretically, several turbulent modes have been identified which could provide a 'pinch' like up-gradient transport mechanism. Up-gradient transport is necessary in a source-free region to obtain a peaked density profile, which is often observed in the core of tokamaks. In an opaque Scrape-Off Layer, without the existence of such a pinch, the pedestal density structure would eventually disappear. Not only would this limit our ability to fuel a fusion reactor, it would also affect the pedestal stability. In this paper we will give an overview of the state of the art on what sets the pedestal density structure based on experimental observations as well as theoretical models. We will complement these studies with new results from Alcator C-Mod and DIII-D where the opaqueness is increased to reach values similar to those in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Digital active control law synthesis for aeroservoelastic systems

This paper presents a formulation for synthesis of digital active control laws for aeroservoelastic systems, which are typically modeled by large order equations in order to accurately represent the rigid and flexible body modes, unsteady aerodynamic forces, actuator dynamics, and gust spectra. The control law is expected to satisfy multiple design requirements on the dynamic loads, responses, actuator deflection and rate limitations, as well as maintain certain stability margins, yet should be simple enough to be implemented by an onboard digital microprocessor. The synthesis procedure minimizes a linear quadratic Gaussian type cost function, by updating selected free parameters of the control law, while satisfying a set of inequality constraints on the design loads, responses and stability margins. A stable classical control law or an estimator based full or reduced order control law can be used as an initial design starting point. The gradients of the cost function and the constraints, with respect to the digital control law design variables are derived analytically, to facilitate rapid convergence. Selected design responses can be treated as constraints instead of lumping them into the cost function, in order to satisfy individual root-mean-square load and response limitations. Constraints are also imposed on the minimum singular value requirements for stability robustness improvement.

Mukhopadhyay, Vivekananda↗

Testing the light:nutrient hypothesis: Insights into biofilm structure and function using metatranscriptomics

Aquatic biofilms are hotspots of biogeochemical activity due to concentrated microbial biomass (Battin, Kaplan, Newbold, & Hansen, 2003). However, biofilms are often considered a single entity when their role in biogeochemical transformations is assessed, even though these biofilms harbour functionally diverse microbial communities (Battin, Besemer, Bengtsson, Romani, & Packmann, 2016; Veach, Stegen, Brown, Dodds, & Jumpponen, 2016). Often overlooked are the biotic interactions among biofilm components that can affect ecosystem-scale processes such as primary production and nutrient cycling. These interactions are likely to be especially important under resource limitation. Light is a primary resource mediating algal photosynthesis and both phototrophic and heterotrophic production due to bacterial reliance on C-rich algal exudates (Cole, 1982). However, current understanding of function–structure linkages in streams has yet to unravel the relative degree of these microbial feedbacks under resource availability gradients. In this issue of Molecular Ecology, Bengtsson, Wagner, Schwab, Urich, and Battin (2018) studied stream biofilm responses to light availability to understand its impact across three domains of life. By integrating biogeochemical rate estimation and metatranscriptomics within a microcosm experiment, they were able to link primary production and nutrient uptake rates to algal and bacterial metabolic processes and specify what taxa contributed to gene expression. Under low light, diatoms and cyanobacteria upregulated photosynthetic machinery and diatom-specific chloroplast rRNA suggesting heightened transcriptional activity under light limitation to maintain phototrophic energy demands. Under high light, heterotrophic bacteria upregulated mRNAs related to phosphorous (P) metabolism while biofilm P uptake increased indicating high bacterial-specific P demand when algal biomass was high. Together, these results indicate that biogeochemical function is mediated by complex microbial interactions across trophic levels.

59 BASIC BIOLOGICAL SCIENCES↗

Inertial currents in isotropic plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasmas, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, MHD plasma. Solutions are developed by taking the MHD limit ot two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗

Inertial Currents in Isotropic Plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasma, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, MED plasma. Solutions are developed by taking the MHD limit of two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗

Designing reinforcement learning algorithms for building HVAC control: From experimental observation to simulation comparisons

Advanced supervisory-level control with reinforcement learning (RL) is regarded as a promising solution for HVAC systems to minimize energy consumption while maintaining thermal comfort and indoor air quality. However, most RL applications were conducted in the simulation environment rather than real-world HVAC systems. This paper developed a value-based RL controller termed Deep Q-Network (DQN) for a typical central HVAC system and evaluated its performance in a building test facility. By comparing DQN with a rule-based controller, the study not only demonstrated the cases where DQN could properly maintain indoor comfort but also discussed possible reasons why DQN failed in some other situations. Recognizing the limitations of value-based RL algorithms from the experimental tests, a simulation study was conducted to compare DQN with an alternative RL approach, an actor–critic algorithm termed Deep Deterministic Policy Gradient (DDPG). In scenarios with a relatively large action space, DDPG outperformed DQN by requiring fewer computational resources and achieving better thermal comfort, lower energy consumption, and more stable control actions. The findings suggest that the ability of DDPG to handle continuous control variables more effectively allows for faster convergence in training and more precise control in practice, which enhances the overall efficiency and reliability of the HVAC system.

Guo, Fangzhou↗

Coupling of dust acoustic and shear mode through velocity shear in a strongly coupled dusty plasma

In the strongly coupled limit, the generalized hydrodynamic model shows that a dusty plasma, acquiring significant rigidity, is able to support a “shear” like mode. It is being demonstrated here that in presence of velocity shear gradient, this shear like mode gets coupled with the dust acoustic mode which is generated by the compressibility effect of the dust fluid due to the finite temperatures of the dust, electron, and ion fluids. In the local analysis, the dispersion relation shows that velocity shear gradient not only couples the two modes but is also responsible for the instabilities of that coupled mode which is confirmed by nonlocal analysis with numerical techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of inhomogeneous surface disorder on the superheating field of superconducting RF cavities

Recent advances in surface treatments of Niobium superconducting radio frequency (SRF) cavities have led to substantially increased Q-factors and maximum surface field. This poses theoretical challenges to identify the mechanisms responsible for such performance enhancements. We report theoretical results for the effects of inhomogeneous surface disorder on the superheating field --- the surface magnetic field above which the Meissner state is globally unstable. We find that inhomogeneous disorder, such as that introduced by infusion of Nitrogen into the surface layers of Niobium SRF cavities, can increase the superheating field above the maximum for superconductors in the clean limit or with homogeneously distributed disorder. Homogeneous disorder increases the penetration of screening current, but also suppresses the maximum supercurrent. Inhomogeneous disorder in the form of an impurity diffusion layer biases this trade-off by increasing the penetration of the screening currents into cleaner regions with larger critical currents, thus limiting the suppression of the screening current to a thin dirty region close to the surface. Our results suggest that the impurity diffusion layers play a role in enhancing the maximum accelerating gradient of Nitrogen treated Niobium SRF cavities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗