Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stokes parameters”

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 451 records · Page 25

Higher Ion Selectivity with Lower Energy Usage Promoted by Electro-osmotic Flow in the Transport through Conical Nanopores

The tradeoff between selectivity and throughput presents fundamental challenges to improve desalination and charge storage, salinity gradient-based energy harvesting, memory device/circuit development, and so forth. The well-known ion-current rectification and the recently resolved time-dependent transport hysteresis in conical nanopores or asymmetric nanointerfaces offer new opportunities for the selective transport of matter. This report shows that electro-osmotic flow (EOF) is an overlooked factor that increases ion selectivity while maintaining enhanced transport throughput in rectified nanoscale electrokinetic transport. The increased selectivity originates primarily from the suppression of anion flux by the fluid flow in the opposite direction under the applied electrical field. By solving the Poisson and Nernst–Planck (PNP) equations without and with coupled Navier–Stokes (PNP–NS) equations, the EOF effects on cation and anion transport are unequivocally revealed in asymmetric nanopipettes. The flux of cations and anions as well as the transference number and flow velocity are elucidated using the models and boundary conditions validated by previous experiments. A dimensionless parameter, radius over the Debye length, reveals optimal ion selectivity and energy cost at intermediate ion concentrations and nanopore sizes, up to hundreds of millimolars and tens of nanometers under time-dependent potential stimulus. Furthermore, the fundamental insights into EOF at nanointerfaces suggest new routes/strategies for better separation, analysis, and energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strategies for estimating the marine geoid from altimeter data

In processing altimeter data from a spacecraft borne altimeter to estimate the fine structure of the marine geoid, a problem is encountered. In order to describe the geoid fine structure, a large number of parameters must be employed and it is not possible to simultaneously estimate all of them. Unless the parameterization exhibits good orthogonality in the data, serious aliasing results. From simulation studies it has been found that amongst several competing parameterizations, the mean free air gravity anomaly model (i.e., Stokes' formula) exhibited promising geoid recovery characteristics. Using covariance analysis techniques, this report provides quantitative measures of the orthogonality properties associated with the above mentioned parameterization. It has been determined that a 5 deg x 5 deg area mean free air gravity anomaly can be estimated with an uncertainty of 1 mgal (40 cm undulation) provided that all free air gravity anomalies within a spherical radius of 10 arc degrees are simultaneously estimated.

Argentiero, P.↗

Numerical simulation of a controlled boundary layer

The problem of interest is the boundary layer over a flat plate. The three standard laminar flow control (LFC) techniques are pressure gradient, suction, and heating. The parameters used to describe the amount of control in the context of the boundary layer equations are introduced. The numerical method required to find the mean flow, the linear eigenvalues of the Orr-Sommerfeld equation, and the full, nonlinear, 3-D solution of the Navier-Stokes equations are outlined. A secondary instability exists for the parallel boundary subject to uniform pressure gradient, suction, or heating. Selective control of the spanwise mode reduces the secondary instability in the parallel boundary layer at low Reynolds number.

Zang, Thomas A.↗

A Newton/upwind method and numerical study of shock wave/boundary layer interactions

The objective of the paper is two-fold. First, an upwind/central differencing method for solving the steady Navier-Stokes equations is described. The symmetric line relation method is used to solve the resulting algebraic system to achieve high computational efficiency. The grid spacings used in the calculations are determined from the triple-deck theory, in terms of Mach and Reynolds numbers and other flow parameters. Thus the accuracy of the numerical solutions is improved by comparing them with experimental, analytical, and other computational results. Secondly, the shock wave/boundary layer interactions are studied numerically, with special attention given to the flow separation. The concept of free interaction is confirmed. Although the separated region varies with Mach and Reynolds numbers, it is found that the transverse velocity component behind the incident shock, which has not been identified heretofore, is also an important parameter. A small change of this quantity is sufficient to eliminate the flow separation entirely.

Liou, Meng-Sing↗

Navier-Stokes computations of a viscous optimized waverider

The performance of a Mach 6 viscous optimized waverider was calculated using the 3-D Navier-Stokes equations. The Mach 6 viscous optimized waverider was generated using MAXWARP, a code developed at the University of Maryland. The computations were performed using CFL3D, an implicit upwind-biased finite-volume algorithm developed at NASA Langley. Results show that good agreement was found between the calculated performance by MAXWARP and results from the Mach 6 Navier-Stokes computation. Furthermore, off-design performance of the Mach 6 optimized waverider was computed at Mach 4 and 8. The performance at these Mach numbers compared well with the performance of the viscous optimized waveriders specifically designed for these Mach numbers. Finally, contours of different flow parameters in the cross-flow plane were examined for the three calculations. The results indicate that the flow gradients are relatively small within the captured flow, and the variation itself is well behaved; thus, making the waverider configuration a promising choice for an engine/airframe design, especially for cruise-type applications.

Takashima, Naruhisa↗

Grid and aerodynamic sensitivity analyses of airplane components

An algorithm is developed to obtain the grid sensitivity with respect to design parameters for aerodynamic optimization. The procedure is advocating a novel (geometrical) parameterization using spline functions such as NURBS (Non-Uniform Rational B-Splines) for defining the wing-section geometry. An interactive algebraic grid generation technique, known as Two-Boundary Grid Generation (TBGG) is employed to generate C-type grids around wing-sections. The grid sensitivity of the domain with respect to geometric design parameters has been obtained by direct differentiation of the grid equations. A hybrid approach is proposed for more geometrically complex configurations such as a wing or fuselage. The aerodynamic sensitivity coefficients are obtained by direct differentiation of the compressible two-dimensional thin-layer Navier-Stokes equations. An optimization package has been introduced into the algorithm in order to optimize the wing-section surface. Results demonstrate a substantially improved design due to maximized lift/drag ratio of the wing-section.

Sadrehaghighi, Ideen↗

Grid Sensitivity and Aerodynamic Optimization of Generic Airfoils

An algorithm is developed to obtain the grid sensitivity with respect to design parameters for aerodynamic optimization. The procedure is advocating a novel (geometrical) parameterization using spline functions such as NURBS (Non-Uniform Rational B- Splines) for defining the airfoil geometry. An interactive algebraic grid generation technique is employed to generate C-type grids around airfoils. The grid sensitivity of the domain with respect to geometric design parameters has been obtained by direct differentiation of the grid equations. A hybrid approach is proposed for more geometrically complex configurations such as a wing or fuselage. The aerodynamic sensitivity coefficients are obtained by direct differentiation of the compressible two-dimensional thin-layer Navier-Stokes equations. An optimization package has been introduced into the algorithm in order to optimize the airfoil surface. Results demonstrate a substantially improved design due to maximized lift/drag ratio of the airfoil.

Sadrehaghighi, Ideen↗

High-Order Hyperbolic Residual-Distribution Schemes on Arbitrary Triangular Grids

In this paper, we construct high-order hyperbolic residual-distribution schemes for general advection-diffusion problems on arbitrary triangular grids. We demonstrate that the second-order accuracy of the hyperbolic schemes can be greatly improved by requiring the scheme to preserve exact quadratic solutions. We also show that the improved second-order scheme can be easily extended to third-order by further requiring the exactness for cubic solutions. We construct these schemes based on the LDA and the SUPG methodology formulated in the framework of the residual-distribution method. For both second- and third-order-schemes, we construct a fully implicit solver by the exact residual Jacobian of the second-order scheme, and demonstrate rapid convergence of 10-15 iterations to reduce the residuals by 10 orders of magnitude. We demonstrate also that these schemes can be constructed based on a separate treatment of the advective and diffusive terms, which paves the way for the construction of hyperbolic residual-distribution schemes for the compressible Navier-Stokes equations. Numerical results show that these schemes produce exceptionally accurate and smooth solution gradients on highly skewed and anisotropic triangular grids, including curved boundary problems, using linear elements. We also present Fourier analysis performed on the constructed linear system and show that an under-relaxation parameter is needed for stabilization of Gauss-Seidel relaxation.

Mazaheri, Alireza↗

COPTRAN - A method of optimum communication systems design

Single set of mathematical expressions describes system cost and probability of error of data transmission in terms of four basic parameters in the link equation. A Lagrange multiplier sets up equations whose solutions yield the optimum values for system design considerations and weight and cost values.

Brinkman, K. L.↗

Viscous starting jets

The transient motion which is produced when a viscous incompressible fluid is forced from an initial state of rest is studied. The equations for unsteady particle paths, written in terms of similarity variables, are analyzed as a quasi-autonomous system with the Reynolds number treated as a parameter. By finding and classifying critical points in the system's phase portrait, the flow structure is examined. It is shown that: (1) bifurcations in the phase portrait occur at specific values of the Reynolds number of the flow in question, and (2) the exact solutions of the Stokes equations for the low-Reynolds-number limit contain two critical Reynolds numbers and three distinct states of motion which culminate in the onset of a vortex roll-up.

Cantwell, Brian J.↗

A solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method.

Rosenfeld, Moshe↗

Supersonic flow computations for an ASTOVL aircraft configuration

A multi-zonal, unified, space/time marching method has been used to solve the Euler and Reynolds-averaged Navier-Stokes equations for supersonic flow past an Advanced Short Take-Off and Vertical Landing (ASTOVL) aircraft configuration. Lift, drag and pitching moment values obtained from the computations compare well with wind-tunnel measurements. The calculations demonstrate the capability of the method used to accurately predict design parameters such as lift, drag and pitching moment for very complex aircraft configurations.

Ramakrishnan, S. V.↗

A novel potential/viscous flow coupling technique for computing helicopter flow fields

Because of the complexity of helicopter flow field, a zonal method of analysis of computational aerodynamics is required. Here, a new procedure for coupling potential and viscous flow is proposed. An overlapping, velocity coupling technique is to be developed with the unique feature that the potential flow surface singularity strengths are obtained directly from the Navier-Stokes at a smoother inner fluid boundary. The closed-loop iteration method proceeds until the velocity field is converged. This coupling should provide the means of more accurate viscous computations of the near-body and rotor flow fields with resultant improved analysis of such important performance parameters as helicopter fuselage drag and rotor airloads.

Summa, J. Michael↗

Assessment of RANS-based Transition Models based on Experimental Data of the Common Research Model with Natural Laminar Flow

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework often rely upon transition correlations that were derived from a limited number of low-speed experiments and these models often fail to account for all of the relevant transition mechanisms and/or the variation in those mechanisms with respect to changes in the significant flow parameters. Available data from a recent experiment on the Common Research Model with Natural Laminar Flow (CRM-NLF) in the National Transonic Facility at the NASA Langley Research Center are used to assess the current transition modeling capability in NASA's OVERFLOW 2.3b solver for a swept wing configuration with nonzero taper and transonic cruise conditions. Specifically, the OVERFLOW solutions are used to evaluate the accuracy and robustness of the transport-equation-based transition models. Results highlight that the Spalart-Allmaras-based amplification factor transport (AFT-2017b) equation model and Menter’s shear-stress transport equation (SST2003)-based Langtry-Menter transition models (either with or without the modeling of crossflow transition) significantly underpredict the reported extent of laminar flow region over the entire span of the wing, irrespective of which instability mechanism(s) is expected to dominate the onset of the transition process. We show that the transition correlations underlying these models fail to account for the stabilizing effect of compressibility on the Tollmien-Schlichting transition, which is likely to be a major contributor to the underprediction of the laminar flow region on the CRM-NLF. The SST-2003-based Langtry-Menter model also appears to inaccurately predict the chordwise pressure variation along the majority of the wing span at all the flow conditions studied herein, due to how the turbulence intensity levels were enforced in the computations and how that was interfering with the functioning of the underlying turbulence model within the boundary layer. The AFT and Langtry-Menter models appear to be sensitive to the level of the freestream turbulence intensity, but the degree of sensitivity varies across the models.

CFD modeling↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

A title-gap flow model for use in aerodynamic loads assessment of space shuttle thermal protection system: Parallel gap faces

The problem of predicting aerodynamic loads on the insulating tiles of the space shuttle thermal protection system (TPS) is discussed and seen to require a method for predicting pressure and mass flux in the gaps between tiles. A mathematical model of the tile-gap flow is developed, based upon a slow viscous (Stokes) flow analysis, and is verified against experimental data. The tile-gap pressure field is derived from a solution of the two-dimensional Laplace equation; the mass-flux vector is then calculated from the pressure gradient. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The flow model shows tile-gap mass flux to be very sensitive to the gap width indicating a need for coupling the TPS flow and tile displacement calculation. Analytical and experimental work to improve TPS flow predictions and a possible shuttle TPS hardware modification are recommended.

Dwoyer, D. L.↗

A survey of simulation and diagnostic techniques for hypersonic nonequilibrium flows

The possible means of simulating nonequilibrium reacting flows in hypersonic environments, and the required diagnostic techniques, are surveyed in two categories: bulk flow behavior and determination of chemical rate parameters. Flow visualization of shock shapes for validation of computational-fluid dynamic calculations is proposed. The facilities and the operating conditions necessary to produce the required nonequilibrium conditions, the suitable optical techniques, and their sensitivity requirements, are surveyed. Shock-tubes, shock-tunnels, and ballistic ranges in a wide range of sizes and strengths are found to be useful for this purpose, but severe sensitivity requirements are indicated for the optical instruments, which can be met only by using highly-collimated laser sources. Likewise, for the determination of chemical parameters, this paper summarizes the quantities that need to be determined, required facilities and their operating conditions, and the suitable diagnostic techniques and their performance requirements. Shock tubes of various strengths are found to be useful for this purpose. Vacuum ultraviolet absorption and fluorescence spectroscopy and coherent anti-Stokes Raman spectroscopy are found to be the techniques best suited for the measurements of the chemical data.

Sharma, Surendra P.↗

The use of a Navier-Stokes code in the wing design process

The feasibility was determined of incorporating the Navier-Stokes computational code, CFL3D, into the supersonic wing design process. The approach taken is of two steps. The first step was to calibrate CFL3D against existing experimental data sets obtained on thin sharp edged delta wings. The experimental data identified six flow types which are dependent on the similarity parameters of Mach number and angle of attack normal to the leading edge. The calibration showed CFL3D capable of simulating these various separated and attached flow conditions. The second step was to use CFL3D to study the initial formation of leading edge separation over delta wings at supersonic speeds. This consisted of examining solutions obtained on a 65 deg delta wing at Mach number of 1.6 with varying cross sectional shapes. Reynolds number was held constant at 1000000 and the Baldwin-Lomax turbulence model was used. The study showed that through the use of leading edge radius and/or camber, the onset of leading edge separation can be delayed to a higher angle of attack than observed on a flat sharp edged wing. Based on the geometries studied, three wind tunnel models are being designed to verify these results.

Mcmillin, S. Naomi↗