Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex geometry”

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 1,135 records · Page 63

An alternative to unstructured grids for computing gas dynamic flows around arbitrarily complex two-dimensional bodies

In this paper we describe an approach for dealing with arbitrary complex, two dimensional geometries, the so-called cartesian boundary method. Conceptually, the cartesian boundary method is quite simple. Solid bodies blank out areas of a background, cartesian mesh, and the resultant cut cells are singled out for special attention. However, there are several obstacles that must be overcome in order to achieve a practical scheme. We present a general strategy that overcomes these obstacles, together with some details of our successful conversion of an adaptive mesh algorithm from a body-fitted code to a cartesian boundary code.

Quirk, James J.↗

A hybrid structured-unstructured grid method for unsteady turbomachinery flow computations

A hybrid grid technique for the solution of 2D, unsteady flows is developed. This technique is capable of handling complex, multiple component geometries in relative motion, such as those encountered in turbomachinery. The numerical approach utilizes a mixed structured-unstructured zonal grid topology along with modeling equations and solution methods that are most appropriate in the individual domains, therefore combining the advantages of both structured and unstructured grid techniques.

Mathur, Sanjay R.↗

Dynamic and Thermal Turbulent Time Scale Modelling for Homogeneous Shear Flows

A new turbulence model, based upon dynamic and thermal turbulent time scale transport equations, is developed and applied to homogeneous shear flows with constant velocity and temperature gradients. The new model comprises transport equations for k, the turbulent kinetic energy; tau, the dynamic time scale; k(sub theta), the fluctuating temperature variance; and tau(sub theta), the thermal time scale. It offers conceptually parallel modeling of the dynamic and thermal turbulence at the two equation level, and eliminates the customary prescription of an empirical turbulent Prandtl number, Pr(sub t), thus permitting a more generalized prediction capability for turbulent heat transfer in complex flows and geometries. The new model also incorporates constitutive relations, based upon invariant theory, that allow the effects of nonequilibrium to modify the primary coefficients for the turbulent shear stress and heat flux. Predictions of the new model, along with those from two other similar models, are compared with experimental data for decaying homogeneous dynamic and thermal turbulence, homogeneous turbulence with constant temperature gradient, and homogeneous turbulence with constant temperature gradient and constant velocity gradient. The new model offers improvement in agreement with the data for most cases considered in this work, although it was no better than the other models for several cases where all the models performed poorly.

Schwab, John R.↗

Numerical Study of Boundary Layer Interaction with Shocks: Method Improvement and Test Computation

The objective is the development of a high-order and high-resolution method for the direct numerical simulation of shock turbulent-boundary-layer interaction. Details concerning the spatial discretization of the convective terms can be found in Adams and Shariff (1995). The computer code based on this method as introduced in Adams (1994) was formulated in Cartesian coordinates and thus has been limited to simple rectangular domains. For more general two-dimensional geometries, as a compression corner, an extension to generalized coordinates is necessary. To keep the requirements or limitations for grid generation low, the extended formulation should allow for non-orthogonal grids. Still, for simplicity and cost efficiency, periodicity can be assumed in one cross-flow direction. For easy vectorization, the compact-ENO coupling algorithm as used in Adams (1994) treated whole planes normal to the derivative direction with the ENO scheme whenever at least one point of this plane satisfied the detection criterion. This is apparently too restrictive for more general geometries and more complex shock patterns. Here we introduce a localized compact-ENO coupling algorithm, which is efficient as long as the overall number of grid points treated by the ENO scheme is small compared to the total number of grid points. Validation and test computations with the final code are performed to assess the efficiency and suitability of the computer code for the problems of interest. We define a set of parameters where a direct numerical simulation of a turbulent boundary layer along a compression corner with reasonably fine resolution is affordable.

Adams, N. A.↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

Analysis of New Composite Architectures

Efficient and accurate specialty finite elements methods to analyze textile composites were developed and are described. Textile composites present unique challenges to the analyst because of the large, complex 'microstructure'. The geometry of the microstructure is difficult to model and it introduces unusual free surface effects. The size of the microstructure complicates the use of traditional homogenization methods. The methods developed constitute considerable progress in addressing the modeling difficulties. The details of the methods and attended results obtained therefrom, are described in the various chapters included in Part 1 of the report. Specific conclusions and computer codes generated are included in Part 2 of the report.

Whitcomb, John D.↗

EUROPA2: Plan Database Services for Planning and Scheduling Applications

NASA missions require solving a wide variety of planning and scheduling problems with temporal constraints; simple resources such as robotic arms, communications antennae and cameras; complex replenishable resources such as memory, power and fuel; and complex constraints on geometry, heat and lighting angles. Planners and schedulers that solve these problems are used in ground tools as well as onboard systems. The diversity of planning problems and applications of planners and schedulers precludes a one-size fits all solution. However, many of the underlying technologies are common across planning domains and applications. We describe CAPR, a formalism for planning that is general enough to cover a wide variety of planning and scheduling domains of interest to NASA. We then describe EUROPA(sub 2), a software framework implementing CAPR. EUROPA(sub 2) provides efficient, customizable Plan Database Services that enable the integration of CAPR into a wide variety of applications. We describe the design of EUROPA(sub 2) from the perspective of both modeling, customization and application integration to different classes of NASA missions.

Bedrax-Weiss, Tania↗

Aerodynamic Optimization of Rocket Control Surface Geometry Using Cartesian Methods and CAD Geometry

Aerodynamic design is an iterative process involving geometry manipulation and complex computational analysis subject to physical constraints and aerodynamic objectives. A design cycle consists of first establishing the performance of a baseline design, which is usually created with low-fidelity engineering tools, and then progressively optimizing the design to maximize its performance. Optimization techniques have evolved from relying exclusively on designer intuition and insight in traditional trial and error methods, to sophisticated local and global search methods. Recent attempts at automating the search through a large design space with formal optimization methods include both database driven and direct evaluation schemes. Databases are being used in conjunction with surrogate and neural network models as a basis on which to run optimization algorithms. Optimization algorithms are also being driven by the direct evaluation of objectives and constraints using high-fidelity simulations. Surrogate methods use data points obtained from simulations, and possibly gradients evaluated at the data points, to create mathematical approximations of a database. Neural network models work in a similar fashion, using a number of high-fidelity database calculations as training iterations to create a database model. Optimal designs are obtained by coupling an optimization algorithm to the database model. Evaluation of the current best design then gives either a new local optima and/or increases the fidelity of the approximation model for the next iteration. Surrogate methods have also been developed that iterate on the selection of data points to decrease the uncertainty of the approximation model prior to searching for an optimal design. The database approximation models for each of these cases, however, become computationally expensive with increase in dimensionality. Thus the method of using optimization algorithms to search a database model becomes problematic as the number of design variables is increased.

Nelson, Andrea↗

Sub-Grid Modeling of Electrokinetic Effects in Micro Flows

Advances in micro-fabrication processes have generated tremendous interests in miniaturizing chemical and biomedical analyses into integrated microsystems (Lab-on-Chip devices). To successfully design and operate the micro fluidics system, it is essential to understand the fundamental fluid flow phenomena when channel sizes are shrink to micron or even nano dimensions. One important phenomenon is the electro kinetic effect in micro/nano channels due to the existence of the electrical double layer (EDL) near a solid-liquid interface. Not only EDL is responsible for electro-osmosis pumping when an electric field parallel to the surface is imposed, EDL also causes extra flow resistance (the electro-viscous effect) and flow anomaly (such as early transition from laminar to turbulent flow) observed in pressure-driven microchannel flows. Modeling and simulation of electro-kinetic effects on micro flows poses significant numerical challenge due to the fact that the sizes of the double layer (10 nm up to microns) are very thin compared to channel width (can be up to 100 s of m). Since the typical thickness of the double layer is extremely small compared to the channel width, it would be computationally very costly to capture the velocity profile inside the double layer by placing sufficient number of grid cells in the layer to resolve the velocity changes, especially in complex, 3-d geometries. Existing approaches using "slip" wall velocity and augmented double layer are difficult to use when the flow geometry is complicated, e.g. flow in a T-junction, X-junction, etc. In order to overcome the difficulties arising from those two approaches, we have developed a sub-grid integration method to properly account for the physics of the double layer. The integration approach can be used on simple or complicated flow geometries. Resolution of the double layer is not needed in this approach, and the effects of the double layer can be accounted for at the same time. With this approach, the numeric grid size can be much larger than the thickness of double layer. Presented in this report are a description of the approach, methodology for implementation and several validation simulations for micro flows.

Chen, C. P.↗

Electromagnetic Radial Forces in a Hybrid Eight-Stator-Pole, Six-Rotor-Pole Bearingless Switched-Reluctance Motor

Analysis and experimental measurement of the electromagnet force loads on the hybrid rotor in a novel bearingless switched-reluctance motor (BSRM) have been performed. A BSRM has the combined characteristics of a switched-reluctance motor and a magnetic bearing. The BSRM has an eight-pole stator and a six-pole hybrid rotor, which is composed of circular and scalloped lamination segments. The hybrid rotor is levitated using only one set of stator poles. A second set of stator poles imparts torque to the scalloped portion of the rotor, which is driven in a traditional switched reluctance manner by a processor. Analysis was done for nonrotating rotor poles that were oriented to achieve maximum and minimum radial force loads on the rotor. The objective is to assess whether simple one-dimensional magnetic circuit analysis is sufficient for preliminary evaluation of this machine, which may exhibit strong three-dimensional electromagnetic field behavior. Two magnetic circuit geometries, approximating the complex topology of the magnetic fields in and around the hybrid rotor, were employed in formulating the electromagnetic radial force equations. Reasonable agreement between the experimental results and the theoretical predictions was obtained with typical magnetic bearing derating factors applied to the predictions.

Morrison, Carlos R.↗

Increased Fidelity in Prediction Methods For Landing Gear Noise

An aeroacoustic prediction scheme has been developed for landing gear noise. The method is designed to handle the complex landing gear geometry of current and future aircraft. The gear is represented by a collection of subassemblies and simple components that are modeled using acoustic elements. These acoustic elements are generic, but generate noise representative of the physical components on a landing gear. The method sums the noise radiation from each component of the undercarriage in isolation accounting for interference with adjacent components through an estimate of the local upstream and downstream flows and turbulence intensities. The acoustic calculations are made in the code LGMAP, which computes the sound pressure levels at various observer locations. The method can calculate the noise from the undercarriage in isolation or installed on an aircraft for both main and nose landing gear. Comparisons with wind tunnel and flight data are used to initially calibrate the method, then it may be used to predict the noise of any landing gear. In this paper, noise predictions are compared with wind tunnel data for model landing gears of various scales and levels of fidelity, as well as with flight data on fullscale undercarriages. The present agreement between the calculations and measurements suggests the method has promise for future application in the prediction of airframe noise.

Lopes, Leonard V.↗

Electromagnetic Forces in a Hybrid Magnetic-Bearing Switched-Reluctance Motor

Analysis and experimental measurement of the electromagnetic force loads on the hybrid rotor in a novel hybrid magnetic-bearing switched-reluctance motor (MBSRM) have been performed. A MBSRM has the combined characteristics of a switched-reluctance motor and a magnetic bearing. The MBSRM discussed in this report has an eight-pole stator and a six-pole hybrid rotor, which is composed of circular and scalloped lamination segments. The hybrid rotor is levitated using only one set of four stator poles, while a second set of four stator poles imparts torque to the scalloped portion of the rotor, which is driven in a traditional switched reluctance manner by a processor. Static torque and radial force analysis were done for rotor poles that were oriented to achieve maximum and minimum radial force loads on the rotor. The objective is to assess whether simple one-dimensional magnetic circuit analysis is sufficient for preliminary evaluation of this machine, which may exhibit strong three-dimensional electromagnetic field behavior. Two magnetic circuit geometries, approximating the complex topology of the magnetic fields in and around the hybrid rotor, were employed in formulating the electromagnetic radial force equations. Reasonable agreement between the experimental and the theoretical radial force loads predictions was obtained with typical magnetic bearing derating factors applied to the predictions.

Morrison, Carlos R.↗

A Domain-Decomposed Multilevel Method for Adaptively Refined Cartesian Grids with Embedded Boundaries

Preliminary verification and validation of an efficient Euler solver for adaptively refined Cartesian meshes with embedded boundaries is presented. The parallel, multilevel method makes use of a new on-the-fly parallel domain decomposition strategy based upon the use of space-filling curves, and automatically generates a sequence of coarse meshes for processing by the multigrid smoother. The coarse mesh generation algorithm produces grids which completely cover the computational domain at every level in the mesh hierarchy. A series of examples on realistically complex three-dimensional configurations demonstrate that this new coarsening algorithm reliably achieves mesh coarsening ratios in excess of 7 on adaptively refined meshes. Numerical investigations of the scheme's local truncation error demonstrate an achieved order of accuracy between 1.82 and 1.88. Convergence results for the multigrid scheme are presented for both subsonic and transonic test cases and demonstrate W-cycle multigrid convergence rates between 0.84 and 0.94. Preliminary parallel scalability tests on both simple wing and complex complete aircraft geometries shows a computational speedup of 52 on 64 processors using the run-time mesh partitioner.

Aftosmis, M. J.↗

The Multiple Use Plug Hybrid for NanoSats (MUPHyN) Miniature Thruster

The Multiple Use Plug Hybrid (for) Nanosats is a prototype thruster is being developed to fill a niche application for NanoSat-scale spacecraft propulsion. When fully developed, the MUPHyN thruster will provide an effective and low-risk propulsive capability that could enable multiple NanoSats to be independently re-positioned after deployment from a parent launch vehicle. Because the environmentally benign, chemically-stable propellants are mixed only within the combustion chamber after ignition and the flow rate of the fuel is determined by a pyrolysis mechanism that is nearly independent of pressure or fuel grain defects, the system is inherently safe and can be piggy-backed near a secondary payload with little or no overall mission risk increase to the primary payload. The MUPHyN thruster uses safe-handling and inexpensive nitrous oxide (N2O) and acrylonitrile-butadiene-styrene (ABS) as propellants. Fused Deposition Modeling (FDM), a direct digital manufacturing process, is used to fabricate short-form-factor solid fuel grains with multiple helical combustion ports from ABS thermoplastic. This manufacturing process allows for the rapid development and manufacture of complex fuel grain geometries that are not possible to extrude or cast using conventional methods. This technology enables the construction of fuel grains with length-to-diameter ratios appropriate for incorporation into CubeSats while maintaining high surface areas and regression rates that allow the system to maintain a near optimal oxidizer to fuel ratio. The MUPHyN system provides attitude control torques by using secondary-injection thrust vectoring on a truncated aerospike nozzle. This configuration allows large impulse delta V burns and small impulse attitude control firings to be performed with the same system. To ensure survivability during extend duration burns, the MUPHyN incorporates a novel regenerative cooling design where the N2O oxidizer flows through a cooling path embedded in the aerospike nozzle before being injected into the combustion chamber near the nozzle base.

Eilers, Shannon D.↗

The Evolving Polarized Jet of Black Hole Candidate Swift J1745-26

Swift J1745−26 is an X-ray binary towards the Galactic Centre that was detected when it went into outburst in September 2012. This source is thought to be one of a growing number of sources that display "failed outbursts", in which the self-absorbed radio jets of the transient source are never fully quenched and the thermal emission from the geometrically-thin inner accretion disk never fully dominates the X-ray flux. We present multifrequency data from the Very Large Array, Australia Telescope Compact Array and Karoo Array Telescope (KAT- 7) radio arrays, spanning the entire period of the outburst. Our rich data set exposes radio emission that displays a high level of large scale variability compared to the X-ray emission and deviations from the standard radio-X-ray correlation that are indicative of an unstable jet and confirm the outburst's transition from the canonical hard state to an intermediate state. We also observe steepening of the spectral index and an increase of the linear polarization to a large fraction (is approx. equal to 50%) of the total flux, as well as a rotation of the electric vector position angle. These are consistent with a transformation from a self-absorbed compact jet to optically-thin ejecta - the first time such a discrete ejection has been observed in a failed outburst - and may imply a complex magnetic field geometry.

black↗

Rapid Fabrication Techniques for Liquid Rocket Channel Wall Nozzles

The functions of a regeneratively-cooled nozzle are to (1) expand combustion gases to increase exhaust gas velocity while, (2) maintaining adequate wall temperatures to prevent structural failure, and (3) transfer heat from the hot gases to the coolant fluid to promote injector performance and stability. Regeneratively-cooled nozzles are grouped into two categories: tube-wall nozzles and channel wall nozzles. A channel wall nozzle is designed with an internal liner containing a series of integral coolant channels that are closed out with an external jacket. Manifolds are attached at each end of the nozzle to distribute coolant to and away from the channels. A variety of manufacturing techniques have been explored for channel wall nozzles, including state of the art laser-welded closeouts and pressure-assisted braze closeouts. This paper discusses techniques that NASA MSFC is evaluating for rapid fabrication of channel wall nozzles that address liner fabrication, slotting techniques and liner closeout techniques. Techniques being evaluated for liner fabrication include large-scale additive manufacturing of freeform-deposition structures to create the liner blanks. Abrasive water jet milling is being evaluated for cutting the complex coolant channel geometries. Techniques being considered for rapid closeout of the slotted liners include freeform deposition, explosive bonding and Cold Spray. Each of these techniques, development work and results are discussed in further detail in this paper.

Gradl, Paul R.↗

Computational Icing Risk Analysis of the D8 "Double Bubble" Aircraft

A computational icing risk analysis utilizing LEWICE3D was performed for the D8 Double Bubble aircraft. A variety of discrete drop sizes spanning the Appendix C and O regimes were simulated. For computational efficiency a 50-bin global discretization was produced and projected onto the distributions of interest, eliminating redundant simulations. The trajectory and impingement characteristics for discrete drop diameters were analyzed to help understand the behavior of the water drops in the presence of a complex flow field. The collection efficiency results for the discrete drop diameters were then weighted by their contributions to the total water content of six different continuous distributions and subsequently superposed to approximate these curves. Results indicate that significant variation in impingement exists as a function of drop diameter for complex wing body geometries, and that current discretization practices may be insufficient to accurately predict water collection on certain regions of the aircraft. Results also indicate that the Appendix O distributions, specifically those with considerable water content at large drops, generates water collection patterns that are markedly different from distributions representative of Appendix C.

icing↗

Computational Icing Risk Analysis of the D8 "Double Bubble" Aircraft

A computational icing risk analysis utilizing LEWICE3D was performed for the D8 Double Bubble aircraft. A variety of discrete drop sizes spanning the Appendix C and O regimes were simulated. For computational efficiency a 50-bin global discretization was produced and projected onto the distributions of interest, eliminating redundant simulations. The trajectory and impingement characteristics for discrete drop diameters were analyzed to help understand the behavior of the water drops in the presence of a complex flow field. The collection efficiency results for the discrete drop diameters were then weighted by their contributions to the total water content of six different continuous distributions and subsequently superposed to approximate these curves. Results indicate that significant variation in impingement exists as a function of drop diameter for complex wing body geometries, and that current discretization practices may be insufficient to accurately predict water collection on certain regions of the aircraft. Results also indicate that the Appendix O distributions, specifically those with considerable water content at large drops, generates water collection patterns that are markedly different from distributions representative of Appendix C.

Porter, Christopher E.↗