Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Generation”

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 757 records · Page 42

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

Flexible Microsensor Array for the Monitoring and Control of Plant Growth System

Testing for plant experiments in space has begun to explore active nutrient delivery concepts in which water and nutrients are replenished on a continuous basis for long-term growth. The goal of this study is to develop a novel microsensor array to provide information on the dissolved oxygen environment in the plant root zone for the optimum control of hydroponics and solid substrate plant cultivation systems in the space environment. Miniaturized polarographic dissolved oxygen sensors have been designed and fabricated on a flexible Kapton (trademark) (polyimide) substrate. Two capabilities of the new microsensor array were explored. First, measurements of dissolved oxygen in the plant root zone in hydroponics and solid substrate culture systems were made. The microsensor array was fabricated on a flexible substrate, and then cut out into a mesh type to make a suspended array that could be placed either in a hydroponics system or in a solid substrate cultivation system to measure the oxygen environments. Second, the in situ self-diagnostic and self-calibration capability (two-point for oxygen) was adopted by dynamically controlling the microenvironment in close proximity to the microsensors. With a built-in generating electrode that surrounds the microsensor, two kinds of microenvironments (oxygen-saturated and oxygen-depleted phases) could be established by water electrolysis depending on the polarity of the generating electrode. The unique features of the new microsensor array (small size, multiple sensors, flexibility and self-diagnosis) can have exceptional benefits for the study and optimization of plant cultivation systems in both terrestrial and microgravity environments. The in situ self-diagnostic and self-calibration features of the microsensor array will also enable continuous verification of the operability during entire plant growth cycles. This concept of automated control of a novel chemical monitoring system will minimize crew time required for maintenance, as well as reduce volume, mass, and power consumption by eliminating bulky diagnosis systems including calibrant (fluid and gas) reservoir and flow system hardware.

Kim, Chang-Soo↗

Computational Non-Destructive Evaluation Improving Ultrasonic Interrogation of Complex Geometry Composite Parts

Finite element simulation was employed in modeling the ultrasound (UT) pressure pulse propagation through a coupled liquid-composite medium to reproduce experimental data. From the simulation point of view, the proposed approach is challenging when there is a large simulation domain. For example, it is shown that a sub-micron wavelength of an ultrasonic wave requires a mesh size of several microns and this in turn requires significant computational resources, as well as special care in modeling. Some of the simulation results are presented considering that such modeling should reproduce experimental data for a healthy and faulty composite structure with complex geometry. Many possible experimental setups are simulated to demonstrate the non-destructive testing technique. This setup includes the generation of pressure pulse propagating through the tested composite plate and possible scattering by discontinuities (area of different impedance) that may be present in the panel. This scattered pulse together with the baseline pressure pulse generates a signature on the probe element which can be used to locate the position of defects in the structures.

modeling↗

Generation of spiral bevel gears with conjugate tooth surfaces and tooth contact analysis

A new method for generation of spiral bevel gears is proposed. The main features of this method are as follows: (1) the gear tooth surfaces are conjugated and can transform rotation with zero transmission errors; (2) the tooth bearing contact is localized; (3) the center of the instantaneous contact ellipse moves in a plane that has a fixed orientation; (4) the contact normal performs in the process of meshing a parallel motion; (5) the motion of the contact ellipse provides improved conditions of lubrication; and (6) the gears can be manufactured by use of Gleason's equipment.

Litvin, Faydor L.↗

Generation Method Improves Spiral Bevel Gears

Conjugate tooth surfaces reduce kinematic errors. Method for generation of spiral bevel gears provides conjugated gear-tooth surfaces and improved bearing contact. Conjugated surfaces preferable because maintain constant gear ratios during tooth-contact cycles. Changing gear ratios of nonconjugated surfaces give rise to kinematical errors in transfer of motion from driving gears to driven gears. Errors major sources of noise in power transmissions. Computer program developed to simulate cutting and meshing processes for pinions and gears to minimize kinematical errors.

Litvin, Faydor L.↗

Numerical techniques for large cosmological N-body simulations

Techniques for carrying out large N-body simulations of the gravitational evolution of clustering in the fundamental cube of an infinite periodic universe are described and compared. The accuracy of the forces derived from several commonly used particle mesh schemes is examined, showing how submesh resolution can be achieved by including short-range forces between particles by direct summation techniques. The time integration of the equations of motion is discussed, and the accuracy of the codes for various choices of 'time' variable and time step is tested by considering energy conservation as well as by direct analysis of particle trajectories. Methods for generating initial particle positions and velocities corresponding to a growing mode representation of a specified power spectrum of linear density fluctuations are described. The effects of force resolution are studied and different simulation schemes are compared. An algorithm is implemented for generating initial conditions by varying the number of particles, the initial amplitude of density fluctuations, and the initial peculiar velocity field.

Efstathiou, G.↗

Raster-Based Approach to Solar Pressure Modeling

An algorithm has been developed to take advantage of the graphics processing hardware in modern computers to efficiently compute high-fidelity solar pressure forces and torques on spacecraft, taking into account the possibility of self-shading due to the articulation of spacecraft components such as solar arrays. The process is easily extended to compute other results that depend on three-dimensional attitude analysis, such as solar array power generation or free molecular flow drag. The impact of photons upon a spacecraft introduces small forces and moments. The magnitude and direction of the forces depend on the material properties of the spacecraft components being illuminated. The parts of the components being lit depends on the orientation of the craft with respect to the Sun, as well as the gimbal angles for any significant moving external parts (solar arrays, typically). Some components may shield others from the Sun. The purpose of this innovation is to enable high-fidelity computation of solar pressure and power generation effects of illuminated portions of spacecraft, taking self-shading from spacecraft attitude and movable components into account. The key idea in this innovation is to compute results dependent upon complicated geometry by using an image to break the problem into thousands or millions of sub-problems with simple geometry, and then the results from the simpler problems are combined to give high-fidelity results for the full geometry. This process is performed by constructing a 3D model of a spacecraft using an appropriate computer language (OpenGL), and running that model on a modern computer's 3D accelerated video processor. This quickly and accurately generates a view of the model (as shown on a computer screen) that takes rotation and articulation of spacecraft components into account. When this view is interpreted as the spacecraft as seen by the Sun, then only the portions of the craft visible in the view are illuminated. The view as shown on the computer screen is composed of up to millions of pixels. Each of those pixels is associated with a small illuminated area of the spacecraft. For each pixel, it is possible to compute its position, angle (surface normal) from the view direction, and the spacecraft material (and therefore, optical coefficients) associated with that area. With this information, the area associated with each pixel can be modeled as a simple flat plate for calculating solar pressure. The vector sum of these individual flat plate models is a high-fidelity approximation of the solar pressure forces and torques on the whole vehicle. In addition to using optical coefficients associated with each spacecraft material to calculate solar pressure, a power generation coefficient is added for computing solar array power generation from the sum of the illuminated areas. Similarly, other area-based calculations, such as free molecular flow drag, are also enabled. Because the model rendering is separated from other calculations, it is relatively easy to add a new model to explore a new vehicle or mission configuration. Adding a new model is performed by adding OpenGL code, but a future version might read a mesh file exported from a computer-aided design (CAD) system to enable very rapid turnaround for new designs

Wright, Theodore W. II↗

LS-DYNA User-Defined Internal Ballistic Modeling

LS-DYNA Explicit Finite Element Analysis software and CADPROG Internal Ballistic Analysis code are integrated as a single modeling tool to analyze pyro-mechanical devices. Coupling these two codes allows the ballistic calculations to be performed without modeling the propellant in an FEA mesh. At each LS-DYNA time step, CADPROG is run analytically using the kinematic data (stroke, velocity, chamber volume) fed back from the LS-DYNA simulation. The calculated pressure is then applied to a piston or any applicable surface mesh in LS-DYNA model interactively. All external loads such as friction, damping loads, assisting/resisting loads, locking/unlocking loads are modeled in LS-DYNA and coupled with pressure loads from CADPROG to create a fully defined equation of motion Both codes are written in Fortran 77. CADPROG is integrated in the LS-DYNA source code as a user-defined subroutine and recompiled to generate a custom executable.

pyrotechnics↗

Direct Simulation of Flow over A Wall-Mounted Hump

Simulating separated flows at high Reynolds numbers using Reynolds Averaged Navier-Stokes (RANS) modeled equations remains a challenge in aeronautics. The main hindrance to progress stems from the lack of extended de- tailed data pertinent to the root cause of failure of RANS models, as well as the little progress in RANS modeling innovations in the past several decades. The goal of the current effort is to generate data for separated flow at a Reynolds numbers where conventional models are challenged. We use Direct Numerical Simulations to model turbulent flow over the wall-mounted hump configuration to investigate the physics of flow separation and boundary layer recovery, as well as provide data relevant to the modeling community. A chord-based Reynolds number of Re(sub c) = 47,500 is considered with a turbulent inflow profile of Re(sub θ) = 1,400 (θ /c = 3%). We use FDL3DI, a code that solves the compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard re- cycling/rescaling method of generating turbulent boundary layers as inflow to the computational domain. Two differ- ent configurations of the upper-wall are analyzed for two sets of boundary conditions (slip and no-slip). The results are compared with the available higher Re(sub c) (= 936,000, Re(sub θ) = 7,200, θ /c = 0.77%) experiment for major flow features. The simulated lower Rec allows for DNS-like mesh resolutions, and adequately wide spans. The results from these simulations show earlier separation and delayed reattachment compared to Re(sub c) = 936,000, and significantly higher skin friction in the forebody of the hump. We also find that the upper-wall shape and boundary condition influence pressure distribution over the hump, whereas skin friction is only influenced by the boundary condition.

Wall-Mounted↗

LS-DYNA User-Defined Internal Ballistic Modeling

- LS-DYNA Explicit Finite Element Analysis software and CADPROG Internal Ballistic Analysis code are integrated as a single modeling tool to analyze pyro-mechanical devices. - Combining these two codes allows the ballistic calculations to be performed without modeling the propellant in an FEA mesh. At each LS-DYNA time step, CADPROG is analytically run using the kinematic data (stroke, velocity, chamber volume) fed back from the LS-DYNA simulation. Calculated pressure is then applied to a piston or any applicable surface mesh in LS-DYNA model interactively. -All external loads such as friction, damping loads, assisting/resisting loads, locking/unlocking loads are modeled in LS-DYNA and coupled with pressure loads from CADPROG to create a fully defined equation of motion -Both codes are written in Fortran 77. CADPROG is inserted in LS-DYNA source code as a user-defined subroutine and re-compiled to generate a custom executable.

pyrotechnics↗

Eight cm technology thruster development

A structural integrated ion thruster with 8-cm beam diameter (SIT-8) was developed for attitude control and stationkeeping of synchronous satellites. As optimized, the system demonstrates a thrust T=1.14 mlb (not corrected for beam V sub B = 1200 V (I sub sp = 2200 sec) total propellant utilization efficiency nu sub u = 59.8% (is approximately 72% without auxiliary pulse-igniter electrode), and electrical efficiency n sub E 61.9%. The thruster incorporates a wire-mesh anode and tantalum cover surfaces to control discharge chamber flake formation and employs an auxiliary pulse-igniter electrode for hollow-cathode ignition. When the SIT-8 is integrated with the compatible SIT-5 propellant tankage, the system envelope is 35 cm long by 13 cm flange bolt circle with a mass of 9.8 kg including 6.8 kg of mercury propellant. Two thrust vectoring systems which generate beam deflections in two orthogonal directions were also developed under the program and tested with the 8-cm thruster. One system vectors the beam over + or - 10 degrees by gimbaling of the entire thruster (not including tankage), while the other system vectors the beam over + or - 7 degrees by translating the accel electrode relative to the screen electrode.

Hyman, J., Jr.↗

Three dimensional spline-generated coordinate transformations for grids around wing-body configurations

A direct algebraic method was developed and applied to generate three dimensional grids around wing-body configurations. The method used is a generalized transfinite interpolation method which generates the desired coordinate transformation using geometric data only on the boundaries of the domain of interest. The geometric data that can be specified includes not only coordinates on the boundaries but also out-of-surface parametric derivatives that give a very precise control over the transformation in the vicinity of the surface. In addition to this, the method gives good control over the stretching of the mesh between different boundaries.

Eriksson, L. E.↗

Aeroacoustic Computations of a Generic Low Boom Concept in Landing Configuration: Part 1 - Aerodynamic Simulations

Computational results are presented for a 15%-scale, full-span, generic low-boom concept aircraft model simulated as installed in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel and in free-air conditions. The simulations were conducted with the lattice Boltzmann solver PowerFLOW® to capture the time-accurate characteristics of the flow. The aerodynamic behavior of the aircraft was investigated with various combinations of flap and landing gear deployments. To replicate the flow speed used during tests of the model, a Mach number of approximately 0.20 was chosen while varying the angle of attack through a limited range. Analyses were performed on local flow quantities, global forces, as well as flow distributions between configurations. Aerodynamic quantities were shown to be sensitive to mesh resolution levels, especially with flap and landing gear systems retracted. When comparing wind tunnel to free-air simulations, the aerodynamic quantities sampled were found to be in good agreement, suggesting that tunnel blockage and model support effects did not adversely impact the aerodynamic behavior of the model during the wind tunnel test. The results of this study improve our understanding of the highly complex, vortical flow generated by configurations with swept, low aspect-ratio wings at speeds encountered during landing operations.

low boom↗

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗

Heat transfer predictions for two turbine nozzle geometries at high Reynolds and Mach numbers

Predictions of turbine vane and endwall heat transfer and pressure distributions are compared with experimental measurements for two vane geometries. The differences in geometries were due to differences in the hub profile, and both geometries were derived from the design of a high rim speed turbine (HRST). The experiments were conducted in the Isentropic Light Piston Facility (ILPF) at Pyestock at a Reynolds number of 5.3 x 10(exp 6), a Mach number of 1.2, and a wall-to-gas temperature ratio of 0.66. Predictions are given for two different steady-state three-dimensional Navier-Stokes computational analyses. C-type meshes were used, and algebraic models were employed to calculate the turbulent eddy viscosity. The effects of different turbulence modeling assumptions on the predicted results are examined. Comparisons are also given between predicted and measured total pressure distributions behind the vane. The combination of realistic engine geometries and flow conditions proved to be quite demanding in terms of the convergence of the CFD solutions. An appropriate method of grid generation, which resulted in consistently converged CFD solutions, was identified.

Boyle, R. J.↗

Unstructured Grid Development for the Space Launch System Liftoff and Transition Lineloads Computational Analysis

Production of aerodynamic lineloads databases for the Space Launch System (SLS) vehicle at Liftoff and Transition (LOT) conditions has required the development of a Computational Fluid Dynamics (CFD) workflow capable of producing high-quality solutions for this unique phase of flight. Aerodynamic considerations included a wide range of flow angles (from 0°up to 90°total angle of attack), resulting leeside separation, and interaction effects between the three bodies of the integrated SLS vehicle, as well as the nearby launch tower. Computational mesh development for similar problems at the NASA Langley Research Center (such as for the Constellation/Ares launch vehicles) has primarily relied on in-house tools such as VGRID/POSTGRID, with grids designed for NASA-developed and maintained flow solvers such as USM3D and FUN3D. The workflow for such problems has evolved over the development of the various SLS configurations to incorporate new tools such as the Heldenpatch/Heldenmesh grid generator (Helden Aerospace) and CREATE-AV Kestrel (US Department of Defense) flow solver. This paper describes efforts to benchmark a grid generation approach for LOT problems using Heldenpatch/Heldenmesh and Kestrel, verified against prior best practices from VGRID/POSTGRID. Parameters studied include surface grid density, first-layer viscous cell height, and volume grid growth rate parameters. Resulting solutions are compared based on total force and moment values, sectional line loads, and surface pressures, all validated against existing wind tunnel aerodynamic databases where available for the SLS Block 1B Cargo configuration.

Space launch systems↗

Development and Validation of a Slosh-Based Ullage Collapse Model

Ullage collapse modeling is an essential tool for the success of spacecraft powered by cryogenic propellants. The degree of ullage collapse can have a major impact on the design of cryogenic propellant tanks as well as the mission’s concept of operations. Additional commodities (e.g., helium gas for pressure control) and their quantities are dependent on the amount of heat transfer and phase change of propellants expected during flight. Currently, high-fidelity computational fluid dynamics (CFD) with heat transfer and phase change is the most accurate method for predicting ullage collapse. CFD simulations of this kind, however, can become very costly for problems requiring large domains, long simulation times, or significant liquid motion. Consequently, there is a need for a reduced order model to provide fast responses which can be used in engineering design decisions. The proposed numerical model is a nodal approach, augmented by outputs from CFD analyses that are generated rapidly and reliably. These outputs include transient values for liquid-vapor interface area, interface velocity, and liquid position using a volume of fluid (VoF) method. There are several benefits of combining a nodal approach with slosh predictions from CFD simulations: leveraging of the Loci/STREAM-VoF solver that has been validated in high-g and low-g environments, rapid modeling of VoF with isothermal propellant simplifications and relatively coarse meshes, and improved model predictions from the nodal tool. The nodal approach, informed by the slosh predictions from CFD, assumes uniform conditions in the liquid and gas domains. Modifications to account for temperature stratification in the liquid and gas may be implemented as a future improvement. Thermodynamic changes within the tank are solved at user-defined timesteps. First, impacts due to mass entering or leaving the control volume are captured. Second, heat transfer impacts are accounted for by tracking the volume of the gas which has been contacted by liquid in that timestep and rebalancing the internal energy. Finally, phase change is approximated using the area and velocity of the interface. The tank conditions are re-computed by conserving mass and energy. Calibration of the model to several test flights and experiments are required to baseline the predictive power and uncertainty. The final product will enable designers and engineers to create quick-turnaround predictions for ullage collapse and inform design feasibility.

Computational Fluid Dynamics↗

New Developments in the Method of Space-Time Conservation Element and Solution Element-Applications to Two-Dimensional Time-Marching Problems

A new numerical discretization method for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods, i.e., finite difference, finite volume, finite element, and spectral methods. It is motivated by several important physical/numerical considerations and designed to avoid several key limitations of the above traditional methods. As a result of the above considerations, a set of key principles for the design of numerical schemes was put forth in a previous report. These principles were used to construct several numerical schemes that model a 1-D time-dependent convection-diffusion equation. These schemes were then extended to solve the time-dependent Euler and Navier-Stokes equations of a perfect gas. It was shown that the above schemes compared favorably with the traditional schemes in simplicity, generality, and accuracy. In this report, the 2-D versions of the above schemes, except the Navier-Stokes solver, are constructed using the same set of design principles. Their constructions are simplified greatly by the use of a nontraditional space-time mesh. Its use results in the simplest stencil possible, i.e., a tetrahedron in a 3-D space-time with a vertex at the upper time level and other three at the lower time level. Because of the similarity in their design, each of the present 2-D solvers virtually shares with its 1-D counterpart the same fundamental characteristics. Moreover, it is shown that the present Euler solver is capable of generating highly accurate solutions for a famous 2-D shock reflection problem. Specifically, both the incident and the reflected shocks can be resolved by a single data point without the presence of numerical oscillations near the discontinuity.

Chang, Sin-Chung↗