Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D velocity models”

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 19 records

Calculation of fully three-dimensional separated flow with an unsteady viscous-inviscid interaction method

Three previous papers have shown that the viscous inviscid numerical methods were capable of calculating separated turbulent flows. The 'Semi-inverse method' and the models of the first author allow the computation of massive separations, stalled flows, and shock wave-boundary layer interactions, in two-dimensional or quasi-three-dimensional conditions, from low speeds to supersonic speeds. The 'Semi-implicit' method for time-consistent coupling allows for the computation of time-accurate transonic separated flow and buffer computations in two-dimensions. The present paper shows that the viscous-inviscid interaction approach is also able to compute the fully three-dimensional flow separation phenomena. The method is based on a thin-layer approximation of the theory of 'Defect-Formulation' that provides the viscous-inviscid splitting of the Navier-Stokes equations. A parametric analytical modelling of the 3D-turbulent velocity profiles is involved. Numerically, the 3D-velocity profiles are discretized in the normal z-direction and driven by parametric integral 3D-equations in direct or inverse modes in the x-direction. The viscous-inviscid coupling is fully 3D and solved the time-consistency problem with an extension of the 'Semi-implicit' method previously suggested in two-dimensions. A 3D inviscid subroutine with TSP approximation is used. Results are obtained for transonic steady flows over wings with shock-induced transonic separation. The method provides results for highly three-dimensional flow separations, such as induced by a three-dimensional through at the wall of a flat plate. The 3D viscous-inviscid coupling and the 3D model of the velocity field provide three-dimensional instantaneous skin-friction lines whose patterns exhibit the same complex topology as Navier-Stokes solvers, with foci, nodes, and saddle-points.

Leballeur, J. C.↗

Retrieval of wind temperature and pressure from single Doppler radar and a numerical model

A 4D data assimilation algorithm to obtain 3D wind and thermodynamic fields from radial velocities and a numerical model is derived. It is a blend of the continuous updating technique of Charney et al. (1969) and a diagnostic pressure and temperature retrieval technique suggested by Gal-Chen (1978). This model is tested only against model-generated data which are viewed as 'data observed from a real atmosphere'. A thermal bubble developing in a dry, neutral environment is simulated. The thermal recovery technique is found to be a necessary procedure to obtain successful data assimilation results. The present assimilation method is found to be capable of reducing the 'observational error' and making the assimilation run converge toward the control run. Without data insertion, errors in the nonobserved wind component are large.

Liou, Yu-Chieng↗

Probabilistic Path Planning of Montgolfier Balloons in Strong, Uncertain Wind Fields

Lighter-than-air vehicles such as hot-air balloons have been proposed for exploring Saturn s moon Titan, as well as other bodies with significant atmospheres. For these vehicles to navigate effectively, it is critical to incorporate the effects of surrounding wind fields, especially as these winds will likely be strong relative to the control authority of the vehicle. Predictive models of these wind fields are available, and previous research has considered problems of planning paths subject to these predicted forces. However, such previous work has considered the wind fields as known a priori, whereas in practical applications, the actual wind vector field is not known exactly and may deviate significantly from the wind velocities estimated by the model. A probabilistic 3D path-planning algorithm was developed for balloons to use uncertain wind models to generate time-efficient paths. The nominal goal of the algorithm is to determine what altitude and what horizontal actuation, if any is available on the vehicle, to use to reach a particular goal location in the least expected time, utilizing advantageous winds. The solution also enables one to quickly evaluate the expected time-to-goal from any other location and to avoid regions of large uncertainty. This method is designed for balloons in wind fields but may be generalized for any buoyant vehicle operating in a vector field. To prepare the planning problem, the uncertainty in the wind field is modeled. Then, the problem of reaching a particular goal location is formulated as a Markov decision process (MDP) using a discretized space approach. Solving the MDP provides a policy of what actuation option (how much buoyancy change and, if applicable, horizontal actuation) should be selected at any given location to minimize the expected time-to-goal. The results provide expected time-to-goal values from any given location on the globe in addition to the action policy. This stochastic approach can also provide insights not accessible by deterministic methods; for example, one can evaluate variability and risk associated with different scenarios, rather than only viewing the expected outcome.

Wolf, Michael↗

Multiple-Relaxation-Time Lattice Boltzmann Models in 3D

This article provides a concise exposition of the multiple-relaxation-time lattice Boltzmann equation, with examples of fifteen-velocity and nineteen-velocity models in three dimensions. Simulation of a diagonally lid-driven cavity flow in three dimensions at Re=500 and 2000 is performed. The results clearly demonstrate the superior numerical stability of the multiple-relaxation-time lattice Boltzmann equation over the popular lattice Bhatnagar-Gross-Krook equation.

dHumieres, Dominique↗

3D Modeling of Solar-Type Stars to Characterize Stellar Jitter

Detection of Earth-mass planets requires measurements of radial velocity with extreme precision. To capture the tiny disturbances caused by a planet's motion, it is necessary to understand and characterize the host star’s turbulent dynamics in order to apply proper filtering to the observational data. We take advantage of current computational and technological capabilities to develop 3D realistic models of the stellar subsurface convection and atmospheres and thereby estimate the photospheric jitter. We have identified an initial set of target stars, obtained initial conditions using the MESA code, and obtained initial 3D radiative models of the stellar surfaces and atmospheres with a spatial resolution of 50km. We present initial 3D radiative hydrodynamic model results of the planet-hosting star HD209458.

SMD↗

Background and Pickup Ion Velocity Distribution Dynamics in Titan's Plasma Environment: 3D Hybrid Simulation and Comparison with CAPS T9 Observations

In this report we discuss the ion velocity distribution dynamics from the 3D hybrid simulation. In our model the background, pickup, and ionospheric ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. The current simulation shows that mass loading by pickup ions H(+); H2(+), CH4(+) and N2(+) is stronger than in the previous simulations when O+ ions are introduced into the background plasma. In our hybrid simulations we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M = 28 amu ions that were generated inside the ionosphere. The moon is considered as a weakly conducting body. Special attention will be paid to comparing the simulated pickup ion velocity distribution with CAPS T9 observations. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of the Alfve n wing-like structures. The simulation also shows that the ring-like velocity distribution for pickup ions relaxes to a Maxwellian core and a shell-like halo.

Lipatov, A. S.↗

Turbo-Design: Open-Source Radial Equilibrium Turbomachinery Solver: Part I - Turbines

Advances in 3D Geometrical Designs and Cooling have played a significant role in improving the efficiency of turbomachinery. However, these advancements must be effectively translated back to the modeler. Machine learning can facilitate this transition. Specifically, machine learning–based loss models can bridge the gap between 3D and 1D designs, enabling modelers not only to predict velocity triangles but also to extract additional geometric features. Currently, the design tools used at NASA have not been updated to support such integration—until now. TurboDesign is an open-source, Python-based framework that replaces TD2 (LEW-11029-1) and AXOD2 (LEW-16323-1), both of which are radial equilibrium solvers for axial turbines. The goal of this update is to enable the integration of machine learning loss models into radial equilibrium equations. Additionally, TurboDesign is designed to support radial machines. This paper presents the governing equations, the assumptions underlying the code, the integration of legacy loss models, an example of machine learning model integration, and a validation comparison with CFD. All code, tutorials, and documentation are available at: https://www.github.com/nasa/turbo-design

Radial Equilibrium↗

LDV flowfield measurements on a straight and swept wing with a simulated ice accretion

Simulated glaze ice accretion effects on a flowfield are presently studied for the case of a 3D semispan wing with NACA 0012 airfoil section on a rectangular untwisted planform. The model was tested at both zero-sweep and 30-deg sweep setting with four-beam/two-color LDV flowfield mapping. The comparison of LDV-measure velocity profiles with 3D Navier-Stokes predictions revealed correct trends, but with several differences that are attributable to the turbulence models and grid resolution used in the computations.

Bragg, M. B.↗

Rotor-Fuselage Interaction: Analysis and Validation with Experiment

The problem of rotor-fuselage aerodynamic interaction has to be considered in industry applications from various aspects. First, in order to increase helicopter speed and reduce operational costs, rotorcraft tend to be more and more compact, with a main rotor closer to the fuselage surface. This creates significant perturbations both on the main rotor and on the fuselage, including steady and unsteady effects due to blade and wake passage and perturbed inflow at the rotor disk. Furthermore,the main rotor wake affects the tail boom, empennage and anti-torque system. This has important consequences for helicopter control and vibrations at low speeds and also on tail rotor acoustics (main rotor wake-tail rotor interactions). This report describes the US Army-France MOD cooperative work on this problem from both the theoretical and experimental aspects. Using experimental 3D velocity field and fuselage surface pressure measurements, three codes that model the interactions of a helicopter rotor with a fuselage are compared. These comparisons demonstrate some of the strengths and weaknesses of current models for the combined rotor-fuselage analysis.

Berry, John D.↗

3D MHD Models of Active Region Loops

Present imaging and spectroscopic observations of active region loops allow to determine many physical parameters of the coronal loops, such as the density, temperature, velocity of flows in loops, and the magnetic field. However, due to projection effects many of these parameters remain ambiguous. Three dimensional imaging in EUV by the STEREO spacecraft will help to resolve the projection ambiguities, and the observations could be used to setup 3D MHD models of active region loops to study the dynamics and stability of active regions. Here the results of 3D MHD models of active region loops are presented, and the progress towards more realistic 3D MHD models of active regions. In particular the effects of impulsive events on the excitation of active region loop oscillations, and the generation, propagations and reflection of EIT waves are shown. It is shown how 3D MHD models together with 3D EUV observations can be used as a diagnostic tool for active region loop physical parameters, and to advance the science of the sources of solar coronal activity.

Ofman, Leon↗

Vortex/surface interaction

This paper considers the interaction of a vortex generated upstream in a flow field with a downstream aerodynamic surface that possesses a large chord. The flow is assumed to be steady, incompressible, inviscid and irrotational, and the surface to be semiinfinite. The vortex is considered to be a straight vortex filament. To lowest order the problem is modeled using potential theory, where the 3D Laplace's equation for the velocity potential on the surface is solved exactly. The closed-form equation for pressure distribution obtained from this theory is found to have a square root singularity at the leading-edge. It also converges, as x goes to infinity, to the solution of the 2D point-vortex/infinite plane problem. The pressure coefficient presents an anti-symmetric behavior, near the leading-edge and a symmetric behavior as x goes to infinity.

Bodstein, G. C. R.↗

Development of an algorithm to model an aircraft equipped with a generic CDTI display

A model of human pilot performance of a tracking task using a generic Cockpit Display of Traffic Information (CDTI) display is developed from experimental data. The tracking task is to use CDTI in tracking a leading aircraft at a nominal separation of three nautical miles over a prescribed trajectory in space. The analysis of the data resulting from a factorial design of experiments reveals that the tracking task performance depends on the pilot and his experience at performing the task. Performance was not strongly affected by the type of control system used (velocity vector control wheel steering versus 3D automatic flight path guidance and control). The model that is developed and verified results in state trajectories whose difference from the experimental state trajectories is small compared to the variation due to the pilot and experience factors.

Driscoll, W. C.↗

Review of planar laser-induced fluorescence measurements in the NASA Langley 31-Inch Mach 10 Air Tunnel

In 2003 a mobile planar laser-induced fluorescence (PLIF) system was constructed at NASA Langley Research Center to visualize and quantify high speed flow fields containing fluorescent species such as nitric oxide (NO). Thereafter, a series of test entries utilizing NO PLIF were conducted in the NASA Langley Research Center’s 31 Inch Mach 10 air wind tunnel, which is a hypersonic blowdown wind tunnel located in Hampton, Virginia. This paper briefly reviews laser-based measurements previously performed in hypersonic blowdown wind tunnel facilities and then details NO PLIF measurement campaigns performed at this facility in the 20 years to date. Measurements were performed to develop 2D and 3D visualization capabilities and quantitatively measure velocity and temperature. Applications described herein include historical models (e.g., Apollo, X-33 and X-43), flight vehicles under development at the time (e.g., Mars Science Lab (MSL), Orion, and LOFTID), as well as fundamental research (e.g., flat plates with and without discrete roughness).

PLIF↗

Explicit 3D continuum fracture modeling with smooth particle hydrodynamics

Impact phenomena shaped our solar system. As usual for most solar system processes, the scales are far different than we can address directly in the laboratory. Impact velocities are often much higher than we can achieve, sizes are often vastly larger, and most impacts take place in an environment where the only gravitational force is the mutual pull of the impactors. The Smooth Particle Hydrodynamics (SPH) technique has been applied in the past to the simulations of giant impacts. In these simulations, the colliding objects were so massive (at least a sizeable fraction of the Earth's mass) that material strength was negligible compared to gravity. This assumption can no longer be made when the bodies are much smaller. To this end, we have developed a 3D SPH code that includes a strength model to which we have added a von Mises yielding relation for stresses beyond the Hugoniot Elastic Limit. At the lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution. Following Grady and Kipp and Melosh et al., we introduce a state variable D ('damage'), 0 less than D less than 1, which expresses the local reduction in strength due to crack growth under tensile loading. Unfortunately for the hydrodynamics, Grady and Kipp's model predicts which fragments are the most probable ones and not the ones that are really formed. This means, for example, that if a given laboratory experiment is modeled, the fragment distribution obtained from the Grady-Kipp theory would be equivalent to a ensemble average over many realizations of the experiment. On the other hand, the hydrodynamics itself is explicit and evolves not an ensemble average but very specific fragments. Hence, there is a clear incompatibility with the deterministic nature of the hydrodynamics equations and the statistical approach of the Grady-Kipp dynamical fracture model. We remedy these shortcomings by making the incipient flaw distribution explicit, i.e., particles carry activation strains which are distributed at random with a probability of occurrence given by the Weibull distribution. If the local principal axis strain exceeds this limit, damage starts to grow. By growing explicit cracks together with statistical cracks (damage) at the sub-particle scale, we ensure that material strength and fragmentation is independent of model resolution. We tested our scheme by simulating laboratory impact experiments on basalt spheres.

Benz, W.↗

Resolution of Lava Tubes with Ground Penetrating Radar: The TubeX project

Remote sensing surveys of the Moon and Mars show evidence of lava tubes, which are potential safe havens for human crews and their equipment. Ground penetrating radar (GPR) can be used to map tubes because the void/rock interface at tube ceilings and floors strongly reflects radar pulses. We have tested the capacity of GPR to sense lava tube geometry at Lava Beds National Monument (LBNM) in California, USA. GPR and detailed LiDAR data are presented for two tubes: Skull Cave, with a few meters of overburden, diameter ~10-20 m, and a rubbly floor; and Valentine Cave, with similarly thin overburden, diameter ~1-3 meters, and a flatter smoother floor. On both caves GPR clearly resolves the ceiling and permits good estimates of the cave width as validated with LiDAR data. Where GPR fails, the primary cause is inferred to be strong out-of-plane effects due to complex 3D geometries. Recovery of the floor position requires migrating the GPR data with a 2D velocity model, as signal velocity is faster in void space. We find that floor position is recoverable in caves whose voids are taller than the radar wavelength (~3 m in this study). Forward modeling assuming planetary parameters suggests GPR should be similarly successful on the Moon or Mars.

S. Esmaeili↗

Time Dependence of Collision Probabilities During Satellite Conjunctions

The NASA Conjunction Assessment Risk Analysis (CARA) team has recently implemented updated software to calculate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm can employ complex dynamical models for orbital motion, and account for the effects of non-linear trajectories as well as both position and velocity uncertainties. This “3D P (sub c)” method entails computing a 3-dimensional numerical integral for each estimated probability. Our analysis indicates that the 3D method provides several new insights over the traditional “2D P (sub c)” method, even when approximating the orbital motion using the relatively simple Keplerian two-body dynamical model. First, the formulation provides the means to estimate variations in the time derivative of the collision probability, or the probability rate, R (sub c). For close-proximity satellites, such as those orbiting in formations or clusters, R (sub c) variations can show multiple peaks that repeat or blend with one another, providing insight into the ongoing temporal distribution of risk. For single, isolated conjunctions, R (sub c) analysis provides the means to identify and bound the times of peak collision risk. Additionally, analysis of multiple actual archived conjunctions demonstrates that the commonly used “2D P (sub c)” approximation can occasionally provide inaccurate estimates. These include cases in which the 2D method yields negligibly small probabilities (e.g., P (sub c)) is greater than 10 (sup -10)), but the 3D estimates are sufficiently large to prompt increased monitoring or collision mitigation (e.g., P (sub c) is greater than or equal to 10 (sup -5)). Finally, the archive analysis indicates that a relatively efficient calculation can be used to identify which conjunctions will have negligibly small probabilities. This small-P (sub c) screening test can significantly speed the overall risk analysis computation for large numbers of conjunctions.

Hall, Doyle T.↗

The Evolution of Oblique Impact Flow Fields Using Maxwell's Z Model

Oblique impacts are the norm rather than the exception for impact craters on planetary surfaces. This work focuses on the excavation of experimental oblique impact craters using the NASA Ames Vertical Gun Range (AVGR). Three-dimensional particle image velocimetry (3D PIV) is used to obtain quantitative data on ejection positions, three dimensional velocities and angles. These data are then used to constrain Maxwell's Z Model and follow the subsurface evolution of the excavation-stage flow-field center during oblique impacts.

Anderson, J. L. B.↗

ROSAT observations of the Orion-Eridanus superbubble

The Orion-Eridanus bubble is a superbubble created by the Orion OB1 stellar association. The X-ray emission from the interior of this bubble forms a prominent feature of the soft X-ray diffuse background extending across approximately 20 x 35 deg. The Rosat all sky survey maps of this object show various soft X-ray shadows cast by neutral foreground clouds against the bright emission of the superbubble interior. The velocities of some of these absorbing clouds were identified and distance limits were established, providing a 3D model of the bubble. The data are consistent with a model for the bubble in which the structure is formed by stellar winds and a recent supernova. Absorption features associated with the shell of the local bubble were identified which enable constraints to be placed on the distance and direction and to obtain measurement of the thermal pressure of the hot gas in this direction.

Burrows, David N.↗