SEARCH · Search NASA

Results for “dual numbers”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Simultaneous fermion and exciton condensations from a model Hamiltonian

Fermion-exciton condensation in which both fermion-pair (i.e., superconductivity) and exciton condensations occur simultaneously in a single coherent quantum state has recently been conjectured to exist. Here, we capture the fermion-exciton condensation through a model Hamiltonian that can recreate the physics of this new class of highly correlated condensation phenomena. We demonstrate that the Hamiltonian generates the large-eigenvalue signatures of fermion-pair and exciton condensations for a series of states with increasing particle numbers. The results confirm that the dual-condensate wave function arises from the entanglement of fermion-pair and exciton wave functions, which we previously predicted in the thermodynamic limit. Furthermore, this model Hamiltonian—generalizing well-known model Hamiltonians for either superconductivity or exciton condensation—can explore a wide variety of condensation behavior. It provides significant insights into the required forces for generating a fermion-exciton condensate, which will likely be invaluable for realizing such condensations in realistic materials with applications from superconductors to excitonic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1982

Galileo spacecraft power distribution and autonomous fault recovery

There is a trend in current spacecraft design to achieve greater fault tolerance through the implemenation of on-board software dedicated to detecting and isolating failures. A combination of hardware and software is utilized in the Galileo power system for autonomous fault recovery. Galileo is a dual-spun spacecraft designed to carry a number of scientific instruments into a series of orbits around the planet Jupiter. In addition to its self-contained scientific payload, it will also carry a probe system which will be separated from the spacecraft some 150 days prior to Jupiter encounter. The Galileo spacecraft is scheduled to be launched in 1985. Attention is given to the power system, the fault protection requirements, and the power fault recovery implementation.

Detwiler, R. C.↗

NASA NTRS1992

Space Transfer Concepts and Analyses for Exploration Missions. Technical Directive 12: Beamed Power Systems Study

Parametric models were constructed for Earth-based laser powered electric orbit transfer from low Earth orbit to geosynchronous orbit. These models were used to carry out performance, cost/benefit, and sensitivity analyses of laser-powered transfer systems including end-to-end life cycle cost analyses for complete systems. Comparisons with conventional orbit transfer systems were made indicating large potential cost savings for laser-powered transfer. Approximate optimization was done to determine best parameter values for the systems. Orbit transfer flights simulations were conducted to explore effects of parameters not practical to model with a spreadsheet. The simulations considered view factors that determine when power can be transferred from ground stations to an orbit transfer vehicle and conducted sensitivity analyses for numbers of ground stations, Isp including dual-Isp transfers, and plane change profiles. Optimal steering laws were used for simultaneous altitude and plane change. Viewing geometry and low-thrust orbit raising were simultaneously simulated. A very preliminary investigation of relay mirrors was made.

Eder, D.↗

NASA NTRS2002

Unsteady Computations of a Jet in a Crossflow with Ground Effect

A numerical study of a jet in crossflow with ground effect is conducted using OVERFLOW with dual time-stepping and low Mach number preconditioning. The results of the numerical study are compared to an experiment to show that the numerical methods are capable of capturing the dominant features of the flow field as well as the unsteadiness associated with the ground vortex.

Pandya, Shishir↗

NASA NTRS1980

Technologies Involved in Configuring an Advanced Earth-to-Orbit Transport for Low Structural Mass

The current space shuttle is expected to adequately meet Government and industry needs for the transport of cargo to and from orbit well into the 1990's. However, continual study of potential follow-on shuttle systems is necessary and desirable in order to complement ongoing research in materials, structures, propulsion, aerodynamics, and other related areas. By studying alternate systems well in advance, it will be possible to explore the various technologies and develop those for which there is the greatest apparent payoff. In this paper a single-stage Earth-to-orbit transport designed for delivery of approximately 29,500 kg (65,000 lb) payload will be described. The vehicle, which takes off vertically and lands horizontally, is 60 m (197 feet) long and weighs approximately 1.8 Gg (4 M lb) at liftoff. In the interest of weight reduction, a simple body of revolution is utilized for the main body shell. In this design the main propulsion tanks serve as a primary load-carrying structure. Further, in order to minimize structural mass, the cargo bay is located between two of the main propellant tanks. The cargo volume, at 396 cu m (14,000 cu feet), exceeds that provided by the shuttle; but the bay itself is nonconforming in shape - being approximately 10 m (32 feet) in diameter by 5 m (17 feet) long. Dual-fuel propulsion is employed, since a number of studies have shown that (though lowering performance) the operation of hydrocarbon (RP) engines in parallel with LOX/LH2 engines results in a net reduction in the vehicle's physical size and structural mass. Other weight-saving features entail the extensive use of honeycomb sandwiches, advanced materials, and advanced fabrication techniques. The vehicle presented is utilized only as a means to study and identify various technologies needed in order to develop a low mass Earth-to-orbit transportation system for the future. The conclusion of this study is that vehicle geometry and structural/materials technology are critical to the development of efficient single-stage Earth-to-orbit transports.

MacConochie, Ian O.↗

DOE OSTI2023

Systems and methods for dual comb spectroscopy

A frequency-measurement method uses a dual frequency-comb spectrometer as an optical wavemeter to measure the frequency of a reference laser that is used to frequency-stabilize the spectrometer. The method includes measuring a walking rate of center bursts in a sequence of interferograms recorded by the spectrometer, determining a number of teeth in each of a plurality of Nyquist windows formed by the dual frequency-comb spectrometer, and determining a Nyquist number of the one Nyquist window covering the laser frequency. The reference laser frequency can then be determined from the number of teeth in each Nyquist window, the Nyquist number, and the comb spacing of either one of the two frequency combs of the dual frequency-comb spectrometer. The reference laser frequency does not need to be measured with a separate wavemeter, or calibrated with respect to a known atomic or molecular transition.

Rieker, Gregory B.↗

NASA NTRS1986

Radar investigation of asteroids

The number of radar detected asteroids has climbed from 6 to 40 (27 mainbelt plus 13 near-Earth). The dual-circular-polarization radar sample now comprises more than 1% of the numbered asteroids. Radar results for mainbelt asteroids furnish the first available information on the nature of these objects at macroscopic scales. At least one object (2 Pallas) and probably many others are extraordinarily smooth at centimeter-to-meter scales but are extremely rough at some scale between several meters and many kilometers. Pallas has essentially no small-scale structure within the uppermost several meters of the regolith, but the rms slope of this regolith exceeds 20 deg., much larger than typical lunar values (approx. 7 deg.). The origin of these slopes could be the hypervelocity impact cratering process, whose manifestations are likely to be different on low-gravity, low-radius-of-curvature objects from those on the terrestrial planets. The range of mainbelt asteroid radar albedoes is very broad and implies big variations in regolith porosity or metal concentration, or both. The highest albedo estimate, for 16 Psyche, is consistent with a surface having porosities typical of lunar soil and a composition nearly completely metallic. Therefore, Psyche might be the collisionally stripped core of a differentiated small plant, and might resemble mineralogically the parent bodies of iron meteorites.

Ostro, S. J.↗

DOE OSTI2021

Topological or rotational non-Abelian gauge fields from Einstein-Skyrme holography

We report analytically known states at non-zero temperature which may serve as a powerful tool to reveal common topological and thermodynamic properties of systems ranging from the QCD phase diagram to topological phase transitions in condensed matter materials. In the holographically dual gravity theory, these are analytic solutions to a five-dimensional non-linear-sigma (Skyrme) model dynamically coupled to Einstein gravity. This theory is shown to be holographically dual to $ \mathcal{N} $ = 4 Super-Yang-Mills theory coupled to an SU(2)-current. All solutions are fully backreacted asymptotically Anti-de Sitter (AdS) black branes or holes. One family of global AdS black hole solutions contains non-Abelian gauge field configurations with positive integer Chern numbers and finite energy density. Larger Chern numbers increase the Hawking-Page transition temperature. In the holographically dual field theory this indicates a significant effect on the deconfinement phase transition. Black holes with one Hawking temperature can have distinct Chern numbers, potentially enabling topological transitions. A second family of analytic solutions, rotating black branes, is also provided. These rotating solutions induce states with propagating charge density waves in the dual field theory. We compute the Hawking temperature, entropy density, angular velocity and free energy for these black holes/branes. These correspond to thermodynamic data in the dual field theory. For these states the energy-momentum tensor, (non-)conserved current, and topological charge are interpreted.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1997

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗

NASA NTRS2003

Implementation of Preconditioned Dual-Time Procedures in OVERFLOW

Preconditioning methods have become the method of choice for the solution of flowfields involving the simultaneous presence of low Mach and transonic regions. It is well known that these methods are important for insuring accurate numerical discretization as well as convergence efficiency over various operating conditions such as low Mach number, low Reynolds number and high Strouhal numbers. For unsteady problems, the preconditioning is introduced within a dual-time framework wherein the physical time-derivatives are used to march the unsteady equations and the preconditioned time-derivatives are used for purposes of numerical discretization and iterative solution. In this paper, we describe the implementation of the preconditioned dual-time methodology in the OVERFLOW code. To demonstrate the performance of the method, we employ both simple and practical unsteady flowfields, including vortex propagation in a low Mach number flow, flowfield of an impulsively started plate (Stokes' first problem) arid a cylindrical jet in a low Mach number crossflow with ground effect. All the results demonstrate that the preconditioning algorithm is responsible for improvements to both numerical accuracy and convergence efficiency and, thereby, enables low Mach number unsteady computations to be performed at a fraction of the cost of traditional time-marching methods.

Pandya, Shishir A.↗

NASA NTRS1941

The Selection of Propellers for High Thrust at Low Airspeed

Problem of improving thrust at low speeds is primarily one of reducing angle of attack of operation of sections to improve L/D or reducing blade helix angle. An analysis, based on recent propeller data, is presented for determining improvements in thrust or efficiency which could be obtained by increased number of blades, increased blade width, increased diameter, dual rotation, and two-speed gearing. All methods were found very effective, particularly two-speed gearing.

PROPELLERS - CURTISS↗

NASA NTRSDate not supplied

How Can We Efficiently Build A Spacecraft That Has Longevity? Experiences From the GSFC Perspective to Inform Lunar Exploration and Science Orbiter’s Architecture

Recent successes in NASA planetary science missions has shown that long-lived missions can yield significant science return. These missions, despite their longevity, were not planned to operate beyond their deign life. For example, the Lunar Reconnaissance Orbiter has a design life of 2-3 years, and yet we are entering its 14th year on orbit. Spacecraft longevity in practice has been related to: 1) mission class; 2) did we test it long enough and resolve all the anomalies to be beyond the early failure curve; 3) consumables budgets, and 4) how components are de-signed and tested. Mission class has been perceived as one of the primary ways to drive longevity. But higher mission class is a significant cost and mission driver. Parts selection only from the limited military standard parts and extensive parts qualification adds to development time and cost. Largely redundant (often erroneously interpreted as fully redundant) adds to launch mass and increases testing complexity (which may reduce the amount of testing in the nominal con-figuration). Lower Risk Posture (reflected in a higher mission class) drives significant additional processes and quality assurance analyses. Higher mission class also drives significantly greater sparing and life testing costs. This discussion focuses on whether this is indeed the best way to achieve longevity efficiently or whether there are more efficient ways to achieve longevity. Empirical evidence shows that lower mission classes that implement effective risk reduction and selective redundancy generally results in long life at a lower Spacecraft bus cost. By evaluating redundancy careful-ly, spacecraft cost can be lowered which can enable a more capable payload. Risk of the mission is highly dependent on the complexity of the mission and is only loosely dependent on Class of Mission. High complexity missions can have many single point failures and require the development of new technology, while lower class mission can have lower risk by baselining or incorporating: (selective) redundancy, fault-tolerant design, design for minimum risk, ability to reset, and/or design for graceful degradation. The team met with Goddard Space Flight Center Space Systems Mission Operations leaders to discuss which avionics have proven in flight to be the most reliable and which have had lifetime issues. The paper proposes a list of components to focus on for redundancy for a long-lived lunar mission. Reliability numbers are presented for both single string and dual string missions. The history of life-time performance versus planned mission life according to mission class is presented. A mission with well thought out selective redundancy and effective on the ground test program, can be expected to last well beyond its mission lifetime and provided enhanced return for the community. This approach minimizes project expenditures that do not retire significant risk and allows the project to focus risk mitigation efforts on those risks that will have a significant likelihood of threatening mission success. This is aided by keeping the amount of technology maturation and mission complexity low, while focusing effort on ensuring all component stress-ing parameters well within the bounds of their capabilities.

Lessons Learned↗

DOE OSTI2024

Uncertainty Quantification and Sensitivity Analysis of Non-Nuclear Advanced Controls Testbed Reactor Mockup

The research presented in this report describes our progress in applying stochastic methods and uncertainty quantification, parametric study, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated with Griffin, with the goal of developing a reduced order (surrogate) model which can be rapidly sampled while perturbing multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums, which are controlled by a hybrid proportional, integral and derivative (PID) controller, actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time and $\beta$ from a transient Griffin simulation executed with the improved quasi-static method to provide the kinetic parameters) as inputs to functions which control the CD rotation angle. Using a number of stochastic method approaches, we developed a dual purpose training-surrogate model of the NTP system using polynomial regression. The trained model can be rapidly sampled while simultaneously perturbing various input parameters of the model, such as coefficients on the PID control, or temperature (directly affect the neutron cross section). The surrogate model delivers accurate results orders-of-magnitude faster (minutes, not days) than the base model. Once the base model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs and their effect on the output. For example, coefficients used in the PID control system may vary due to some physical interference, or there may be uncertainty in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters and operational boundaries can be determined. The goal of this work is to support development of an advanced control system to operate CDs in a functioning NTP system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS2020

Cramer-Rao Lower Bound for SoOp-R-Based Root-Zone Soil Moisture Remote Sensing

Signals of opportunity (SoOp) reflectometry (SoOp-R) is a maturing field for geophysical remote sensing as evidenced by the growing number of airborne and spaceborne experiments. As this approach receives more attention, it is worth analyzing SoOp-R’s capabilities to retrieve subsurface soil moisture (SM) by leveraging communication and navigation satellite transmitters. In this research, the CRLB is used to identify the effects of variable SoOp-R parameters on the best achievable estimation error for root-zone soil moisture (RZSM). This study investigates the use of multiple frequency, polarization, and incidence angle measurement configurations on a two-layered dielectric profile. The results also detail the effects of variable SM conditions on the capability of SoOp-R systems to predict subsurface SM. The most prevalent observation is the importance of using at least two frequencies to limit uncertainties from subsurface SM estimates. If at least two frequencies are used, the CRLB of a profile is retrievable within the root-zone depending on the surface SM content as well as the number of independent measurements of the profile. For a depth of 30 cm, it is observed that a CRLB corresponding to 4% RZSM estimation accuracy is achievable with as few as 2 dual-frequency-based SoOp-R measurements. For this depth, increasing number of measurements provided by polarization and incidence angle allow for sensing of increasingly wet SM profile structures. This study, overall, details a methodology by which SoOp-R receiver system can be designed to achieve a desired CRLB using a trade-off study between the available measurements and SM profile.

Signals of Opportunity (SoOp)↗

DOE OSTI2023

Anion-tethered Single Lithium-ion Conducting Polyelectrolytes through UV-induced Free Radical Polymerization for Improved Morphological Stability of Lithium Metal Anodes

Single Li + ion conducting polyelectrolytes (SICs), which feature covalently tethered counter-anions along their backbone, have the potential to mitigate dendrite formation by reducing concentration polarization and preventing salt depletion. However, due to their low ionic conductivity and complicated synthetic procedure, the successful validation of these claimed advantages in lithium metal (Li 0 ) anode batteries remains limited. In this study, we fabricated a SIC electrolyte using a single-step UV polymerization approach. The resulting electrolyte exhibited a high Li + transference number (t + ) of 0.85 and demonstrated good Li + conductivity (6.3×10 -5 S/cm at room temperature), which is comparable to that of a benchmark dual ion conductor (DIC, 9.1×10 -5 S/cm). Benefitting from the high transference number of SIC, it displayed a three-fold higher critical current density (2.4 mA/cm 2 ) compared to DIC (0.8 mA/cm 2 ) by successfully suppressing concentration polarization-induced short-circuiting. Additionally, the t + significantly influenced the deposition behavior of Li0, with SIC yielding a uniform, compact, and mosaic-like morphology, while the low t + DIC resulted in a porous morphology with Li 0 whiskers. Using the SIC electrolyte, Li 0 ||LiFePO 4 cells exhibited stable operation for 4500 cycles with 70.5 % capacity retention at 22 °C.

25 ENERGY STORAGE↗

DOE OSTI2023

Anion-tethered Single Lithium-ion Conducting Polyelectrolytes through UV-induced Free Radical Polymerization for Improved Morphological Stability of Lithium Metal Anodes

Single Li + ion conducting polyelectrolytes (SICs), which feature covalently tethered counter-anions along their backbone, have the potential to mitigate dendrite formation by reducing concentration polarization and preventing salt depletion. However, due to their low ionic conductivity and complicated synthetic procedure, the successful validation of these claimed advantages in lithium metal (Li 0 ) anode batteries remains limited. In this study, we fabricated a SIC electrolyte using a single-step UV polymerization approach. The resulting electrolyte exhibited a high Li + transference number (t + ) of 0.85 and demonstrated good Li + conductivity (6.3×10 -5 S/cm at room temperature), which is comparable to that of a benchmark dual ion conductor (DIC, 9.1×10 -5 S/cm). Benefitting from the high transference number of SIC, it displayed a three-fold higher critical current density (2.4 mA/cm 2 ) compared to DIC (0.8 mA/cm 2 ) by successfully suppressing concentration polarization-induced short-circuiting. Additionally, the t + significantly influenced the deposition behavior of Li 0 , with SIC yielding a uniform, compact, and mosaic-like morphology, while the low t + DIC resulted in a porous morphology with Li 0 whiskers. In conclusion, using the SIC electrolyte, Li 0 ||LiFePO 4 cells exhibited stable operation for 4500 cycles with 70.5 % capacity retention at 22 °C.

25 ENERGY STORAGE↗

NASA NTRS2010

Reflector Surface Error Compensation in Dual-Reflector Antennas

By probing the field on a small subreflector at a minimal number of points, the main reflector surface errors can be obtained and subsequently used to design a phase-correction subreflector that can compensate for main reflector errors. The compensating phase-error profile across the subreflector can be achieved either by a surface deformation or by the use of an array of elements such as patch antennas that can cause a phase shift between the incoming and outgoing fields. The second option is of primary interest here, but the methodology can be applied to either case. The patch array is most easily implemented on a planar surface. Therefore, the example of a flat subreflector and a parabolic main reflector (a Newtonian dual reflector system) is considered in this work. The subreflector is assumed to be a reflector array covered with patch elements. The phase variation on a subreflector can be detected by a small number of receiving patch elements (probes). By probing the phase change at these few selected positions on the subreflector, the phase error over the entire surface can be recovered and used to change the phase of all the patch elements covering the subreflector plane to compensate for main reflector errors. This is accomplished by using a version of sampling theorem on the circular aperture. The sampling is performed on the phase-error function on the circular aperture of the main reflector by a method developed using Zernike polynomials. This method is based upon and extended from a theory previously proposed and applied to reflector aperture integration. This sampling method provides for an exact retrieval of the coefficients of up to certain orders in the expansion of the phase function, from values on a specifically calculated set of points in radial and azimuthal directions in the polar coordinate system, on the circular reflector aperture. The corresponding points on the subreflector are then obtained and, by probing the fields at these points, a set of phase values is determined that is then transferred back to the main reflector aperture for recovering the phase function. Once this function is recovered, the corresponding phase function on the subreflector is calculated and used to compensate for main reflector surface errors. In going back and forth between sub and main reflectors, geometrical (ray) optics is employed, which even though it ignores edge diffraction and other effects, is shown to be accurate for phase recovery.

Jamnejad, Vahraz↗

NASA NTRS2018

Generation of Performance Model for the Aeolian Wind Tunnel (AWT) Rotor at Reduced Pressure

The NASA Jet Propulsion Laboratory (JPL) designed the Mars Helicopter (MH) in collaboration with AeroVironment Inc., NASA Ames Research Center, and NASA Langley Research Center to explore the possibility of a vertical takeoff and landing (VTOL) Unmanned Aerial Vehicle (UAV) for flight on Mars. A 40-inch-diameter Aeolian Wind Tunnel (AWT) rotor, roughly approximating the proposed MH design by JPL, was tested in forward flight at Mars atmospheric pressure at the NASA Ames Planetary Aeolian Laboratory (PAL) in support of MH research efforts. This report describes the generation of the rotor model used to correlate with that experimental effort as reported by Ament and Koning. The 40-inch-diameter rotor was 3D-scanned and transformed into an airfoil deck. The scanned rotor airfoil sections are analyzed using C81 Generator (C81Gen) to generate the sectional aerodynamic coefficients for comprehensive analyses. A mid-fidelity computational fluid dynamics (CFD) simulation using Rotorcraft CFD (RotCFD) is pursued to efficiently estimate rotor hover and forward flight performance. Simulations at two pressures, 7 mbar (approximate Martian atmospheric pressure) and 1018 mbar (1 atmosphere), are performed to gain an understanding of the performance differences and Reynolds number effects observed. Experimental 1-atmosphere thrust for single- and dual-rotor isolated hover cases correlate well with the modeled rotor. Performance results at reduced pressure (7 mbar) show a drastic decrease in lift for equivalent RPMs tested at 1 atmosphere. Although this is primarily due to pressure reduction, Reynolds number effects also contribute to this decrease, as airfoil lift and drag coefficients are affected when compared with 1-atmosphere results. Further, simulated rotor power coefficient shows drastic increases at reduced pressures, attributed to laminar boundary layer separation, as described in Koning et al. for the MH rotor analysis. PAL experimental Martian Surface Wind Tunnel (MARSWIT) results are presented in the paper by Ament and Koning. The very low Reynolds number range is currently not well understood and presents various challenges for both experimentation and simulation.

Koning, Witold J. F.↗