Downloads & Free Reading Options - Results

We introduce a modified and simplified version of the pre-existing fully parallelized three-dimensional Navier--Stokes flow solver known as TPLS. We demonstrate how the simplified version can be used as a pedagogical tool for the study of computational fluid dynamics and parallel computing. TPLS is at its heart a two-phase flow solver, and uses calls to a range of external libraries to accelerate its performance. However, in the present context we narrow the focus of the study to basic hydrodynamics and parallel computing techniques, and the code is therefore simplified and modified to simulate pressure-driven single-phase flow in a channel, using only relatively simple Fortran 90 code with MPI parallelization, but no calls to any other external libraries. The modified code is analysed in order to both validate its accuracy and investigate its scalability up to 1000 CPU cores. Simulations are performed for several benchmark cases in pressure-driven channel flow, including a turbulent simulation, wherein the turbulence is incorporated via the large-eddy simulation technique. The work may be of use to advanced undergraduate and graduate students as an introductory study in computational fluid dynamics, while also providing insight for those interested in more general aspects of high-performance computing.

A computational study was conducted to better understand experimental results obtained from wind tunnel tests of a Mach 4 waverider model and a comparative reference configuration. The experimental results showed that the performance of the reference configuration was slightly better than that of the waverider model. These results contradict waverider design theory, which suggests that a waverider optimized for maximum lift-to-drag should provide better performance than any other non-waverider configuration at a given design point, especially at hypersonic speeds. The computational results showed that the predicted surface pressure values and the integrated lift and drag coefficients from the pressure distributions were much lower for the reference model than for the flat-top model, due to the reference model bottom surface having a slight expansion. The lift-to-drag ratios for the flat-top model were higher due to a relatively low drag for the same amount of lift. These results indicate that the performance advantage of the reference model was due to the shape of the bottom surface and not due to the flat top surface. The results also showed that the reference model exhibited the same shock attachment characteristics as the waverider because the planform shapes were identical. CFD predictions show that the planform shape gives the waverider an advantage in performance over conventional hypersonic vehicles and that altering the bottom surface of a waverider does not cause significant performance degradation.

This article is from Journal of Human Kinetics , volume 29 . Abstract Nowadays the underwater gliding after the starts and the turns plays a major role in the overall swimming performance. Hence, minimizing hydrodynamic drag during the underwater phases should be a main aim during swimming. Indeed, there are several postures that swimmers can assume during the underwater gliding, although experimental results were not conclusive concerning the best body position to accomplish this aim. Therefore, the purpose of this study was to analyse the effect in hydrodynamic drag forces of using different body positions during gliding through computational fluid dynamics (CFD) methodology. For this purpose, two-dimensional models of the human body in steady flow conditions were studied. Two-dimensional virtual models had been created: (i) a prone position with the arms extended at the front of the body; (ii) a prone position with the arms placed alongside the trunk; (iii) a lateral position with the arms extended at the front and; (iv) a dorsal position with the arms extended at the front. The drag forces were computed between speeds of 1.6 m/s and 2 m/s in a two-dimensional Fluent® analysis. The positions with the arms extended at the front presented lower drag values than the position with the arms aside the trunk. The lateral position was the one in which the drag was lower and seems to be the one that should be adopted during the gliding after starts and turns.

This article is from Journal of Cardiovascular Magnetic Resonance , volume 16 . Abstract None

This article is from Journal of Cardiovascular Magnetic Resonance , volume 16 . Abstract None

This article is from BioMedical Engineering OnLine , volume 13 . Abstract Background: Left pulmonary artery sling (LPAS) is a rare but severe congenital anomaly, in which the stenoses are formed in the trachea and/or main bronchi. Multi-detector computed tomography (MDCT) provides useful anatomical images, but does not offer functional information. The objective of the present study is to quantitatively analyze the airflow in the trachea and main bronchi of LPAS subjects through computational fluid dynamics (CFD) simulation. Methods: Five subjects (four LPAS patients, one normal control) aging 6-19 months are analyzed. The geometric model of the trachea and the two main bronchi is extracted from the MDCT images. The inlet velocity is determined based on the body weight and the inlet area. Both the geometric model and personalized inflow conditions are imported into CFD software, ANSYS. The pressure drop, mass flow ratio through two bronchi, wall pressure, flow velocity and wall shear stress (WSS) are obtained, and compared to the normal control. Results: Due to the tracheal and/or bronchial stenosis, the pressure drop for the LPAS patients ranges 78.9 - 914.5 Pa, much higher than for the normal control (0.7 Pa). The mass flow ratio through the two bronchi does not correlate with the sectional area ratio if the anomalous left pulmonary artery compresses the trachea or bronchi. It is suggested that the C-shaped trachea plays an important role on facilitating the air flow into the left bronchus with the inertia force. For LPAS subjects, the distributions of velocities, wall pressure and WSS are less regular than for the normal control. At the stenotic site, high velocity, low wall pressure and high WSS are observed. Conclusions: Using geometric models extracted from CT images and the patient-specified inlet boundary conditions, CFD simulation can provide vital quantitative flow information for LPAS. Due to the stenosis, high pressure drops, inconsistent distributions of velocities, wall pressure and WSS are observed. The C-shaped trachea may facilitate a larger flow of air into the left bronchus under the inertial force, and decrease the ventilation of the right lung. Quantitative and personalized information may help understand the mechanism of LPAS and the correlations between stenosis and dyspnea, and facilitate the structural and functional assessment of LPAS.

Computational Fluid Dynamics (CFD) Modeling of Heat Transfer in a Polymeric-Membrane using Finite Volume Method 作者： Ahsan, M. 1 Hussain, A. 1 作者单位： 1. NUST, SCME, Islamabad, Pakistan. 提交时间： 2017-11-02 摘要: The efficiency, robustness and reliability of recent numerical methods for finding solutions to flow problems have given rise to the implementation of computational fluid dynamics (CFD) as a broadly used analysis method for engineering problems like membrane separation system. The CFD modeling in this study observes steady and unsteady (transient) heat flux and temperature profiles in a polymeric (cellulose acetate) membrane. This study is novel due to the implementation of user defined scalar (LIDS) diffusion equation by using user-defined functions (UDFs) infinite volume method (FVM). Some details of the FVM used by the solver are carefully discussed when implementing terms in the governing equation and boundary conditions (BC). The contours of temperature due to high-temperature gradient are reported for steady and unsteady problems. Computational fluid dynamics finite volume method heat transfer cellulose acetate based membrane numerical-simulation turbulence promoters narrow channel 分类： 物理学 >> 普通物理:统计和量子力学,量子信息等 引用： ChinaXiv:201711.00015 (或此版本 ChinaXiv:201711.00015V1 ) doi:10.12074/201711.00015 CSTR:32003.36.ChinaXiv.201711.00015.V1 推荐引用方式： Ahsan, M.,Hussain, A..(2017).Computational Fluid Dynamics (CFD) Modeling of Heat Transfer in a Polymeric-Membrane using Finite Volume Method.热科学学报.[ChinaXiv:201711.00015] 版本历史 [V1] 2017-11-02 14:52:11 ChinaXiv:201711.00015V1 下载全文

A systematic approach is presented for developing entropy stable (SS) formulations of any order for the Navier-Stokes equations. These SS formulations discretely conserve mass, momentum, energy and satisfy a mathematical entropy inequality. They are valid for smooth as well as discontinuous flows provided sufficient dissipation is added at shocks and discontinuities. Entropy stable formulations exist for all diagonal norm, summation-by-parts (SBP) operators, including all centered finite-difference operators, Legendre collocation finite-element operators, and certain finite-volume operators. Examples are presented using various entropy stable formulations that demonstrate the current state-of-the-art of these schemes.

Experimental techniques to measure rotorcraft aerodynamic performance are widely used. However, most of them are either unable to capture interference effects from bodies, or require an extremely large computational budget. The objective of the present research is to develop an XV-15 Tilt Rotor Research Aircraft rotor model for investigation of wind tunnel wall interference using a novel Computational Fluid Dynamics (CFD) solver for rotorcraft, RotCFD. In RotCFD, a mid-fidelity URANS solver is used with an incompressible flow model and a realizable k-e turbulence model. The rotor is, however, not modeled using a computationally expensive, unsteady viscous body-fitted grid, but is instead modeled using a blade element model with a momentum source approach. Various flight modes of the XV-15 isolated rotor, including hover, tilt and airplane mode, have been simulated and correlated to existing experimental and theoretical data. The rotor model is subsequently used for wind tunnel wall interference simulations in the National Full-Scale Aerodynamics Complex (NFAC) at NASA Ames Research Center in California. The results from the validation of the isolated rotor performance showed good correlation with experimental and theoretical data. The results were on par with known theoretical analyses. In RotCFD the setup, grid generation and running of cases is faster than many CFD codes, which makes it a useful engineering tool. Performance predictions need not be as accurate as high-fidelity CFD codes, as long as wall effects can be properly simulated. For both test sections of the NFAC wall interference was examined by simulating the XV-15 rotor in the test section of the wind tunnel and with an identical grid but extended boundaries in free field. Both cases were also examined with an isolated rotor or with the rotor mounted on the modeled geometry of the Tiltrotor Test Rig (TTR). A 'quasi linear trim' was used to trim the thrust for the rotor to compare the power as a unique variable. Power differences between free field and wind tunnel cases were found from -7 % to 0 % in the 80- by 120-Foot Wind Tunnel test section and -1.6 % to 4.8 % in the 40- by 80-Foot Wind Tunnel, depending on the TTR orientation, tunnel velocity and blade setting. The TTR will be used in 2016 to test the Bell 609 rotor in a similar fashion to the research in this report.

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

An experimental and analytical program was conducted to investigate heat transfer and pressure losses in rotating multipass passages with configurations and dimensions typical of modern turbine blades. The objective of this program is the development and verification of improved analysis methods that will form the basis for a design system that will produce turbine components with improved durability. As part of this overall program, a technique is developed for computational fluid dynamics. The specific objectives were to: select a baseline CFD computer code, assess the limitations of the baseline code, modify the baseline code for rotational effects, verify the modified code against benchmark experiments in the literature, and to identify shortcomings in the code as revealed by the verification. The Pratt and Whitney 3D-TEACH CFD code was selected as the vehicle for this program. The code was modified to account for rotating internal flows, and these modifications were evaluated for flow characteristics of those expected in the application. Results can make a useful contribution to blade internal cooling.

This work involves the coordination of necessary resources, facilities, and special personnel to provide a workshop to promote the exchange of CFD technology between industry, universities, and government. Critical flow problems have been isolated and simulation of these is being done.

SAGE is a user-friendly, highly efficient, two-dimensional self-adaptive grid code based on Nakahashi and Deiwert's variational principles method. Grid points are redistributed into regions of high flowfield gradients while maintaining smoothness and orthogonality of the grid. Efficiency is obtained by splitting the adaption into 2 directions and applying one-sided torsion control, thus producing a 1-D elliptic system that can be solved as a set of tridiagonal equations.

The methodology is described for converting a large, long-running applications code that executed on a single processor of a CRAY-2 supercomputer to a version that executed efficiently on multiple processors. Although the conversion of every application is different, a discussion of the types of modification used to achieve gigaflop performance is included to assist others in the parallelization of applications for CRAY computers, especially those that were developed for other computers. An existing application, from the discipline of computational fluid dynamics, that had utilized over 2000 hrs of CPU time on CRAY-2 during the previous year was chosen as a test case to study the effectiveness of multitasking on a CRAY-2. The nature of dominant calculations within the application indicated that a sustained computational rate of 1 billion floating-point operations per second, or 1 gigaflop, might be achieved. The code was first analyzed and modified for optimal performance on a single processor in a batch environment. After optimal performance on a single CPU was achieved, the code was modified to use multiple processors in a dedicated environment. The results of these two efforts were merged into a single code that had a sustained computational rate of over 1 gigaflop on a CRAY-2. Timings and analysis of performance are given for both single- and multiple-processor runs.

Thermal Protection System (TPS) Cavity Heating is predicted using Computational Fluid Dynamics (CFD) on unstructured grids for both simplified cavities and actual cavity geometries. Validation was performed using comparisons to wind tunnel experimental results and CFD predictions using structured grids. Full-scale predictions were made for simplified and actual geometry configurations on the Space Shuttle Orbiter in a mission support timeframe.

This paper covers the development of an integrated nonlinear dynamic simulation for a variable cycle turbofan engine and nozzle that can be integrated with an overall vehicle Aero-Propulso-Servo-Elastic (APSE) model. A previously developed variable cycle turbofan engine model is used for this study and is enhanced here to include variable guide vanes allowing for operation across the supersonic flight regime. The primary focus of this study is to improve the fidelity of the model's thrust response by replacing the simple choked flow equation convergent-divergent nozzle model with a MacCormack method based quasi-1D model. The dynamic response of the nozzle model using the MacCormack method is verified by comparing it against a model of the nozzle using the conservation element/solution element method. A methodology is also presented for the integration of the MacCormack nozzle model with the variable cycle engine.

A system and method for generating fluid flow parameter data for use in aerodynamic heating analysis. Computational fluid dynamics data is generated for a number of points in an area on a surface to be analyzed. Sub-areas corresponding to areas of the surface for which an aerodynamic heating analysis is to be performed are identified. A computer system automatically determines a sub-set of the number of points corresponding to each of the number of sub-areas and determines a value for each of the number of sub-areas using the data for the sub-set of points corresponding to each of the number of sub-areas. The value is determined as an average of the data for the sub-set of points corresponding to each of the number of sub-areas. The resulting parameter values then may be used to perform an aerodynamic heating analysis.

In phase II, additional inducer rotations are simulated in order to understand the root cause of the flowliner crack problem. CFD results confirmed that there is a strong unsteady interaction between the backflow regions caused by the LPFTP inducer and secondary flow regions in the bellows cavity through the flowliner slots. It is observed that the swirl on the duct side of the downstream flowliner is stronger than on the duct side of the upstream flowliner. Due to this swirl, there are more significant unsteady flow interactions through the downstream slots than those observed in the upstream slots. Averaged values of the local velocity at the slots were provided to the NESC-ITA flow physics acoustics team to guide them in designing the acoustics experiment. A parametric study was performed to compare the flow field in the flowliner area when one upstream slot and one corresponding downstream slot were enlarged. No significant differences were observed between the flow field obtained from the enlarged slot configuration when compared with the original configuration. More cases must be analyzed with various enlarged slot configurations to generalize this observation. The flow through the A1 test stand and the flow through the orbiter fuel feedline manifold were simulated without the LPFTP. It was observed that incoming flow to the flowliner and inducer was more uniform in the A1 test stand then in the orbiter manifold. Additionally, each engine LPFTP in the orbiter receives significantly different velocity distributions. Because of the differences observed in the computed results, it is not possible for the A1 test stand to represent the three different engine feedlines simultaneously.

This slide presentation reviews the use of computational fluid dynamics in performing analysis of the space shuttle with particular reference to the return to flight analysis and other shuttle problems. Slides show a comparison of pressure coefficient with the shuttle ascent configuration between the wind tunnel test and the computed values. the evolution of the grid system for the space shuttle launch vehicle (SSLv) from the early 80's to one in 2004, the grid configuration of the bipod ramp redesign from the original design to the current configuration, charts with the computations showing solid rocket booster surface pressures from wind tunnel data, calculated over two grid systems (i.e., the original 14 grid system, and the enhanced 113 grid system), and the computed flight orbiter wing loads are compared with strain gage data on STS-50 during flight. The loss of STS-107 initiated an unprecedented review of all external environments. The current SSLV grid system of 600+ grids, 1.8 Million surface points and 95+ million volume points is shown. The inflight entry analyses is shown, and the use of Overset CFD as a key part to many external tank redesign and debris assessments is discussed. The work that still remains to be accomplished for future shuttle flights is discussed.

Computational fluid dynamics (CFD) and other advanced modeling and simulation (M&S) methods are increasingly relied on for predictive performance, reliability and safety of engineering systems. Analysts, designers, decision makers, and project managers, who must depend on simulation, need practical techniques and methods for assessing simulation credibility. The AIAA Guide for Verification and Validation of Computational Fluid Dynamics Simulations (AIAA G-077-1998 (2002)), originally published in 1998, was the first engineering standards document available to the engineering community for verification and validation (V&V) of simulations. Much progress has been made in these areas since 1998. The AIAA Committee on Standards for CFD is currently updating this Guide to incorporate in it the important developments that have taken place in V&V concepts, methods, and practices, particularly with regard to the broader context of predictive capability and uncertainty quantification (UQ) methods and approaches. This paper will provide an overview of the changes and extensions currently underway to update the AIAA Guide. Specifically, a framework for predictive capability will be described for incorporating a wide range of error and uncertainty sources identified during the modeling, verification, and validation processes, with the goal of estimating the total prediction uncertainty of the simulation. The Guide's goal is to provide a foundation for understanding and addressing major issues and concepts in predictive CFD. However, this Guide will not recommend specific approaches in these areas as the field is rapidly evolving. It is hoped that the guidelines provided in this paper, and explained in more detail in the Guide, will aid in the research, development, and use of CFD in engineering decision-making.

In aerospace computational fluid dynamics (CFD) calculations, the Delaunay triangulation of suitable quadrilateral meshes can lead to unsuitable triangulated meshes. Here, we present case studies which illustrate the limitations of using structured grid generation methods which produce points in a curvilinear coordinate system for subsequent triangulations for CFD applications. We discuss conditions under which meshes of quadrilateral elements may not produce a Delaunay triangulation suitable for CFD calculations, particularly with regard to high aspect ratio, skewed quadrilateral elements.

Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This paper describes an approach to quantify the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft without the use of test data. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional "validation by test only" mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions. Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computational Fluid Dynamics can be used to verify these requirements; however, the model must be validated by test data. This research includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available and open source solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT, STARCCM+, and OPENFOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid Dynamics model using the methodology found in "Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations". This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System /spacecraft system. Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. For the flow regime being analyzed (turbulent, three-dimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

Design of Experiments (DOE) was applied to the LAS geometric parameter study to efficiently identify and rank primary contributors to integrated drag over the vehicles ascent trajectory in an order of magnitude fewer CFD configurations thereby reducing computational resources and solution time. SME s were able to gain a better understanding on the underlying flowphysics of different geometric parameter configurations through the identification of interaction effects. An interaction effect, which describes how the effect of one factor changes with respect to the levels of other factors, is often the key to product optimization. A DOE approach emphasizes a sequential approach to learning through successive experimentation to continuously build on previous knowledge. These studies represent a starting point for expanded experimental activities that will eventually cover the entire design space of the vehicle and flight trajectory.

Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This proposal describes an approach to validate the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft. The research described here is absolutely cutting edge. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional"validation by test only'' mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions. Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computationaf Fluid Dynamics can be used to veritY these requirements; however, the model must be validated by test data. The proposed research project includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT and OPEN FOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid . . . Dynamics model using the methodology found in "Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations". This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System /spacecraft system. Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. To date, the author is the only person to look at the uncertainty in the entire computational domain. For the flow regime being analyzed (turbulent, threedimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

There have been few discussions on using Computational Fluid Dynamics (CFD) without experimental validation. Pairing experimental data, uncertainty analysis, and analytical predictions provides a comprehensive approach to verification and is the current state of the art. With pressed budgets, collecting experimental data is rare or non-existent. This paper investigates and proposes a method to perform CFD uncertainty analysis only from computational data. The method uses current CFD uncertainty techniques coupled with the Student-T distribution to predict the heat transfer coefficient over a at plate. The inputs to the CFD model are varied from a specified tolerance or bias error and the difference in the results are used to estimate the uncertainty. The variation in each input is ranked from least to greatest to determine the order of importance. The results are compared to heat transfer correlations and conclusions drawn about the feasibility of using CFD without experimental data. The results provide a tactic to analytically estimate the uncertainty in a CFD model when experimental data is unavailable

Automation of flow field zoning in two dimensions is an important step towards reducing the difficulty of three-dimensional grid generation in computational fluid dynamics. Using a knowledge-based approach makes sense, but problems arise which are caused by aspects of zoning involving perception, lack of expert consensus, and design processes. These obstacles are overcome by means of a simple shape and configuration language, a tunable zoning archetype, and a method of assembling plans from selected, predefined subplans. A demonstration system for knowledge-based two-dimensional flow field zoning has been successfully implemented and tested on representative aerodynamic configurations. The results show that this approach can produce flow field zonings that are acceptable to experts with differing evaluation criteria.

In the photograph, the contour lines indicate temperature at Mach 19 around a generic vehicle similar to the X-30. Because of high pressure and skin friction, the temperature is highest on the surface of the model. The vivid blue color indicates temperature 18 times hotter than the atmosphere, requiring active cooling from within the aircraft to control temperature on the planes surface. Pink color on the outer rings is much cooler than the blue near the body, but is still 2 times hotter than the environment in which the plane travels.

The process of remote laser welding is simulated in this study to identify the keyhole-induced porosity generation mechanisms and keyhole. Three processes are simulated and discussed: laser power levels, laser-beam shaping configurations, and laser keyhole process. The simulation finding reveals that pore development is caused by strong melt flow behind the keyhole. As verification, the equivalent experimental test is also carried out. According to the findings, a welding speed with a high level helps to keep the keyholes released and prevents the flow of strong melt; a big advanced leaning-angle also provides inactive molten pool flow, making it difficult for bubbles to float to the backside of the molten pool. The conclusions of this study offer crucial insight into the method of porosity of aluminum (Al) alloys laser welding, as well as advice on how to avoid keyhole-induced porosity. It is also obtained that a smaller laser beam with constant power raises the velocity, welding pool depth, and liquid metal temperature.

The Fluid Dynamics Panel of AGARD arranged a Symposium on Applications of Computational Fluid Dynamics in Aeronautics, on 7 to 10 April 1986 in Aix-en-Provence, France. The purpose of the Symposium was to provide an assessment of the status of CFD in aerodynamic design and analysis, with an emphasis on emerging applications of advanced computational techniques to complex configurations. Sessions were devoted specifically to grid generation, methods for inviscid flows, calculations of viscous-inviscid interactions, and methods for solving the Navier-Stokes equations. The 31 papers presented at the meeting are published in AGARD Conference Proceedings CP-412 and are listed in the Appendix of this report. A brief synopsis of each paper and some general conclusions and recommendations are given.

With the recent interest in novel control effectors there is a need to determine the stability and control derivatives of new aircraft configurations early in the design process. These derivatives are central to most control law design methods and would allow the determination of closed-loop control performance of the vehicle. Early determination of the static and dynamic behavior of an aircraft may permit significant improvement in configuration weight, cost, stealth, and performance through multidisciplinary design. The classical method of determining static stability and control derivatives - constructing and testing wind tunnel models - is expensive and requires a long lead time for the resultant data. Wind tunnel tests are also limited to the preselected control effectors of the model. To overcome these shortcomings, computational fluid dynamics (CFD) solvers are augmented via automatic differentiation, to directly calculate the stability and control derivatives. The CFD forces and moments are differentiated with respect to angle of attack, angle of sideslip, and aircraft shape parameters to form these derivatives. A subset of static stability and control derivatives of a tailless aircraft concept have been computed by two differentiated inviscid CFD codes and verified for accuracy with central finite-difference approximations and favorable comparisons to a simulation database.

This document in large part is based on the Masters Thesis of Cole Stephens. The document encompasses a variety of technical and practical issues involved when using the STARS codes for Aeroservoelastic analysis of vehicles. The document covers in great detail a number of technical issues and step-by-step details involved in the simulation of a system where aerodynamics, structures and controls are tightly coupled. Comparisons are made to a benchmark experimental program conducted at NASA Langley. One of the significant advantages of the methodology detailed is that as a result of the technique used to accelerate the CFD-based simulation, a systems model is produced which is very useful for developing the control law strategy, and subsequent high-speed simulations.

The field of computational fluid dynamics (CFD) has advanced to the point where it can now be used for the purpose of fluid dynamics physics education. Because of the tremendous wealth of information available from numerical simulation, certain fundamental concepts can be efficiently communicated using an interactive graphical interrogation of the appropriate numerical simulation data base. In other situations, a large amount of aerodynamic information can be communicated to the student by interactive use of simple CFD tools on a workstation or even in a personal computer environment. The emphasis in this presentation is to discuss ideas for how this process might be implemented. Specific examples, taken from previous publications, will be used to highlight the presentation.

A CFD analysis of a low pressure nuclear rocket concept is presented with the use of an advanced chemical kinetics, Navier-Stokes code. The computations describe the flow field in detail, including gas dynamic, thermodynamic and chemical properties, as well as global performance quantities such as specific impulse. Computational studies of several rocket nozzle shapes are conducted in an attempt to maximize hydrogen recombination. These Navier-Stokes calculations, which include real gas and viscous effects, predict lower performance values than have been reported heretofore.

An investigation was conducted from October 1, 1990 to May 31, 1994 on the development of methodologies to improve the designs (more specifically, the shape) of aerodynamic surfaces of coupling optimization algorithms (OA) with Computational Fluid Dynamics (CFD) algorithms via sensitivity analyses (SA). The study produced several promising methodologies and their proof-of-concept cases, which have been reported in the open literature.

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

A combined geometric approach for computational fluid dynamics is presented for the analysis of unsteady flow about mechanisms in which its components are in moderate relative motion. For a CFD analysis, the total dynamics problem involves the dynamics of the aspects of geometry modeling, grid generation, and flow modeling. The interrelationships between these three aspects allow for a more natural formulation of the problem and the sharing of information which can be advantageous to the computation of the dynamics. The approach is applied to planar geometries with the use of an efficient multi-block, structured grid generation method to compute unsteady, two-dimensional and axisymmetric flow. The applications presented include the computation of the unsteady, inviscid flow about a hinged-flap with flap deflections and a high-speed inlet with centerbody motion as part of the unstart / restart operation.

Computational fluid dynamics (CFD) is necessary in the design processes of all current aerospace plane programs. Single-stage-to-orbit (STTO) aerospace planes with air-breathing supersonic combustion are going to be largely designed by means of CFD. The challenge of the aerospace plane design is to provide credible CFD results to work from, to assess the risk associated with the use of those results, and to certify CFD codes that produce credible results. To establish the credibility of CFD results used in design, the following topics are discussed: CFD validation vis-a-vis measurable fluid dynamics (MFD) validation; responsibility for credibility; credibility requirement; and a guide for establishing credibility. Quantification of CFD uncertainties helps to assess success risk and safety risks, and the development of CFD as a design tool requires code certification. This challenge is managed by designing the designers to use CFD effectively, by ensuring quality control, and by balancing the design process. For designing the designers, the following topics are discussed: how CFD design technology is developed; the reasons Japanese companies, by and large, produce goods of higher quality than the U.S. counterparts; teamwork as a new way of doing business; and how ideas, quality, and teaming can be brought together. Quality control for reducing the loss imparted to the society begins with the quality of the CFD results used in the design process, and balancing the design process means using a judicious balance of CFD and MFD.

Applications of computational aerodynamics to aeronautical research, design, and analysis have increased rapidly over the past decade, and these applications offer significant benefits to aeroelasticians. The past developments are traced by means of a number of specific examples, and the trends are projected over the next several years. The crucial factors that limit the present capabilities for unsteady analyses are identified; they include computer speed and memory, algorithm and solution methods, grid generation, turbulence modeling, vortex modeling, data processing, and coupling of the aerodynamic and structural dynamic analyses. The prospects for overcoming these limitations are presented, and many improvements appear to be readily attainable. If so, a complete and reliable numerical simulation of the unsteady, transonic viscous flow around a realistic fighter aircraft configuration could become possible within the next decade. The possibilities of using artificial intelligence concepts to hasten the achievement of this goal are also discussed.

Approximations applicable to the radiating, reacting, and conducting stagnation region of a hypervelocity vehicle were incorporated into a method for rapidly obtaining approximate solutions. This solution utilizes a coordinate system based upon the origin of the radiative losses and includes in a phenomenologically correct manner the effects of chemical and thermal nonequilibrium, and nonequilibrium, nongray radiative transfer. Results were presented which demonstrate the usefulness of the method and indicate which radiation parameters require further study and definition. Excellent comparisons were obtained with published results for the Fire2 data. An axisymmetric nonequilibrium inverse method was modified and extended and used to investigate and compare various vibration dissociation chemistry coupling models and radiative heat transfer approximations. The similarities, differences, and consequences of using these models in the Aero-assist Orbital Transfer Vehicles flight regime will be discussed.

Since the Columbia accident, the threat to the Shuttle launch vehicle from debris during the liftoff timeframe has been assessed by the Liftoff Debris Team at NASA/MSFC. In addition to engineering methods of analysis, CFD-generated flow fields during the liftoff timeframe have been used in conjunction with 3-DOF debris transport methods to predict the motion of liftoff debris. Early models made use of a quasi-steady flow field approximation with the vehicle positioned at a fixed location relative to the ground; however, a moving overset mesh capability has recently been developed for the Loci/CHEM CFD software which enables higher-fidelity simulation of the Shuttle transient plume startup and liftoff environment. The present work details the simulation of the launch pad and mobile launch platform (MLP) with truncated solid rocket boosters (SRBs) moving in a prescribed liftoff trajectory derived from Shuttle flight measurements. Using Loci/CHEM, time-accurate RANS and hybrid RANS/LES simulations were performed for the timeframe T0+0 to T0+3.5 seconds, which consists of SRB startup to a vehicle altitude of approximately 90 feet above the MLP. Analysis of the transient flowfield focuses on the evolution of the SRB plumes in the MLP plume holes and the flame trench, impingement on the flame deflector, and especially impingment on the MLP deck resulting in upward flow which is a transport mechanism for debris. The results show excellent qualitative agreement with the visual record from past Shuttle flights, and comparisons to pressure measurements in the flame trench and on the MLP provide confidence in these simulation capabilities.

The Human Powered Centrifuge (HPC) is a facility that is planned to be installed on board the International Space Station (ISS) to enable crew exercises under the artificial gravity conditions. The HPC equipment includes a "bicycle" for long-term exercises of a crewmember that provides power for rotation of HPC at a speed of 30 rpm. The crewmember exercising vigorously on the centrifuge generates the amount of carbon dioxide of about two times higher than a crewmember in ordinary conditions. The goal of the study is to analyze the airflow and carbon dioxide distribution within Pressurized Multipurpose Module (PMM) cabin when HPC is operating. A full unsteady formulation is used for airflow and CO2 transport CFD-based modeling with the so-called sliding mesh concept when the HPC equipment with the adjacent Bay 4 cabin volume is considered in the rotating reference frame while the rest of the cabin volume is considered in the stationary reference frame. The rotating part of the computational domain includes also a human body model. Localized effects of carbon dioxide dispersion are examined. Strong influence of the rotating HPC equipment on the CO2 distribution detected is discussed.

This paper covers the development of an integrated nonlinear dynamic simulation for a variable cycle turbofan engine and nozzle that can be integrated with an overall vehicle Aero-Propulso-Servo-Elastic (APSE) model. A previously developed variable cycle turbofan engine model is used for this study and is enhanced here to include variable guide vanes allowing for operation across the supersonic flight regime. The primary focus of this study is to improve the fidelity of the model's thrust response by replacing the simple choked flow equation convergent-divergent nozzle model with a MacCormack method based quasi-1D model. The dynamic response of the nozzle model using the MacCormack method is verified by comparing it against a model of the nozzle using the conservation element/solution element method. A methodology is also presented for the integration of the MacCormack nozzle model with the variable cycle engine.

The Air Vehicles Unstructured Solver (AVUS) is a three-dimensional (3-D) flow simulator running on unstructured grids. Embedded in the simulator is a sparse matrix package for solving coupled partial differential equations. Unstructured grid solvers usually do not exhibit optimal computer memory access patterns. The purpose of this report is to detail the memory access patterns and provide a roadmap for improving it.

Combined Computational Fluid Dynamics Reynolds Averaged Navier-Stokes and LES model calculations were completed for comparisons with wind tunnel measurements of flow along the afterbody and in the base region of a circular cylinder aligned to a Mach 2.5 freestream. This cylinder is representative of a missile at zero degrees angle-of-attack. Comparison data includes afterbody and base pressure, high quality non-intrusive Laser Doppler Velocimetry and Pressure Sensitive Paint measurements, and turbulent kinetic energy. The flow in the base/separated region of a tactical missile is very complex even though the geometry is quite simple. There are several reasons this is the case. First, the external flow over the missile body is supersonic while the flow in the base region is mixed supersonic/subsonic flow. Second, this flow is characterized by shock wave/boundary layer interactions at the lip of the base region. Third, the flow is turbulent in the base/separated region and highly transient in nature. There has been a long history of attempts to accurately describe this flow going back to the 1950s where integral models were used and continuing to the 2000s where sophisticated computation fluid dynamic models were used. However, none of these models have produced results that provide accurate descriptions of the flowfield.

This report discusses the continuing efforts to port the F3D computational fluid dynamics code to RISC-based SMPs. Originally, this program was optimized for Cray vector supercomputers such as the Cray C90. Previous attempts to run this code on 561 Power Challenges, and Convex Exemplars, as well as systems from SUN and Digital Equipment, demonstrated a level of performance that was so low as to be utterly useless (in many cases, it became necessary to kill the job before the first time step had completed). After making a concerted effort to port the program to an 561 Power Challenge (R8000 processor), acceptable levels of performance were finally achieved (Pressel 1997). Using this version of the code as the starting point, an effort was made to produce a program that ran efficiently on both systems from 561 and Convex. Unfortunately, a number of limitations with the Convex Exemplar were discovered that limited the success of this effort.

Computational fluid dynamics approaches were used to compute the flow fields and aerodynamic forces and moments of a slender-body finned projectile. Steady-state numerical results have been obtained for a series of cases, with Mach numbers ranging from 1.5 to 5.0 and at angles of attack from 0 deg to 5 deg. Full three-dimensional computations were performed using a cubic k-epsilon turbulence model. In general, the computed aerodynamic coefficients compared well with available semi-empirical data.

In Europe, as in the US, the use of the computer to solve complex problems in fluid dynamics is a burgeoning field of endeavor. This report briefly describes some of the current European efforts in computational fluid dynamics. The citations are intended as 'leads' to the individuals and institutions involved in such activities, and recent references are given to provide guidance to the latest published information.

The research reported herein addresses the reduced basis method, which is based on the combination of a basic iteractive method and a projection method. The convergence is analyzed and some error bounds are established. The relationship between the reduced basis method and the preconditioned conjugate gradient method is discussed. Also included is a practical implementation of the reduced basis method when a pseudoresidual-based Krylov space is chosen. The final section of this report dwells on the weight selection procedures of hybrid difference methods for the linear convective equation. The procedures are based on the ability of hybrid difference methods to conserve the discrete weighted energy. Keywords: Reduced basis; Linear solver; Contraction number; Projection method; Sparse symmetric linear systems; Hybrid difference methods.

Greater computational power is needed for solving Computational Fluid Dynamics (CFD) problems of interest in engineering design. Parallel architecture computers offer the promise of providing orders of magnitude greater computational power. In this paper we quantify that promise by considering an explicit CFD method and analyze the potential parallelism for three different parallel computer architectures. The use of an explicit method gives us a bast case analysis from the point of view of parallelism, and allows us to uncover potential problems in exploiting significant parallelism. The analysis is validated against experiments on three representative parallel computers. The results allow us to predict the performance of different parallel architectures. In particular, our results show that distributed memory parallel processors offer greater potential speedup. We discuss the import of our model for the development of parallel CFD algorithms and parallel computers. We also discuss our experiences in converting our model code to run on the three different parallel computers.

This article is from Pulmonary Circulation , volume 2 . Abstract Previous clinical studies in pulmonary arterial hypertension (PAH) have concentrated predominantly on distal pulmonary vascular resistance, its contribution to the disease process, and response to therapy. However, it is well known that biomechanical factors such as shear stress have an impact on endothelial health and dysfunction in other parts of the vasculature. This study tested the hypothesis that wall shear stress is reduced in the proximal pulmonary arteries of PAH patients with the belief that reduced shear stress may contribute to pulmonary endothelial cell dysfunction and as a result, PAH progression. A combined MRI and computational fluid dynamics (CFD) approach was used to construct subject-specific pulmonary artery models and quantify flow features and wall shear stress (WSS) in five PAH patients with moderate-to-severe disease and five age- and sex-matched controls. Three-dimensional model reconstruction showed PAH patients have significantly larger main, right, and left pulmonary artery diameters (3.5 ± 0.4 vs. 2.7 ± 0.1 cm, P = 0.01; 2.5 ± 0.4 vs. 1.9 ± 0.2 cm, P = 0.04; and 2.6 ± 0.4 vs. 2.0 ± 0.2 cm, P = 0.01, respectively), and lower cardiac output (3.7 ± 1.2 vs. 5.8 ± 0.6 L/min, P = 0.02.). CFD showed significantly lower time-averaged central pulmonary artery WSS in PAH patients compared to controls (4.3 ± 2.8 vs. 20.5 ± 4.0 dynes/cm2, P = 0.0004). Distal WSS was not significantly different. A novel method of measuring WSS was utilized to demonstrate for the first time that WSS is altered in some patients with PAH. Using computational modeling in patient-specific models, WSS was found to be significantly lower in the proximal pulmonary arteries of PAH patients compared to controls. Reduced WSS in proximal pulmonary arteries may play a role in the pathogenesis and progression of PAH. This data may serve as a basis for future in vitro studies of, for example, effects of WSS on gene expression.