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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effects of Dead‐End Fractures on Non‐Fickian Transport in Three‐Dimensional Discrete Fracture Networks

Abstract Understanding mechanistic causes of non‐Fickian transport in fractured media is important for many hydrogeologic processes and subsurface applications. This study elucidates the effects of dead‐end fractures on non‐Fickian transport in three‐dimensional (3D) fracture networks. Although dead‐end fractures have been identified as low‐velocity regions that could delay solute transport, the direct relation between dead‐end fractures and non‐Fickian transport has been elusive. We systematically generate a large number of 3D discrete fracture networks with different fracture length distributions and fracture densities. We then identify dead‐end fractures using a novel graph‐based method. The effect of dead‐end fractures on solute residence time maximizes at the critical fracture density of the percolation threshold, leading to strong late‐time tailing. As fracture density increases beyond the percolation threshold, the network connectivity increases, and dead‐end fractures diminish. Consequently, the increase in network connectivity leads to a reduction in the degree of late‐time tailing. We also show that dead‐end fractures can inform about main transport paths, such as the mean tortuosity of particle trajectories. This study advances our mechanistic understanding of solute transport in 3D fracture networks.

58 GEOSCIENCES↗
DOE OSTI2022

Effects of local order parameter dependent transport coefficient in diblock copolymers under applied electric fields

In this work, we present an approach for constructing thermodynamically consistent time-dependent models relevant to thin films of diblock copolymers in applied electric fields. The approach is based on the principles of linear irreversible thermodynamics, and, in this work, it is applied to study the effects of electric fields on thin films of incompressible diblock copolymers. Enforcement of local incompressibility constraint at all times leads to a local order parameter dependent transport coefficient in the model for the diblock copolymers. The dependence of the transport coefficient on the local order parameter is used to relate it with the diffusion constant of Rouse chains and leads to sensitivity of the model to initial conditions. In addition, transient behavior is found to be affected when compared with an ad hoc model assuming a constant transport coefficient. Numerical results such as electric field induced alignment of lamellae domains due to the field are found to be in qualitative agreement with experiments. This approach opens up a systematic way of developing kinetic models for simulating effects of electrolytes added to thin films containing diblock copolymers in the presence of applied electric fields.

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

Fracture Network Influence on Rock Damage and Gas Transport following an Underground Explosion

Simulations of rock damage and gas transport following underground explosions that omit preexisting fracture networks in the subsurface cannot fully characterize the influence of geo-structural variability on gas transport. Previous studies do not consider the impact that fracture network structure and variability have on gas seepage. In this study, we develop a sequentially coupled, axi-symmetric model to look at the damage pattern and resulting gas breakthrough curves following an underground explosion given different fracture network realizations. We simulate 0.327 and 0.164 kT chemical explosives with burial depths of 100 m for 90 stochastically generated fracture networks. Gases quickly reach the surface in 30% of the higher yield simulations and 5% of the lower yield simulations. The fast breakthrough can be attributed to the formation of connected pathways between fractures to the surface. The formation of a connected damage pathway to the surface is not clearly correlated with the fracture intensity (P32) in our simulations. Breakthrough curves with slower transport are highly variable depending on the fracture network sample. The variability in the breakthrough behavior indicates that ignoring the influence of fracture networks on rock damage, which strongly influences the hydraulic properties following an underground explosion, will likely lead to a large underestimation of the uncertainty in the gas transport to the surface. This work highlights the need for incorporation of fracture networks into models for accurately predicting gas seepage following underground explosions.

58 GEOSCIENCES↗
NASA NTRS1979

Numerical simulation of ozone production, transport and distribution with a global atmospheric general circulation model

The production, transport, and distribution of ozone are simulated for a January with a global atmospheric general circulation model. In this model, the ozone influences the radiational heating and the photochemical ozone production and destruction, the radiational heating influences the atmospheric circulation, and the circulation redistributes the ozone. The model has simulated the synoptic and time-averaged observed large-scale fields of temperature, mass, and velocity in the troposphere and stratosphere, although with deficiencies which are seen in the time-averaged O3 distribution. Transport of ozone in the summer and winter hemispheres, in different seasons, and in various latitudes is discussed, and it is concluded that the divergence of the ozone transport maintains it below its photochemical equilibrium concentration in the tropics and subtropics, and the convergence of the ozone transport maintains it above its photochemical equilibrium concentration in the middle and high latitudes of both hemispheres. Thus, both the atmospheric motions and the O3 photochemistry determine the O3 sources and sinks.

Schlesinger, M. E.↗
NASA NTRS1981

Numerical modeling of diffusive-convective physical vapor transport in cylindrical vertical ampoules

Diffusive-convective physical vapor transport (PVT) in cylindrical, vertical ampoules of aspect ratio (length/radius) between 0.5 and 10 was modeled numerically. The transport of a crystal forming component through an inert component that undergoes zero net transport was considered. Systems were treated in which: (1) with unequal molecular weight of the components and with temperature gradients typically employed in PVT, convective flow arises dominantly from solutal density gradients; and (2) with equal molecular weight of the components, convective flow can arise only from thermal expansion. It was found that, due to the diffusion-induced horizontal density gradients, buoyancy-driven convective flows are superimposed on the diffusive-advective fluxes without threshold. Net recirculation sets in adjacent to the growing interface, in contrast to the corresponding monocomponent situation where marginally stable convective modes fill the whole fluid space. Depending on the orientation of the main transport direction with respect to gravity, convection can either reduce or enhance the diffusion-induced radial concentration gradients. Significant enhancement of the net transport rate was found to occur only when the whole vapor space between source and growing crystal is filled by a convective recirculation roll. Solutal and thermal convection results are similar; yet for quantitative discussions, thermal and solutal Rayleigh numbers are not interchangeable in contrast to convective situations that lack net mass transport across the fluid space.

Markham, B. L.↗
NASA NTRS1984

A numerical study of the vertical transport of momentum in a tropical rainband

The vertical transport of horizontal momentum in a convective tropical rainband is studied using a two-dimensional cloud ensemble model. Twelve simulations are made under the same large-scale conditions. The vertical transports of v momentum (parallel to the rainband) are essentially the same in all of the simulations, even though the structure of the clouds is different in each of the runs. The magnitude of the v-momentum transport by clouds is fairly large. It takes only half of a day to smooth out the tropical low-level easterly jet parallel to the rainband if no other processes are operating. The vertical transports of u momentum (perpendicular to the rainband) are quite different in all of the simulations. This difference can be explained by the dissimilarities in the distributions of horizontal momentum associated with various cloud configurations. The simulated vertical transports of horizontal momentum are compared with those computed with the Schneider and Lindzen scheme. The results suggest that their scheme is basically correct and usable if some improvements are made.

Soong, S.-T.↗
NASA NTRS1989

Solution of the hydrodynamic device model using high-order non-oscillatory shock capturing algorithms

A micron n+ - n - n+ silicon diode is simulated via the hydrodynamic model for carrier transport. The numerical algorithms employed are for the non-steady case, and a limiting process is used to reach steady state. The simulation employs shock capturing algorithms, and indeed shocks, or very rapid transition regimes, are observed in the transient case for the coupled system, consisting of the potential equation and the conservation equations describing charge, momentum, and energy transfer for the electron carriers. These algorithms, termed essentially non-oscillatory, were successfully applied in other contexts to model the flow in gas dynamics, magnetohydrodynamics, and other physical situations involving the conservation laws in fluid mechanics. The method here is first order in time, but the use of small time steps allows for good accuracy. Runge-Kutta methods allow one to achieve higher accuracy in time if desired. The spatial accuracy is of high order in regions of smoothness.

Fatemi, Emad↗
NASA NTRS1991

Solution of the hydrodynamic device model using high-order non-oscillatory shock capturing algorithms

A micron n+ - n - n+ silicon diode is simulated via the hydrodynamic model for carrier transport. The numerical algorithms employed are for the non-steady case, and a limiting process is used to reach steady state. The simulation employs shock capturing algorithms, and indeed shocks, or very rapid transition regimes, are observed in the transient case for the coupled system, consisting of the potential equation and the conservation equations describing charge, momentum, and energy transfer for the electron carriers. These algorithms, termed essentially nonoscillatory, were successfully applied in other contexts to model the flow in gas dynamics, magnetohydrodynamics, and other physical situations involving the conservation laws in fluid mechanics. The method here is first order in time, but the use of small time steps allows for good accuracy. Runge-Kutta methods allow one to achieve higher accuracy in time if desired. The spatial accuracy is of high order in regions of smoothness.

Fatemi, Emad↗
NASA NTRS1991

Current collection by conducting bodies in the low earth orbit

Collection of electrons by a conducting cylinder in a flowing magnetized plasma is studied by means of numerical simulations. Transport of electrons across magnetic field lines in regions of large electric fields in the sheath structure of the cylinder is demonstrated. The potential structures, current flow in the plasma and plasma distributions around the cylinder are presented with and without a release of a neutral gas plume near the cylinder. The simulation with the plume shows current enhancement and formation of a double layer, which is found to be highly dynamic. The simulation with the neutral gas plume simulates the plasma processes occurring in a plasma contactor.

Singh, Nagendra↗
NASA NTRS1993

Numerical modeling of physical vapor transport under microgravity conditions: Effect of thermal creep and stress

One of the most promising applications of microgravity (micro-g) environments is the manufacture of exotic and high-quality crystals in closed cylindrical ampoules using physical vapor transport (PVT) processes. The quality enhancements are believed to be due to the absence of buoyant convection in the weightless environment - resulting in diffusion-limited transport of the vapor. In a typical experiment, solid-phase sample material is initially contained at one end of the ampoule. The sample is made to sublime into the vapor phase and deposit onto the opposite end by maintaining the source at an elevated temperature with respect to the deposit. Identification of the physical factors governing both the rates and uniformity of crystal growth, and the optimization of the micro-g technology, will require an accurate modeling of the vapor transport within the ampoule. Previous micro-g modeling efforts have approached the problem from a 'classical' convective/diffusion formulation, in which convection is driven by the action of buoyancy on thermal and solutal density differences. The general conclusion of these works have been that in low gravity environments the effect of buoyancy on vapor transport is negligible, and vapor transport occurs in a diffusion-limited mode. However, it has been recently recognized than in the non-isothermal (and often low total pressure) conditions encountered in ampoules, the commonly-assumed no-slip boundary condition to the differential equations governing fluid motion can be grossly unrepresentative of the actual situation. Specifically, the temperature gradients can give rise to thermal creep flows at the ampoule side walls. In addition, temperature gradients in the vapor itself can, through the action of thermal stress, lead to bulk fluid convection.

Mackowski, Daniel W.↗
NASA NTRS1998

Numerical Modeling of Physical Vapor Transport in Contactless Crystal Growth Geometry

Growth from the vapor under conditions of limited contact with the walls of the growth ampoule is beneficial for the quality of the growing crystal due to reduced stress and contamination which may be caused by interactions with the growth container. The technique may be of a particular interest for studies on crystal growth under microgravity conditions: elimination of some factors affecting the crystal quality may make interpretation of space-conducted processes more conclusive and meaningful. For that reason, and as a part of our continuing studies on 'contactless' growth technique, we have developed a computational model of crystal growth process in such system. The theoretical model was built, and simulations were performed using the commercial computational fluid dynamics code, (CFD) ACE. The code uses an implicit finite volume formulation with a gray discrete ordinate method radiation model which accounts for the diffuse absorption and reflection of radiation throughout the furnace. The three-dimensional model computes the heat transfer through the crystal, quartz, and gas both inside and outside the ampoule, and mass transport from the source to the crystal and the sink. The heat transport mechanisms by conduction, natural convection, and radiation, and mass transport by diffusion and convection are modeled simultaneously and include the heat of the phase transition at the solid-vapor interfaces. As the thermal boundary condition, temperature profile along the walls of the furnace is used. For different thermal profiles and furnace and ampoule dimensions, the crystal growth rate and development of the crystal-vapor and source-vapor interfaces (change of the interface shape and location with time) are obtained. Super/under-saturation in the ampoule is determined and critical factors determining the 'contactless' growth conditions are identified and discussed. The relative importance of the ampoule dimensions and geometry, the furnace dimensions and its temperature, and the properties of the grown material are analyzed. The results of the simulations are compared with related experimental results on growth of CdTe, CdZnTe, ZnTe, PbTe, and PbSnTe crystals by this technique.

Palosz, W.↗
NASA NTRS1998

Temperature Dependence of Diffusivities in Liquid Elements (LMD)

This research was to advance the understanding of diffusion mechanisms in liquid metals and alloys through accurate diffusivity measurements over a wide range of temperatures, including the proximity of the materials melting points. Specifically, it was driven towards developing a methodology (and subsequent flight hardware) to enable several diffusion coefficient measurements (i.e., at several different temperatures) to be performed using a single sample. The Liquid Metal Diffusion (LMD) was funded as a Flight Definition Project in February 1993 in response to NRA 91-OSSA-20 (Microgravity Science and Applications Division). The Science Concept Review for LAID was held during April 1994. In January 1995 we were informed that we had failed this review and the project was change to ground-based activities only. A new proposal was submitted for the next NRA addressing the panels concerns. As part of NASA's Risk Mitigation program, a scaled-down version of the hardware was funded in July of 1995 for a flight opportunity utilizing experiment on the Microgravity Isolation Mount. This experiment was to determine the self-diffusivity of indium at 185 C. The LMD was transferred to the Mir Space Station in STS-81 and returned on STS-84 (January - May 1997). Three, out of five, self-diffusion data sets were returned. A description of this experiment/hardware is included below. This summary is only intended to give the reader an overview of the results obtained for the tasks outlined in the original proposal. Research that was not published is explained in more detail. At the end of this report is a list of refereed publications and invited talks that were given as a result of this work. The reader is directed to these for further details. Attachment: Real-time diffusivity measurements in liquids at several temperatures with one sample, On the insensitivity of liquid diffusivity measurements to deviations from 1D transport, and Numerical simulations of the convective contamination of diffusivity measurements in liquids.

Banish, R. Michael↗
NASA NTRS1999

Towards a Strategy of Remote Sensing of Aerosol Direct and Indirect Radiative Forcing

The role of aerosol forcing remains one of the largest uncertainties in estimating man's impact on the global climate system. One school of thought suggests that remote sensing by satellite sensors will provide the data necessary to narrow these uncertainties. Much effort has gone into the development of new satellite sensors specifically designed to retrieve aerosol loading and some information about the sizes of the aerosols. These next generation remote sensing instruments (EOS-MODIS, EOS-MISR, POLDER, ATSR) will provide unprecedented accuracy in the retrieval of aerosol loading. Although the new generation of sensors has excellent accuracy compared to the heritage instruments of the past, they still have measurement limitations. In clean, pristine regions the absolute magnitude of the uncertainty in the aerosol retrieval becomes comparable in magnitude to the signal itself. If much of the aerosol forcing is occurring at very low magnitudes of aerosol concentrations, satellite remote sensing will miss it. This study attempts to quantify remote sensing limitations due to the accuracy limits of the retrieval algorithms. We use a combination of numerical aerosol transport models, ground-based AERONET data and ISCCP cloud climatology to determine how much of the forcing occurs in regions too clean to determine from satellite retrievals. This study is not an intercomparison of global transport models. It is not an estimation of global aerosol forcing. This study is an exercise to determine whether satellite remote sensing can live up to our high expectations. The results show that while remote sensing will be able to measure a majority of the direct forcing, it will miss roughly half of the indirect forcing.

Remer, L. A.↗
NASA NTRS1999

Satellite Remote Sensing of Aerosol Forcing

The role of aerosol forcing remains one of the largest uncertainties in estimating man's impact on the global climate system. One school of thought suggests that remote sensing by satellite sensors will provide the data necessary to narrow these uncertainties. While satellite measurements of direct aerosol forcing appear to be straightforward, satellite measurements of aerosol indirect forcing will be more complicated. Pioneering studies identified indirect aerosol forcing using AVHRR data in the biomass burning regions of Brazil. We have expanded this analysis with AVHRR to include an additional year of data and assimilated water vapor fields. The results show similar latitudinal dependence as reported by Kaufman and Fraser, but by using water vapor observations we conclude that latitude is not a proxy for water vapor and the strength of the indirect effect is not correlated to water vapor amounts. In addition to the AVHRR study we have identified indirect aerosol forcing in Brazil at much smaller spatial scales using the MODIS Airborne Simulator. The strength of the indirect effect appears to be related to cloud type and cloud dynamics. There is a suggestion that some of the cloud dynamics may be influenced by smoke destabilization of the atmospheric column. Finally, this study attempts to quantify remote sensing limitations due to the accuracy limits of the retrieval algorithms. We use a combination of numerical aerosol transport models, ground-based AERONET data and ISCCP cloud climatology to determine how much of the forcing occurs in regions too clean to determine from satellite retrievals.

Remer, Lorraine↗
NASA NTRS2000

Uncertainties in Remote Sensing of Aerosol Radiative Forcing

The role of aerosol forcing remains one of the largest uncertainties in estimating man's impact on the global climate system. One school of thought suggests that remote sensing by satellite sensors will provide the data necessary to narrow these uncertainties. Much effort has gone into the development of new satellite sensors specifically designed to retrieve aerosol loading and some information about the sizes of the aerosols. These next generation remote sensing instruments (EOS-MODIS, EOS-MISR, POLDER, ATSR) will provide unprecedented accuracy in the retrieval of aerosol loading. Although the new generation of sensors has excellent accuracy compared to the heritage instruments of the past, they still have measurement limitations. In clean, pristine regions the absolute magnitude of the uncertainty in the aerosol retrieval becomes comparable in magnitude to the signal itself. If much of the aerosol forcing is occurring at very low magnitudes of aerosol concentrations, satellite remote sensing will miss it. This study attempts to quantify remote sensing limitations due to the accuracy limits of the retrieval algorithms. We use a combination of numerical aerosol transport models, ground-based AERONET data and ISCCP cloud climatology to determine how much of the forcing occurs in regions too clean to determine from satellite retrievals. This study is not an intercomparison of global transport models. It is not an estimation of global aerosol forcing. This study is an exercise to determine whether satellite remote sensing can live up to our high expectations.

Remer, L.↗