Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smoothed boundary method”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

On three dimensional object recognition and pose-determination: An abstraction based approach

A method of computing reliable Gaussian and mean curvature sign-map descriptors from the polynomial approximation of surfaces was demonstrated. Such descriptors which are invariant under perspective variation are suitable for hypothesis generation. A means for determining the pose of constructed geometric forms whose algebraic surface descriptors are nonlinear in terms of their orienting parameters was developed. This was done by means of linear functions which are capable of approximating nonlinear forms and determining their parameters. It was shown that biquadratic surfaces are suitable companion linear forms for cylindrical approximation and parameter estimation. The estimates provided the initial parametric approximations necessary for a nonlinear regression stage to fine tune the estimates by fitting the actual nonlinear form to the data. A hypothesis-based split-merge algorithm for extraction and pose determination of cylinders and planes which merge smoothly into other surfaces was developed. It was shown that all split-merge algorithms are hypothesis-based. A finite-state algorithm for the extraction of the boundaries of run-length regions was developed. The computation takes advantage of the run list topology and boundary direction constraints implicit in the run-length encoding.

Quek, Kok How Francis↗
NASA NTRSDate not supplied

Development of a C-ELS Specimen-Based Numerical Benchmark for Mode II Delamination and Assessment of Two VCCT-Based Propagation Strategies

A finite element (FE) benchmark example inspired by the calibrated end-loaded split (C-ELS) specimen is developed and used to assess the performance of delamination propagation capabilities based on linear elastic fracture mechanics (LEFM). The C-ELS specimen has the advantage of a longer region of stable delamination propagation compared to the existing mode II benchmark case. The new benchmark example may therefore provide a better assessment tool by enabling more stable crack growth in regions further away from the boundary conditions or load application. First, a benchmark result is created manually using two-dimensional finite element models of the C-ELS specimen with different delamination lengths. Second, the performance of the virtual crack closure technique (VCCT) delamination propagation capabilities in the Abaqus/Standard®1 FE code and the recently developed Progressive Release eXplicit-VCCT (PRX-VCCT) method are assessed by comparing the results to the benchmark case. Two examples with different starter delamination lengths are studied. A shorter starter length is chosen to create a scenario with unstable delamination propagation. A longer delamination encourages stable delamination propagation. Detailed results from three-dimensional analyses with aligned and misaligned meshes and two levels of mesh refinement are provided. In general, good agreement can be achieved between the results obtained from the quasi-static propagation analysis and the benchmark analysis. Numerical artifacts including anomalous unreleased nodes in the crack wake and zig-zag crack fronts occur for propagation analyses using Abaqus/Standard VCCT. In comparison, continuous, smooth, delamination fronts are observed for PRX-VCCT. The use of the benchmark case to assess different VCCT-based propagation strategies illustrates the value of establishing benchmark cases.

Composite Materials↗
NASA NTRS2021

Micromechanics-based Modeling of Laminated SiC/SiC Ceramic Matrix Composites

The behavior and response of ceramic matrix composites (CMCs), in particular silicon carbide fiber reinforced silicon carbide matrix (SiC/SiC), is affected by many factors such as variation of fiber volume fraction, residual stresses resulting from processing of the composites at high temperature, random microstructures, and the presence of matrix flaws (e.g., voids, pores, cracks etc.) as well as general material nonlinearity and heterogeneity that occurs randomly in a composite. Residual stresses arising from the phase change of constituents are evaluated in this paper and it is shown that they do influence composite strength and need to be properly accounted for. Additionally, the microstructures (location of fiber centers, coating thickness etc.) of advanced CMCs are usually disordered (or random) and fiber diameter and strength typically have a distribution. They rarely resemble the ordered fiber packing (square, rectangular, or hexagonal) that is generally assumed in micromechanics-based models with periodic boundary conditions for computational expediency. These issues raise the question of how should one model such systems effectively? Can an ordered hexagonal packed repeating unit cell (RUC) accurately represent the random microstructure behavior? How many fibers need to be included to enable accurate representation? Clearly, the number of fibers within an RUC must be limited to insure a balance between accuracy and efficiency. NASA’s in-house micromechanics-based code MAC/GMC provides a framework to analyze such RUCs for the overall composite behavior and the FEAMAC computer code provides linkage of MAC/GMC to the commercial FEA code, ABAQUS. The appropriate level of discretization of the RUC as well as the analysis method employed, i.e., Generalized Method of Cells (GMC) or High Fidelity Generalized Method of Cells (HFGMC), is investigated in this paper in the context of a unidirectional as well as a cross-ply laminated CMC. Results including effective composite properties, proportional limit stress (an important design parameter) and fatigue are shown utilizing both GMC as well as HFGMC. Finally, a few multiscale analyses are performed on smooth bar test coupons as well as test coupons with features such as open-hole and double notches using FEAMAC. Best practices and guidance are provided to take these phenomena into account and keep a proper balance between fidelity (accuracy) and efficiency. Following these guidelines can account for important physics of the problem and provide significant advantages when performing large multiscale composite structural analyses.

Ceramic Matrix Composites↗
NASA NTRS2017

Design of an Axisymmetric Afterbody Test Case for CFD Validation

As identified in the CFD Vision 2030 Study commissioned by NASA, validation of advanced RANS models and scale-resolving methods for computing turbulent flow fields must be supported by continuous improvements in fundamental, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, the underpinnings of a new test platform referred to herein as the NASA Axisymmetric Afterbody are presented. The devised body-of-revolution is a modular platform consisting of a forebody section and afterbody section, allowing for a range of flow behaviors to be studied on interchangeable afterbody geometries. A body-of-revolution offers advantages in shape definition and fabrication, in avoiding direct contact with wind tunnel sidewalls, and in tail-sting integration to facilitate access to higher Reynolds number tunnels. The current work is focused on validation of smooth-body turbulent flow separation, for which a six-parameter body has been developed. A priori RANS computations are reported for a risk-reduction test configuration in order to demonstrate critical variation among turbulence model results for a given afterbody, ranging from barely-attached to mild separated flow. RANS studies of the effects of forebody nose (with/without) and wind tunnel boundary (slip/no-slip) on the selected afterbody are presented. Representative modeling issues that can be explored with this configuration are the effect of higher Reynolds number on separation behavior, flow physics of the progression from attached to increasingly-separated afterbody flows, and the effect of embedded longitudinal vortices on turbulence structure.

Disotell, Kevin J.↗
DOE OSTI2023

Machine learning unifies flexibility and efficiency of spinodal structure generation for stochastic biomaterial design

Abstract Porous biomaterials design for bone repair is still largely limited to regular structures (e.g. rod-based lattices), due to their easy parameterization and high controllability. The capability of designing stochastic structure can redefine the boundary of our explorable structure–property space for synthesizing next-generation biomaterials. We hereby propose a convolutional neural network (CNN) approach for efficient generation and design of spinodal structure—an intriguing structure with stochastic yet interconnected, smooth, and constant pore channel conducive to bio-transport. Our CNN-based approach simultaneously possesses the tremendous flexibility of physics-based model in generating various spinodal structures (e.g. periodic, anisotropic, gradient, and arbitrarily large ones) and comparable computational efficiency to mathematical approximation model. We thus successfully design spinodal bone structures with target anisotropic elasticity via high-throughput screening, and directly generate large spinodal orthopedic implants with desired gradient porosity. This work significantly advances stochastic biomaterials development by offering an optimal solution to spinodal structure generation and design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Preliminary large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids

Large-eddy simulation (LES) has matured to the point where application to complex flows is desirable. The extension to higher Reynolds numbers leads to an impractical number of grid points with existing structured-grid methods. Furthermore, most real world flows are rather difficult to represent geometrically with structured grids. Unstructured-grid methods offer a release from both of these constraints. However, just as it took many years for structured-grid methods to be well understood and reliable tools for LES, unstructured-grid methods must be carefully studied before we can expect them to attain their full potential. In the past two years, important building blocks have been put into place making possible a careful study of LES on unstructured grids. The first building block was an efficient mesh generator which allowed the placement of points according to smooth variation of physical length scales. This variation of length scales is in all three directions independently, which allows a large reduction in points when compared to structured-grid methods, which can only vary length scales in one direction at a time. The second building block was the development of a dynamic model appropriate for unstructured grids. The principle obstacle was the development of an unstructured-grid filtering operator. In the past year, some of the new filters developed by Jansen have been implemented into a highly parallelized finite element code based on the Galerkin/least-squares finite element method. We have chosen the NACA 4412 airfoil at maximum lift as the first simulation for a variety of reasons. First, it is a problem of significant interest since it would be the first LES of an aircraft component. Second, this flow has been the subject of three experimental studies. The third reason for considering this flow is the variety of flow features which provide an important test of the dynamic model. Only the dynamic model can be expected to perform satisfactorily in this variety of situations: from the laminar regions where it must not modify the flow at all to the turbulent boundary layers and wake where it must represent a wide variety of subgrid-scale structures. The flow configuration we have chosen is that of Wadcock (1987) at Reynolds number based on chord Re(sub c) = u(sub infinity)c/v = 1.64 x 10(exp 6), Mach number M = 0.2, and 12 deg angle of attack.

Jansen, Kenneth↗
NASA NTRS1984

Real Surface Effects in Elastohydrodynamic Lubrication

Following a general discussion of real effects in EHL, consideration is given in detail to the role played by stochastic roughness superimposed on the nominal solid boundaries. Particular attention is given to the full-film EHL regime where incipient asperity contact bears a negligible fraction of the load. Flow, from which an averaged Reynolds equation can be formed, is nonetheless modified by the amplitude and texture of the roughness patterns. In describing these effects by means of certain well-defined flow factors, it is found that only two extra parameters are needed - the rms surface height and the ratio of the correlation lengths in the two principal roughness directions. In a lowest order perturbation expansion of these flow factors in powers of the ratio of rms roughness to nominal film thickness, no other properties of the roughness appear. In most cases of practical interest, the factors describing Poiseuille flow are separable into the sum of two single-surface flow factors which means that a combination of a single equivalent rough surface versus an ideal smooth surface can always be found. For flow entrained by slip velocity the factors separate instead into a difference of single-surface factors and it becomes significant which of the two surfaces carries the equivalent roughness. Results are discussed of some applications of the averaged Reynolds equation based on the flow factor method to the EHL line contact problem. Finally, the partial EHL regime is considered where comparable load fractions are carried by the hydrodynamic film and by incipient mechanical contact. An extension of the method into this regime by combining it with asperity contact models appears most encouraging.

John H Tripp↗
NASA NTRS1960

Investigation by Schlieren Technique of Methods of Fixing Fully Turbulent Flow on Models at Supersonic Speeds

An investigation has been made of the effectiveness of two-dimensional wire and three-dimensional grandular-type roughnesses for fixing transition on a 27 deg cone at a free-stream Mach number of 2.20. One height of each type of roughness was investigated over a free-stream Reynolds number per foot range of about 1.0 x 10 (exp 6) to 7.0 x 10 (exp 6). The tests were made at an angle of attack of 0 deg and with essentially zero heat transfer. The results indicate that, within the Reynolds number per foot range investigated, both two-dimensional wire and three-dimensional granular types of roughness caused large forward movements of transition at constant Reynolds number per foot; however, the forward movement with increase in Reynolds number per foot was very slow for the two-dimensional wire roughness and transition never reached the roughness location as it did with the three-dimensional granular roughness. An extrapolated value of the roughness Reynolds number (based on roughness height and flow conditions at the outer edge of the roughness) for the first appearance of forward movement from the smooth-cone case for the three-dimensional granular roughness was in good agreement with the value of about 600 for the first appearance of turbulent bursts obtained on a similar cone by the hot-wire technique. The range of values of roughness Reynolds number for which a fully turbulent boundary layer existed in the vicinity of the roughness element was about 4,500 to 6,500.

Jackson, Mary W.↗
NASA NTRS1994

Experimental studies of transpiration cooling with shock interaction in hypersonic flow, part B

This report describes the result of experimental studies conducted to examine the effects of the impingement of an oblique shock on the flowfield and surface characteristics of a transpiration-cooled wall in turbulent hypersonic flow. The principal objective of this work was to determine whether the interaction between the oblique shock and the low-momentum region of the transpiration-cooled boundary layer created a highly distorted flowfield and resulted in a significant reduction in the cooling effectiveness of the transpiration-cooled surface. As a part of this program, we also sought to determine the effectiveness of transpiration cooling with nitrogen and helium injectants for a wide range of blowing rates under constant-pressure conditions in the absence of shock interaction. This experimental program was conducted in the Calspan 48-Inch Shock Tunnel at nominal Mach numbers of 6 and 8, for a Reynolds number of 7.5 x 10(exp 6). For these test conditions, we obtained fully turbulent boundary layers upstream of the interaction regions over the transpiration-cooled segment of the flat plate. The experimental program was conducted in two phases. In the first phase, we examined the effects of mass-addition level and coolant properties on the cooling effectiveness of transpiration-cooled surfaces in the absence of shock interaction. In the second phase of the program, we examined the effects of oblique shock impingement on the flowfield and surface characteristics of a transpiration-cooled surface. The studies were conducted for a range of shock strengths with nitrogen and helium coolants to examine how the distribution of heat transfer and pressure and the characteristics of the flowfield in the interaction region varied with shock strength and the level of mass addition from the transpiration-cooled section of the model. The effects of the distribution of the blowing rate along the interaction regions were also examined for a range of blowing rates through the transpiration-cooled panels. The regions of shockwave/boundary layer interaction examined in these studies were induced by oblique shocks generated with a sharp, flat plate, inclined to the freestream at angles of 5 degrees, 7.5 degrees, and 10 degrees. It was found that, in the absence of an incident shock, transpiration cooling was a very effective method for reducing both the heat transfer and the skin friction loads on the surface. The helium coolant was found to be significantly more effective than nitrogen, because of its low molecular weight and high specific heat. The studies of shock-wave/transpiration-cooled surface interaction demonstrated that the interaction region between the incident shock and the low-momentum transpiration-cooled boundary layer did not result in a significant increase in the size of attached or separated interaction regions, and did not result in significant flowfield distortions above the interaction region. The increase in heating downstream of the shock-impingement point could easily be reduced to the values without shock impingement by a relatively small increase in the transpiration cooling in this region. Surprisingly, this increase in cooling rate did not result in a significant increase in size of the region ahead of the incident shock or create a significantly enlarged interaction region with a resultant increase in the distortion level in the inviscid flow. Thus, transpiration cooling appears to be a very effective technique to cool the internal surfaces of scramjet engines, where shocks in the engine would induce large local increases in wall heating and create viscous/inviscid interactions that could significantly disturb the smooth flow through the combustor. However, if hydrogen is used as the coolant, burning upstream of shock impingement might result in localized hot spots. Clearly, further research is needed in this area.

Holden, Michael S.↗
NASA NTRS2004

Three-Dimensional High-Order Spectral Volume Method for Solving Maxwell's Equations on Unstructured Grids

A three-dimensional, high-order, conservative, and efficient discontinuous spectral volume (SV) method for the solutions of Maxwell's equations on unstructured grids is presented. The concept of discontinuous 2nd high-order loca1 representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) method, but instead of using a Galerkin finite-element formulation, the SV method is based on a finite-volume approach to attain a simpler formulation. Conventional unstructured finite-volume methods require data reconstruction based on the least-squares formulation using neighboring cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In the SV method, one starts with a relatively coarse grid of triangles or tetrahedra, called spectral volumes (SVs), and partition each SV into a number of structured subcells, called control volumes (CVs), that support a polynomial expansion of a desired degree of precision. The unknowns are cell averages over CVs. If all the SVs are partitioned in a geometrically similar manner, the reconstruction becomes universal as a weighted sum of unknowns, and only a few universal coefficients need to be stored for the surface integrals over CV faces. Since the solution is discontinuous across the SV boundaries, a Riemann solver is thus necessary to maintain conservation. In the paper, multi-parameter and symmetric SV partitions, up to quartic for triangle and cubic for tetrahedron, are first presented. The corresponding weight coefficients for CV face integrals in terms of CV cell averages for each partition are analytically determined. These discretization formulas are then applied to the integral form of the Maxwell equations. All numerical procedures for outer boundary, material interface, zonal interface, and interior SV face are unified with a single characteristic formulation. The load balancing in a massive parallel computing environment is therefore easier to achieve. A parameter is introduced in the Riemann solver to control the strength of the smoothing term. Important aspects of the data structure and its effects to communication and the optimum use of cache memory are discussed. Results will be presented for plane TE and TM waves incident on a perfectly conducting cylinder for up to fifth order of accuracy, and a plane wave incident on a perfectly conducting sphere for up to fourth order of accuracy. Comparisons are made with exact solutions for these cases.

Liu, Yen↗
NASA NTRS2005

The Effects of the Uncertainty of Departures on Multi-Center Traffic Management Advisor (TMA) Scheduling

The Multi-center Traffic Management Advisor (McTMA) provides a platform for regional or national traffic flow management, by allowing long-range cooperative time-based metering to constrained resources, such as airports or air traffic control center boundaries. Part of the demand for resources is made up of proposed departures, whose actual departure time is difficult to predict. For this reason, McTMA does not schedule the departures in advance, but rather relies on traffic managers to input their requested departure time. Because this happens only a short while before the aircraft's actual departure, McTMA is unable to accurately predict the amount of delay airborne aircraft will need to take in order to accommodate the departures. The proportion of demand which is made up by such proposed departures increases as the horizon over which metering occurs gets larger. This study provides an initial analysis of the severity of this problem in a 400-500 nautical mile metering horizon and discusses potential solutions to accommodate these departures. The challenge is to smoothly incorporate departures with the airborne stream while not excessively delaying the departures.' In particular, three solutions are reviewed: (1) scheduling the departures at their proposed departure time; (2) not scheduling the departures in advance; and (3) scheduling the departures at some time in the future based on an estimated error in their proposed time. The first solution is to have McTMA to automatically schedule the departures at their proposed departure times. Since the proposed departure times are indicated in their flight times in advance, this method is the simplest, but studies have shown that these proposed times are often incorrect2 The second option is the current practice, which avoids these inaccuracies by only scheduling aircraft when a confirmed prediction of departure time is obtained from the tower of the departure airport. Lastly, McTMA can schedule the departures at a predicted departure time based on statistical data of past departure time performance. It has been found that departures usually have a wheels-up time after their indicated proposed departure time, as shown in Figure 1. Hence, the departures were scheduled at a time in the future based on the mean error in proposed departure times for their airport.

Thipphavong, Jane↗
NASA NTRS2000

Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-capturing Methods

In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Yee, H. C.↗
NASA NTRS1996

Independent Pixel and Two Dimensional Estimates of LANDSAT-Derived Cloud Field Albedo

A theoretical study has been conducted on the effects of cloud horizontal inhomogeneity on cloud albedo bias. A two-dimensional (2D) version of the Spherical Harmonic Discrete Ordinate Method (SHDOM) is used to estimate the albedo bias of the plane parallel (PP-IPA) and independent pixel (IPA-2D) approximations for a wide range of 2D cloud fields obtained from LANDSAT. They include single layer trade cumulus, open and closed cell broken stratocumulus, and solid stratocumulus boundary layer cloud fields over ocean. Findings are presented on a variety of averaging scales and are summarized as a function of cloud fraction, mean cloud optical depth, cloud aspect ratio, standard deviation of optical depth, and the gamma function parameter Y (a measure of the width of the optical depth distribution). Biases are found to be small for small cloud fraction or mean optical depth, where the cloud fields under study behave linearly. They are large (up to 0.20 for PP-IPA bias, -0.12 for IPA-2D bias) for large v. On a scene average basis PP-IPA bias can reach 0.30, while IPA-2D bias reaches its largest magnitude at -0.07. Biases due to horizontal transport (IPA-2D) are much smaller than PP-IPA biases but account for 20% RMS of the bias overall. Limitations of this work include the particular cloud field set used, assumptions of conservative scattering, constant cloud droplet size, no gas absorption or surface reflectance, and restriction to 2D radiative transport. The LANDSAT data used may also be affected by radiative smoothing.

Chambers, L. H.↗
NASA NTRS1995

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗
NASA NTRS2007

Building Operations Efficiencies into NASA's Ares I Crew Launch Vehicle Design

The U.S. Vision for Space Exploration guides the National Aeronautics and Space Administration's (NASA's) challenging missions that expand humanity's boundaries and open new routes to the space frontier. With the Agency's commitment to complete the International Space Station (ISS) and to retire the venerable Space Shuttle by 2010, the NASA Administrator commissioned the Exploration Systems Architecture Study (ESAS) in 2005 to analyze options for safe, simple, cost-efficient launch solutions that could deliver human-rated space transportation capabilities in a timely manner within fixed budget guidelines. The Exploration Launch Projects (ELP) Office, chartered by the Constellation Program in October 2005, has been conducting systems engineering studies and business planning to successively refine the design configurations and better align vehicle concepts with customer and stakeholder requirements, such as significantly reduced life-cycle costs. As the Agency begins the process of replacing the Shuttle with a new generation of spacecraft destined for missions beyond low-Earth orbit to the Moon and Mars, NASA is designing the follow-on crew and cargo launch systems for maximum operational efficiencies. To sustain the long-term exploration of space, it is imperative to reduce the $4 billion NASA typically spends on space transportation each year. This paper gives toplevel information about how the follow-on Ares I Crew Launch Vehicle (CLV) is being designed for improved safety and reliability, coupled with reduced operations costs. These methods include carefully developing operational requirements; conducting operability design and analysis; using the latest information technology tools to design and simulate the vehicle; and developing a learning culture across the workforce to ensure a smooth transition between Space Shuttle operations and Ares vehicle development.

Dumbacher, Daniel L.↗
NASA NTRS1994

Optimized design of a hypersonic nozzle

Conventional procedures for designing nozzles involve the design of an inviscid contour (using the method of characteristics) that is corrected with a displacement thickness calculated from boundary-layer theory. However, nozzles designed using this classical procedure have been shown to exhibit poor flow quality at Mach numbers characteristic of hypersonic applications. The nozzle to be designed will be a part of the NASA HYPULSE facility which is being used for hypervelocity flight research. Thus, the flow quality of the nozzle is a critical question that needs to be addressed. Design of nozzles for hypersonic applications requires a proper assessment of the effects of the thick boundary layer on the inviscid flowfield. Since the flow field is largely supersonic, the parabolized form of the Navier-Stokes (PNS) equations can be used. The requirement of a uniform flow at the exit plane of the nozzle can be used to define an objective function as part of an optimization procedure. The design procedure used in this study involves the coupling of a nonlinear (least-squares) optimization algorithm with an efficient, explicit PNS solver. The thick boundary layers growing on the walls of the nozzle limit the extent of the usable core region (region with uniform flow) for testing models (especially rectangular). In order to maximize the region of uniform flow, it was decided to have the exit plane of this nozzle to be (nearly) rectangular. Thus, an additional constraint on the nozzle shape resulted, namely the nozzle will have a shape transitioning from a circular one at the inlet to that of a rectangle at the exit. In order to provide for a smooth shape transition, the cross sectional contour of the nozzle is defined by a superellipse. The nozzle is taken to be a meter in length. The axial variations of the major and minor radii of the superellipse are governed by cubic splines. The design parameters are the coefficients of the splines associated with the local nozzle wall slopes. Extensive calculations have been made (with a three-dimensional Euler code) to understand the effects of various parameters such as location of the knot points of the spline function, different ways of characterizing the uniformity of the flow in the exit plane, as well as the effect of constraining the area of the nozzle to be invariant. Turbulent flow (measurements indicate that the flow at the nozzle inlet is turbulent) calculations are now being performed (with the inviscidly designed nozzle contours) to assess the flow quality.

Krishnamurthy, Ramesh↗