Prediction of drag reduction with a viscoelastic model.
Drag reduction prediction in turbulent flow with viscoelastic model, discussing solutions exhibiting elastic properties in laminar flow
SEARCH · Search NASA
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.
Drag reduction prediction in turbulent flow with viscoelastic model, discussing solutions exhibiting elastic properties in laminar flow
There is yet to be a well-established drag model for non-spherical particles required in a particle-laden flow that could cover a wide range of sphericities. This talk will explore the development of a general drag model for non-spherical particles by applying deep learning using available experimental data available in the literature. The integration of several raw experimental measurements from different sources and research directions allows the training of robust Artificial Intelligence and Machine Learning (ML) models. Neural networks are an ML approach inspired by the inner biological workings of the brain. This work aims to develop a Deep Neural Network (DNN) that predicts drag coefficient values with the ability to adapt appropriately to unseen data. Given the limited number of data points available and the variance found within the data collected from various sources, challenges may arise when looking to train the model. Our study tests and implements various model regularization techniques and assesses different loss and activation functions for the proposed DNN. The proposed model considers a broader range of features other than sphericity and Reynold number. These features include density ratio, solid volume fraction, lengthwise and crosswise sphericity, and more. Furthermore, we present the features that play a significant role in predicting different drag coefficients through feature importance. Within the investigated parameter ranges in this study, the following conclusions can be achieved and summarized below: • An improved drag coefficient model can be developed by considering more features such as, aspect ratio, lengthwise sphericity, crosswise sphericity, and density ratio. • DNN model can predict better results compared to traditional methods using MAE metric. • The proposed model addresses data challenges such as limited data and extreme data points through expanded feature-set and regularization. • Three major features that mostly affect the drag coefficient were identified from a feature importance analysis.
Rough surfaces in contact with a flow of fluid exhibit alternating no-slip and free shear boundary conditions at the solid–liquid and air–liquid interfaces, respectively, thereby potentially offering drag reduction benefits. The balance between the dynamic pressure in the flow and the restoring capillary pressure in the interasperity spaces determines the stability of the Cassie state of wettability and is a function of the relative extent of no-slip and free shear regions per unit surface area. In the present study, using a fractal representation of rough surface topography, an analytical model is developed to quantify the stability of the Cassie state of wettability as well as drag reduction and the friction factor for laminar flow in a rectangular channel between nonwetting multiscale rough surfaces. A systematic study is conducted to quantify the effects of fractal parameters of the surfaces and the flow Reynolds number on drag reduction and the friction factor. The studies are used to develop friction factor curves extending the classical Moody diagram to hydrophobic and superhydrophobic surfaces. On the basis of the studies, regime maps are derived for estimating the extent of drag reduction offered by hydrophobic and superhydrophobic surfaces, revealing that superhydrophobic surfaces do not always offer the best drag reduction performance. The application of the fractal model to practical topographies of nonwetting surfaces of copper, aluminum, and zinc oxide fabricated via electrodeposition and etching is also discussed.
The two body problem was analyzed with a specific drag model. The model treats drag as a force proportional to the vector velocity and inversely proportional to the distance to the center of attraction. The solution is expressed in terms of known functions and is of a simple and compact form. The time of flight is expressed as a quadrature in the true anomaly. The results are: (1) development of a vector differential equation which allows analysis of an infinite number of gravitational and drag models; and (2) obtaining the solution of a linear differential equation using the inverse method of laplace transforms.
The present tests illustrate how the aerodynamic drag of a flying boat hull may be reduced by following closely the form of a low drag aerodynamic body and the manner in which the extent of the aerodynamic refinement is limited by poorer hydrodynamic performance. This limit is not sharply defined but is first evidenced by an abnormal flow of water over certain parts of the form accompanied by a sharp increase in resistance. In the case of models 74-A and 75, the resistance (sticking) occurs only at certain combinations of speed, load, and trim and can be avoided by proper control of the trim at high water speeds. Model 75 has higher water resistance at very high speeds than does model 74-A. With constant speed propellers and high takeoff speeds, it appears that the form of model 75 would give slightly better takeoff performance. Model 74-A, however, has lower aerodynamic drag than does model 75 for the same volume of hull.
We examine the accuracy of density prediction by the first principles model Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM) developed by the National Center for Atmospheric Research and compare it to the accuracy of three empirical models: Jacchia 71, the Naval Research Laboratory Mass Spectrometer Incoherent Scatter Extended 2000 (NRLMSIS), Jacchia 1971, and Jacchia-Bowman 2008. Comparisons are made for three large storms: the October 2003 storm, the March 2013 storm, and the March 2015 storm. To evaluate the accuracy of these models we use tracking data for nine space objects in low Earth orbit. Additionally, we evaluate the accuracy of the TIEGCM and NRLMSIS with data from high precision accelerometers on the Challenging Minisatellite Payload (CHAMP) and Gravity field and Circulation Explorer (GOCE) satellites. The goal is to assess the use of a first principles model as a potential tool for forecasting satellite drag during large magnetic storms. For the storms considered, we found the TIEGCM, JB2008, and NRLMSIS models to be substantially more accurate than the Jacchia 71 model. The accuracies of the TIEGCM and JB2008 models were similar, but overall, the TIEGCM was more accurate. We found smaller differences for TIEGCM versus CHAMP than for NRLMIS for the Halloween Storm, and smaller differences than results published for JB2008 and the assimilative model HASDM. The empirical models are at present more practical for operational purposes, but the TIEGCM, developed as a research model, with a greater focus on operational use offers the potential for improved utility during stressing conditions.
Procedures and data for estimating drag were used to calculate the drag characteristics of the Model 25 airplane. Based on cruise flight test data obtained on the Model 25, these methods generally predicted the total drag characteristics within current acceptable and reasonable engineering accuracy.
Drag and longitudinal trim at low lift of the North American YF-100A airplane at Mach numbers from 0.76 to 1.77 as determined from the flight test of a 0.11-scale rocket model are presented herein. Also included are some longitudinal stability and some qualitative pitch-damping data. The subsonic external-drag-coefficient level was about 0.012, and the supersonic level was about 0.043. The drag rise occurred at a Mach number of 0.95. The longitudinal trim change at low lift consisted basically of a mild nose-up tendency at a Mach number of 0.90. An indication of wing flutter was present at Mach numbers from 0.95 to 1.11. However, the full-scale airplane wing has approximately twice the scaled first-bending frequency as the model tested and, hence, will probably be free of this type of flutter. The aerodynamic-center location was 71 percent behind the leading edge of the mean aerodynamic chord at a Mach number of 1.03 and 62 percent at a Mach number of 1.74. Qualitative measurement of damping in pitch indicates that at low lift coefficients damping will be low at a Mach number of 1.03.
Due to the challenging design requirements, elbows are often unavoidable in duct configuration, and these 90-degree bends introduce swirl, velocity variations, recirculation, and secondary flow. These disturbances make it difficult for nuclear facilities to meet particle sampling standards. A series of numerical analyses are conducted to track aerosols in a nuclear facility duct having a 90-degree elbow with the assistance of computational fluid dynamics (CFD). A turbulence model, a continuity, and a momentum, a discrete phase model, and species transport equations are solved simultaneously to track aerosols in the duct. The effect of turbulence models, turbulent dispersion models, droplet drag model, aerosol amount, aerosol spray configuration, guide vanes, and mixers are investigated. Simulation results are analyzed per relevant testing codes such as DOE-HDBK-1169, ASME AG1, ISO 14644-3, ACGIH, and ANSI/HPS N13.1.
An analytic solution to the two-body problem with a specific drag model is obtained. The model treats drag as a force proportional to the vector velocity and inversely proportional to the square of the distance to the center of attraction. The solution is expressed in terms of known functions and is of a simple and compact form. The time-of-flight is expressed as a quadrature in the 'true anomaly'.
Results of flight tests of an NACA RM-2 drag research model equipped with a pilot's canopy having a vee windshield are presented for a Mach number range from 0.75 to 1.43. Comparison is made with test results of a similar canopy having a flat windshield. The vee-windshield canopy produced lower drag-coefficient values than the flat-windshield canopy for Mach numbers from 0.85 to about 1.2. From M - 1.2 to 1.4 both canopies produced the same drag coefficient.
An experimental study of low Reynolds number drag on laboratory models of dendrite fragments has been conducted. The terminal velocities of the dendrites undergoing free fall along their axis of symmetry were measured in a large Stokes flow facility. Corrections for wall interference give nearly linear drag vs Reynolds number curves. Corrections for both wall interference and inertia effects show that the dendrite Stokes settling velocities are always less than that of a sphere of equal mass and volume. In the Stokes limit, the settling speed ratio is found to correlate well with primary dendrite arm aspect ratio and a second dimensionless shape paremeter which serves as a measure of the fractal-like nature of the dendrite models. These results can be used to estimate equiaxed grain velocities and distance of travel in metal castings. The drag measurements may be used in numerical codes to calculate the movement of grains in a convecting melt in an effort to determine macrosegregation patterns caused by the sink/float mechanism.
An experimental study of low Reynolds number drag on laboratory models of dendrite fragments has been conducted. The terminal velocities of the dendrites undergoing free fall along their axis of symmetry were measured in a large Stokes flow facility. Corrections for wall interference give nearly linear drag vs Reynolds number curves. Corrections for both wall interference and inertia effects show that the dendrite Stokes settling velocities are always less than that of a sphere of equal mass and volume. In the Stokes limit, the settling speed ratio is found to correlate well with the primary dendrite arm aspect ratio and a second dimensionless shape parameter which serves as a measure of the fractal-like nature of the dendrite models. These results can be used to estimate equiaxed grain velocities and distance of travel in metal castings. The drag measurements may be used in numerical codes to calculate the movement of grains in a convecting melt in an effort to determine macrosegregation patterns caused by the sink/float mechanism.
Here, this study applies machine learning-based approaches to develop a drag force model for irregular-shaped particles in incompressible flows. The particle-laden flows are studied through an in-house particle-resolved direct numerical simulation (PR-DNS) at low-intermediate Reynolds numbers (Re). We utilize the PR-DNS to obtain drag force coefficients and flow fields of single particles. A variational auto-encoder model is used to obtain latent vectors to represent the geometrical features of the particles, and artificial neural networks (ANN) are developed to predict drag force coefficients and flow fields of the single particles. This study applies a pairwise interaction extended point-particle (PIEP) model to obtain the coefficients assuming the flow fields of neighboring particles can be linearly superposed. The PIEP method shows significant improvement on prediction for a few neighbored particles. In addition, the results reveal R2 scores of 0.56-0.62 and errors of 9.1-10.0 % for the dilute, polydispersed systems with a volume fraction of 0.5 %.
A free-flight investigation over a Mach number range from 0.6 to 2.0 has been conducted to determine the longitudinal aerodynamic characteristics and effect of rocket jet on zero-lift drag of 1/5-scale models of two ballistic-type missiles, the Hermes A-3A and A-3B. Models of both types of missiles exhibited very nearly linear normal forces and pitching moments over the angle-of-attack range of 8 deg to -4 deg and Mach number range tested. The centers of pressure for both missiles were not appreciably affected by Mach number over the subsonic range; however, between a Mach number of 1.02 and 1.50 the center of pressure for the A-3A model moved forward 0.34 caliber with increasing Mach number. At a trim angle-of-attack of approximately 30 deg, the A-3A model indicated a total drag coefficient 30% higher than the power-off zero-lift drag over the subsonic Mach number range and 10% higher over the supersonic range. Under the conditions of the present test, and excluding the effect of the jet on base drag, there was no indicated effect of the propulsive jet on the total drag of the A-3A model. The propulsive jet operating at a jet pressure ratio p(sub j)/p(sub o) of 0.8 caused approximately 100% increase in base drag over the Mach number range M = 0.6 to 1.0. This increase in base drag amounts to 15% of the total drag. An underexpanded jet operating at jet pressure ratios corresponding approximately to those of the full-scale missile caused a 22% reduction in base drag at M = 1.55 (p(sub j)/p(sub o) = 1.76) but indicated no change at M = 1.30 (p(sub j)/p(sub o) = 1.43). At M = 1.1 and p(sub j)/p(sub o) = 1.55, the jet caused a 50% increase in base drag.
Heavy-duty vehicles (HDVs) are a significant source of fuel consumption and greenhouse gas emissions, prompting solutions such as HDV platooning to mitigate these negative impacts through air drag reduction. The intervehicle distance in an HDV platoon needs to be carefully selected, such that the platoon-level energy efficiency and safety considerations can be well balanced. Underlying this problem lies in accurately modeling the relationship between HDV air drag coefficient and intervehicle distance. Through comprehensive evaluation and comparison, we analyze five control-oriented HDV air drag coefficient models, including the polynomial model, rational polynomial model, rational model, semi-quadratic model, and ridge model. Leveraging Scipy Curve-Fit toolbox and our previously compiled air drag coefficient datasets, we optimally identify the parameters inside each model. The calibrated models are then thoroughly evaluated via five complementary metrics. The comparison results reveal that the semi-quadratic model has the highest overall performance, while the widely adopted rational model only exhibits suboptimal performance.
Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.
Large, global dust storms that periodically form in the Martian atmosphere pose a threat to NASA missions due to a lack of understanding of how they will affect the heat loading and mass loss of the entry vehicle’s thermal protection system (TPS). Previous studies have laid out viable methodologies for estimating the amount of surface recession from a significant dust storm, but rely on underlying modeling assumptions that are difficult to validate. The current work utilizes the DUST library, which is a Lagrangian particle solver under development at NASA Ames Research Center. DUST includes variable-fidelity fluid-particle coupling and has been tested for a wide range of problem types. The DUST library is used in conjunction with the US3D Naiver-Stokes solver. Using these tools, a sensitivity study is performed to determine the influence of the underlying dust modeling (drag coefficient, Nusselt number, and size distribution) on the surface erosion through a hypersonic Martian entry trajectory. The surface erosion was found to be insensitive to the drag and heating models, but was affected significantly by the size distributions parameters. The total recession was comparable to the thermochemical ablation predicted by material response codes, but the lack of a particle erosion model for PICA limits the predictive capability of the current work.