Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Methods in transport”

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 937 records · Page 52

Monte Carlo calculations of high energy nucleon meson cascades and applications to galactic cosmic ray transport

Results obtained using a recently developed calculational method for determining the nucleon-meson cascade induced in thick materials by high-energy nucleons and charged pions are presented. The calculational method uses the intranuclear-cascade-evaporation model to treat nonelastic collisions by particles with energies approximately or smaller than GeV and an extrapolation model at higher energies. The following configurations are considered: (1) 19.2-GeV/c protons incident on iron; (2) 30.3-GeV/c protons incident on iron; (3) solar and galactic protons incident on the moon, and (4) galactic protons incident on tissue. For the first three configurations, experimental results are available and comparisons between the experimental and calculated results are given.

Armstrong, T. W.↗

Computer program to perform cost and weight analysis of transport aircraft. Volume 2: Technical volume

An improved method for estimating aircraft weight and cost using a unique and fundamental approach was developed. The results of this study were integrated into a comprehensive digital computer program, which is intended for use at the preliminary design stage of aircraft development. The program provides a means of computing absolute values for weight and cost, and enables the user to perform trade studies with a sensitivity to detail design and overall structural arrangement. Both batch and interactive graphics modes of program operation are available.

Source record↗

Radiation shielding methods

System aids in charged-particle radiation transport analysis and shielding design. System can be adapted for shielding design and analysis in any charged-particle radiation environment.

Davis, H. S.↗

Combined assimilation of hydrography and TOPEX data into an Indian Ocean GCM using the adjoint method

To study seasonal circulation and meridional heat transport of the Indian Ocean by synthesizing dynamics with data, climatological monthly temperatures and salinities, surface heat and freshewater fluxes, and wind stresses, together with monthly ensembles of three years (93-95) of TOPEX-derived surface geostrophic velocity anomalies, are assimilated into an Indian Ocean GCM.

Hydrography Indian Ocean TOPEX meridional heat tra↗

Relationship between Weather, Traffic and Delay Based on Empirical Methods

The steady rise in demand for air transportation over the years has put much emphasis on the need for sophisticated air traffic flow management (TFM) within the National Airspace System (NAS). The NAS refers to hardware, software and people, including runways, radars, networks, FAA, airlines, etc., involved in air traffic management (ATM) in the US. One of the metrics that has been used to assess the performance of NAS is the actual delays provided through FAA's Air Traffic Operations Network (OPSNET). The OPSNET delay data includes those reportable delays, i.e. delays of 15 minutes or more experienced by Instrument Flight Rule (IFR) flights, submitted by the FAA facilities. These OPSNET delays are caused by the application of TFM initiatives in response to, for instance, weather conditions, increased traffic volume, equipment outages, airline operations, and runway conditions. TFM initiatives such as, ground stops, ground delay programs, rerouting, airborne holding, and miles-in-trail restrictions, are actions which are needed to control the air traffic demand to mitigate the demand-capacity imbalance due to the reduction in capacity. Consequently, TFM initiatives result in NAS delays. Of all the causes, weather has been identified as the most important causal factor for NAS delays. Therefore, in order to accurately assess the NAS performance, it has become necessary to create a baseline for NAS performance and establish a model which characterizes the relation between weather and NAS delays.

Sridhar, Banavar↗

Methods, systems, and apparatuses for calculating global fluence for neutron and photon monte carlo transport using expected value estimators

Global fluence estimators may be calculated on accelerators and processors for neutron and photon Monte Carlo transport. Monte Carlo random walk simulation may be performed on the processors and the calculation of a Volumetric-Ray-Casting (VRC) estimator may be offloaded to the accelerators. The VRC estimator may modify an expected-value estimator to extend a pseudo-particle ray along the direction of the emitted particle from source and collision event through not only the event volume, but also through all volumes that describe the problem geometry. Additionally, many pseudo-particle rays may be sampled per event, rather than just a single pseudo-particle ray per event, in order to provide more complete angular coverage.

Sweezy, Jeremy Ed↗

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

spopt: a python package for solving spatial optimization problems in PySAL

Spatial optimization is a major spatial analytical tool in management and planning, the significance of which cannot be overstated. Spatial optimization models play an important role in designing and managing effective and efficient service systems such as transportation, education, public health, environmental protection, and commercial investment among others. To this end, spopt (spatial optimization) is under active development for the inclusion of newly proposed models and methods for regionalization, facility location, and transportation-oriented solutions (Feng et al., 2021). Spopt is a submodule in the open-source spatial analysis library PySAL (Python Spatial Analysis Library) founded by Dr. Sergio J. Rey and Dr. Luc Anselin in 2005 (Rey et al., 2015, 2021; Rey & Anselin, 2007). The goal of developing spopt is to provide management and decision-making support to all relevant practitioners and to further promote the appropriate and meaningful application of spatial optimization models in practice.

97 MATHEMATICS AND COMPUTING↗

Extending TOUGH + HYDRATE with a parallel particle transport simulator: numerical investigation of sand production during gas production from hydrate deposits

A new parallel code for simulating particle transport in porous media is integrated with the TOUGH + HYDRATE simulator to investigate sand production associated with gas production from unconsolidated gas hydrate-bearing sediments (HBS). Here, the parallel coupled simulator is named THMPT and uses the integral finite difference method to describe the Darcian and non-Darcian flow of fluids and heat transport, the finite element method to describe the associated geomechanical changes, and the discrete element method to track the trajectory of individual sand particles within the HBS. The THMPT simulator is written in Fortran, incorporates multiple optimized algorithms, and can comprehensively address the coupled flow, thermal, chemical, geomechanical, and particle transport processes that characterize the system behaviors during gas production from HBS. The simulator can capture all processes involved in sand particle transport in porous media, including sand detachment, collision, clogging (i.e., bridging), and migration. A benchmark case study of sand production in the course of depressurization-induced gas production from a representative HBS reveals various distinct microscopic particle migration mechanisms and the adverse impact of sand particle detachment, transport, and clogging. The numerical investigation also examines the effect of bottomhole pressure on mitigating sand production. The simulation results indicate that sand clogging near the wellbore significantly reduces permeability, decreasing gas production by at least 50%. Lastly, the efficiency of gravel packing in mitigating sand production is numerically evaluated, revealing that the structure of the porous media appears to profoundly influence the macroscopic motion behavior of sand particles and sand clogging characteristics.

discrete element method↗

Random element method for numerical modeling of diffusional processes

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Ghoniem, A. F.↗

LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report

In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.

07 ISOTOPE AND RADIATION SOURCES↗

Space Transportation Systems Using Tethers

The groundwork has been laid for tether space transportation systems. NASA has developed tether technology for space applications since the 1960's. Important recent milestones include retrieval of a tether in space (TSS-1, 1992), successful deployment of a 20-km-long tether in space (SEDS-1, 1993), and operation of an electrodynamic tether with tether current driven in both directions-power and thrust modes (PMG, 1993). Various types of tethers and systems can be used for space transportation. Short electrodynamic tethers can use solar power to 'push' against a planetary magnetic field to achieve propulsion without the expenditure of propellant. The planned Propulsive Small Expendable Deployer System (ProSEDS) experiment will demonstrate electrodynamic tether thrust during its flight in early 2OOO. Utilizing completely different physical principles, long non-conducting tethers can exchange momentum between two masses in orbit to place one body into a transfer orbit for lunar and planetary missions. Recently completed system studies of this concept indicate that it would be a relatively low-cost in-space asset with long-term multimission capability. Both methods of using tethers for space transportation and propulsion are described in the paper.

Johnson, Les↗

Modeling of nonequilibrium space plasma flows

Godunov-type numerical solution of the 20 moment plasma transport equations. One of the centerpieces of our proposal was the development of a higher order Godunov-type numerical scheme to solve the gyration dominated 20 moment transport equations. In the first step we explored some fundamental analytic properties of the 20 moment transport equations for a low b plasma, including the eigenvectors and eigenvalues of propagating disturbances. The eigenvalues correspond to wave speeds, while the eigenvectors characterize the transported physical quantities. In this paper we also explored the physically meaningful parameter range of the normalized heat flow components. In the second step a new Godunov scheme type numerical method was developed to solve the coupled set of 20 moment transport equations for a quasineutral single-ion plasma. The numerical method and the first results were presented at several national and international meetings and a paper describing the method has been published in the Journal of Computational Physics. To our knowledge this is the first numerical method which is capable of producing stable time-dependent solutions to the full 20 (or 16) moment set of transport equations, including the full heat flow equation. Previous attempts resulted in unstable (oscillating) solutions of the heat flow equations. Our group invested over two man-years into the development and implementation of the new method. The present model solves the 20 moment transport equations for an ion species and thermal electrons in 8 domain extending from a collision dominated to a collisionless region (200 km to 12,000 km). This model has been applied to study O+ acceleration due to Joule heating in the lower ionosphere.

Gombosi, Tamas↗

Measurement of flow properties coupled to experimental and numerical analyses of dense, granular flows for solar thermal energy storage

Granular flows of sintered bauxite proppants were researched along an inclined plane for solar thermal energy storage applications. Granular flow properties needed to drive numeric granular models were measured for improved numerical model predictions for Carbobead CP 50/140, 40/100, and 30/60 particles. Particle shape and size distributions were determined by coupling optical microscopy to an in-house image processing algorithm. The impulse excitation technique was used to measure elastic and shear moduli, and compute Poisson’s ratio. The coefficient of static sliding friction was measured using the slip-stick method, and the static rolling friction was determined from measured shear on particles positioned between two hot-pressed plates. The coefficient of restitution was measured by dropping particles on a surface and determining the kinetic energy before and after impact with the surface using high resolution particle tracking velocimetry. Particle size did not significantly impact the coefficients of restitution and static rolling friction, however, particle shape distribution resulted in a large variation in measurements. An inclined flow experiment was performed to characterize granular flows of Carbobead CP 30/60 particles using particle image velocimetry. Numerical models of the experiment using discrete element method were generated with the measured mechanical properties as inputs for comparison with experimental results. A constant directional torque rolling friction model best predicted bulk granular flow behavior. Furthermore, good agreement between the model and experiment was achieved at ambient, steady state conditions, with average velocity differences <10%.

14 SOLAR ENERGY↗

Plans to improve the experimental limit in the comparison of the east-west and west-east one-way light propagation times on the rotating earth

The preceding paper describes the results so far (interrupted in the Spring of 1989 because of lack of funds) of an experiment comparing the one-way light propagation times on the surface of the rotating Earth. For the 20 Km path length component in the East-West direction the predicted difference between the opposite sense propagation times would be 160 ps, if the approximately 360 Km/s surface speed of the Earth gives effective light speeds of 3 x 10(exp 8) m/s +/- 360 m/s. This could lead to a prediction of the difference between the clock transport and the light pulse synchronization methods described in the preceding paper: delta(T) = 0.5 (160) = 80 ps. The current upper bound of approximately 100 ps for delta(T) is limited by poorly understood systematic errors. The most important seems to be intensity-dependent time delays in the remote light pulse avalanche photo-diode detector. This will be replaced by a continuously operating circular scan streak camera having single photon sensitivity and a time resolution of approximately 5 ps. (This camera has recently been developed by the Xian Institute of Optics and Precision Mechanics in the P.R.C.). Better isolation from shocks and vibration for the Sigma-Tau hydrogen maser during transport will be provided. It is hoped that delta(T) is less than 20 ps can be achieved.

Alley, C. O.↗