Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ASPECT RATIO”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Initial transport and turbulence analysis and gyrokinetic simulation validation in NSTX-U L-mode plasmas

Transport analysis, ion-scale turbulence measurements, and initial linear and nonlinear gyrokinetic simulations are reported for a transport validation study based on low aspect ratio NSTX-U L-mode discharges. The relatively long, stationary L-modes enabled by the upgraded centerstack provide a more ideal target for transport validation studies that were not available during NSTX operation. Transport analysis shows that anomalous electron transport dominates energy loss while ion thermal transport is well described by neoclassical theory. Linear gyrokinetic GYRO analysis predicts that ion temperature gradient (ITG) modes are unstable around normalized radii $\rho$=0.6-0.8, although $E\timesB$ shearing rates are larger than the linear growth rates over much of that region. Deeper in the core ($\rho$=0.4-0.6), electromagnetic microtearing modes (MTM) are unstable as a consequence of the relatively high beta and collisionality in these particular discharges. Consistent with the linear analysis, local, nonlinear ion-scale GYRO simulations predict strong ITG transport at $\rho$=0.76, whereas electromagnetic MTM transport is important at $\rho$=0.47. The prediction of ion-scale turbulence is consistent with 2D beam emission spectroscopy (BES) that measures the presence of broadband ion-scale fluctuations. Interestingly, the BES measurements also indicate the presence of bi-modal poloidal phase velocity propagation that could be indicative of two different turbulence types. However, in the region between ($\rho$=0.56, 0.66), ion-scale simulations are strongly suppressed by the locally large $E\timesB$ shear. Instead, electron temperature gradient (ETG) turbulence simulations predict substantial transport, illustrating electron-scale contributions can be important in low aspect ratio L-modes, similar to recent analysis at conventional aspect ratio. However, agreement within experimental uncertainties has not been demonstrated, which requires additional simulations to test parametric sensitivities. The potential need to include profile-variation effects (due to the relatively large value of $\rho_*$=$\rho_i$/a at low aspect ratio), including electromagnetic and possibly multi-scale effects, is also discussed.

Validation↗
NASA NTRS1973

A study of the effects of aeroelastic divergence on the wing structure of an oblique-wing supersonic transport configuration

The aerodynamic characteristics of transport aircraft with oblique wing flying at supersonic speeds are discussed. Aeroelastic divergence of the forward swept portion of the wing is analyzed. The effect of aspect ratio as a method for avoiding aeroelastic divergence is examined. A relatively low aspect ratio appears necessary for an oblique wing when constructed of conventional aluminum alloy materials. The aspect ratio may be increased by increasing the wing thickness ratio and by utilizing materials with higher moduli of elasticity and rigidity.

Source record↗
NASA NTRS2015

The Effects of Sweeping Jet Actuator Parameters on Flow Separation Control

A parametric experimental study was performed with sweeping jet actuators (fluidic oscillators) to determine their effectiveness in controlling flow separation on an adverse pressure gradient ramp. Actuator parameters that were investigated include blowing coefficients, operation mode, pitch and spreading angles, streamwise location, aspect ratio, and scale. Surface pressure measurements and surface oil flow visualization were used to characterize the effects of these parameters on the actuator performance. 2D Particle Image Velocimetry measurements of the flow field over the ramp and hot-wire measurements of the actuator's jet flow were also obtained for selective cases. In addition, the sweeping jet actuators were compared to other well-known flow control techniques such as micro-vortex generators, steady blowing, and steady vortex-generating jets. The results confirm that the sweeping jet actuators are more effective than steady blowing and steady vortex-generating jets. The results also suggest that an actuator with a larger spreading angle placed closer to the location where the flow separates provides better performance. For the cases tested, an actuator with an aspect ratio, which is the width/depth of the actuator throat, of 2 was found to be optimal. For a fixed momentum coefficient, decreasing the aspect ratio to 1 produced weaker vortices while increasing the aspect ratio to 4 reduced coverage area. Although scaling down the actuator (based on the throat dimensions) from 0.25 inch x 0.125 inch to 0.15 inch x 0.075 inch resulted in similar flow control performance, scaling down the actuator further to 0.075 inch x 0.0375 inch reduced the actuator efficiency by reducing the coverage area and the amount of mixing in the near-wall region. The results of this study provide insight that can be used to design and select the optimal sweeping jet actuator configuration for flow control applications.

Koklu, Mehti↗
NASA NTRS1955

Aerodynamic Characteristics of Various Configurations of a Model of a 45 Degree Swept-wing Airplane at a Mach Number of 2.01

An investigation has been conducted at the Langley 4- by 4-foot supersonic pressure tunnel at a Mach nmber of 2.01 to determine the aerodynamic characteristics of several configurations of a model of a 45 deg swept-wing airplane. The basic configuratin had a wing with 45 deg sweepback at the quarter-chord line, aspect ration 3.2, taper ration 0.468, NACA 65A005.5 sections just outboard of the inlet and NACA 65A003.7 sections at the tip. The wing was mounted slightly above the body center line and an all-movable horizantal tail was located slightly below the extended chord line of the wing. Tre design incorporated twin wing-root supersonic inlets ducted to a single exit at the base of the fuselage. The configurations investigated included an extended nose length, a bumped-fuselage afterbody, an inlet droop, an lncreased wing aspect ratio, and a revised canopy shape. Configurations employing the wing of increased aspect ratio of 3.7, which constituted the bulk of the tests, produced about a 10-percent increase in lift and in longitudinal stability as compared with the basic wing of aspect ratio 3.2. There was a slight but masurable increase in minimum drag and maximum lift-drag ratio.

Spearman, Leroy M↗
DOE OSTI2013

EFFECTS OF RESISTIVITY ON MAGNETIZED CORE-COLLAPSE SUPERNOVAE

We studied the role of turbulent resistivity in the core-collapse of a strongly magnetized massive star, carrying out two-dimensional resistive-MHD simulations. Three cases with different initial strengths of magnetic field and rotation are investigated: (1) a strongly magnetized rotating core, (2) a moderately magnetized rotating core, and (3) a very strongly magnetized non-rotating core. In each case, one ideal-MHD model and two resistive-MHD models are computed. As a result of these computations, each model shows an eruption of matter assisted by magnetic acceleration (and also by centrifugal acceleration in the rotating cases). We found that resistivity attenuates the explosion in cases 1 and 2, while it enhances the explosion in case 3. We also found that in the rotating cases, the main mechanisms for the amplification of a magnetic field in the post-bounce phase are an outward advection of the magnetic field and a twisting of poloidal magnetic field lines by differential rotation, which are somewhat dampened down with the presence of resistivity. Although magnetorotational instability seems to occur in the rotating models, it plays only a minor role in magnetic field amplification. Another impact of resistivity is that on the aspect ratio. In the rotating cases, a large aspect ratio of the ejected matter, >2.5, attained in an ideal-MHD model is reduced to some extent in a resistive model. These results indicate that resistivity possibly plays an important role in the dynamics of strongly magnetized supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Investigation of precipitate refinement in Mg alloys by an analytical composite failure model

An analytical model is developed to simulate precipitate refinement in second phase strengthened magnesium alloys. The model is developed based on determination of the stress fields inside elliptical precipitates embedded in a rate dependent inelastic matrix. The stress fields are utilized to determine the failure mode that governs the refinement behavior. Using an AZ31 Mg alloy as an example, the effects the applied load, aspect ratio and orientation of the particle is studied on the macroscopic failure of a single α-Mg17Al12 precipitate. Additionally, a temperature dependent version of the corresponding constitutive law is used to incorporate the effects of temperature. In plane strain compression, an extensional failure mode always fragments the precipitates. The critical strain rate at which the precipitates start to fail strongly depends on the orientation of the precipitate with respect to loading direction. The results show that the higher the aspect ratio is, the easier the precipitate fractures. Precipitate shape is another factor influencing the failure response. In contrast to elliptical precipitates with high aspect ratio, spherical precipitates are strongly resistant to sectioning. In pure shear loading, in addition to the extensional mode of precipitate failure, a shearing mode may get activated depending on orientation and aspect ratio of the precipitate. The effect of temperature in relation to strain rate was also verified for plane strain compression and pure shear loading cases.

Tabei, Ali↗
NASA NTRS2005

Aperture-coupled thin-membrane microstrip array antenna for beam scanning application

A microstrip array using aperture-slot-coupling technique with very thin membranes has been developed at the L-band frequency for beam scanning application. This technology-demonstration array with 4 x 2 elements achieved a relatively wide bandwidth of 100 MHz (8%) and +/-45(deg) beam scan. Very narrow coupling slots were used with each having an aspect ratio of 160 (conventional slot aspect ratio is between 10 to 30) for coupling through very thin membrane (0.05mm thickness). This thin-membrane aperture-coupling technique allows the array antenna elements to be more easily integrated with transmit/receive amplifier (T/R) and phase shifter modules. This paper addresses only the radiator portion of the array. The array and active components will be presented in a separate pape.

array antenna↗
DOE OSTI2014

Limitation of parallel flow in double diffusive convection: Two- and three-dimensional transitions in a horizontal porous domain

Two-dimensional (2D) and three-dimensional (3D) numerical simulations of double diffusion natural convection in an elongated enclosure filled with a binary fluid saturating a porous medium are carried out in the present work. The Boussinesq approximation is made in the formulation of the problem, and Neumann boundary conditions for temperature and concentration are adopted, respectively, on vertical and horizontal walls of the cavity. The used numerical method is based on the control volume approach, with the third order quadratic upstream interpolation scheme in approximating the advection terms. A semi implicit method algorithm is used to handle the velocity-pressure coupling. To avoid the excessively high computer time inherent to the solution of 3D natural convection problems, full approximation storage with full multigrid method is used to solve the problem. A wide range of the controlling parameters (Rayleigh-Darcy number Ra, lateral aspect ratio Ay, Lewis number Le, and the buoyancy ration N) is investigated. We clearly show that increasing the depth of the cavity (i.e., the lateral aspect ratio) has an important effect on the flow patterns. The 2D perfect parallel flows obtained for small lateral aspect ratio are drastically destabilized by increasing the cavity lateral dimension. This yields a 3D fluid motion with a much more complex flow pattern and the usually considered 2D parallel flow model cannot be applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2008

Fabricating Nanodots using Lift-Off of a Nanopore Template

A process for fabricating a planar array of dots having characteristic dimensions of the order of several nanometers to several hundred nanometers involves the formation and use of a thin alumina nanopore template on a semiconductor substrate. The dot material is deposited in the nanopores, then the template is lifted off the substrate after the dots have been formed. This process is expected to be a basis for development of other, similar nanofabrication processes for relatively inexpensive mass production of nanometerscale optical, optoelectronic, electronic, and magnetic devices. Alumina nanopore templates are self-organized structures that result from anodization of aluminum under appropriate conditions. Alumina nanopore templates have been regarded as attractive for use in fabricating the devices mentioned above, but prior efforts to use alumina nanopore templates for this purpose have not been successful. One reason for the lack of success is that the aspect ratios (ratios between depth and diameter) of the pores have been too large: large aspect ratios can result in blockage of deposition and/or can prevent successful lift-off. The development of the present process was motivated partly by a requirement to reduce aspect ratios to values (of the order of 10) for which there is little or no blockage of deposition and attempts at lift-off are more likely to be successful. The fabrication process is outlined.

Yang, Eui-Hyeok↗
NASA NTRS1994

Oscillatory/chaotic thermocapillary flow induced by radiant heating

The objective of this paper is to conduct ground-based experiments to measure the onset conditions of oscillatory Marangoni flow in laser-heated silicone oil in a cylindrical container. For a single fluid, experimental data are presented using the aspect ratio and the dynamic Bond number. It is found that for a fixed aspect ratio, there seems to be an asymptotic limit of the dynamic Bond number beyond which no onset of flow oscillation could occur. Experimental results also suggested that there could be a lower limit of the aspect ratio below which there is no onset of oscillatory flow.

Hsieh, Kwang-Chung↗
NASA NTRS2006

Improved Oxygen-Beam Texturing of Glucose-Monitoring Optics

An improved method has been devised for using directed, hyperthermal beams of oxygen atoms and ions to impart desired textures to the tips of polymethylmethacrylate [PMMA] optical fibers to be used in monitoring the glucose content of blood. The improved method incorporates, but goes beyond, the method described in Texturing Blood-Glucose- Monitoring Optics Using Oxygen Beams (LEW-17642-1), NASA Tech Briefs, Vol. 29, No. 4 (April 2005), page 11a. The basic principle of operation of such a glucose-monitoring sensor is as follows: The textured surface of the optical fiber is coated with chemicals that interact with glucose in such a manner as to change the reflectance of the surface. Light is sent down the optical fiber and is reflected from, the textured surface. The resulting change in reflectance of the light is measured as an indication of the concentration of glucose. The required texture on the ends of the optical fibers is a landscape of microscopic cones or pillars having high aspect ratios (microscopic structures being taller than they are wide). The average distance between hills must be no more than about 5 mso that blood cells (which are wider) cannot enter the valleys between the hills, where they would interfere with optical sensing of glucose in the blood plasma. On the other hand, the plasma is required to enter the valleys, and high aspect ratio structures are needed to maximize the surface area in contact with the plasma, thereby making it possible to obtain a given level of optical glucose-measurement sensitivity with a relatively small volume of blood. There is an additional requirement that the hills be wide enough that a sufficient amount of light can propagate into them and, after reflection, can propagate out of them. The method described in the cited prior article produces a texture comprising cones and pillars that conform to the average-distance and aspect-ratio requirements. However, a significant fraction of the cones and pillars are so narrow that not enough light can propagate along them. The improved method makes it possible to form wider cones and pillars while still satisfying the average-distance and aspect-ratio requirements. In the improved method, as in the previously reported method, multiple optical fibers are first bundled together for simultaneous texturing of their distal tips. However, prior to texturing by exposure to an oxygen beam, the tips are first coated by vapor deposition of a thin, sparse layer of aluminum: The exposure to the aluminum vapor must be short enough (typically of the order of seconds) so that the aluminum nucleates into islands separated by uncoated areas. The coated tips are textured by exposure to a directed beam of hyperthermal (kinetic energy >1 eV) oxygen atoms and/or ions in a vacuum chamber, as in the previously reported method. The aluminum islands partially shield the underlying PMMA from oxidation and erosion by the beam, so that the cones or pillars remaining after texturing are wider than they would otherwise be. To some extent, the dimensions of the hills and the distances between them can be tailored through choice of the thickness of the aluminum coat and/or the oxygen-beam fluence. The figure illustrates an example of texturing of the tip of a PMMA optical fiber without and with prior aluminum coating.

Banks, Bruce A.↗
NASA NTRS1975

Preliminary analysis of the span-distributed-load concept for cargo aircraft design

A simplified computer analysis of the span-distributed-load airplane (in which payload is placed within the wing structure) has shown that the span-distributed-load concept has high potential for application to future air cargo transport design. Significant increases in payload fraction over current wide-bodied freighters are shown for gross weights in excess of 0.5 Gg (1,000,000 lb). A cruise-matching calculation shows that the trend toward higher aspect ratio improves overall efficiency; that is, less thrust and fuel are required. The optimal aspect ratio probably is not determined by structural limitations. Terminal-area constraints and increasing design-payload density, however, tend to limit aspect ratio.

Whitehead, A. H., Jr.↗
DOE OSTI2015

Boundary conditions for the solution of the three-dimensional Poisson equation in open metallic enclosures

Numerical solution of the Poisson equation in metallic enclosures, open at one or more ends, is important in many practical situations, such as high power microwave or photo-cathode devices. It requires imposition of a suitable boundary condition at the open end. In this paper, methods for solving the Poisson equation are investigated for various charge densities and aspect ratios of the open ends. It is found that a mixture of second order and third order local asymptotic boundary conditions is best suited for large aspect ratios, while a proposed non-local matching method, based on the solution of the Laplace equation, scores well when the aspect ratio is near unity for all charge density variations, including ones where the centre of charge is close to an open end or the charge density is non-localized. The two methods complement each other and can be used in electrostatic calculations where the computational domain needs to be terminated at the open boundaries of the metallic enclosure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗
NASA NTRS2017

Global Statistics of Microphysical Properties of Cloud-Top Ice Crystals

Ice crystals in clouds are highly complex. Their sizes, macroscale shape (i.e., habit), mesoscale shape (i.e., aspect ratio of components) and microscale shape (i.e., surface roughness) determine optical properties and affect physical properties such as fall speeds, growth rates and aggregation efficiency. Our current understanding on the formation and evolution of ice crystals under various conditions can be considered poor. Commonly, ice crystal size and shape are related to ambient temperature and humidity, but global observational statistics on the variation of ice crystal size and particularly shape have not been available. Here we show results of a project aiming to infer ice crystal size, shape and scattering properties from a combination of MODIS measurements and POLDER-PARASOL multi-angle polarimetry. The shape retrieval procedure infers the mean aspect ratios of components of ice crystals and the mean microscale surface roughness levels, which are quantifiable parameters that mostly affect the scattering properties, in contrast to a habit. We present global statistics on the variation of ice effective radius, component aspect ratio, microscale surface roughness and scattering asymmetry parameter as a function of cloud top temperature, latitude, location, cloud type, season, etc. Generally, with increasing height, sizes decrease, roughness increases, asymmetry parameters decrease and aspect ratios increase towards unity. Some systematic differences are observed for clouds warmer and colder than the homogeneous freezing level. Uncertainties in the retrievals will be discussed. These statistics can be used as observational targets for modeling efforts and to better constrain other satellite remote sensing applications and their uncertainties.

Van Diedenhoven, Bastiaan↗
NASA NTRS2005

Texturing Blood-Glucose-Monitoring Optics Using Oxygen Beams

A method has been invented for utilizing directed, hyperthermal oxygen atoms and ions for texturing tips of polymeric optical fibers or other polymeric optical components for use in optical measurement of concentration of glucose in blood. The required texture of the sensory surface of such a component amounts to a landscape of microscopic hills having high aspect ratios (hills taller than they are wide), with an average distance between hills of no more than about 5 m. This limit on the average distance between hills is chosen so that blood cells (which are wider) cannot enter the valleys between the hills, where they could obstruct optical sensing of glucose in the blood plasma. On the other hand, the plasma is required to enter the valleys, and a high aspect ratio is intended to maximize the hillside and valley surface area in contact with the plasma, thereby making it possible to obtain a given level of optical glucose- measurement sensitivity with a relatively small volume of blood. The present method of texturing by use of directed, hyperthermal (particle energy >1 eV) oxygen atoms and ions stands in contrast to a prior method of texturing by use of thermal monatomic oxygen characterized by a temperature of the order of 0.5 eV. The prior method yields low-aspect- ratio (approximately hemispherical) craters that are tens of microns wide . too wide to exclude blood cells. The figure schematically depicts parts of a typical apparatus for texturing according to the present method. One or more polymeric optical components to be textured (e.g., multiple optical fibers bundled together for simultaneous processing) are mounted in a vacuum chamber facing a suitable ion- or atom-accelerating device capable of generating a beam of oxygen atoms and/or ions having kinetic energies >1 eV. Typically, such a device includes a heated cathode, in which case it is desirable to interpose a water-cooled thermal-radiation shield to prevent melting of the polymeric component(s) to be textured. In operation, the chamber is evacuated to a pressure .10.5 torr (less than or equal to approximately 1.3 mPa), then the beam is turned on. The resulting texture is characterized by approximately conical hills having aspect ratios greater than 1. In experiments, it was demonstrated that separations between adjacent hills can be made .1 m and that the separations and heights of the hills can be varied by varying the fluence of monatomic oxygen and/or oxygen ions.

Banks, Bruce↗
DOE OSTI2022

Drift kinetic theory of the NTM magnetic islands in a finite beta general geometry tokamak plasma

In (Imada et al 2019 Nucl. Fusion 59 046016 and references therein) a new 4D drift kinetic nonlinear theory, valid in the limit of a low beta, small inverse aspect ratio, circular cross section, toroidal geometry, to describe the plasma response to the neoclassical tearing mode (NTM) magnetic perturbation is derived. In (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) this theory is reduced in a low collisionality limit, which allows a dimensionality reduction to a 3D problem to efficiently resolve the collisional dissipation layer in the vicinity of the trapped-passing boundary. (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) adopts an improved model for the magnetic drift frequency, which reduces the threshold magnetic island half-width from $8.73\rho_{b i}$, where $\rho_{b i}$ is the trapped ion banana orbit width, to $1.46 \rho_{b i}$, making it in closer agreement with experimental observations for the large aspect ratio tokamak equilibrium. In the present paper, the theory is extended to a high beta, arbitrary tokamak geometry to capture the plasma shaping effects on the NTM threshold physics with the focus on the non-zero triangularity discharges that are known to have a strong impact on the plasma MHD stability. First, it is found that the higher triangularity plasma is more prone to NTMs which is in agreement with the 2/1 tearing mode onset relative frequency measurements in DIII-D. Second, the NTM threshold dependence on the tokamak inverse aspect ratio obtained in (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) is refined and extended to a finite aspect ratio limit. Third, the NTM threshold dependence on poloidal beta is obtained and benchmarked against the EAST threshold island width measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Pore-scale and Continuum Simulations of Solute Transport Micromodel Benchmark Experiments

Four sets of micromodel nonreactive solute transport experiments were conducted with flow velocity, grain diameter, pore-aspect ratio, and flow focusing heterogeneity as the variables. The data sets were offered to pore-scale modeling groups to test their simulators. Each set consisted of two learning experiments, for which all results was made available, and a challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing. Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice-Boltzmann (LB) approach, and one employed a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used these experiments to appropriately discretize the grid representations. The continuum model use published non-linear relations between transverse dispersion coefficients and Peclet numbers to compute the required dispersivity input values. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values and, resulting in less dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models needed up to several days on supercomputers to resolve the more complex problems.

benchmark↗