Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sum reduction”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-T-08: A Beam Source Model for Monte Carlo Simulations of a Double-Scattering Proton Beam Delivery System Using Beam Current Modulation

Purpose: To reconstruct phase-space information upstream of patient specific collimators for Monte Carlo simulations using only radiotherapy planning system data. Methods: The proton energies are calculated based on residual ranges, e.g., sum of prescribed ranges in a patient and SSD. The Kapchinskij and Vladimirskij (KV) distribution was applied to sample proton’s x-y positions and momentum direction and the beam shape was assumed to be a circle. Free parameters, e.g., the initial energy spread and the emittance of KV distribution were estimated from the benchmarking with commissioning data in a commercial treatment planning system for an operational proton therapy center. The number of histories, which defines the height of individual pristine Bragg peaks (BP) of Spread-out Bragg peak (SOBP), are weighted based on beam current modulation and a correction factor is applied to take into account the fluence reduction as the residual range decreases due to the rotation of the range modulator wheel. The timedependent behaviors, e.g., the changes of the residual range and histories per a pristine BP, are realized by utilizing TOPAS (Tool for Particle Simulation). Results: Benchmarking simulations for selected SOBPs ranging 7.5 cm to 15.5 cm matched within 2 mm in range and up to 5 mm in SOBP width against measurement data in water phantom. We found this model tends to underestimate entrance dose by about 5 % in comparison to measurement. This was attributed to the situation that the energy distribution used in the model was limited in its granularity at the limit of single energy spectrum for the narrow angle modulator steps used in the proximal pull back region of the SOBPs. Conclusion: Within these limitations the source modeling method proved itself an acceptable alternative of a full treatment head simulation when the machine geometry and materials information are not available.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1982

Dipole-field sums and Lorentz factors for orthorhombic lattices, and implications for polarizable molecules

A method for computing the Lorentz tensor components in single crystals via rapidly convergent sums of Bessels functions is developed using the relationship between dipole-field sums and the tensor components. The Lorentz factors for simple, body-centered, and base-centered orthorhombic lattices are computed using this method, and the derivative Lorentz factors for simple orthorhombic lattices are also determined. Both the Lorentz factors and their derivatives are shown to be very sensitive to a lattice structure. The equivalent of the Clausius-Mossotti relation for general orthorhombic lattices is derived using the Lorentz-factor formalism, and the permanent molecular dipole moment is related to crystal polarization for the case of a ferroelectric of polarizable point dipoles. It is concluded that the polarization enhancement due to self-polarization familiar from classical theory may actually be a reduction in consequences of negative Lorentz factors in one or two lattice directions for noncubic crystals.

Purvis, C. K.↗
NASA NTRS1994

Hydrochloric acid and the chlorine budget of the lower statosphere

Concentrations of HCl measurements in the lower stratosphere in 1993 by the ALIAS instrument on the ER-2 aircraft reveal that only 40% of inorganic chlorine (Cl(y), inferred from in situ measurements of organic chlorinated sources gases) is present as HCl, significantly lower than model predictions. Although the sum of measured HCl, ClO and ClONO2, the latter inferred from measurements of ClO and NO2 equals Cl(y) to within the incertainty of measurement, it is systematically less than Cl(y) by 30-50%. This discrepancy suggests that concentrations of ClONO2 may exceed those of HC; near 20 km altitude, consistent with a slower photolysis rate for ClONO2 than calculated using recommended cross sections. Comparison of profiles of HCl measured during 1992 and 1193 at mid-latitudes by balloon (BLISS and MARKIV), space shuttle (ATMOS), and satellite (HALOE) instruments with the aircraft data reveal an apparent pressure dependence to the HCl to Cl(y) ratio, consistent with a factor of 3-10 reduction in the photolysis rate for ClONO2 at ER-2 altitudes. However, the diurnal variation of ClO is well-simulated by models using the recommended photolysis rate, and simulations measurements of ClONO2 and HCl at mid-latitudes by ATMOS and MARKIV report HCl (HCL+ ClONO2) ratios greater than or equal to 50%. Premliminary measurements by ALIAS in the southern hemisphere report HCl/Cl(y) values of about 75%.

Webster, C. R.↗
NASA NTRS1981

A method for obtaining reduced-order control laws for high-order systems using optimization techniques

A method of synthesizing reduced-order optimal feedback control laws for a high-order system is developed. A nonlinear programming algorithm is employed to search for the control law design variables that minimize a performance index defined by a weighted sum of mean-square steady-state responses and control inputs. An analogy with the linear quadractic Gaussian solution is utilized to select a set of design variables and their initial values. To improve the stability margins of the system, an input-noise adjustment procedure is used in the design algorithm. The method is applied to the synthesis of an active flutter-suppression control law for a wind tunnel model of an aeroelastic wing. The reduced-order controller is compared with the corresponding full-order controller and found to provide nearly optimal performance. The performance of the present method appeared to be superior to that of two other control law order-reduction methods. It is concluded that by using the present algorithm, nearly optimal low-order control laws with good stability margins can be synthesized.

Mukhopadhyay, V.↗
DOE OSTI2026

Selective Chemical Looping Combustion of Terminal Alkynes in Mixtures with Alkenes

The selective combustion of terminal alkynes in mixtures with alkenes is demonstrated during anaerobic reduction half-cycles on bulk bismuth oxide (Bi 2 O 3 ) as an approach to remove alkynes, which act as inhibitors in olefin polymerization. Bi 2 O 3 combusts phenylacetylene in styrene, 3-methylphenylacetylene in 3-methylstyrene, propyne in propylene, 1-hexyne in 1-hexene, and 1-octyne in 1-octene, with alkyne combustion selectivities exceeding 96%. Near unity reaction orders for hydrocarbon consumption during reduction half-cycles are consistent with combustion pathways initiated by rate-determining initial C–H activation, which drive selective alkyne combustion through intrinsic differences in the first-order rate constants for alkyne and alkene combustion rather than preferential adsorption of alkynes on Bi 2 O 3 surfaces. Computational assessments of initial C–H activation pathways for alkynes and alkenes on (010) α-Bi 2 O 3 surfaces using density functional theory illustrate that heterolytic transition states which form proton-carbanion pairs on Bi–O sites kinetically favor the activation of alkynes rather than alkenes due to differences in C–H bond acidity, and the barrier for heterolytic C–H activation is dictated in part by the sum of the molecular deprotonation energy and the energy to bind an R – carbanion to a Bi site in its transition-state geometry. Finally, these heterolytic reactivity channels during selective chemical looping combustion present novel routes for purifying olefin gas streams containing alkyne impurities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The combined effects of Mg 2+ and Sr 2+ incorporation during CaCO 3 precipitation and crystal growth

Calcite (CaCO 3 ) composition and properties are defined by the chemical environment in which CaCO 3 forms. However, a complete understanding of the relationship between aqueous chemistry during calcite precipitation and resulting chemical and physical CaCO 3 properties remains elusive; therefore, we present an investigation into the coupled effects of divalent cations Sr 2+ and Mg 2+ on CaCO 3 precipitation and subsequent crystal growth. Through chemical analysis of the aqueous phases and microscopy of the resulting calcite phases in compliment with density functional theory calculations, we elucidate the relationship between crystal growth and the resulting composition (elemental and isotopic) of calcite. The results of this experimental and modeling work suggest that Mg 2+ and Sr 2+ have cation-specific impacts that inhibit calcite crystal growth, including: (1) Sr 2+ incorporates more readily into calcite than Mg 2+ (D Sr > D Mg ), and increasing [Sr 2+ ] t or [Mg 2+ ] t increases D Sr ; (2) the inclusion of Mg 2+ into structure leads to a reduction in the calcite unit cell volume, whereas Sr 2+ leads to an expansion; (3) the inclusion of both Mg 2+ and Sr 2+ results in a distribution of unit cell impacts based on the relative positions of the Sr 2+ and Mg 2+ in the lattice. These experiments were conducted at saturation indices of CaCO 3 of ~4.1, favoring rapid precipitation. This rapid precipitation resulted in observed Sr isotope fractionation confirming Sr isotopic fractionation is dependent upon the precipitation rate. We further note that the precipitation and growth of calcite favors the incorporation of the lighter 86 Sr isotope over the heavier 87 Sr isotope, regardless of the initial solution conditions, and the degree of fractionation increases with D Sr . In sum, these results demonstrate the impact of solution environment to influence the incorporation behavior and crystal growth behavior of calcite. These factors are important to understand in order to effectively use geochemical signatures resulting from calcite precipitation or dissolution to gain specific information.

58 GEOSCIENCES↗
NASA NTRS2011

Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration

A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.

Scott, James R.↗
DOE OSTI2023

Higher-Order Corrections to Optimisers based on Newton's Method

The Newton, Gauss–Newton and Levenberg–Marquardt methods all use the first derivative of a vector function (the Jacobian) to minimise its sum of squares. When the Jacobian matrix is ill-conditioned, the function varies much faster in some directions than others and the space of possible improvement in sum of squares becomes a long narrow ellipsoid in the linear model. This means that even a small amount of nonlinearity in the problem parameters can cause a proposed point far down the long axis of the ellipsoid to fall outside of the actual curved valley of improved values, even though it is quite nearby. This paper presents a differential equation that ‘follows’ these valleys, based on the technique of geodesic acceleration, which itself provides a 2 nd order improvement to the Levenberg–Marquardt iteration step. Higher derivatives of this equation are computed that allow n th order improvements to the optimisation methods to be derived. These higher-order accelerated methods up to 4 th order are tested numerically and shown to provide substantial reduction of both number of steps and computation time.

43 PARTICLE ACCELERATORS↗
NASA NTRS2005

Series-Connected Buck Boost Regulators

A series-connected buck boost regulator (SCBBR) is an electronic circuit that bucks a power-supply voltage to a lower regulated value or boosts it to a higher regulated value. The concept of the SCBBR is a generalization of the concept of the SCBR, which was reported in "Series-Connected Boost Regulators" (LEW-15918), NASA Tech Briefs, Vol. 23, No. 7 (July 1997), page 42. Relative to prior DC-voltage-regulator concepts, the SCBBR concept can yield significant reductions in weight and increases in power-conversion efficiency in many applications in which input/output voltage ratios are relatively small and isolation is not required, as solar-array regulation or battery charging with DC-bus regulation. Usually, a DC voltage regulator is designed to include a DC-to-DC converter to reduce its power loss, size, and weight. Advances in components, increases in operating frequencies, and improved circuit topologies have led to continual increases in efficiency and/or decreases in the sizes and weights of DC voltage regulators. The primary source of inefficiency in the DC-to-DC converter portion of a voltage regulator is the conduction loss and, especially at high frequencies, the switching loss. Although improved components and topology can reduce the switching loss, the reduction is limited by the fact that the converter generally switches all the power being regulated. Like the SCBR concept, the SCBBR concept involves a circuit configuration in which only a fraction of the power is switched, so that the switching loss is reduced by an amount that is largely independent of the specific components and circuit topology used. In an SCBBR, the amount of power switched by the DC-to-DC converter is only the amount needed to make up the difference between the input and output bus voltage. The remaining majority of the power passes through the converter without being switched. The weight and power loss of a DC-to-DC converter are determined primarily by the amount of power processed. In the SCBBR, the unswitched majority of the power is passed through with very little power loss, and little if any increase in the sizes of the converter components is needed to enable the components to handle the unswitched power. As a result, the power-conversion efficiency of the regulator can be very high, as shown in the example of Figure 1. A basic SCBBR includes a DC-to-DC converter (see Figure 2). The switches and primary winding of a transformer in the converter is connected across the input bus, while the secondary winding and switches are connected in series with the output bus, so that the output voltage is the sum of the input voltage and the secondary voltage of the converter. In the breadboard SCBBR, the input voltage applied to the primary winding is switched by use of metal oxide/semiconductor field-effect transistors (MOSFETs) in a full bridge circuit; the secondary winding is center-tapped, with two MOSFET switches and diode rectifiers connected in opposed series in each leg. The sets of opposed switches and rectifiers are what enable operation in either a boost or a buck mode. In the boost mode, input voltage and current, and the output voltage and current are all positive; that is, the secondary voltage is added to the input voltage and the net output voltage can be regulated at a value equal or greater than the input voltage. In the buck mode, input voltage is still positive and the current still flows in the same direction in the secondary, but the switches are controlled such that some power flows from the secondary to the primary. The voltage across the secondary and the current into the primary are reversed. The result is that the output voltage is lower than the input voltage, and some power is recirculated from the converter secondary back to the input. Quantitatively, the advantage of an SCBBR is a direct function of the regulation range required. If, for example, a regulation range of 20 percent is required for a 500-W supply, th it suffices to design the DC-to-DC converter in the SCBBR for a power rating of only 100 W. The switching loss and size are much smaller than those of a conventional regulator that must be rated for switching of all 500 W. The reduction in size and the increase in efficiency are not directly proportional to switched-power ratio of 5:1 because the additional switches contribute some conduction loss and the input and output filters must be larger than those typically required for a 100-W converter. Nevertheless, the power loss and the size can be much smaller than those of a 500-W converter.

Birchenough, Arthur G.↗
NASA NTRS2012

Optimal Force Control of Vibro-Impact Systems for Autonomous Drilling Applications

The need to maintain optimal energy efficiency is critical during the drilling operations performed on future and current planetary rover missions (see figure). Specifically, this innovation seeks to solve the following problem. Given a spring-loaded percussive drill driven by a voice-coil motor, one needs to determine the optimal input voltage waveform (periodic function) and the optimal hammering period that minimizes the dissipated energy, while ensuring that the hammer-to-rock impacts are made with sufficient (user-defined) impact velocity (or impact energy). To solve this problem, it was first observed that when voice-coil-actuated percussive drills are driven at high power, it is of paramount importance to ensure that the electrical current of the device remains in phase with the velocity of the hammer. Otherwise, negative work is performed and the drill experiences a loss of performance (i.e., reduced impact energy) and an increase in Joule heating (i.e., reduction in energy efficiency). This observation has motivated many drilling products to incorporate the standard bang-bang control approach for driving their percussive drills. However, the bang-bang control approach is significantly less efficient than the optimal energy-efficient control approach solved herein. To obtain this solution, the standard tools of classical optimal control theory were applied. It is worth noting that these tools inherently require the solution of a two-point boundary value problem (TPBVP), i.e., a system of differential equations where half the equations have unknown boundary conditions. Typically, the TPBVP is impossible to solve analytically for high-dimensional dynamic systems. However, for the case of the spring-loaded vibro-impactor, this approach yields the exact optimal control solution as the sum of four analytic functions whose coefficients are determined using a simple, easy-to-implement algorithm. Once the optimal control waveform is determined, it can be used optimally in the context of both open-loop and closed-loop control modes (using standard realtime control hardware).

Aldrich, Jack B.↗
NASA NTRS2023

Investigating Hydroclimatic Impacts of the 168–158 BCE Volcanic Quartet and Their Relevance to the Nile River Basin and Egyptian History

The Ptolemaic era (305–30 BCE) is an important period of Ancient Egyptian history known for its material and scientific advances, but also intermittent political and social unrest in the form of (sometimes widespread) revolts against the Ptolemaic elites. While the role of environmental pressures has long been overlooked in this period of Egyptian history, ice-core-based volcanic histories have identified the period as experiencing multiple notable eruptions, and a repeated temporal association between explosive volcanism and revolt has recently been noted. Here we analyze the global and regional (Nile River basin) hydroclimatic response to a unique historical sequence of four large and closely timed volcanic eruptions (first a tropical one, followed by three extratropical northern hemispheric events) between 168 and 158 BCE, a particularly troubled period in Ptolemaic history for which we now provide a more detailed hydroclimatic context. The NASA (National Aeronautics and Space Administration) GISS (Goddard Institute for Space Studies) ModelE2.1 Earth system model simulates a strong radiative response with a radiative forcing (top of atmosphere) of −7.5 W m−2 (following the first eruption) and −2.5 W m−2 (after each of the three remaining eruptions) at a global scale. Associated with this, we observe a global surface cooling of the order of 1.5 ∘C following the first (tropical) eruption, with the following three extratropical eruptions extending the cooling period for more than 15 years. Consequently, this series of eruptions is observed to constrain the northward migration of the inter-tropical convergence zone (ITCZ) during the Northern Hemisphere summer monsoon season, and major monsoon zones (African, South Asian, and East Asian) were seen to experience a suppression of rainfall of >1 mm d−1 during the monsoon (JJAS) season averaged for 2 years after each eruption. A substantial suppression of the Indian and North African summer monsoon (over the Nile River headwater region) was seen to strongly affect the modeled river flow in the catchment and discharge at river mouth. River mass flow over the basin was observed to decrease by 29 % and 38 % relative to an unperturbed (non-volcanic) annual mean flow in the first and second year, respectively, after the first (i.e., tropical) eruption. A moderate decrease ranging between 5 % and 18 % was observed after the third and fourth (extratropical) eruptions. These results indicate, in sum, that the first eruption likely produced a strong hydroclimate response, with the following extratropical eruptions prolonging this. These results also support the recently hypothesized association between ice-core-based signals of explosive volcanism and hydroclimatic variability during the Ptolemaic era, including the suppression of the agriculturally critical Nile summer flooding.

hydroclimate impacts↗
DOE OSTI2018

Cold rolling texture evolution of TiB2 particle reinforced Al-based composites by Neutron Diffraction and EBSD analysis

Highlights: • Texture development of in-situ TiB{sub 2} particle reinforced Al matrix composites after cold rolling were characterized. • The sum of rolling texture volume fractions has been significantly reduced by the weak texture in particle deformation zone. • Recrystallization nuclei with the Cube and Rotated-Cube orientations were observed in the particle deformation zones. • Small TiB{sub 2} particles within Al grains kept the parallel relationship with surrounding matrix during the deformation. - Abstract: Texture development of particle-reinforced metal matrix materials (P-MMCs) has not been understood well, especially the effects of second-phase particles during the deformation. In this work, we have investigated the texture evolution of TiB{sub 2}/Al composites during cold rolling by Neutron Diffraction and EBSD analysis. Two forms of TiB{sub 2} particles influenced the texture evolution in different ways. The particle-clusters aggregated along grain boundaries led to a reduction of texture volume fraction; recrystallization nuclei with weak texture appeared around these particle-clusters. Small non-shearable TiB{sub 2} particles dispersed inside Al matrix contributed to strong Copper texture component at high strains; the individual particles rotated together with Al matrix during the deformation. These revealed effects of particles on the texture development in TiB{sub 2}/Al composites provide a reference for the texture control of P-MMCs.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Optimizing and Diversifying Electric Vehicle Driving Range for U.S. Drivers

Properly determining the driving range is critical for accurately predicting the sales and social benefits of battery electric vehicles (BEVs). This study proposes a framework for optimizing the driving range by minimizing the sum of battery price, electricity cost, and range limitation cost referred to as the "range-related cost" as a measurement of range anxiety. The objective function is linked to policy-relevant parameters, including battery cost and price markup, battery utilization, charging infrastructure availability, vehicle efficiency, electricity and gasoline prices, household vehicle ownership, daily driving patterns, discount rate, and perceived vehicle lifetime. Qualitative discussion of the framework and its empirical application to a sample (N=36664) representing new car drivers in the United States is included. The quantitative results strongly suggest that ranges of less than 100 miles are likely to be more popular in the BEV market for a long period of time. The average optimal range among U.S. drivers is found to be largely inelastic. Still, battery cost reduction significantly drives BEV demand toward longer ranges, whereas improvement in the charging infrastructure is found to significantly drive BEV demand toward shorter ranges. In conclusion, the bias of a single-range assumption and the effects of range optimization and diversification in reducing such biases are both found to be significant.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2014

Toward the development of intrafraction tumor deformation tracking using a dynamic multi-leaf collimator

Purpose: Intrafraction deformation limits targeting accuracy in radiotherapy. Studies show tumor deformation of over 10 mm for both single tumor deformation and system deformation (due to differential motion between primary tumors and involved lymph nodes). Such deformation cannot be adapted to with current radiotherapy methods. The objective of this study was to develop and experimentally investigate the ability of a dynamic multi-leaf collimator (DMLC) tracking system to account for tumor deformation. Methods: To compensate for tumor deformation, the DMLC tracking strategy is to warp the planned beam aperture directly to conform to the new tumor shape based on real time tumor deformation input. Two deformable phantoms that correspond to a single tumor and a tumor system were developed. The planar deformations derived from the phantom images in beam's eye view were used to guide the aperture warping. An in-house deformable image registration software was developed to automatically trigger the registration once new target image was acquired and send the computed deformation to the DMLC tracking software. Because the registration speed is not fast enough to implement the experiment in real-time manner, the phantom deformation only proceeded to the next position until registration of the current deformation position was completed. The deformation tracking accuracy was evaluated by a geometric target coverage metric defined as the sum of the area incorrectly outside and inside the ideal aperture. The individual contributions from the deformable registration algorithm and the finite leaf width to the tracking uncertainty were analyzed. Clinical proof-of-principle experiment of deformation tracking using previously acquired MR images of a lung cancer patient was implemented to represent the MRI-Linac environment. Intensity-modulated radiation therapy (IMRT) treatment delivered with enabled deformation tracking was simulated and demonstrated. Results: The first experimental investigation of adapting to tumor deformation has been performed using simple deformable phantoms. For the single tumor deformation, the A{sub u}+A{sub o} was reduced over 56% when deformation was larger than 2 mm. Overall, the total improvement was 82%. For the tumor system deformation, the A{sub u}+A{sub o} reductions were all above 75% and the total A{sub u}+A{sub o} improvement was 86%. Similar coverage improvement was also found in simulating deformation tracking during IMRT delivery. The deformable image registration algorithm was identified as the dominant contributor to the tracking error rather than the finite leaf width. The discrepancy between the warped beam shape and the ideal beam shape due to the deformable registration was observed to be partially compensated during leaf fitting due to the finite leaf width. The clinical proof-of-principle experiment demonstrated the feasibility of intrafraction deformable tracking for clinical scenarios. Conclusions: For the first time, we developed and demonstrated an experimental system that is capable of adapting the MLC aperture to account for tumor deformation. This work provides a potentially widely available management method to effectively account for intrafractional tumor deformation. This proof-of-principle study is the first experimental step toward the development of an image-guided radiotherapy system to treat deforming tumors in real-time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Nonadiabatic excited-state molecular dynamics: On-the-fly limiting of essential excited states

The simulation of nonadiabatic dynamics in extended molecular systems involving hundreds of atoms and large densities of states is particularly challenging. Nonadiabatic coupling terms (NACTs) represent a significant numerical bottleneck in surface hopping approaches. Rather than using unreliable NACT cutting schemes, here we develop “on-the-fly” state limiting methods to eliminate states that are no longer essential for the non-radiative relaxation dynamics as a trajectory proceeds. Here, we propose a state number criteria and an energy-based state limit. The latter is more physically relevant by requiring a user-imposed energy threshold. For this purpose, we introduce a local kinetic energy gauge by summing contributions from atoms within the spatial localization of the electronic wavefunction to define the energy available for upward hops. The proposed state limiting schemes are implemented within the nonadiabatic excited-state molecular dynamics framework to simulate photoinduced relaxation in poly-phenylene vinylene (PPV) and branched poly-phenylene ethynylene (PPE) oligomers for benchmark evaluation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Antigen‐binding fragments with improved crystal lattice packing and enhanced conformational flexibility at the elbow region as crystallization chaperones

It has been shown previously that a set of three modifications—termed S1, Crystal Kappa, and elbow—act synergistically to improve the crystallizability of an antigen-binding fragment (Fab) framework. Here, we prepared a phage-displayed library and performed crystallization screenings to identify additional substitutions—located near the heavy-chain elbow region—which cooperate with the S1, Crystal Kappa, and elbow modifications to increase expression and improve crystallizability of the Fab framework even further. One substitution (K141Q) supports the signature Crystal Kappa-mediated Fab:Fab crystal lattice packing interaction. Another substitution (E172G) improves the compatibility of the elbow modification with the Fab framework by alleviating some of the strain incurred by the shortened and bulkier elbow linker region. A third substitution (F170W) generates a split-Fab conformation, resulting in a powerful crystal lattice packing interaction comprising the biological interaction interface between the variable heavy and light chain domains. In sum, we have used K141Q, E172G, and F170W substitutions—which complement the S1, Crystal Kappa, and elbow modifications—to generate a set of highly crystallizable Fab frameworks that can be used as chaperones to enable facile elucidation of Fab:antigen complex structures by x-ray crystallography.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Instability Analysis of a Low-Density Gas Jet Injected into a High-Density Gas

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas were performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. It was found that the presence of variable density within the shear layer resulted in an increase in the temporal amplification rate of the disturbances and an increase in the range of unstable frequencies, accompanied by a reduction in the phase velocities of the disturbances. Also, the temporal growth rates of the disturbances were increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity. The spatio-temporal stability analysis was performed to determine the nature of the absolute instability of the jet. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the jet s absolute instability were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be. Experiments were carried out to observe the qualitative differences between a round low-density gas jet injected into a high-density gas (helium jet injected into air) and a round constant density jet (air jet injected into air). Flow visualizations and velocity measurements in the near-injector region of the helium jet show more mixing and spreading of the helium jet than the air jet. The vortex structures develop and contribute to the jet spreading causing the helium jet to oscillate.

Lawson, Anthony Layiwola↗
NASA NTRS1993

Interannual Sea Level Variations in the Tropical Indian Ocean from Geosat and Shallow Water Simulations

Sea level variations of the Indian Ocean north of 20 deg S are analyzed from Geosat satellite altimeter data over April 1985-September 1989. These variations are compared and interpreted with numerical simulations derived from a reduced gravity model forced by FSU observed winds over the same period. After decomposition into complex empirical orthogonal functions, the low-frequency anomalies are described by the first two modes for observations as well as for simulations. The sums of the two modes contain 34% and 40% of the observed and simulated variances, respectively. Averaged over the basin, the observed and simulated sea level changes are correlated by 0.92 over 1985-1988. The strongest change happens during the El Ninio 1986-1987: between winter 1986 and summer 1987 the basin-averaged sea level rises by approx. 1 cm. These low-frequency variations can partly be explained by changes in the Sverdrup circulation. The southern tropical Indian Ocean between 1O deg and 20 deg S is the domain where those changes are strongest: the averaged sea level rises by approx. 4.5 cm between winter 1986 and winter 1987. There, the signal propagates southwestward across the basin at a speed similar to free Rossby waves. Sensitivity of observed anomalies is examined over 1987-1988, with different orbit ephemeris, tropospheric corrections, and error reduction processes. The uncertainty of the basin-averaged sea level estimates is mostly due to the way the orbit error is reduced and reaches approx. 1 cm. Nonetheless, spatial correlation is good between the various observations and better than between observations and simulations. Sensitivity of simulated anomalies to the wind uncertainty, examined with Former Soviet Union (FSU) and European Center for Medium-Range Weather Forecasting (ECMWF) forcings over 1985-1988, shows that the variance of the simulations driven by ECMWF is 52% smaller, as FSU winds are stronger than ECMWF. Results show that the wind strength also affects the dynamic response of the ocean: anomalies propagate westward across the basin more than twice as fast with FSU than with ECMWF. It is found that the discrepancy is larger between ECMWF and FSU simulations than between observations and FSU simulations.

Perigaud, Claire↗