Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “topological gradient”

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.

51 records · Page 3

Sizing and Layout Design of an Aeroelastic Wingbox Through Nested Optimization

The goals of this work are to 1) develop an optimization algorithm that can simultaneously handle a large number of sizing variables and topological layout variables for an aeroelastic wingbox optimization problem and 2) utilize this algorithm to ascertain the benefits of curvilinear wingbox components. The algorithm used here is a nested optimization, where the outer level optimizes the rib and skin stiffener layouts with a surrogate-based optimizer, and the inner level sizes all of the components via gradient-based optimization. Two optimizations are performed: one restricted to straight rib and stiffener components only, the other allowing curved members. A moderate 1.18% structural mass reduction is obtained through the use of curvilinear members.

Stanford, Bret K.↗

Topology of three-dimensional, variable density flows

This paper is concerned with the interpretation of unsteady, variable-density flow fields. The topology of the flow is determined by finding critical points and identifying the character of local solution trajectories. The time evolution of the flow is studied by following the paths of the critical points in the three-dimensional space of invariants of the local deformations tensor. The methodology can be applied to any smooth vector field and its associated gradient tensor including the vorticity and pressure gradient fields. This approach provides a framework for describing the geometry of complex flow patterns. Concisely summarizing that geometry in the space of invariants of the local gradient tensor may be a useful way of gaining insight into time-dependent processes described by large computational data bases. Applications to the descriptions of a flickering diffusion flame and a compressible wake are discussed.

Cantwell, Brian↗

A study of the turbulence structures of wall-bounded shear flows

This project extends the study of the structure of wall-bounded flows using the topological properties of eddying motions as developed by Chong et al. (1990), Soria et al. (1992, 1994), and as recently extended by Blackburn et al. (1996) and Chacin et al. (1996). In these works, regions of flow which are focal in nature are identified by being enclosed by an isosurface of a positive small value of the discriminant of the velocity gradient tensor. These regions resemble the attached vortex loops suggested first by Theodorsen (1955). Such loops are incorporated in the attached eddy model versions of Perry & Chong (1982), Perry et al. (1986), and Perry & Marusic (1995), which are extensions of a model first formulated by Townsend (1976). The DNS data of wall bounded flows studied here are from the zero pressure gradient flow of Spalart (1988) and the boundary layer with separation and reattachment of Na & Moin (1996). The flow structures are examined from the viewpoint of the attached eddy hypothesis.

Chong, M. S.↗

Thermoplastic Space Point Design (TSPD) Tall Tower Lunar Thermal Analysis

The Thermoplastics Development for Exploration Applications (TDEA) project seeks to advance NASA’s thermoplastic composites capabilities by developing structurally efficient joining solutions for large-scale space structures to support NASA’s future exploration missions. A key objective is to develop and understand advanced thermoplastic joining techniques relevant to space environments and applicable to unitized and/or reconfigurable composite structures. The TDEA project uses NASA Langley’s Tall Lunar Tower (TLT) concept as a target application for the design and development of in-space assembled thermoplastic composites. TLT is a 50-meter-tall tower constructed on the lunar surface near the south pole. The tower’s height at this location on the moon allows for extended durations of incident solar flux which can be used to generate power via photovoltaic arrays. The thermal environment present on lunar south pole represents a challenging environment for any space-based asset. This study provides an overview of the analysis techniques used to determine the temperature extremes, gradients, and timelines of a tall composite tower over a one-year duration at the Shackleton Connecting Ridge. The developed thermal model includes tower and lunar surface elements with anisotropic and temperature-dependent material properties and incorporates a solar vector orbit with time-varying solar flux. Study topics include: a comparison between the effects of simplified lunar surface topologies on thermal results, prediction of surface temperature values versus Lunar Reconnaissance Orbiter (LRO) flight data, and tower structural-thermal deformation.

Lunar↗

Structure of coronal neutral sheets

A qualitative model for the structure of the neutral sheet lying along the axis of coronal streamers is developed. The internal topology of the sheet is that of extremely thin magnetic tongues greatly distended outward by the solar wind expansion inside the sheet. Due to finite conductivity effects, expansion is taking place across the field lines but is retarded relative to the external flow by the reverse jxB force. The sheet thickness is determined by three considerations: the electrical conductivity that specifies the magnitude of the gradients in field strength, the expansion velocity that stretches the field lines outward decreasing the sheet thickness, and finally, the lateral pressure balance that limits the approach of the oppositely directed external field toward the neutral plane. If sigma is the electrical conductivity, the sheet thickness is shown to be proportional to sigma-1/3. For an electron conductivity evaluated perpendicular to the internal field in the sheet, the thickness is of the order of 100 km in the inner corona and 10,000 km at 1 AU. Microturbulence and instabilities are expected to yield dimensions greater than these theoretical values since these effects tend to reduce the effective conductivity.

Pneuman, G. W.↗

On the Use of CAD and Cartesian Methods for Aerodynamic Optimization

The objective for this paper is to present the development of an optimization capability for Curt3D, a Cartesian inviscid-flow analysis package. We present the construction of a new optimization framework and we focus on the following issues: 1) Component-based geometry parameterization approach using parametric-CAD models and CAPRI. A novel geometry server is introduced that addresses the issue of parallel efficiency while only sparingly consuming CAD resources; 2) The use of genetic and gradient-based algorithms for three-dimensional aerodynamic design problems. The influence of noise on the optimization methods is studied. Our goal is to create a responsive and automated framework that efficiently identifies design modifications that result in substantial performance improvements. In addition, we examine the architectural issues associated with the deployment of a CAD-based approach in a heterogeneous parallel computing environment that contains both CAD workstations and dedicated compute engines. We demonstrate the effectiveness of the framework for a design problem that features topology changes and complex geometry.

Nemec, M.↗

Aerodynamic Shape Optimization Based on Free-form Deformation

This paper presents a free-form deformation technique suitable for aerodynamic shape optimization. Because the proposed technique is independent of grid topology, we can treat structured and unstructured computational fluid dynamics grids in the same manner. The proposed technique is an alternative shape parameterization technique to a trivariate volume technique. It retains the flexibility and freedom of trivariate volumes for CFD shape optimization, but it uses a bivariate surface representation. This reduces the number of design variables by an order of magnitude, and it provides much better control for surface shape changes. The proposed technique is simple, compact, and efficient. The analytical sensitivity derivatives are independent of the design variables and are easily computed for use in a gradient-based optimization. The paper includes the complete formulation and aerodynamics shape optimization results.

Samareh, Jamshid A.↗

Deformation of the heliospheric current sheet

A kinematic analysis of the previously neglected effect of velocity inhomogeneity on the topology of the heliospheric current sheet in a radially flowing solar wind shows how the originally smooth current sheet becomes 'ruffled'. In the highly idealized case of a totally uniform, radial, steady solar wind, the shape of the current sheet is independent of distance from the sun. However, the real solar wind is inhomogeneous; the velocity varies from point to point along the current sheet, causing a distortion in the current sheet of progressively greater amplitude with increasing distance from the sun. This is true even for purely radial flow. Significant and observable distortion is produced by relatively small gradients in velocity; thus to predict or understand the shape of the heliospheric current sheet it is essential to know the solar wind in which the current sheet is embedded. Examples are given of mild velocity gradients which demonstrate the principles, the magnitude, and the character of the effect; deformation of the actual heliospheric current sheet in the highly variable solar wind is expected to be of far greater amplitude and complexity than in the simplified, tutorial examples. A new expression for the inclination of the current sheet as a function of velocity inhomogeneity and distance from the sun that is easily applied to the interpretation of solar wind data is also derived.

Suess, S. T.↗

Stress studies in EFG

Stress distributions were calculated for a creep law to predict a rate of plastic deformation. The expected reduction in stresses is obtained. Improved schemes for calculating growth system temperature distributions were evaluated. Temperature field modeling examined the possibility of using horizontal temperature gradients to influence stress distribution in ribbon. The defect structure of 10 cm wide ribbon grown in the cartridge system was examined. A new feature is identified from an examination of cross sectional micrographs. It consists of high density dislocation bands extending through the ribbon thickness. A four point bending apparatus was constructed for high temperature study of the creep response of silicon, to be used to generate defects for comparison with as grown defects in ribbon. The feasibility of laser interferometric techniques for sheet residual stress distribution measurement is examined. The mathematical formalism for calculating residual stress from changes in surface topology caused by an applied stress in a rectangular specimen was developed, and the system for laser interferometric measurement to obtain surface topology data was tested on CZ silicon.

Source record↗

HLPW-4: Wall-Modeled LES and Lattice-Boltzmann Technology Focus Group Workshop Summary

A summary of the nine submissions to the Wall-Modeled LES and Lattice-Boltzmann(WMLESLB) Technical Focus Group (TFG) at the 4th High lift Prediction Workshop is provided. The focus of this TFG was to assess the current capabilities of WMLES and LB methods on a complex high-lift configuration across a wide range of angles-of-attack. Analysis of the submitted data suggests that>250𝑀spatial degrees of freedom are needed to accurately predict pitching moments at high angles-of-attack due to large pressure gradients present on the outboard slat and main element for 𝛼 >17◦(corrected for free-air). While some scatter is reported in pitching moment coefficient at the low-angles of attack (𝛼 <11◦) - excellent agreement is observed between submissions near the 𝐶𝐿,max state. Objective superiority of WMLES methods over RANS can be seen in terms of lack of excess outboard separation; a majority of the WMLES and LB submissions predict wedge-shaped separation patterns consistent with the experimental oil flow. The in-tunnel simulations show excellent agreement with the experiment in terms of a) integrated loads, b) surface flow-topology, and c) mechanism for the onset of inboard stall. Further evidence is provided to demonstrate both qualitative and quantitative superiority of the WMLES submissions over RANS.

TTT↗

The Dynamics of Miscible Interfaces: A Space Flight Experiment

Experiments as well as accompanying simulations are described that serve in preparation of a space flight experiment to study the dynamics of miscible interfaces. The investigation specifically addresses the importance of both nonsolenoidal effects as well as nonconventional Korteweg stresses in flows that give rise to steep but finite concentration gradients. The investigation focuses on the flow in which a less viscous fluid displaces one of higher viscosity and different density within a narrow capillary tube. The fluids are miscible in all proportions. An intruding finger forms that occupies a fraction of the total tube diameter. Depending on the flow conditions, as expressed by the Peclet number, a dimensionless viscosity ratio, and a gravity parameter, this fraction can vary between approximately 0.9 and 0.2. For large Pe values, a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the specific flow conditions, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip. For small Pe values, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). Comparisons between experiments and numerical simulations based on the 'conventional' assumption of solenoidal velocity fields and without Korteweg stresses yield poor agreement as far as the Pe value is concerned that distinguishes these two regimes. As one possibility, we attribute this lack of agreement to the disregard of these terms. An attempt is made to use scaling arguments in order to evaluate the importance of the Korteweg stresses and of the assumption of solenoidality. While these effects should be strongest in absolute terms when steep concentration fronts exist, i.e., at large Pe, they may be relatively most important at lower values of Pe. We subsequently compare these conventional simulations to more complete simulations that account for nonvanishing divergence as well as Korteweg stresses. While the exact value of the relevant stress coefficients are not known, ballpark numbers do exist, and their use in the simulations indicates that these stresses may indeed be important. We plan to evaluate these issues in detail by means of comparing a space experiment with corresponding simulations, in order to extract more accurate Korteweg stress coefficients, and to confirm or deny the importance of such stresses.

Maxworthy, Tony↗

HLPW-4/GMGW-3: Wall-Modeled LES and Lattice-Boltzmann Technology Focus Group Workshop Summary

A summary of the nine submissions to the Wall-Modeled LES and Lattice-Boltzmann(WMLESLB) Technical Focus Group (TFG) at the 4th High lift Prediction Workshop is provided. The focus of this TFG was to assess the current capabilities of WMLES and Lattice-Boltzmann methods on a complex high-lift configuration across a wide range of angles of attack. Preliminary analysis of the submitted data suggests that>250𝑀spatial degrees of freedom are needed to accurately predict pitching moments at high angles of attack due to large pressure gradients present on the outboard slat and main element for𝛼 >17◦(corrected for free-air).While some scatter is reported in pitching moments at the low-angles of attack (𝛼 <11◦) - likely caused by differences in flap separation possibly due to low Reynolds number effects - excellent agreement is observed between the submissions near the𝐶𝐿,max state. Objective superiority over steady state RANS can be seen in terms of lack of excess outboard separation; a majority of the good quality WMLES and LB submissions predict wedge-shaped separation patterns consistent with the experimental oil flow. Differences in the onset of stall mechanism in the free-air configuration for𝛼 >20◦is reported with two distinct topologies observed. Topology A is characterized by the onset of corner-flow separation which progressively grows to produce a pitch break in free-air with an angle of attack increase from𝛼=20.55◦to𝛼=21.47◦+𝜖where𝜖is a small perturbation (varying between submissions and likely to be within±0.3◦). Topology Bis characterized by boundary layer weakness emanating from the inboard side of the wing-pylon juncture substantially larger than any weakness in the wing-body juncture (due to corner-flow);submissions in this category do not show any tendency for occurrence of a pitch break for the free-air configuration within the vicinity of𝛼=21.47◦. The in-tunnel simulations submitted by3 participants using different discretizations, grids and closure models show excellent agreement with the experiment in terms of a) integrated loads, b) surface flow-topology, and c) mechanism for the onset of inboard stall. Further evidence is provided to demonstrate both qualitative and quantitative superiority of all 3 WMLES submissions over a single steady state RANS submission to the workshop.

TTT↗

Force balance in the magnetospheres of Jupiter and Saturn

Spacecraft measurements of the plasma populations and magnetic fields near Jupiter and Saturn have revealed that large magnetospheres surround both planets. Magnetic field measurements have indicated closed field line topologies in the dayside magnetospheres of both planets while plasma instruments have shown these regions to be populated by both hot and cold plasma components convected azimuthally in the sense of planetary rotation. By using published data from the Voyager Plasma Science (PLS), Low Energy Charged Particle (LECP), and Magnetometer (MAG) instruments, it is possible to investigate the validity of the time stationary MHD momentum equation in the middle magnetospheres of Jupiter and Saturn. At Saturn, the hot plasma population is negligible in the dynamic sense and the centrifugal force of the cold rotating plasma appears to balance the Lorentz force. At Jupiter, the centrifugal force balances about 25 percent of the Lorentz force. The remaining inward Lorentz force is balanced by pessure gradients in the hot, high-beta plasma of the Jovian magnetodisk.

Mcnutt, R. L., Jr.↗

Pioneer 10 and 11 Data Analysis

This report finishes the work of NASA Grant NAS2-153, which supported data analysis for the UCSD instruments on Pioneers 10 and 11. The data analyzed under this grant span 22 years of interplanetary measurements in the inner and outer heliosphere. The UCSD instruments made their mark in cosmic ray research based upon their high energy thresholds, directional responses, and reliable data streams. one of their primary scientific objectives concerns the size, configuration, and time behavior of the heliosphere. The size scale is inferred from the radial intensity gradient, which is measured between the two spacecraft and extrapolated to interstellar intensity levels at the cosmic ray modulation boundary. This boundary still eludes us, and its position, motion, and the best method of extrapolation are ongoing problems. Current projections place the boundary beyond 100 AU, which may be beyond the termination shock, and raises the question of possible modulation in the heliosheath. Probably our only hope of seeing this region in the immediate future rides on the possibility that the boundary will move inward. Our instruments have recorded many Forbush, or transient, decreases in the outer heliosphere. These observations led us to a model that attributes many of the decreases to solar wind stream-stream interactions, and relates the cosmic ray variations to the locally observed magnetic field magnitude. As the cosmic ray variations in this model result only from topological changes in the modulation integral, the model is a tool for studying the possibility that the 11 year cosmic ray modulation cycle can be accounted for by a superposition of Forbush decreases. The cosmic ray angular distribution function is measurable, given a good telemetry rate, by the UCSD Cerenkov detector which counts particles of energy greater than 500 MeV/n. We obtained statistically significant samples from 1 to 9 AU, at 13 AU, and at 34 AU. The anisotropy tends to be a few tenths of a per cent at all radial distances. A quasiperiodic variation in the east-west anisotropy with period of about 50 days remains unexplained.

Fillius, Walker↗

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗