Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamics methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 343 records · Page 19

Oxgen-burning hydrodynamics. 1: Steady shell burning

With new hydrodynamic techniques, the relatively fast evolutionary stages of a star prior to core collapse may be explicitly computed in two spatial dimensions, with a treatment of the microphysics (e.g., nuclear reactions, equation of state, neutrino cooling) which is comparable to typical one-dimensional simulations. The nature of shell oxygen burning in a massive star, prior to core collapse, is used as a first example; it is of particular interest because it is (1) the region in which Ni-56 will be produced by the supernova shock, (2) the region of the 'mass cut', which will separate the collapsed core from the ejected mantle, (3) the site of much of the explosive nucleosynthesis, and (4) a suggested source of symmetry breaking to drive mixing instabilities which were observed in SN 1987A. The nature of the shell burning affects the size of the core which will collapse. The method is illustrated on this test case, and the character of the convection is examined.

Arnett, David↗

Hydrodynamic water impact

The hydrodynamic impact of a falling body upon a viscous incompressible fluid was investigated by numerically solving the equations of motion. Initially the mathematical model simulated the axisymmetric impact of a rigid right circular cylinder upon the initially quiescent free surface of a fluid. A compressible air layer exists between the falling cylinder and the liquid free surface. The mathematical model was developed by applying the Navier-Stokes equations to the incompressible air layer and the incompressible fluid. Assuming the flow to be one dimensional within the air layer, the average velocity, pressure and density distributions were calculated. The liquid free surface was allowed to deform as the air pressure acting on it increases. For the liquid the normalized equations were expressed in two-dimensional cylindrical coordinates. The governing equations for the air layer and the liquid were expressed in finite difference form and solved numerically. For the liquid a modified version of the Marker-and-Cell method was used. The mathematical model has been reexamined and a new approach has recently been initiated. Essentially this consists of examining the impact of an inclined plate onto a quiesent water surface with the equations now formulated in cartesian coordinates.

Kettleborough, C. F.↗

Method to Increase Performance of Foil Bearings Through Passive Thermal Management

This invention is a new approach to designing foil bearings to increase their load capacity and improve their reliability through passive thermal management. In the present case, the bearing is designed in such a way as to prevent the carryover of lubricant from the exit of one sector to the inlet of the ensuing sector of the foil bearing. When such passive thermal management techniques are used, bearing load capacity is improved by multiples, and reliability is enhanced when compared to current foil bearings. This concept has recently been tested and validated, and shows that load capacity performance of foil bearings can be improved by a factor of two at relatively low speeds with potentially greater relative improvements at higher speeds. Such improvements in performance with respect to speed are typical of foil bearings. Additionally, operation of these newly conceived bearings shows much more reliability and repeatable performance. This trait can be exploited in machine design to enhance safety, reliability, and overall performance. Finally, lower frictional torque has been demonstrated when operating at lower (non-load capacity) loads, thus providing another improvement above the current state of the art. The objective of the invention is to incorporate features into a foil bearing that both enhance passive thermal management and temperature control, while at the same time improve the hydrodynamic (load capacity) performance of the foil bearing. Foil bearings are unique antifriction devices that can utilize the working fluid of a machine as a lubricant (typically air for turbines and motors, liquids for pumps), and as a coolant to remove excess energy due to frictional heating. The current state of the art of foil bearings utilizes forced cooling of the bearing and shaft, which represents poor efficiency and poor reliability. This invention embodies features that utilize the bearing geometry in such a manner as to both support load and provide an inherent and passive cooling mechanism. This cooling mechanism functions in such a way as to prevent used (higher temperature) lubricant from being carried over from the exit of one sector into the entry of the next sector of the foil bearing. The disclosed innovation is an improved foil bearing design that reduces or eliminates the need for force cooling of the bearing, while at the same time improving the load capacity of the bearing by at least a factor of two. These improvements are due to the elimination of lubricant carryover from the trailing edge of one sector into the leading edge of the next, and the mixing of used lubricant with the surrounding ambient fluid.

Bruckner, Robert↗

High-Resolution Spectroscopy of Mars: Recent Results and Implications for Atmospheric Evolution

It is believed that Earth, Venus, and Mars were formed by the same rocky and icy planetesimals, which resembled meteorites and comets in their composition, respectively. These planets are thus expected to have initially had the same chemical and isotope composition. Scaling the mass of the terrestrial ocean by the planetary mass ratio, the expected initial H2O abundance on Mars is a layer of about 1 km thick. Scaling the abundance of CO2 on Venus, the expected initial CO2 abundance on Mars is 15 bars. Evidently, significant parts of the initial H2O and CO2 abundances have been lost. Intense meteorite impact erosion and hydrodynamic escape of hydrogen (which could drag to escape more heavy species) were dominant loss processes in the first 0.8 Byr. Later, atmospheric sputtering by O+ ions resulted in the dissociation of CO2 and massive losses of O, C, and H. Formation of carbonates also reduced CO2 to its present abundance which currently exists in the atmosphere, on the polar caps, and is absorbed by regolith. Water loss is currently due to thermal escape of H and nonthermal escape of O, both formed by photodissociation of H2O. All loss processes resulted in fractionation of the H, O, and C isotopes. Therefore, the current isotope ratios in H2O and CO2 are clues to the history of volatiles on Mars. There are three tools to study H2O and CO2 isotopes in the martian atmosphere: (i) mass spectrometry from landing probes, (ii) analyses of Mars' gases trapped in the SNC meteorites which were ejected from Mars, and (iii) high-resolution spectroscopy of the H2O andCO2 bands. Method (i) is the best but is the most expensive. Mass spectrometers to be used should be designed for high-precision isotope measurements. Method (ii) makes it possible to reach an uncertainty +/- 0.1%. However, the obtained results are affected by some uncontrolled interactions: isotope fractionations of (1) trapped gases and (2) those released in pyrolysis, (3) contribution of the impactor, isotope exchanges (4) in the terrestrial environment and (5) with the host rock during pyrolysis. Therefore, the spectroscopic data are of great interest, though their formal accuracy is lower. High-resolution spectroscopy is also a tool to study the current atmosphere of Mars by mapping of some photochemically important species and searching for some minor constituents and their variations. Additional information is contained in the original extended abstract.

Krasnopolsky, V. A.↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

A provisional RR Lyrae distance scale

Hydrodynamic pulsation models are matched with observations of globular cluster RRc stars to derive masses, luminosities and temperatures as functions of two observables: pulsation period, P1, and Fourier phase parameter, phi sub 31. We find that mean RRc masses and luminosities increase, and mean temperatures fall, with decreasing cluster metallicity. The Oosterhoff dichotomy is explained as mainly a temperature effect, while hints of a new dichotomy (in RRc mass and luminosity) are found among certain Oo II clusters. A provisional RR Lyrae distance scale emerges for the RRc stars in the form of a relation between luminosity and the two observables, P1 and phi sub 31. Finally, the scheme we propose passes three independent tests: (1) it reproduces the observed hierarchy of relative luminosity among a large sample of RRc stars in Omega Centauri; (2) it yields mean RRc masses which are fully consistent with the RRd masses derived for the clusters M68 and M15; and (3) it gives an LMC distance modulus in agreement with that obtained by other methods.

Simon, Norman R.↗

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

Electric-field-driven Phenomena for Manipulating Particles in Micro-Devices

Compared to other available methods, ac dielectrophoresis is particularly well-suited for the manipulation of minute particles in micro- and nano-fluidics. The essential advantage of this technique is that an ac field at a sufficiently high frequency suppresses unwanted electric effects in a liquid. To date very little has been achieved towards understanding the micro-scale field-and shear driven behavior of a suspension in that, the concepts currently favored for the design and operation of dielectrophoretic micro-devices adopt the approach used for macro-scale electric filters. This strategy considers the trend of the field-induced particle motions by computing the spatial distribution of the field strength over a channel as if it were filled only with a liquid and then evaluating the direction of the dielectrophoretic force, exerted on a single particle placed in the liquid. However, the exposure of suspended particles to a field generates not only the dielectrophoretic force acting on each of these particles, but also the dipolar interactions of the particles due to their polarization. Furthermore, the field-driven motion of the particles is accompanied by their hydrodynamic interactions. We present the results of our experimental and theoretical studies which indicate that, under certain conditions, these long-range electrical and hydrodynamic interparticle interactions drastically affect the suspension behavior in a micro-channel due to its small dimensions.

Khusid, Boris↗

Method for Reducing Pumping Damage to Blood

Methods are provided for minimizing damage to blood in a blood pump wherein the blood pump comprises a plurality of pump components that may affect blood damage such as clearance between pump blades and housing, number of impeller blades, rounded or flat blade edges, variations in entrance angles of blades, impeller length, and the like. The process comprises selecting a plurality of pump components believed to affect blood damage such as those listed herein before. Construction variations for each of the plurality of pump components are then selected. The pump components and variations are preferably listed in a matrix for easy visual comparison of test results. Blood is circulated through a pump configuration to test each variation of each pump component. After each test, total blood damage is determined for the blood pump. Preferably each pump component variation is tested at least three times to provide statistical results and check consistency of results. The least hemolytic variation for each pump component is preferably selected as an optimized component. If no statistical difference as to blood damage is produced for a variation of a pump component, then the variation that provides preferred hydrodynamic performance is selected. To compare the variation of pump components such as impeller and stator blade geometries, the preferred embodiment of the invention uses a stereolithography technique for realizing complex shapes within a short time period.

Bozeman, Richard J., Jr.↗

Neutral–neutral synthesis of organic molecules in cometary comae

Remote and in situ observations of cometary gases have revealed the presence of a wealth of complex organic molecules, including carbon chains, alcohols, imines, and the amino acid glycine. Such chemical complexity in cometary material implies that impacts by comets could have supplied reagents for prebiotic chemistry to young planetary surfaces. However, the assumption that some of the molecules observed in cometary comae at millimetre wavelengths originate from ices stored inside the nucleus has not yet been proven. In fact, the comae of moderately-active comets reach sufficient densities within a few thousand kilometres of the nucleus for an active (solar radiation-driven) photochemistry to ensue. Here, we present results from our latest chemical-hydrodynamic models incorporating an updated reaction network, and show that the commonly-observed HC3N (cyanoacetylene) and NH2CHO (formamide) molecules can be efficiently produced in cometary comae as a result of two-body, neutral–neutral, gas-phase reactions involving well-known coma species. In the presence of a near-nucleus distributed source of CN (similar to that observed by the Rosetta spacecraft at comet 67P), we find that sufficient HC3N and NH2CHO can be synthesized to match the abundances of these molecules observed previously in Oort cloud comets. The precise coma origins of these (and other) complex organic molecules can be verified through radio interferometric mapping observations, for example using the Atacama Large Millimeter/submillimeter Array.

Comets↗

Unshrouded Centrifugal Turbopump Impeller

The ratio of rocket engine thrust to weight is a limiting constraint in placing more payload into orbit at a low cost. A key component of an engine's overall weight is the turbopump weight, Reducing the turbopump weight can result in significant engine weight reduction and hence, increased delivered payload. There are two main types of pumps: centrifugal and axial pumps. These types of pumps can be further sub-divided into those with shrouds and those without shrouds (unshrouded pumps). Centrifugal pumps can achieve the same pump discharge pressure as an axial pump and it requires fewer pump stages and lower pump weight than an axial pump. Also, with unshrouded centrifugal pumps (impeller), the number of stages and weight can be further reduced. However. there are several issues with regard to using an unshrouded impeller: 1) there is a pump performance penalty due to the front open face recirculation flow; 2) there is a potential pump axial thrust problem from the unbalanced front open face and the back shroud face; and, 3) since test data is very linu'ted for this configuration, there is uncertainty in the magnitude and phase of the rotordynamics forces due to the front impeller passage. The purpose of the paper is to discuss the design of an unshrouded impeller and to examine the design's hydrodynamic performance, axial thrust, and rotordynamics performance. The design methodology will also be discussed. This work will help provide some guidelines for unshrouded impeller design. In particular, the paper will discuss the design of three unshrouded impellers - one with 5 full and 5 partial blades (5+5). one with 6+6 blades and one with 8+8 blades. One of these designs will be selected for actual fabrication and flow test. Computational fluid dynamics (CFD) is used to help design and optimize the unshrouded impeller. The relative pump performance penalty is assessed by comparing the CFD results of the unshrouded impeller with the equivalent shrouded impeller for a particular design. Limited unshrouded - versus - shrouded impeller data from the J-2 pump is used to anchor the CFD. Since no detailed impeller blade force data is available, axial thrust and rotordynamic force predictions are based on the CFD model. For the axial thrust, the impeller front flow passage axial force is integrated from the CFD results and compared to the equivalent shrouded impeller axial force. For the rotordynamics forces, the fluid reaction forces are computed from unsteady flow CFD results using a moving boundary method; the rotor- shaft is moved at several whirl-to-speed frequency ratios to extract the rotordynamics coefficients.

Prueger, George↗

Enthalpy-Based Thermal Evolution of Loops: II. Improvements to the Model

This paper further develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) originally proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating. It has typically been applied to impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. It is found that while the evolution of the loop temperature is rather insensitive to the details of the model, accurate tracking of the density requires the inclusion of our new features. In particular, we are able to now obtain highly over-dense loops in the late cooling phase and decreases to the coronal density arising due to stratification. The 0D results are compared to a 1D hydro code (Hydrad). The agreement is acceptable, with the exception of the flare case where some versions of Hydrad can give significantly lower densities. This is attributed to the method used to model the chromosphere in a flare. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function and (b) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Cargill, P. J.↗

Microgravity-Enhanced Stem Cell Selection

Stem cells, both embryonic and adult, promise to revolutionize the practice of medicine in the future. In order to realize this potential, a number of hurdles must be overcome. Most importantly, the signaling mechanisms necessary to control the differentiation of stem cells into tissues of interest remain to be elucidated, and much of the present research on stem cells is focused on this goal. Nevertheless, it will also be essential to achieve large-scale expansion and, in many cases, assemble cells in 3D as transplantable tissues. To this end, microgravity analog bioreactors can play a significant role. Microgravity bioreactors were originally conceived as a tool to study the cellular responses to microgravity. However, the technology can address some of the shortcomings of conventional cell culture systems; namely, the deficiency of mass transport in static culture and high mechanical shear forces in stirred systems. Unexpectedly, the conditions created in the vessel were ideal for 3D cell culture. Recently, investigators have demonstrated the capability of the microgravity bioreactors to expand hematopoietic stem cells compared to static culture, and facilitate the differentiation of umbilical cord stem cells into 3D liver aggregates. Stem cells are capable of differentiating into functional cells. However, there are no reliable methods to induce the stem cells to form specific cells or to gain enough cells for transplantation, which limits their application in clinical therapy. The aim of this study is to select the best experimental setup to reach high proliferation levels by culturing these cells in a microgravity-based bioreactor. In typical cell culture, the cells sediment to the bottom surface of their container and propagate as a one-cell-layer sheet. Prevention of such sedimentation affords the freedom for self-assembly and the propagation of 3D tissue arrays. Suspension of cells is easily achievable using stirred technologies. Unfortunately, in conventional bioreactors, stirring invokes deleterious forces that disrupt cell aggregation and results in cell death. First-generation rotating bioreactors provided rotation on the horizontal axis, which resulted in the suspension of cells without stirring, thus providing a suitable environment to propagate cells without sedimentation to a surface. The rotating wall bioreactors did not provide a way to remove air bubbles that were causing shear and disrupting 3D cultures. Johnson Space Center successfully engineered the hydrofocusing bioreactor (HFB) that resolved the problem of removing the air bubbles from the fluid medium of NASA's rotating-wall space bioreactors. The HFB uses the principle of hydrodynamic focusing that simultaneously produces a low-shear fluid culture environment and a variable hydrofocusing force that can control the movement, location, and removal of suspended cells, tissues, and air bubbles from the bioreactor. The HFB is a rotating, domeshaped cell culture vessel with a centrally located sampling port and an internal viscous spinner. The vessel and spinner can rotate at different speeds either in the same or opposite directions. Rotation of the vessel and viscous interaction at the spinner generate a hydrofocusing force. Adjusting the differential rotation rate between vessel and spinner controls the magnitude of the force.

Claudio, Pier Paolo↗

The Surface Water and Ocean Topography Satellite Mission - An Assessment of Swath Altimetry Measurements of River Hydrodynamics

The Surface Water and Ocean Topography (SWOT) satellite mission, scheduled for launch in 2020 with development commencing in 2015, will provide a step-change improvement in the measurement of terrestrial surface water storage and dynamics. In particular, it will provide the first, routine two-dimensional measurements of water surface elevations, which will allow for the estimation of river and floodplain flows via the water surface slope. In this paper, we characterize the measurements which may be obtained from SWOT and illustrate how they may be used to derive estimates of river discharge. In particular, we show (i) the spatia-temporal sampling scheme of SWOT, (ii) the errors which maybe expected in swath altimetry measurements of the terrestrial surface water, and (iii) the impacts such errors may have on estimates of water surface slope and river discharge, We illustrate this through a "virtual mission" study for a approximately 300 km reach of the central Amazon river, using a hydraulic model to provide water surface elevations according to the SWOT spatia-temporal sampling scheme (orbit with 78 degree inclination, 22 day repeat and 140 km swath width) to which errors were added based on a two-dimension height error spectrum derived from the SWOT design requirements. Water surface elevation measurements for the Amazon mainstem as may be observed by SWOT were thereby obtained. Using these measurements, estimates of river slope and discharge were derived and compared to those which may be obtained without error, and those obtained directly from the hydraulic model. It was found that discharge can be reproduced highly accurately from the water height, without knowledge of the detailed channel bathymetry using a modified Manning's equation, if friction, depth, width and slope are known. Increasing reach length was found to be an effective method to reduce systematic height error in SWOT measurements.

Wilson, Matthew D.↗

VLTI-AMBER Velocity-Resolved Aperture-Synthesis Imaging of Eta Carinae with a Spectral Resolution of 12 000: Studies of the Primary Star Wind and Innermost Wind-Wind Collision Zone

The mass loss from massive stars is not understood well. Eta Carinae is a unique object for studying the massive stellar wind during the luminous blue variable phase. It is also an eccentric binary with a period of 5.54 yr. The nature of both stars is uncertain, although we know from X-ray studies that there is a wind-wind collision whose properties change with orbital phase. Aims. We want to investigate the structure and kinematics of Car's primary star wind and wind-wind collision zone with a high spatial resolution of approx.6 mas (approx.14 au) and high spectral resolution of R = 12 000. Methods. Observations of Car were carried out with the ESO Very Large Telescope Interferometer (VLTI) and the AMBER instrument between approximately five and seven months before the August 2014 periastron passage. Velocity-resolved aperture-synthesis images were reconstructed from the spectrally dispersed interferograms. Interferometric studies can provide information on the binary orbit, the primary wind, and the wind collision. Results. We present velocity-resolved aperture-synthesis images reconstructed in more than 100 di erent spectral channels distributed across the Br(gamma) 2.166 micron emission line. The intensity distribution of the images strongly depends on wavelength. At wavelengths corresponding to radial velocities of approximately -140 to -376 km/s measured relative to line center, the intensity distribution has a fan-shaped structure. At the velocity of -277 km/s, the position angle of the symmetry axis of the fan is 126. The fan-shaped structure extends approximately 8.0 mas (approx.18:8 au) to the southeast and 5.8 mas (approx.13:6 au) to the northwest, measured along the symmetry axis at the 16% intensity contour. The shape of the intensity distributions suggests that the obtained images are the first direct images of the innermost wind-wind collision zone. Therefore, the observations provide velocity-dependent image structures that can be used to test three-dimensional hydrodynamical, radiative transfer models of the massive interacting winds of Eta Car.

Velocity-resolved aperture-synthesis images↗

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗

VLTI-MATISSE Chromatic Aperture-Synthesis Imaging of η Carinae’s Stellar Wind Across the Brα Line: Periastron Passage Observations in February 2020

Context: Eta Carinae is a highly eccentric, massive binary system (semimajor axis∼15.5 au) with powerful stellar winds and a phase-dependent wind-wind collision (WWC) zone. The primary star, ηCar A, is a luminous blue variable (LBV); the secondary, η Car B, is a Wolf-Rayet or O star with a faster but less dense wind. Aperture-synthesis imaging allows us to study the mass loss from the enigmatic LBV η Car. Understanding LBVs is a crucial step toward improving our knowledge about massive stars and their evolution. Aims: Our aim is to study the intensity distribution and kinematics of η Car’s WWC zone. Methods: Using the VLTI-MATISSE mid-infrared interferometry instrument, we perform Brα imaging of η Car’s distorted wind. Results: We present the first VLTI-MATISSE aperture-synthesis images of η Car A’s stellar wind in several spectral channels distributed across the Brα 4.052 μm line (spectral resolving power R∼960). Our observations were performed close to periastron passage in February 2020 (orbital phase∼14.0022). The reconstructed iso-velocity images show the dependence of the primary stellar wind on wavelength or line-of-sight (LOS) velocity with a spatial resolution of 6 mas (∼14 au). The radius of the faintest outer wind regions is∼26 mas (∼60 au). At several negative LOS velocities, the primary stellar wind is less extended to the northwest than in other directions. This asymmetry is most likely caused by the WWC. Therefore, we see both the velocity field of the undisturbed primary wind and the WWC cavity. In continuum spectral channels, the primary star wind is more compact than in line channels. A fit of the observed continuum visibilities with the visibilities of a stellar wind CMFGEN model (CMFGEN is an atmosphere code developed to model the spectra of a variety of objects) provides a full width at half maximum fit diameter of the primary stellar wind of 2.84±0.06 mas (6.54±0.14 au). We compare the derived intensity distributions with the CMFGEN stellar wind model and hydrodynamic WWC models.

G Weigelt↗