Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical 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 343 records · Page 19

Space Transport and Fatigue

Space travel presents numerous challenges to fatigue management. Fatigue in space results from many of the same challenges that shiftworkers face on Earth, including inadequate sleep, circadian misalignment, and elevated workload. However, spaceflight also presents unique challenges including circadian misalignment due to non-24- hour light-dark cycles, erratic work schedules due to operational constraints, a microgravity sleep environment, constraints to habitat design, and an isolated and confined living environment. These challenges can be mitigated with fatigue management strategies and countermeasures specifically targeted for space travel. Those strategies can include light interventions, advanced scheduling tools, improved habitat design, and pharmacological countermeasures. What is needed is a future-facing approach to most effectively manage the challenges that may arise as advancements in space travel continue to develop at a rapid pace.

sleep deprivation and performance↗

Derivation of a three dimensional numerical water quality model for estuary and continental shelf application

A derivation is given for a three dimensional mass transport equation which is appropriate for numerical modeling of estuary and continental shelf water quality variations for both the time dependent and steady state cases. A finite difference approximation to the derived equation is presented and a solution scheme for the resulting equations outlined. Preliminary results are obtained using the model for the extremely simple problems which have analytical solutions. The numerical model, as presented, will provide a scheme to study water quality problems in coastal waters for both steady state and time dependent cases.

Spaulding, M.↗

Simulation of lithium transport using the BOUT++ framework

A numerical model that calculates the collisional interactions between the lithium atoms from a lithium pellet and the background plasmas has been upgraded. The ion density ($N_t$), electron temperature ($T_e$), ion temperature ($T_i$) and parallel ion velocity ($V_{∥, i}$) are used to characterize the background plasmas. The lithium atom density ($N^{a}_{Li}$) and parallel velocity ($V_{∥,a}$) of lithium atoms evolve with time. For each lithium ion, the density ($N_{Li^{n+}}$), temperature ($T_{Li^{n+}}$) and parallel velocity ($V_{∥, Li^{n+}}$) are self-consistently calculated. A C-mod lower single null equilibrium is used to generate the grid for the BOUT++ simulation. The lithium atoms can be fully ionized to $Li^{3+}$ in ~2 μs. The rapid radial and poloidal expansion of the lithium ions are found in the simulation. After the collision interaction process, the electron temperature rapidly decreases at the pellet location; then, it rapidly poloidally expands, and the temperature at the pellet location starts to recover. The electron pressure increases at the pellet location despite the decrease in electron temperature because of the extra electrons from the lithium ionization. The ion pressure profile decreases in the pellet location due to the decrease in ion temperature.

74 ATOMIC AND MOLECULAR PHYSICS↗

1D model of tokamak scrape-off layer and divertor plasmas: Dimensionless parameters and the results of numerical simulations

Solutions of a 1D model of the scrape-off layer and a divertor plasma in a tokamak are analyzed on the basis of dimensionless parameters governing the physics of both plasma and neutral transport and atomic physics processes. Numerical simulations based on fluid plasma and Monte Carlo neutral modeling confirm earlier results obtained with simplified analytic considerations and numerical simulations performed with fluid neutral description and demonstrate bifurcation of plasma parameters. Further, these simulations also confirm that the major physical reasons for the rollover of the plasma flux to the target, which is the key feature of the phenomenon of divertor plasma detachment, are the impurity radiation loss and plasma recombination. It is suggested that simulations of the edge plasma based on dimensionless parameters provide a powerful tool for checking the consistency of the physics incorporated in existing complex 2D codes used for edge plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimizing Carbon Capture, Transport, and Storage: Overcoming Challenges with Machine Learning and Cost-Benefit Analysis

Carbon Capture, Utilization, and Storage (CCUS) is a critical strategy for reducing CO₂ emissions and mitigating climate change. However, its widespread deployment faces numerous challenges across the capture, transport, and storage phases. These challenges include the technical complexity of predicting subsurface behaviors during CO₂ injection, ensuring long-term storage integrity, optimizing transportation networks, and balancing the economic and environmental trade-offs. Addressing these issues requires an integrated approach combining advanced subsurface modeling with system-level analyses to assess costs, risks, and benefits. This presentation provides an overview of studies conducted by the National Energy Technology Laboratory (NETL) to tackle these challenges. NETL’s efforts encompass cutting-edge research in subsurface fluid behavior machine learning predictions, alongside the development of innovative tools for system optimization and economic evaluation. By bridging technical expertise and strategic analysis, NETL aims to advance the deployment of CCUS technologies to support global decarbonization efforts. Presented at the Carnegie Mellon University CEE IESS Student Seminar October 4, 2024.

Shih, Chung Yan↗

Specialization Restricts the Evolutionary Paths Available to Yeast Sugar Transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence–function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ~400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

59 BASIC BIOLOGICAL SCIENCES↗

Least-squares finite element method for fluid dynamics

An overview is given of new developments of the least squares finite element method (LSFEM) in fluid dynamics. Special emphasis is placed on the universality of LSFEM; the symmetry and positiveness of the algebraic systems obtained from LSFEM; the accommodation of LSFEM to equal order interpolations for incompressible viscous flows; and the natural numerical dissipation of LSFEM for convective transport problems and high speed compressible flows. The performance of LSFEM is illustrated by numerical examples.

Jiang, Bo-Nan↗

Least-squares finite element method for fluid dynamics

An overview is given of new developments of the least squares finite element method (LSFEM) in fluid dynamics. Special emphasis is placed on the universality of LSFEM; the symmetry and positiveness of the algebraic systems obtained from LSFEM; the accommodation of LSFEM to equal order interpolations for incompressible viscous flows; and the natural numerical dissipation of LSFEM for convective transport problems and high speed compressible flows. The performance of LSFEM is illustrated by numerical examples.

Jiang, Bo-Nan↗

A direct method for the solution of unsteady two-dimensional incompressible Navier-Stokes equations

The unsteady incompressible Navier-Stokes equations are formulated in terms of vorticity and stream function in generalized curvilinear orthogonal coordinates to facilitiate analysis of flow configurations with general geometries. The numerical method developed solves the conservative form of the transport equation using the alternating-direction implicit method, whereas the stream-function equation is solved by direct block Gaussian elimination. The method is applied to a model problem of flow over a back-step in a doubly infinite channel, using clustered conformal coordinates. One-dimensional stretching functions, dependent on the Reynolds number and the asymptotic behavior of the flow, are used to provide suitable grid distribution in the separation and reattachment regions, as well as in the inflow and outflow regions. The optimum grid distribution selected attempts to honor the multiple length scales of the separated-flow model problem. The asymptotic behavior of the finite-differenced transport equation near infinity is examined and the numerical method is carefully developed so as to lead to spatially second-order accurate wiggle-free solutions, i.e., with minimum dispersive error. Results have been obtained in the entire laminar range for the backstep channel and are in good agreement with the available experimental data for this flow problem.

Ghia, K. N.↗

Assessing the release, transport, and retention of radioactive aerosols from hypothetical breaches in spent fuel storage canisters

Interim dry storage of spent nuclear fuel involves storing the fuel in welded stainless-steel canisters. Under certain conditions, the canisters could be subjected to environments that may promote stress corrosion cracking leading to a risk of breach and release of aerosol-sized particulate from the interior of the canister to the external environment through the crack. Research is currently under way by several laboratories to better understand the formation and propagation of stress corrosion cracks, however little work has been done to quantitatively assess the potential aerosol release. The purpose of the present work is to introduce a reliable generic numerical model for prediction of aerosol transport, deposition, and plugging in leak paths similar to stress corrosion cracks, while accounting for potential plugging from particle deposition. The model is dynamic (changing leak path geometry due to plugging) and it relies on the numerical solution of the aerosol transport equation in one dimension using finite differences. The model’s capabilities were also incorporated into a Graphical User Interface (GUI) that was developed to enhance user accessibility. Model validation efforts presented in this paper compare the model’s predictions with recent experimental data from Sandia National Laboratories (SNL) and results available in literature. We expect this model to improve the accuracy of consequence assessments and reduce the uncertainty of radiological consequence estimations in the remote event of a through-wall breach in dry cask storage systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Behavior of nitric oxide formed by the space shuttle in the mesosphere

A numerical model of eddy diffusive transport of emitted gases from the space shuttle wake, including chemical reactions between the emitted constituents and the ambient atmosphere, has been constructed for 75 km altitude. The numerical methods involve explicit solution of the diffusion equation and Runge-Kutta method for the chemical reactions. The time required to reach background levels of nitric oxide concentration of 7 x 10 to the 7th power molecules/cc has been calculated. This relaxation time depends strongly on atmospheric conditions.

Eberstein, I. J.↗

Guidelines for Computing Longitudinal Dynamic Stability Characteristics of a Subsonic Transport

A systematic study is presented to guide the selection of a numerical solution strategy for URANS computation of a subsonic transport configuration undergoing simulated forced oscillation about its pitch axis. Forced oscillation is central to the prevalent wind tunnel methodology for quantifying aircraft dynamic stability derivatives from force and moment coefficients, which is the ultimate goal for the computational simulations. Extensive computations are performed that lead in key insights of the critical numerical parameters affecting solution convergence. A preliminary linear harmonic analysis is included to demonstrate the potential of extracting dynamic stability derivatives from computational solutions.

Thompson, Joseph R.↗

Magnetic flux transport on the sun

Although most of the magnetic flux observed on the sun originates in the low-latitude sunspot belts, this flux is gradually dispersed over a much wider range of latitudes by supergranular convective motions and meridional circulation. Numerical simulations show how these transport processes interact over the 11-year sunspot cycle to produce a strong 'topknot' polar field, whose existence near sunspot minimum is suggested by the observed strength of the interplanetary magnetic field and by the observed areal extent of polar coronal holes. The required rates of diffusion and flow are consistent with the decay rates of active regions and with the rotational properties of the large-scale solar magnetic field.

Wang, Y.-M.↗

Particle erosion mechanisms and mass redistribution in Saturn's rings

A variety of physical processes can erode the surfaces of planetary ring particles. According to current estimates, the most efficient of these over the bulk of Saturn's rings is hypervelocity impact by 100 microns to one centimeter radius meteoroids. The atoms, molecules, and fragments ejected from ring particles by erosion arc across the rings along elliptical orbits to produce a tenuous halo of solid ejecta and an extensive gaseous atmosphere. Continuous exchange of ejecta between different ring regions can lead to net radial transport of mass and angular momentum. The equations governing this ballistic transport process are presented and discussed. Both numerical and analytic studies of idealized ring systems illustrate that ballistic transport can cause significant mass redistribution in the rings, especially near regions of high density contrast, such as the inner edges of the A and B rings. Ejecta exchanges can also alter local particle sizes and compositions and may produce pulverized regoliths at least several centimeters deep. The meteoroid erosion rate is so high that significant global torques and mass loss are possible on times shorter than a solar system life time.

Durisen, R. H.↗

A numerical simulation of the magnetospheric gate model for the X-ray bursters

A Lagrangian, stellar interior computer code which assumes radiation transport by diffusion is used to numerically simulate an X-ray burst occurring on a rapidly rotating, 1 solar mass neutron star accreting material from a companion in noncontinuous, episodic fashion. The blob of gas from the companion extends from the surface of the neutron star to a radius of 185 km, is optically thick, has mass of 10 to the -11th solar masses, and is initially at rest with respect to the surface. The gas initially has a temperature of about 10 million K and a surface luminosity of 0.1 that of the sun. The infall results in a burst which lasts about 0.1 second and reaches a peak luminosity and effective temperature of 240,000 suns and nine million K. The burst is followed by a phase of oscillations with a period of 0.2 second. It is found that the burst is too cool and rapid to resemble normal bursts, although the time scale is in agreement with some observations.

Starrfield, S.↗

Amesh v1.0

Amesh generates discretized grids for numerical modeling of flow and transport problems in which the formulation is based on the integral finite difference method (IFDM). For example, the output of Amesh can be used directly as input to the TOUGH2 or TOUGH3 numerical simulators. Amesh can generate 1D, 2D, or 3D numerical grids for a given set of locations, i.e., the centers of each discrete sub-domain. In the 2D areal plane Voronoi tessellation is used.

Fuller, Peter↗

The Wiedemann-Franz law in doped Mott insulators without quasiparticles

Many metallic quantum materials display anomalous transport phenomena that defy a Fermi liquid description. In this report we use numerical methods to calculate thermal and charge transport in the doped Hubbard model and observe a crossover separating high- and low-temperature behaviors. Distinct from the behavior at high temperatures, the Lorenz number L becomes weakly doping dependent and less sensitive to parameters at low temperatures. At the lowest numerically accessible temperatures, L roughly approaches the Wiedemann-Franz constant L, even in a doped Mott insulator that lacks well-defined quasiparticles. Decomposing the energy current operator indicates a compensation between kinetic and potential contributions, which may help to clarify the interpretation of transport experiments beyond Boltzmann theory in strongly correlated metals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

n-Propylbenzene Droplet Combustion: Experiments and Numerical Modeling

n-Propylbenzene (PB, C9H12, 432K boiling point) is important as a constituent in surrogates for some transportation fuels but little is understood about its droplet burning characteristics. We present an experimental and numerical study of a fuel droplet burning under conditions where the sole means of gas transport is fuel evaporation at the droplet surface thereby creating a one-dimensional gas transport that is amenable to detailed numerical modeling. The one-dimensional gas transport that results is amenable to detailed numerical modeling. Initial droplet diameters (Do) ranged from 0.3 mm to about 5 mm to better understand the mechanisms that control burning over such a wide range of Do. For Do < 1 mm droplets were burned while anchored to 14 m SiC fibers within a sealed container in free-fall during which time the droplets were ignited by spark discharge and the burning process recorded by video cameras. Larger and free-floating droplets were studied in the spaced-based platform of the orbiting International Space Station for Do > 1.0 mm which provided unlimited times for experimental observation. Measurements of droplet, flame and soot shell diameters were made and are discussed in the presentation. The experimental results showed formation of optically thick soot clouds developing during burning that required a robust image analysis algorithm to enable measurement of droplet diameter. The measured droplet diameters were well predicted by a numerical model that included soot formation, a PB kinetic mechanism comprising 576 species and 27369 reactions, unsteady gas and liquid transport, variable properties and radiative transport. Flame diameters showed more variability compared to the simulations, depending on how the flame was defined (peak gas temperature and OH concentration). Droplets burned slower as Do was increased and the PB sooting propensity increased with Do. Flame extinction was observed for Do > 3.3 mm indicating a strong effect of radiation as Do was increased compared to the smaller droplet diameters examined where extinction was not experimentally observed. For the largest droplets examined the possibility that the observed extinction is indicative of a transition to cool flame burning is still under investigation.

Guo, Songtao↗