Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “domination number”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Wind Tunnel Investigation of the Effects of Surface Porosity and Vertical Tail Placement on Slender Wing Vortex Flow Aerodynamics at Supersonic Speeds

A wind tunnel experiment was conducted in the NASA Langley Research Center (LaRC) Unitary Plan Wind Tunnel (UPWT) to determine the effects of passive surface porosity and vertical tail placement on vortex flow development and interactions about a general research fighter configuration at supersonic speeds. Optical flow measurement and flow visualization techniques were used that featured pressure sensitive paint (PSP), laser vapor screen (LVS), and schlieren, These techniques were combined with conventional electronically-scanned pressure (ESP) and six-component force and moment measurements to quantify and to visualize the effects of flow-through porosity applied to a wing leading edge extension (LEX) and the placement of centerline and twin vertical tails on the vortex-dominated flow field of a 65 cropped delta wing model. Test results were obtained at free-stream Mach numbers of 1.6, 1.8, and 2.1 and a Reynolds number per foot of 2.0 million. LEX porosity promoted a wing vortex-dominated flow field as a result of a diffusion and weakening of the LEX vortex. The redistribution of the vortex-induced suction pressures contributed to large nose-down pitching moment increments but did not significantly affect the vortex-induced lift. The trends associated with LEX porosity were unaffected by vertical tail placement. The centerline tail configuration generally provided more stable rolling moments and yawing moments compared to the twin wing-mounted vertical tails. The strength of a complex system of shock waves between the twin tails was reduced by LEX porosity.

Erickson, Gary E.↗
NASA NTRS2009

Automation of a Versatile Crane (the LSMS) for Lunar Outpost Construction, Maintenance and Inspection

Devices for lifting, translating and precisely placing payloads are critical for efficient Earth-based construction operations. Both recent and past studies have demonstrated that devices with similar functionality will be needed to support lunar outpost operations. Although several designs have been developed for Earth based applications, these devices lack unique design characteristics necessary for transport to and use on the harsh lunar surface. These design characteristics include: a) lightweight components, b) compact packaging for launch, c) automated deployment, d) simple in-field reconfiguration and repair, and e) support for tele-operated or automated operations. Also, because the cost to transport mass to the lunar surface is very high, the number of devices that can be dedicated to surface operations will be limited. Thus, in contrast to Earth-based construction, where many single-purpose devices dominate a construction site, a lunar outpost will require a limited number of versatile devices that provide operational benefit from initial construction through sustained operations. The first generation test-bed of a new high performance device, the Lunar Surface Manipulation System (LSMS) has been designed, built and field tested. The LSMS has many unique features resulting in a mass efficient solution to payload handling on the lunar surface. Typically, the LSMS device mass is estimated at approximately 3% of the mass of the heaviest payload lifted at the tip, or 1.8 % of the mass of the heaviest mass lifted at the elbow or mid-span of the boom for a high performance variant incorporating advanced structural components. Initial operational capabilities of the LSMS were successfully demonstrated during field tests at Moses Lake, Washington using a tele-operated approach. Joint angle sensors have been developed for the LSMS to improve operator situational awareness. These same sensors provide the necessary information to support fully automated operations, greatly expanding the operational versatility of the LSMS. This paper develops the equations describing the forward and inverse relation between LSMS joint angles and Cartesian coordinates of the LSMS tip. These equations allow a variety of schemes to be used to maneuver the LSMS to optimize the maneuver. One such scheme will be described in detail that eliminates undesirable swinging of the payload at the conclusion of a maneuver, even when the payload is suspended from a passive rigid link. The swinging is undesirable when performing precision maneuvers, such as aligning an object for mating or positioning a camera. Use of the equations described here enables automated control of the LSMS greatly improving its operational versatility.

Doggett, William R.↗
DOE OSTI2024

Discoveries in Blast-Driven Turbulence of Astrophysical Relevance

The fluid mixing caused by variable-density instabilities is important in a wide variety of scenarios from ocean mixing and astrophysical phenomena to nuclear fusion techniques and atomic weapons. This thesis explores the mixing resulting from a specific instability known as the Blast Driven Instability (BDI). This work investigates the variable density mixing in an explosively driven environment due to the fluid instabilities at the material interfaces. Specifically, diverging Richtmyer-Meshkov (impulsive-acceleration environment) and Rayleigh-Taylor (variable-acceleration environment) instabilities (present in supernova and inertial confinement fusion) are studied using advanced high-speed diagnostics in carefully designed laboratory experiments. The BDI morphology is presented through a time development of Mie scattering images, and steps through the parameter space (varying density ratio and driver speed), highlighting the development of the structures that form during mixing. A scaling criterion is used to relate the two systems of vastly different spatiotemporal scales. Velocity fields in the BDI have been captured for the first time using the high temporal resolution PIV technique. Subsequent analysis of the dynamics of the instability from the velocity fields illustrates the distribution of kinetic energy, the transition to turbulence, and the characteristic growth of the instability are discussed. This study furthers understanding of how blast-driven instability pertains to supernova and inertial confinement fusion science. The morphology of the BDI has been characterized for the first time. This work steps through the parameter space covered in Mie scattering experiments, and how the different parameters contribute to development of structures and mixing. It also examines a scaling of the Atwood number for expanding predictive capabilities to other experimental conditions and simulations. The first collection of velocity fields acquired for the BDI are recorded, and subsequent analysis evaluating the distribution of kinetic energy throughout space and time for two density ratios from the overall parameter space, as well as the transition to turbulence, estimated from a Reynolds number calculated based on momentum mixing are all presented. This information is useful in advancing the development of models to predict physics of high energy density applications where experiments are not always readily available. This research has successfully demonstrated understanding for the time criteria defining regimes where the shock driven (Richtmyer-Meshkov instability) and the buoyancy driven (Rayleigh-Taylor instability) dominates through a parametric study of density variation (Atwood number) and driver speed (Mach number). All this furthers understanding of how the BDI pertains to SN and ICF.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

(abstract) The Impactor Flux in the Pluto-Charon System

The impactor flux on Pluto and Charon consists of long-period comets from the Oort cloud passing through the region, and intermediate-period comets evolving inward to (or outward from) short-period orbits. The source of the intermediate-period comets is likely an extended disk of remnant planetesimals beyond the orbit of Neptune, commonly referred to as the Kuiper belt, though the Oort cloud may be an additional source of some short-period comets. The Pluto-Charon system actually circulates within the Kuiper belt, and with an aphelion of 49.3 AU, moves through a region where remnant planetesimal orbits would likely be stable over the age of the solar system. The recent discovery of two objects at Kuiper belt distances, 1992 QB 1 and 1993 FW, and the fact that the orbit determined for 1992 QB 1 is likely of low eccentricity and low inclination, have provided strong observational evidence for the Kuiper belt. Weisman et al (1989) showed that, for the expected number of objects in the Kuiper belt necessary to provide the short-period comet flux, cometary cratering on Pluto and Charon is dominated by Kuiper belt comets. We will provide updated estimates of the numbers, based on an improved understanding of Oort cloud and Kuiper belt dynamics since 1989, and on improved estimates for the radii and masses for Pluto and Charon. Because Pluto and Charon actually penetrate the stable region of the Kuiper belt beyond 40 AU, expected cratering rates for the pair are substantially higher than for Triton, which is located at the inner edge of the Kuiper belt. However, the Pluto-Charon system is tidally evolved and this may have led to resurfacing events at various times in its history. Thus, detailed predictions of the integrated crater density or comparisons with Triton are not very meaningful without a knowledge of that history. Differences in the suspected compositions of Pluto and Charon may also lead to different surface rheologies and thus, different retention times for cratered surfaces.

Pluto↗
NASA NTRS2002

Development of an Efficient Binaural Simulation for the Analysis of Structural Acoustic Data

Binaural or "virtual acoustic" representation has been proposed as a method of analyzing acoustic and vibroacoustic data. Unfortunately, this binaural representation can require extensive computer power to apply the Head Related Transfer Functions (HRTFs) to a large number of sources, as with a vibrating structure. This work focuses on reducing the number of real-time computations required in this binaural analysis through the use of Singular Value Decomposition (SVD) and Equivalent Source Reduction (ESR). The SVD method reduces the complexity of the HRTF computations by breaking the HRTFs into dominant singular values (and vectors). The ESR method reduces the number of sources to be analyzed in real-time computation by replacing sources on the scale of a structural wavelength with sources on the scale of an acoustic wavelength. It is shown that the effectiveness of the SVD and ESR methods improves as the complexity of the source increases. In addition, preliminary auralization tests have shown that the results from both the SVD and ESR methods are indistinguishable from the results found with the exhaustive method.

Lalime, Aimee L.↗
NASA NTRS2011

Consequences of the Breakout Model for Particle Acceleration in CMEs and Flares

The largest and most efficient particle accelerators in the solar system are the giant events consisting of a fast coronal mass ejection (CME) and an intense X-class solar flare. Both flares and CMEs can produce l0(exp 32) ergs or more in nonthermal particles. Two general processes are believed to be responsible: particle acceleration at the strong shock ahead of the CME, and reconnection-driven acceleration in the flare current sheet. Although shock acceleration is relatively well understood, the mechanism by which flare reconnection produces nonthermal particles is still an issue of great debate. We address the question of CME/flare particle acceleration in the context of the breakout model using 2.5D MHD simulations with adaptive mesh refinement (AMR). The AMR capability allows us to achieve ultra-high numerical resolution and, thereby, determine the detailed structure and dynamics of the flare reconnection region. Furthermore, we employ newly developed numerical analysis tools for identifying and characterizing magnetic nulls, so that we can quantify accurately the number and location of magnetic islands during reconnection. Our calculations show that flare reconnection is dominated by the formation of magnetic islands. In agreement with many other studies, we find that the number of islands scales with the effective Lundquist number. This result supports the recent work by Drake and co-workers that postulates particle acceleration by magnetic islands. On the other hand, our calculations also show that the flare reconnection region is populated by numerous shocks and other indicators of strong turbulence, which can also accelerate particles. We discuss the implications of our calculations for the flare particle acceleration mechanism and for observational tests of the models.

Antiochos, S. K.↗
DOE OSTI2023

A helical actuator driven by biased SMA: design, model, and experiment

A helical actuator driven by biased shape memory alloy (SMA) patterns embedded into a soft composite ribbon base is presented in this work. Instead of common U-shape SMA wires, a single SMA wire is woven into planar patterns, which enable helical deformation of the composite ribbon from an initially flat geometry. An analytical static model is established for accurate and rapid prediction of the helical reconfiguration arising from the shape memory effect of woven SMA patterns, followed by validation of the static model using the finite element method (FEM). The finite element results are compared with the analytical solutions given by this static model, which show a high agreement. Parametric study of the influences of eight independent design variables on the dependent helical parameters, such as combined curvatures, pitches, and helical angles, is completed. It is found that the helically deformed geometry is mainly dominated by diameters, biased positions, inclined angles, and numbers of skewed segments of the SMA wire. Fabrication and in-situ experimental test of a prototype of such helical actuators qualitatively demonstrate its dramatic three-dimensional (3D) spiral reconfiguration from a two-dimensional (2D) flat ribbon. In conclusion, such SMA patterns will allow more diverse designs of soft actuators for a wider range of robotic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Contributions From Cloud Morphological Changes to the Interannual Shortwave Cloud Feedback Based on MODIS and ISCCP Satellite Observations

The surface temperature-mediated change in cloud properties, referred to as the cloud feedback, continues to dominate the uncertainty in climate projections. A larger number of contemporary global climate models (GCMs) project a higher degree of warming than the previous generation of GCMs. This greater projected warming has been attributed to a less negative cloud feedback in the Southern Ocean. Here, we apply a novel “double decomposition method” that employs the “cloud radiative kernel” and “cloud regime” concepts, to two data sets of satellite observations to decompose the interannual cloud feedback into contributions arising from changes within and shifts between cloud morphologies. Our results show that contributions from the latter to the cloud feedback are large for certain regimes. We then focus on interpreting how both changes within and between cloud morphologies impact the shortwave cloud optical depth feedback over the Southern Ocean in light of additional observations. Results from the former cloud morphological changes reveal the importance of the wind response to warming increases low- and mid-level cloud optical thickness in the same region. Results from the latter cloud morphological changes reveal that a general shift from thick storm-track clouds to thinner oceanic low-level clouds contributes to a positive feedback over the Southern Ocean that is offset by shifts from thinner broken clouds to thicker mid- and low-level clouds. Our novel analysis can be applied to evaluate GCMs and potentially diagnose shortcomings pertaining to their physical parameterizations of particular cloud morphologies.

58 GEOSCIENCES↗
DOE OSTI2023

Critical behavior in a chiral molecular model

Understanding the condensed-phase behavior of chiral molecules is important in biology as well as in a range of technological applications, such as the manufacture of pharmaceuticals. Here, we use molecular dynamics simulations to study a chiral four-site molecular model that exhibits a second-order symmetry-breaking phase transition from a supercritical racemic liquid into subcritical D-rich and L-rich liquids. We determine the infinite-size critical temperature using the fourth-order Binder cumulant, and we show that the finite-size scaling behavior of the order parameter is compatible with the 3D Ising universality class. We also study the spontaneous D-rich to L-rich transition at a slightly subcritical temperature of T = 0.985Tc, and our findings indicate that the free energy barrier for this transformation increases with system size as N2/3, where N is the number of molecules, consistent with a surface-dominated phenomenon. The critical behavior observed herein suggests a mechanism for chirality selection in which a liquid of chiral molecules spontaneously forms a phase enriched in one of the two enantiomers as the temperature is lowered below the critical point. Furthermore, the increasing free energy barrier with system size indicates that fluctuations between the L-rich and D-rich phases are suppressed as the size of the system increases, trapping it in one of the two enantiomerically enriched phases. Such a process could provide the basis for an alternative explanation for the origin of biological homochirality. We also conjecture the possibility of observing nucleation at subcritical temperatures under the action of a suitable chiral external field.

Chemistry↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Observation of Directed Flow of Hypernuclei $^3_Λ\text{H}$ and $^4_Λ\text{H}$ in $\sqrt{s_{\text{NN}}}$ = 3 GeV Au + Au Collisions at RHIC

We report here the first observation of directed flow ($v_1$) of the hypernuclei $^3_Λ\text{H}$ and $^4_Λ\text{H}$ in mid-central Au + Au collisions at $\sqrt{s_{\text{NN}}}$ = 3 GeV at RHIC. These data are taken as part of the beam energy scan program carried out by the STAR experiment. From 165 × 10 6 events in 5%–40% centrality, about 8400 $^3_Λ\text{H}$ and 5200 $^4_Λ\text{H}$ candidates are reconstructed through two- and three-body decay channels. In this work, we observe that these hypernuclei exhibit significant directed flow. Comparing to that of light nuclei, it is found that the midrapidity $v_1$ slopes of $^3_Λ\text{H}$ and $^4_Λ\text{H}$ follow baryon number scaling, implying that the coalescence is the dominant mechanism for these hypernuclei production in the 3 GeV Au + Au collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Low-Frequency Correlated Charge-Noise Measurements Across Multiple Energy Transitions in a Tantalum Transmon

Transmon qubits fabricated with tantalum metal have been shown to possess energy relaxation times greater than 400 𝜇⁢s and, as such, present an attractive platform for high precision, correlated noise studies across multiple higher-energy transitions. Tracking the multilevel fluctuating qudit frequencies with a precision enabled by the high coherence of the device allows us to extract the charge offset and quasiparticle dynamics. We observe qualitatively different charge-offset behavior in the tantalum device than those measured in previous low-frequency charge-noise studies. In particular, we find the charge-offset dynamics are dominated by rare, discrete jumps between a finite number of quasistationary charge configurations, a previously unobserved charge-noise process in superconducting qubits.

3-dimensional systems↗
DOE OSTI2023

Unitarity cuts of the worldsheet

We compute the imaginary parts of genus-one string scattering amplitudes. Following Witten’s i\varepsilon i ε prescription for the integration contour on the moduli space of worldsheets, we give a general algorithm for computing unitarity cuts of the annulus, Möbius strip, and torus topologies exactly in \alpha' α ′ . With the help of tropical analysis, we show how the intricate pattern of thresholds (normal and anomalous) opening up arises from the worldsheet computation. The result is a manifestly-convergent representation of the imaginary parts of amplitudes, which has the analytic form expected from Cutkosky rules in field theory, but bypasses the need for performing laborious sums over the intermediate states. We use this representation to study various physical aspects of string amplitudes, including their behavior in the (s,t) ( s , t ) plane, exponential suppression, decay widths of massive strings, total cross section, and low-energy expansions. We find that planar annulus amplitudes exhibit a version of low-spin dominance: at any finite energy, only a finite number of low partial-wave spins give an appreciable contribution to the imaginary part.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The 100 pc White Dwarf Sample in the SDSS Footprint

In this work, we present follow-up spectroscopy of 711 white dwarfs within 100 pc, and we present a detailed model atmosphere analysis of the 100 pc white dwarf sample in the Sloan Digital Sky Survey footprint. Our spectroscopic follow-up is complete for 83% of the white dwarfs hotter than 6000 K, where the atmospheric composition can be constrained reliably. We identify 1508 DA white dwarfs with pure hydrogen atmospheres. The DA mass distribution has an extremely narrow peak at 0.59 M ⊙ and reveals a shoulder from relatively massive white dwarfs with M = 0.7–0.9 M ⊙. Comparing this distribution with binary population synthesis models, we find that the contribution from single stars that form through mergers cannot explain the overabundance of massive white dwarfs. In addition, the mass distribution of cool DAs shows a near absence of M > 1 M ⊙ white dwarfs. The pile-up of 0.7–0.9 M ⊙ and the disappearance of M > 1 M ⊙ white dwarfs is consistent with the effects of core crystallization. Even though the evolutionary models predict the location of the pile-up correctly, the delay from the latent heat of crystallization by itself is insufficient to create a significant pile-up, and additional cooling delays from related effects like phase separation are necessary. We also discuss the population of infrared-faint (ultracool) white dwarfs and demonstrate for the first time the existence of a well-defined sequence in color and magnitude. Curiously, this sequence is connected to a region in the color–magnitude diagrams where the number of white dwarfs with a helium-dominated atmosphere is low. This suggests that the infrared-faint white dwarfs likely have mixed H/He atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Influence of planetary-scale topography on the diurnal thermal tide during the 1971 Martian dust storm

Data obtained with the Mariner 9 infrared spectroscopy experiment during the Martian Dust Storm of 1971 to 1972 are examined for evidence of topographic influence on the atmospheric thermal structure. Temperature perturbations which are well correlated with the planetary scale topography are found superposed on the large amplitude diurnal thermal tide. Applications of tidal theory to the data indicate that the observed perturbations result from the kinematic interaction of the westward traveling diurnal wave with the large scale components of topography. The dominant mode is excited by the wave-number two topography component and is a vertically evanescent eastward traveling wave with an equivalent depth comparable to the atmospheric scale height. The principle dynamic effect of this mode is the enhancement of the amplitude of the near-surface diurnal wind to over 40m/sec in limited areas near 30 deg south latitude. It appears likely that dust was injected into the atmosphere in these regions during the storm.

Conrath, B. J.↗
NASA NTRS1981

Ion composition and energy distribution during 10 magnetic storms

Data from the plasma composition experiment of ISEE 1 are used to investigate the relative quantities and energy characteristics of H(+), He(++), He(+), and O(+) ions in the near-equatorial magnetosphere during magnetic storm conditions. The ions in the study had energies between 0.1 and 17 keV/e and pitch angles between 45 and 135 deg. The data were obtained during 10 storms, for the most part at or immediately following the peak Dst, covering all major local time sectors and geocentric distances between 2 and 15 earth radii. The ion fluxes are averaged over the spacecraft spin angle and over time for periods ranging from about 20 min close to the earth to more than an hour in most distant regions. The inferred 'isotropic' number densities are characterized by a large to dominant fraction of terrestrial ions throughout the energy range covered. The data are found to be consistent with a terrestrial origin for all of the O(+), most of the He(+), and a large but varying fraction of the H(+), whereas the He(++) and part of the H(+) appear to be of solar wind origin.

Lennartsson, W.↗
NASA NTRS1984

Simulation of blunt-fin-induced shock wave and turbulent boundary-layer interaction

A supersonic flow over a blunt fin mounted on a flat plate is numerically simulated. The fin shock causes the boundary layer to separate and results in a complicated, three-dimensional shockwave and boundary-layer interaction. The computed result is in good agreement with the measured pressure on the fin and the flat plate. The main features, such as peak pressure on the fin leading edge and a double peak pressure on the plate, are closely predicted. The role of the horseshoe vortex is discussed. The vortex leads to the development of high-speed flow, and, hence, low-pressure regions on the fin and the plate. Different thicknesses of the incoming boundary layer have been studied. Varying the thicknesses by an order of magnitude shows that the size of the horseshoe vortex and therefore the spatial extent of the interaction are inviscid-dominated, and are weakly dependent on the Reynolds number. Colored graphics are used to show details of the interaction flow field.

Hung, C.-M.↗