Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully linear model”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Static stability of a three-dimensional space truss

In order to deploy large flexible space structures it is necessary to develop support systems that are strong and lightweight. The most recent example of this aerospace design need is vividly evident in the space station solar array assembly. In order to accommodate both weight limitations and strength performance criteria, ABLE Engineering has developed the Folding Articulating Square Truss (FASTMast) support structure. The FASTMast is a space truss/mechanism hybrid that can provide system support while adhering to stringent packaging demands. However, due to its slender nature and anticipated loading, stability characterization is a critical part of the design process. Furthermore, the dire consequences surely to result from a catastrophic instability quickly provide the motivation for careful examination of this problem. The fundamental components of the space station solar array system are the (1) solar array blanket system, (2) FASTMast support structure, and (3) mast canister assembly. The FASTMast once fully deployed from the canister will provide support to the solar array blankets. A unique feature of this structure is that the system responds linearly within a certain range of operating loads and nonlinearly when that range is exceeded. The source of nonlinear behavior in this case is due to a changing stiffness state resulting from an inability of diagonal members to resist applied loads. The principal objective of this study was to establish the failure modes involving instability of the FASTMast structure. Also of great interest during this effort was to establish a reliable analytical approach capable of effectively predicting critical values at which the mast becomes unstable. Due to the dual nature of structural response inherent to this problem, both linear and nonlinear analyses are required to characterize the mast in terms of stability. The approach employed herein is one that can be considered systematic in nature. The analysis begins with one and two-dimensional failure models of the system and its important components. From knowledge gained through preliminary analyses a foundation is developed for three-dimensional analyses of the FASTMast structure. The three-dimensional finite element (FE) analysis presented here involves a FASTMast system one-tenth the size of the actual flight unit. Although this study does not yield failure analysis results that apply directly to the flight article, it does establish a method by which the full-scale mast can be evaluated.

Shaker, John F.↗
DOE OSTI2026

Ab Initio Polariton Spectra of ZnTPP Molecules Collectively Coupled Inside an Optical Cavity

Exciton-polaritons are quasi-particles formed by the quantum mechanical hybridization of electronic and photonic excitations. Despite extensive investigations, a fundamental understanding of molecular polariton spectra and the polariton delocalization from an ab initio theoretical perspective remains elusive. We simulate experimentally measured linear transmission spectroscopy of many Zinc(II) tetraphenylporphyrin (ZnTPP) molecules collectively coupled to a cavity from first principles. Our theoretical approach incorporates many low-lying electronic excitations in ZnTPP molecules, as well as collective light-matter couplings between ZnTPP and the quantized radiation modes, both of which are shown to be the key to accurately recovering the experimental spectra. We further analyzed to what extent the polariton and dark states are delocalized over many molecules, for the first time, using fully ab initio descriptions of the molecules. We finally investigate the line width as a function of detuning, providing new theoretical insights into the experimentally observed motional narrowing behavior. Our work presents first-ofits- kind theoretical studies on molecular polariton spectra, offering a new perspective on molecular polariton formation in realistic ab initio molecular systems whose rich, many-state nature provides spectral features enabled by the high density of electronic states beyond simple quantum optics models.

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

Low- n global ideal MHD instabilities in the CFETR baseline scenario

This article reports an evaluation on the linear ideal magnetohydrodynamic (MHD) stability of the China Fusion Engineering Test Reactor (CFETR) baseline scenario for various first-wall locations. The initial-value code NIMROD and eigen-value code AEGIS are employed in this analysis. Despite the distinctly different approaches in modeling the scrape off layer region, the dominant growth in each of the low-n (n=1–10) modes are consistent between the two codes. We report the higher-n modes are dominated by ballooning modes and localized in the pedestal region, while the lower-n modes have more prominent external kink components and broader mode profiles. The influence of the plasma-vacuum profile and wall shape are also examined using NIMROD. In the presence of a resistive wall, the low-n ideal MHD instabilities are further studied using AEGIS. For the designed first-wall location, the n=1 resistive wall mode is found to be unstable, which could be fully stabilized by uniform toroidal rotation above the 2.9% core Alfven speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Optimum Planning for Inductively Charged On-Demand Automated Electric Shuttles at Greenville, South Carolina

Wireless charging technology presents perfect fits for autonomous electric vehicles for realizing a fully automated system (vehicle and charger). This paper presents a planning optimization analysis for a pilot project of in-route wireless charging infrastructure serving fixed-route on-demand shared automated electric shuttles (SAESs) at Greenville, South Carolina, USA. A single-objective non-linear integer planning optimization problem is formulated. A comprehensive cost function representing the inductively charged SAESs is developed, considering road construction, power electronics and materials, traction battery, and installation costs. The optimization problem is solved to determine the best combination of the system key design parameters (number and allocations of wireless chargers, charging power level, track length and on-board battery capacity) that show the most cost-effective solution and allow the SAESs realizing charge sustaining operation. The planning platform incorporates representative simulated traffic data (drive cycles and routes) for four SAESs at Greenville project using the Simulation of Urban Mobility (SUMO) tool. These data are fed to a vehicle powertrain model and a wireless charger power model to predict the battery power, energy and state-of-charge (SOC) profiles, which are provided to the search algorithm to assess the design objectives under specific constraints. The results indicated that implementing wireless charging at a few designated stops for fixed-route SAESs with proper design allows the vehicles realizing charge sustaining operation, infinite range and zero recharge downtime, with a significant reduction in the on-board battery (36 % ) and road coverage (69%), at minimum cost.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Less Common Topological Phenomena in Bulk Materials

"Prediction and subsequent discovery of topological insulators is considered to be one of the main results in condensed matter physics in the last decades. Not surprisingly, it has received major attention of both researchers and funding agencies. This attention is well-deserved; yet, one cannot but note that nearly all this research revolves about essentially the same concept: electronic excitations with linear dispersion, covering, of course, such diverse and intriguing phenomena as topological insulator, bulk Dirac states (or Weyl, if not spin degenerate), Mayorana fermions. In this project, we will address, mainly, other topological phenomena, such as topologically nontrivial magnetic patterns (as, for instance, topological Hall and related phenomena). Specifically, we propose three interrelated trusts: (1) Time-reversal symmetry breaking nonrelativistic antiferromagnets, called altermagnets. These are materials that break Kramers degeneracy of electronic bands, despite having zero net magnetization by symmetry and being fully collinear, and not necessarily non-centrosymmetric. The corresponding band structure is very similar to the band structure in non-centrosymmetric spin-orbital materials, but materials that we propose to study are distinctly different, first and foremost in the sense that despite sharing many aspects of their electronic properties with the latter, they break the time-reversal symmetry without either spin-orbit coupling or lack of inversion symmetry. (2) Topologically nontrivial magnetic spirals. The PI has been engaged with the experimental group of Dr. Ghimire at GMU investigating Dirac materials with helical magnetism, based on stacked magnetic Kagome layers, with a generic formula of RMn6Sn6. In particular, Y Mn6Sn6 demonstrates a component of the Hall effect that is naturally interpreted in terms of a topological spin texture, as well as linear magnetoresistance. Our calculation identify Dirac states that are robust with respect to the spiral formation, and let us derive an advance mean-field model explaining the observed phase diagram. This model predicted four distinct phases, with very distinct properties, which have now been seen in neutron experiments. The same compound is known to demonstrate topological Hall effect in a particular magnetic phase, and only at elevated temperature. Based on our understanding of the phase diagram, we have worked out a phenomenological theory of a chiral (skyrmionic) response to an external magnetic field, similar to the nematic response to external strain in Fe-based superconductors, which is possible in a centrosymmetric lattice and without interplanar Dzyaloshinskii-Moriya interaction. This phenomenological theory agrees quantitatively with the experiment. It is in our plans to research other similar materials for this effect. (3) Search for 3D analogues of Fe-based superconductors. We want to investigate materials that can be viewed as 3D analogues of FeSe. Specifically, we want materials that are good metals and host antiferromagnetism, which can be suppressed by pressure and generate an s-wave superconductivity, as in Fe-based superconductors. We have in mind some candidates already. This work will proceed in close collaboration with the experimental group of Prof. Nirmal Ghimire in the same department, whose expertise lies in sample making, magnetometry and transport measurements of materials with complex magnetic structures."

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1979

Marangoni bubble motion in zero gravity

It was shown experimentally that the Marangoni phenomenon is a primary mechanism for the movement of a gas bubble in a nonisothermal liquid in a low gravity environment. A mathematical model consisting of the Navier-Stokes and thermal energy equations, together with the appropriate boundary conditions for both media, is presented. Parameter perturbation theory is used to solve this boundary value problem; the expansion parameter is the Marangoni number. The zeroth, first, and second order approximations for the velocity, temperature and pressure distributions in the liquid and in the bubble, and the deformation and terminal velocity of the bubble are determined. Experimental zero gravity data for a nitrogen bubble in ethylene glycol, ethanol, and silicone oil subjected to a linear temperature gradient were obtained using the NASA Lewis zero gravity drop tower. Comparison of the zeroth order analytical results for the bubble terminal velocity showed good agreement with the experimental measurements. The first and second order solutions for the bubble deformation and bubble terminal velocity are valid for liquids having Prandtl numbers on the order of one, but there is a lack of appropriate data to test the theory fully.

Thompson, R. L.↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2017

Initial CEED API (CEED-MS10)

Key components of CEED software development involve fast finite element operator storage and evaluation, architecture optimizations, performant algorithms for all orders, global kernels for finite element operators, efficient use of the memory sub-system, optimal data locality and motion, enhanced scalability and parallelism, and fast tensor contractions. As ultimate goals of the project, CEED is exploring and identifying the best algorithms for the full range of discretizations and applying algorithmic and software development to support ECP applications’ needs. In this milestone, we initiated the development and implementation of a common CEED API consisting of two main components: low-level and high-level APIs. We accomplished the goal of having the low-level component ready for the release of CEED v1.0 in the next quarter, see CEED-MS13. In this report we also provide a highlight from CEED’s integration with the MARBL app and the corresponding Laghos miniapp. A main theme in this work has been the systematic transition of advanced algorithms, developed and tested in Laghos, to the large-scale multi-physics settings of the MARBL code. The low-level CEED API provides a set of FE kernels and components for writing new low-level kernels. Examples include: vector and sparse linear algebra, element matrix assembly over a batch of elements, partial assembly and action for efficient high-order operators like mass, diffusion, advection, etc. The main goal of the low-level API is to establish the basis for the high-level API. Also, identifying such low-level kernels and providing a reference implementation for them serves as the basis for collaboration with vendors and STs. Another goal for the low-level API is to make it easy for applications with their own discretization infrastructure, who do not want to switch directly to using the high-level CEED API, to still evaluate and use the core operations provided by the low-level API. One of the challenges with high-order methods is that a global sparse matrix is no longer a good representation of a high-order linear operator, both with respect to the FLOPs needed for its evaluation, as well as the memory transfer needed for a matvec. Thus, another important question we tried to address while developing the low-level API is: is it possible to abstract out in a natural and efficient way the different implementations and data structures that arise when dealing with various computational device types: CPUs, GPUs, etc. Fully answering this question is beyond the scope of this particular milestone, since that requires longer time as we develop and expand the devices we support. The future high-level API will be the frontend provided by CEED consisting of high-level abstractions such as meshes, finite element spaces and solutions, bilinear forms, solvers, etc. At this high-level, an important aspect of the API will be to provide a uniform abstraction layer for writing code in a portable manner. To achieve that, we plan to utilize the MPI programming model combined with a device abstraction layer (DAL) for finite element computations that we will develop. The DAL will serve as the bridge between the low- and the high-level APIs, allowing a variety of programming models to be used for the low-level kernel implementations. The software artifacts delivered as part of this milestone include an initial CEED API software implementation (libCEED) provided through the CEED website, http://ceed.exascaleproject.org and the CEED GitHub organization, http://github.com/ceed. In addition to details and results from the above R&D efforts, in this document we are also reporting on other project-wide activities performed in Q1 of FY18, including: the Nek5000 v17.0 MFEM v3.3.2 and PETSc v3.8 releases, collaborations with ECP/ST and SciDAC projects, the inaugural Nek5000 hackathon and various outreach activities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Method for modeling the gradual physical degradation of a porous material

Cementitious and other engineered porous materials encountered in waste disposals may degrade over time due to one or more mechanisms. Physical degradation may take the form of cracking (fracturing) and/or altered (e.g. increased) porosity, depending on the material and underlying degradation mechanism. In most cases, the hydraulic properties of degrading materials are expected to evolve due to physical changes occurring over roughly the pore to decimeter scale, which is conducive to calculating equivalent or effective material properties. The exact morphology of a degrading material in its end-state may or may not be known. In the latter case, the fully-degraded condition can be assumed to be similar to a more-permeable material in the surrounding environment, such as backfill soil. Then the fully-degraded waste form or barrier material is hydraulically neutral with respect to its surroundings, constituting neither a barrier to nor conduit for moisture flow and solute transport. Unless the degradation mechanism is abrupt, a gradual transition between the intact initial and fully-degraded final states is desired. Linear interpolation through time is one method for smoothly blending hydraulic properties between those of an intact matrix and those of a soil or other surrogate for the end-state.

36 MATERIALS SCIENCE↗
NASA NTRS1994

High speed civil transport: Sonic boom softening and aerodynamic optimization

An improvement in sonic boom extrapolation techniques has been the desire of aerospace designers for years. This is because the linear acoustic theory developed in the 60's is incapable of predicting the nonlinear phenomenon of shock wave propagation. On the other hand, CFD techniques are too computationally expensive to employ on sonic boom problems. Therefore, this research focused on the development of a fast and accurate sonic boom extrapolation method that solves the Euler equations for axisymmetric flow. This new technique has brought the sonic boom extrapolation techniques up to the standards of the 90's. Parallel computing is a fast growing subject in the field of computer science because of its promising speed. A new optimizer (IIOWA) for the parallel computing environment has been developed and tested for aerodynamic drag minimization. This is a promising method for CFD optimization making use of the computational resources of workstations, which unlike supercomputers can spend most of their time idle. Finally, the OAW concept is attractive because of its overall theoretical performance. In order to fully understand the concept, a wind-tunnel model was built and is currently being tested at NASA Ames Research Center. The CFD calculations performed under this cooperative agreement helped to identify the problem of the flow separation, and also aided the design by optimizing the wing deflection for roll trim.

Cheung, Samson↗
NASA NTRS2014

A Self-Stabilizing Hybrid-Fault Tolerant Synchronization Protocol

In this report we present a strategy for solving the Byzantine general problem for self-stabilizing a fully connected network from an arbitrary state and in the presence of any number of faults with various severities including any number of arbitrary (Byzantine) faulty nodes. Our solution applies to realizable systems, while allowing for differences in the network elements, provided that the number of arbitrary faults is not more than a third of the network size. The only constraint on the behavior of a node is that the interactions with other nodes are restricted to defined links and interfaces. Our solution does not rely on assumptions about the initial state of the system and no central clock nor centrally generated signal, pulse, or message is used. Nodes are anonymous, i.e., they do not have unique identities. We also present a mechanical verification of a proposed protocol. A bounded model of the protocol is verified using the Symbolic Model Verifier (SMV). The model checking effort is focused on verifying correctness of the bounded model of the protocol as well as confirming claims of determinism and linear convergence with respect to the self-stabilization period. We believe that our proposed solution solves the general case of the clock synchronization problem.

Malekpour, Mahyar R.↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
NASA NTRS2025

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗
DOE OSTI2014

Advances in the Understanding of ELM Suppression by Resonant Magnetic Perturbations (RMPs) in DIII-D and Implications for ITER

Experiments on DIII-D have expanding the operating window for RMP ELM suppression to higher q95 with dominant electron heating and fully non-inductive current drive relevant to advanced modes of ITER operation. Robust ELM suppression has also been obtained with a reduced coil set, mitigating the risk of coil failure in maintaining ELM suppression in ITER. These results significantly expand the operating space and reduce risk for obtaining RMP ELM suppression in ITER. Efforts have also been made to search for 3D cause of ELM suppression. No internal non-axisymmetric structure is detected at the top of the pedestal, indicating that the dominant effect of the RMP is to produce an n=0 transport modification of the profiles. Linear two fluid MHD simulations using M3D-C1 indicate resonant field penetration and significant magnetic stochasticity at the top of the pedestal, consistent with the absence of detectable 3D structure in that region. A profile database was developed to compare the scaling of the pedestal and global confinement with the applied 3D field strength in ELM suppressed and ELM mitigated plasmas. The EPED pedestal model accurately predicts the measured pedestal pressure at the threshold of ELM suppression, increasing confidence in theoretical projections to ITER pedestal conditions. Both the H-factor (H(sub)98y2) and thermal energy confinement time do not degrade substantially with applied RMP fields near the threshold of ELM suppression, enhancing confidence in the compatibility of ITER high performance operation with RMP ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Reactive Landing of Gramicidin S and Ubiquitin Ions onto Activated Self-Assembled Monolayer Surfaces

Using mass-selected ion deposition combined with in situ infrared reflection absorption spectroscopy (IRRAS), we examined the reactive landing of gramicidin S and ubiquitin ions onto activated self-assembled monolayer (SAM) surfaces terminated with N-hydroxysuccinimidyl ester (NHS-SAM) and acyl fluoride (COF-SAM) groups. Doubly protonated gramicidin S, [GS+2H]2+, and two charge states of ubiquitin, [U+5H]5+ and [U+13H]13+, were used as model systems, allowing us to explore the effect of the number of free amino groups and the secondary structure on the efficiency of covalent bond formation between the projectile ion and the surface. For all projectile ions, ion deposition resulted in the depletion of IRRAS bands corresponding to the terminal groups on the SAM and the appearance of several new bands not associated with the deposited species. These new bands were assigned to the C=O stretching vibrations of COOH and COO- groups formed on the surface as a result of ion deposition. The presence of these bands was attributed to an alternative reactive landing pathway that competes with covalent bond formation. This pathway with similar yields for both gramicidin S and ubiquitin ions is analogous to the hydrolysis of the NHS ester bond in solution. The covalent bond formation efficiency increased linearly with the number of free amino groups and was found to be lower for the more compact conformation of ubiquitin compared with the fully unfolded conformation. This observation was attributed to the limited availability of amino groups on the surface of the folded conformation. Our results have provided new insights on the efficiency and mechanism of reactive landing of peptides and proteins onto activated SAMs

Environmental Molecular Sciences Laboratory↗
DOE OSTI2016

Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach

We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a number of existing methods, we discretize rigid bodies using a collection of minimally resolved spherical blobs constrained to move as a rigid body, to arrive at a potentially large linear system of equations for the unknown Lagrange multipliers and rigid-body motions. Here we develop a block-diagonal preconditioner for this linear system and show that a standard Krylov solver converges in a modest number of iterations that is essentially independent of the number of particles. Key to the efficiency of the method is a technique for fast computation of the product of the blob-blob mobility matrix and a vector. For unbounded suspensions, we rely on existing analytical expressions for the Rotne-Prager-Yamakawa tensor combined with a fast multipole method (FMM) to obtain linear scaling in the number of particles. For suspensions sedimented against a single no-slip boundary, we use a direct summation on a graphical processing unit (GPU), which gives quadratic asymptotic scaling with the number of particles. For fully confined domains, such as periodic suspensions or suspensions confined in slit and square channels, we extend a recently developed rigid-body immersed boundary method by B. Kallemov, A. P. S. Bhalla, B. E. Griffith, and A. Donev (Commun. Appl. Math. Comput. Sci. 11 (2016), no. 1, 79-141) to suspensions of freely moving passive or active rigid particles at zero Reynolds number. We demonstrate that the iterative solver for the coupled fluid and rigid-body equations converges in a bounded number of iterations regardless of the system size. In our approach, each iteration only requires a few cycles of a geometric multigrid solver for the Poisson equation, and an application of the block-diagonal preconditioner, leading to linear scaling with the number of particles. We optimize a number of parameters in the iterative solvers and apply our method to a variety of benchmark problems to carefully assess the accuracy of the rigid multiblob approach as a function of the resolution. We also model the dynamics of colloidal particles studied in recent experiments, such as passive boomerangs in a slit channel, as well as a pair of non-Brownian active nanorods sedimented against a wall.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2003

Route Monopolie and Optimal Nonlinear Pricing

To cope with air traffic growth and congested airports, two solutions are apparent on the supply side: 1) use larger aircraft in the hub and spoke system; or 2) develop new routes through secondary airports. An enlarged route system through secondary airports may increase the proportion of route monopolies in the air transport market.The monopoly optimal non linear pricing policy is well known in the case of one dimension (one instrument, one characteristic) but not in the case of several dimensions. This paper explores the robustness of the one dimensional screening model with respect to increasing the number of instruments and the number of characteristics. The objective of this paper is then to link and fill the gap in both literatures. One of the merits of the screening model has been to show that a great varieD" of economic questions (non linear pricing, product line choice, auction design, income taxation, regulation...) could be handled within the same framework.VCe study a case of non linear pricing (2 instruments (2 routes on which the airline pro_ddes customers with services), 2 characteristics (demand of services on these routes) and two values per characteristic (low and high demand of services on these routes)) and we show that none of the conclusions of the one dimensional analysis remain valid. In particular, upward incentive compatibility constraint may be binding at the optimum. As a consequence, they may be distortion at the top of the distribution. In addition to this, we show that the optimal solution often requires a kind of form of bundling, we explain explicitly distortions and show that it is sometimes optimal for the monopolist to only produce one good (instead of two) or to exclude some buyers from the market. Actually, this means that the monopolist cannot fully apply his monopoly power and is better off selling both goods independently.We then define all the possible solutions in the case of a quadratic cost function for a uniform distribution of agent types and explain the implications for airlines in terms of service differentiation.

Tournut, Jacques↗
DOE OSTI2016

Predicting the 5-Year Risk of Biochemical Relapse After Postprostatectomy Radiation Therapy in ≥PT2, pN0 Patients With a Comprehensive Tumor Control Probability Model

Purpose: To fit the individual biochemical recurrence-free survival (bRFS) data from patients treated with postprostatectomy radiation therapy (RT) with a comprehensive tumor control probability (TCP) model. Methods and Materials: Considering pre-RT prostate-specific antigen (PSA) as a surrogate of the number of clonogens, bRFS may be expressed as a function of dose-per-fraction–dependent radiosensitivity (α{sub eff}), the number of clonogens for pre-RT PSA = 1 ng/mL (C), and the fraction of patients who relapse because of clonogens outside the treated volume (K), assumed to depend (linearly or exponentially) on pre-RT PSA and Gleason score (GS). Data from 894 node-negative, ≥pT2, pN0 hormone-naive patients treated with adjuvant (n=331) or salvage (n=563) intent were available: 5-year bRFS data were fitted grouping patients according to GS (<7:392, =7:383, >7:119). Results: The median follow-up time, pre-RT PSA, and dose were 72 months, 0.25 ng/mL, and 66.6 Gy (range 59.4-77.4 Gy), respectively. The best-fit values were 0.23 to 0.26 Gy{sup −1} and 10{sup 7} for α{sub eff} and C for the model considering a linear dependence between K and PSA. Calibration plots showed good agreement between expected and observed incidences (slope: 0.90-0.93) and moderately high discriminative power (area under the curve [AUC]: 0.68-0.69). Cross-validation showed satisfactory results (average AUCs in the training/validation groups: 0.66-0.70). The resulting dose-effect curves strongly depend on pre-RT PSA and GS. bRFS rapidly decreases with PSA: the maximum obtainable bRFS (defined as 95% of the maximum) declined by about 2.7% and 4.5% for each increment of 0.1 ng/mL for GS <7 and ≥7, respectively. Conclusions: Individual data were fitted by a TCP model, and the resulting best-fit parameters were radiobiologically consistent. The model suggests that relapses frequently result from clonogens outside the irradiated volume, supporting the choice of lymph-node irradiation, systemic therapy, or both for specific subgroups (GS <7: PSA >0.8-1.0 ng/mL; GS ≥7: PSA >0.3 ng/mL). Early RT should be preferred over delayed RT; the detrimental effect of PSA increase can never be fully compensated by increasing the dose, especially for patients with GS ≥7.

62 RADIOLOGY AND NUCLEAR MEDICINE↗