Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear potential flow”

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 307 records · Page 17

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗

A generic scenario analysis of end-of-life plastic management: Chemical additives

Plastic growing demand and the increment in global plastics production have raised the number of spent plastics, out of which over 90% are either landfilled or incinerated. Both methods for handling spent plastics are susceptible to releasing toxic substances, damaging air, water, soil, organisms, and public health. Improvements to the existing infrastructure for plastics management are needed to limit chemical additive release and exposure resulting from the end-of-life (EoL) stage. This article analyzes the current plastic waste management infrastructure and identifies chemical additive releases through a material flow analysis. Additionally, we performed a facility-level generic scenario analysis of the current U.S. EoL stage of plastic additives to track and estimate their potential migration, releases, and occupational exposure. Potential scenarios were analyzed through sensitivity analysis to examine the merit of increasing recycling rates, using chemical recycling, and implementing additive extraction post-recycling. Our analyses identified that the current state of plastic EoL management possesses high mass flow intensity toward incineration and landfilling. Although maximizing the plastic recycling rate is a reasonably straightforward goal for enhancing material circularity, the conventional mechanical recycling method requires improvement because major chemical additive release and contamination routes act as obstacles to achieving high-quality plastics for future reuse and should be mitigated through chemical recycling and additive extraction. In conclusion, the potential hazards and risks identified in this research create an opportunity to design a safer closed-loop plastic recycling infrastructure to handle additives strategically and support sustainable materials management efforts to transform the US plastic economy from linear to circular.

54 ENVIRONMENTAL SCIENCES↗

An improved method for the aerodynamic analysis of wing-body-tail configurations in subsonic and supersonic flow. Part 1: Theory and application

A new method has been developed for calculating the pressure distribution and aerodynamic characteristics of wing-body-tail combinations in subsonic and supersonic potential flow. A computer program has been developed to perform the numerical calculations. The configuration surface is subdivided into a large number of panels, each of which contains an aerodynamic singularity distribution. A constant source distribution is used on the body panels, and a vortex distribution having a linear variation in the streamwise direction is used on the wing and tail panels. The normal components of velocity induced at specified control points by each singularity distribution are calculated and make up the coefficients of a system of linear equations relating the strengths of the singularities to the magnitude of the normal velocities. The singularity strengths which satisfy the boundary condition of tangential flow at the control points for a given Mach number and angle of attack are determined by solving this system of equations using an interactive procedure. Once the singularity strengths are known, the pressure coefficients are calculated, and the forces and moments acting on the configuration determined by numerical integration.

Woodward, F. A.↗

The effect of solar array voltage patterns on plasma power losses

The use of high-voltage solar arrays in space is discussed in connection with the draining of array power by currents flowing between exposed surfaces through the surrounding plasma. The possibility of reducing the power loss by arranging solar cell strings in repeated small-area modules to eliminate any large areas at high potentials is investigated. It is found that the difference in power loss between modular and linear patterned high-voltage arrays is fairly small. Although the use of modular patterns can reduce the effective mean potential by about 10%, for the type of configuration being considered there is also a 10% increase in sheath area, leading to only a few percent change in total power loss. It is concluded that plasma power loss should not be a primary consideration in designing the physical arrangement of high-voltage arrays.

Mandell, M. J.↗

High-velocity, multistage, nozzled, ion driven wind generator and method of operation of the same adaptable to mesoscale realization

Gas flows of modest velocities are generated when an organized ion flux in an electric field initiates an ion-driven wind of neutral molecules. When a needle in ambient air is electrically charged to a potential sufficient to produce a corona discharge near its tip, such a gas flow can be utilized downstream of a ring-shaped or other permeable earthed electrode. In view of the potential practical applications of such devices, as they represent blowers with no moving parts, a methodology for increasing their flow velocities includes exploitation of the divergence of electric field lines, avoidance of regions of high curvature on the second electrode, control of atmospheric humidity, and the use of linear arrays of stages, terminating in a converging nozzle. The design becomes particularly advantageous when implemented in mesoscale domains.

Dunn-Rankin, Derek↗

Secondary Instability of Second Modes in Hypersonic Boundary Layers

Second mode disturbances dominate the primary instability stage of transition in a number of hypersonic flow configurations. The highest amplification rates of second mode disturbances are usually associated with 2D (or axisymmetric) perturbations and, therefore, a likely scenario for the onset of the three-dimensionality required for laminar-turbulent transition corresponds to the parametric amplification of 3D secondary instabilities in the presence of 2D, finite amplitude second mode disturbances. The secondary instability of second mode disturbances is studied for selected canonical flow configurations. The basic state for the secondary instability analysis is obtained by tracking the linear and nonlinear evolution of 2D, second mode disturbances using nonlinear parabolized stability equations. Unlike in previous studies, the selection of primary disturbances used for the secondary instability analysis was based on their potential relevance to transition in a low disturbance environment and the effects of nonlinearity on the evolution of primary disturbances was accounted for. Strongly nonlinear effects related to the self-interaction of second mode disturbances lead to an upstream shift in the upper branch neutral location. Secondary instability computations confirm the previously known dominance of subharmonic modes at relatively small primary amplitudes. However, for the Purdue Mach 6 compression cone configuration, it was shown that a strong fundamental secondary instability can exist for a range of initial amplitudes of the most amplified second mode disturbance, indicating that the exclusive focus on subharmonic modes in the previous applications of secondary instability theory to second mode primary instability may not have been fully justified.

Li, Fei↗

Physics-constrained deep learning of nonlinear normal modes of spatiotemporal fluid flow dynamics

In this study, we present a physics-constrained deep learning method to discover and visualize from data the invariant nonlinear normal modes (NNMs) which contain the spatiotemporal dynamics of the fluid flow potentially containing strong nonlinearity. Specifically, we develop a NNM-physics-constrained convolutional autoencoder (NNM-CNN-AE) integrated with a multi-temporal-step dynamics prediction block to learn the nonlinear modal transformation, the NNMs containing the spatiotemporal dynamics of the flow, and reduced-order reconstruction and long-time future-state prediction of the flow fields, simultaneously. In test cases, we apply the developed method to analyze different flow regimes past a cylinder, including laminar flows with low Reynolds number in transient and steady states (RD = 100) and high Reynolds number flow (RD = 1000), respectively. The results indicate that the identified NNMs are able to reveal the nonlinear spatiotemporal dynamics of these flows, and the NNMs-based reduced-order modeling consistently achieves better accuracy with orders of magnitudes smaller errors in construction and prediction of the nonlinear velocity and vorticity fields, compared to the linear proper orthogonal decomposition (POD) method and the Koopman-constrained-CNN-AE using the same number or dimension of modes. We perform an analysis of the modal energy distribution of NNMs and find that compared to POD modes, the few fundamental NNMs capture a very high level of total energy of the flow, which is advantageous for reduced-order modeling and representation of the complex flows. Finally, we discuss the potentials and limitations of the presented method.

Mechanics↗

Aerodynamic Characteristics of a Supersonic Fighter Aircraft Model at Mach 0.40 to 2.47

The aerodynamic characteristics of an advanced twin-engine fighter aircraft designed for supersonic cruise have been studied in the Langley 16-Foot Transonic Tunnel and the Lewis 10- by 10-Foot Supersonic Tunnel. The objective of this investigation was to establish an aerodynamic data base for the configuration with flow-through nacelles and representative inlets. The use of a canard for trim and the effects of fairing over the inlets were assessed. Comparisons between experimental and theoretical results were also made. The theoretical results were determined by using a potential vortex lift code for subsonic speeds and a linear aerodynamic code for supersonic speeds. This investigation was conducted at Mach numbers from 0.40 to 2.47, at angles of attack from 0 deg to about 20 deg, and at inlet capture ratios of about 0.5 to 1.4.

Capone, F. J.↗

The effects of Kelvin-Helmholtz instability on resonance absorption layers in coronal loops

One of the long-standing uncertainties in the wave-resonance theory of coronal heating is the stability of the resonance layer. The wave motions in the resonance layer produce highly localized shear flows which vary sinusoidally in time with the resonance period. This configuration is potentially susceptible to the Kelvin-Helmholtz instability (KHI), which can enhance small-scale structure and turbulent broadening of shear layers on relatively rapid ideal timescales. We have investigated numerically the response of a characteristic velocity profile, derived from resonance absorption models, to finite fluid perturbations comparable to photospheric fluctuations. We find that the KHI primarily should affect long (approximately greater than 6 x 10(exp 4) km) loops where higher velocity flows (M approximately greater than 0.2) exist in resonance layers of order 100 km wide. There, the Kelvin-Helmholtz growth time is comparable to or less than the resonance quarter-period, and the potentially stabilizing magnetic effects are not felt until the instability is well past the linear growth stage. Not only is the resonance layer broadened by the KHI, but also the convective energy transport out of the resonance layer is increased, thus adding to the efficiency of the wave-resonance heating process. In shorter loops, e.g., those in bright points and compact flares, the stabilization due to the magnetic field and the high resonance frequency inhibit the growth of the Kelvin-Helmholtz instability beyond a minimal level.

Karpen, Judith T.↗

Turbulent mixing in supercritical jets: effect of compressibility factor and inflow condition

Fuel injection and turbulent mixing at supercritical pressures determines ignition and combustion in numerous engineering applications. Flow evolution under such conditions is characterized by strong non-linear coupling between dynamics, transport coefficients, and thermodynamics. Experimental studies observe that the jets injected at supercritical pressures exhibit significantly different dynamics from the jets at subcritical conditions, owing to the lack of distinct liquid and gas phases in supercritical state. Thus, the averaged flow quantities such as the potential core length, jet spatial growth rate and velocity decay profiles differ in the two conditions, resulting in different mixed-fluid distributions. In this study, turbulent jet direct numerical simulations (DNS) are performed to examine the variations in flow statistics between injection of Nitrogen (N₂) in Nitrogen (N₂) at both subcritical (perfect-gas) and supercritical conditions. In all cases, isothermal round jets at Reynolds number (Re_{D}), based on jet diameter (D) and jet orifice velocity (U₀), of 5000 are considered. For mixing analyses, a passive scalar transported with the flow is examined.

Sharan, Nek↗

PROGRAM VSAERO: A computer program for calculating the non-linear aerodynamic characteristics of arbitrary configurations: User's manual

VSAERO is a computer program used to predict the nonlinear aerodynamic characteristics of arbitrary three-dimensional configurations in subsonic flow. Nonlinear effects of vortex separation and vortex surface interaction are treated in an iterative wake-shape calculation procedure, while the effects of viscosity are treated in an iterative loop coupling potential-flow and integral boundary-layer calculations. The program employs a surface singularity panel method using quadrilateral panels on which doublet and source singularities are distributed in a piecewise constant form. This user's manual provides a brief overview of the mathematical model, instructions for configuration modeling and a description of the input and output data. A listing of a sample case is included.

Maskew, B.↗

Design and Development of a Particle Flow Control Mechanism for Particle-Based CSP System Applications

The submission is a two-page extended abstract. Within particle-based Concentrated Solar Power Systems (CSP), particles are the heat transfer medium, and therefore, the heat transfer effectiveness of these systems is reliant on their ability to effectively control particle mass flow rate. The vast majority of approaches for particle flow control in high temperature CSP is based on a slide-gate approach where a gate is linearly actuated across a tear-drop shaped opening. The slide-gate mechanism has a sensitivity of 0.354 g/s-mm, which is defined as the mass flow rate over the position of the gate. In the present work, we have designed and tested a new particle flow control mechanism with the potential for faster response times and different sensitivities than a slide gate approach enabling more control of particle flow.

14 SOLAR ENERGY↗

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

Sedimentation and shear-induced dynamics of spheroids in fluids with spatial viscosity variations

A generalized reciprocal theorem is used to relate the force and torque induced on a particle in an inertia-less fluid with small variation in viscosity to integrals involving Stokes flow fields and the spatial dependence of viscosity. These resistivity expressions are analytically evaluated using spheroidal harmonics and then used to obtain the mobility of the spheroid during sedimentation, and in linear flows, of a fluid with linear viscosity stratification. The coupling between the rotational and translational motion induced by stratification rotates the spheroid’s centerline, creating a variety of rotational and translational dynamics dependent upon the particle’s aspect ratio, κ, and the component of the stratification unit vector in the gravity direction, d g . Spheroids with 0.55 ⪅ κ ⪅ 2.0 exhibit the largest variety of settling behaviors. Interestingly, this range covers most microplastics and typical microorganisms. One of the modes include a stable orientation dependent only on κ and d g , but independent of initial orientation, thus allowing for the potential control of settling angles and sedimentation rates. In a simple shear flow, cross-streamline migration occurs due to the stratification-induced force generated on the particle. Similarly, a particle no longer stays at the stagnation point of a uniaxial extensional flow. While fully analytical results are obtained for spheroids, numerical simulations provide a source of validation. These simulations also provide additional insights into the stratification-induced force- and torque-producing mechanisms through the stratification-induced stress, which is not accessed in the reciprocal theorem-based analytical calculations.

Geophysical and Geological Flows: Stratified flows↗

International Space Station Bipropellant Plume Contamination Model Update for Short Thruster Pulse Widths

The International Space Station (ISS) Bipropellant Plume Contamination Model has been a vital tool for characterizing the thruster plume-induced contamination environment at the ISS but was not intended to be used for very short thruster pulse widths. This paper presents an updated model that ensures full start-up and shut-down phases are modeled for each thruster firing, when the majority of liquid phase contaminant mass is released. The updated ISS Bipropellant Plume Contamination Model prevents potential under-prediction of thruster plume-induced contamination due to visiting vehicle proximity operations and provides a way to take advantage of thruster start-up and shut-down data performance data gathered during thruster test programs, if available. The International Space Station (ISS) Bipropellant Plume Contamination Model developed by Boeing Space Environments is a semi-empirical model anchored in flight experiment data and has been a vital tool for characterizing the thruster plume-induced contamination environment at the ISS.[1] The current model utilizes flight experiment data from the Plume Impingement Contamination (PIC) and Shuttle Plume Impingement Flight Experiment (SPIFEX) studies, which include Orbiter 3870 N Primary Reaction Control System (PRCS) and Russian 130 N thrusters operating in pulse mode with 80-100 ms.[2] As the next generation of crew and cargo visiting vehicles are developed and arrive at ISS, minimum pulse widths of vehicle thrusters used for ISS proximity operations have decreased significantly below 80 ms. An update to the ISS Bipropellant Contamination Model is needed to prevent potential under-prediction of thruster plume-induced contamination for these very short pulse widths. Contamination due to thruster plumes occurs in the liquid phase (i.e., unburned or partially burned propellant in the plume). Liquid phase releases primarily occur during thruster start-up and shut-down phases, with the steady state phase contributing a small amount to the total contaminant mass released. The current ISS Bipropellant Plume Contamination Model includes functional dependencies on thruster parameters (thrust, mass flow rate, specific impulse (Isp)), distance to receiver surface, and angle off plume centerline. For modeling purposes, contaminant mass released scales linearly with pulse width. This is a conservative approach for pulse widths greater than 80 ms but could under-predict contamination for very short pulse widths (i.e., with little or no steady state phase). This paper presents a model update to add a functional dependency on commanded pulse width to ensure the portion of the pulse spent in start-up and shut-down phases is appropriately modeled, when the majority of liquid phase contaminant mass is released. The updated model prevents potential under-prediction of thruster plume-induced contamination due to visiting vehicle proximity operations while providing a way to take advantage of thruster start-up and shut-down data gathered during test programs, if available. Developing the plume model for a specific thruster using the ISS Bipropellant Plume Contamination Model must be done in consideration of all available thruster performance data. Example thruster data will be used to illustrate this point and options for implementing the model update with existing ISS Bipropellant Plume Contamination Model code will be discussed.

plume contamination↗

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S=αN/f, whereα,N, andfare the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases withS, asymptotically matching Lentz and Chapman forS< 1. When the solution conserves interior potential vorticity, the momentum transport is proportional toS 2 forS< 1 and is in better agreement with numerical simulations.

Oceanography↗

The effects of canard-wing flow-field interactions on longitudinal stability, effective dihedral and potential deep-stall trim

The literature available on high aspect ratio canard configurations shows them to have some unique stability characteristics. Using a generic canard-wing model, the effects of canard-wing flow-field interactions on stability were investigated in the NASA Langley Twelve-Foot Low-Speed Wind Tunnel. Results for the attached flow regime indicate linear interactions shift the neutral point of a canard configuration forward, but the effect of a canard on a wing can change significantly when the flow over the surface begins to separate, even several degrees below stall. The asymmetry of the canard downwash in a sideslip condition can result in an increment in effective dihedral roughly proportional to canard lift coefficient. At very high angles of attack the presence of a wing can cause an incremental normal force on a canard, contributing to the possibility of a deep-stall trim point. This effect is greater for a high canard and less for a low one.

Muchmore, C. B., Jr.↗

Extending A Correction Method for Unsteady Transonic Aerodynamics to Variable Camber Continuous Trailing Edge Flap

In this study, we have extended our correction method for transonic aerodynamics to baseline and variable camber continuous trailing edge (VCCTEF) geometries of super-critical airfoils at cruise conditions. This correction method modifies Theodorsen theory, which is based on incompressible potential flow over a flat plate. The method is partly based on CFD RANS simulations of oscillating NASA generic transport model (GTM) baseline cruise and VCCTEF airfoil geometries in pitch and shows promise for addressing transonic flutter problems of aeroelasticity. In the earlier study1, the Theodorsen functions were modulated with non-linear functions of Mach number, M, and reduced frequency, k, and thus the correction was demonstrated for the symmetrical NACA0012 airfoil. In the present study, the modulation functions are also allowed to account for the effects of camber and thickness of the airfoil. This is demonstrated by applying the correction method to various symmetric NACA00xx airfoil geometries (NACA001, NACA0002, NACA0012, NACA0015 and NACA0018), the NASA GTM super-critical baseline and VCCTEF airfoil geometries, e.g., the circular airfoil, VCCTEF222. First results indicate that the correction method accounts for a large variation of thickness of airfoils and camber. The present correction method will guide the development of a new state space model for the VCCTEF system and eventually a new transfer function that will be incorporated in a new aeroelastic framework leading to an appropriate transonic flutter model for use in the future aircraft systems in development under the NASA Advanced Air Transportation Technologies (AATT) project.

Kaul, Upender K.↗