Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moving mesh”

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 91 records · Page 5

Cartesian based grid generation/adaptive mesh refinement

Grid adaptation has recently received attention in the computational fluid dynamics (CFD) community as a means to capture the salient features of a flowfield by either moving grid points of a structured or by adding cells in an unstructured manner. An approach based on a background cartesian mesh is investigated from which the geometry is 'cut' out of the mesh. Once the mesh is obtained, a solution on this coarse grid is found, that indicates which cells need to be refined. This process of refining/solving continues until the flow is grid refined in terms of a user specified global parameter (such as drag coefficient etc.). The advantages of this approach are twofold: the generation of the base grid is independent of the topology of the bodies or surfaces around/through which the flow is to be computed, and the resulting grid (in uncut regions) is highly isotropic, so that the truncation error is low. The flow solver (which, along with the grid generation is still under development) uses a completely unstructured data base, and is a finite volume, upwinding scheme. Current and future work will address generating Navier-Stokes suitable grids by using locally aligned and normal face/cell refining. The attached plot shows a simple grid about two turbine blades.

Coirier, William J.↗

Numerical simulation of low gravity draining

A boundary value problem was solved numerically for a liquid that is assumed to be inviscid and incompressible, having a motion that is irrotational and axisymmetric, and having a constant (5 degrees) solid-liquid contact angle. The avoidance of excessive mesh distortion, encountered with strictly Lagrangian or Eulerian kinematics, was achieved by introducing an auxiliary kinematic velocity field along the free surface in order to vary the trajectories used in integrating the ordinary differential equations simulating the moving boundary. The computation of the velocity potential was based upon a nonuniform triangular mesh which was automatically revised to varying depths to accommodate the motion of the free surface. These methods permitted calculation of draining induced axisymmetric slosh through the many (or fractional) finite amplitude oscillations that can occur depending upon the balance of draining, gravitational, and surface tension forces. Velocity fields, evolution of the free surface with time, and liquid residual volumes were computed for three and one half decades of Weber number and for two Bond numbers, tank fill levels, and drain radii. Comparisons with experimental data are very satisfactory.

Bizzell, G. D.↗

On the accuracy of shape sensitivity derivatives

Six methods for calculating shape sensitivity derivatives are compared for a two-material beam problem with a moving interface. It is found that as the finite-element mesh is refined, displacement sensitivity derivatives converge more slowly than the displacement themselves. Five of the methods agree fairly well, but the adjoint variational surface method provides substantially different results. However, the difference is found to reflect convergence from another direction to the same answer, rather than reduced accuracy.

Haftka, Raphael T.↗

Simulation of disk-disk encounters with co-moving polar grids

The two-grid simulation method combining advantages of both polar and Cartesian mesh-codes is described. In addition to the stellar component reacting solely to gravitational forces, the gas component is included with dissipatively colliding particles. This allows fairly realistic simulation of planar encounters where both systems contain star plus gas disks.

Salo, Heikki↗

Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes

High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel Benchmarks (NPB). In this paper, we present some interesting performance results of ow OpenMP parallel implementation on different architectures such as the SGI Origin2000, SGI Altix, and Cray MTA-2.

Feng, Huiyu↗

Finite elements: Theory and application

Recent advances in FEM techniques and applications are discussed in reviews and reports presented at the ICASE/LaRC workshop held in Hampton, VA in July 1986. Topics addressed include FEM approaches for partial differential equations, mixed FEMs, singular FEMs, FEMs for hyperbolic systems, iterative methods for elliptic finite-element equations on general meshes, mathematical aspects of FEMS for incompressible viscous flows, and gradient weighted moving finite elements in two dimensions. Consideration is given to adaptive flux-corrected FEM transport techniques for CFD, mixed and singular finite elements and the field BEM, p and h-p versions of the FEM, transient analysis methods in computational dynamics, and FEMs for integrated flow/thermal/structural analysis.

Dwoyer, D. L.↗

Investigation of Moving Belt Radiator Technology Issues

The development of an advanced spacecraft radiator technology is reported. The moving belt radiator is a thermal radiator concept with the promise of lower specific mass (per kW rejected) than that afforded by existing technologies. The results of a parametric study to estimate radiator mass for future space power systems is presented. It is shown that this technology can be scaled up to 200 MW for higher rejection temperatures. Several aspects of the design concept are discussed, including the dynamics of a large rotating belt in microgravity. The results of a computer code developed to model the belt dynamics are presented. A series of one-g experiments to investigate the dynamics of small belts is described. A comprehensive test program to investigate belt dynamics in microgravity aboard the NASA KC-135 aircraft is discussed. It was found that the desired circular shape can readily be achieved in microgravity. It is also shown that a rotating belt is stable when subjected to simulated attitude control maneuvers. Heat exchanger design is also investigated. Several sealing concepts were examined experimentally, and are discussed. Overall heat transfer coefficients to the rotating belt are presented. Material properties for various belt materials, including screen meshes, are also presented. The results presented in this report indicate that the moving belt radiator concept is technically feasible.

Teagan, W. Peter↗

Generation of spiral bevel gears with zero kinematical errors and computer aided simulation of their meshing and contact

There is proposed a method for generation of Gleason's spiral bevel gears which provides the following properties of meshing and contact: (1) the contact normal keeps its original direction within the neighborhood of the main contact point; (2) the contact ellipse moves along the gear tooth surface; and (3) the kinematical errors caused by Gleason's method of cutting are almost zero. Computer programs for the simulation of meshing and bearing contact are developed.

Litvin, F. L.↗

Research on regional numerical weather prediction

Extension of the predictive power of dynamic weather forecasting to scales below the conventional synoptic or cyclonic scales in the near future is assessed. Lower costs per computation, more powerful computers, and a 100 km mesh over the North American area (with coarser mesh extending beyond it) are noted at present. Doubling the resolution even locally (to 50 km mesh) would entail a 16-fold increase in costs (including vertical resolution and halving the time interval), and constraints on domain size and length of forecast. Boundary conditions would be provided by the surrounding 100 km mesh, and time-varying lateral boundary conditions can be considered to handle moving phenomena. More physical processes to treat, more efficient numerical techniques, and faster computers (improved software and hardware) backing up satellite and radar data could produce further improvements in forecasting in the 1980s. Boundary layer modeling, initialization techniques, and quantitative precipitation forecasting are singled out among key tasks.

Kreitzberg, C. W.↗

A Comparative Study of Four Impedance Eduction Methodologies Using Several Test Liners

A comparative study of four commonly used impedance eduction methods is presented for a range of liner structures and test conditions. Two of the methods are restricted to uniform flow while the other two accommodate both uniform and boundary layer flows. Measurements on five liner structures (a rigid-wall insert, a ceramic tubular liner, a wire mesh liner, a low porosity conventional liner, and a high porosity conventional liner) are obtained using the NASA Langley Grazing Flow Impedance Tube. The educed impedance of each liner is presented for forty-two test conditions (three Mach numbers and fourteen frequencies). In addition, the effects of moving the acoustic source from upstream to downstream and the refractive effects of the mean boundary layer on the wire mesh liner are investigated. The primary conclusions of the study are that: (1) more accurate results are obtained for the upstream source, (2) the uniform flow methods produce nearly identical impedance spectra at and below Mach 0.3 but significant scatter in the educed impedance occurs at the higher Mach number, (3) there is better agreement in educed impedance among the methods for the conventional liners than for the rigid-wall insert, ceramic, or wire mesh liner, and (4) the refractive effects of the mean boundary layer on the educed impedance of the wire mesh liner are generally small except at Mach 0.5.

Watson, Willie R.↗

New Development in NASA's Rodent Research Hardware for Conducting Long Duration Biomedical and Basic Research in Space

Animal models, particularly rodents, are the foundation of pre-clinical research to understand human diseases and evaluate new therapeutics, and play a key role in advancing biomedical discoveries both on Earth and in space. The National Research Councils Decadal survey emphasized the importance of expanding NASAs life sciences research to perform long duration, rodent experiments on the International Space Station (ISS). To accomplish this objective, flight hardware, operations, and science capabilities were developed at NASA Ames Research Center (ARC) to enhance science return for both commercial (CASIS) and government-sponsored rodent research. The Rodent Research program at NASA ARC has pioneered a new research capability on the International Space Station and has progressed toward translating research to the ISS utilizing commercial rockets, collaborating with academia and science industry, while training crewmembers to assist in performing research on orbit. Throughout phases of these missions, our practices, hardware and operations have evolved from tested to developed standards, and we are able to modify and customize our procedure and operations for mission specific requirements. The Rodent Research Habitat is capable of providing a living environment for animals on ISS according to standard animal welfare requirements. Using the cameras in the Habitat, the Rodent Research team has the ability to perform daily health checks on animals, and further analyze the collected videos for behavioral studies. A recent development of the Rodent Research hardware is inclusion of enrichment, to provide the animals the ability to rest and huddle. The Enrichment Hut is designed carefully for adult mice (up to 35 week old) within animal welfare, engineering, and operations constraints. The Hut is made out of the same stainless steel mesh as the cage interior, it has an ingress and an egress to allow animals move freely, and a hinge door to allow crewmembers remove the animals easily. The Rodent Research team has also developed Live Animal Return (LAR) capability, which will be implemented during Rodent Research-5 mission for the first time. The animals will be transported from the Habitat to a Transporter, which will return on the Dragon capsule and splashes down in the Pacific Ocean. Once SpaceX retrieves the Dragon, all powered payloads will be transferred to a SeaVan and transferred to the Long Beach pier. The NASA team then receives the transporter and delivers to a PI-designated laboratory within 120 mile radius of Long Beach. This is a significant improvement allowing researchers to examine animals within 72 hrs. of reentry or to conduct recovery experiments. Together, the hardware improvements and experience that the Rodent Research team has gained working with principal investigators and ISS crew to conduct complex experiments on orbit are expanding capabilities for long duration rodent research on the ISS to achieve both basic science and biomedical objectives.

ISS↗

Optimal expression evaluation for data parallel architectures

A data parallel machine represents an array or other composite data structure by allocating one processor (at least conceptually) per data item. A pointwise operation can be performed between two such arrays in unit time, provided their corresponding elements are allocated in the same processors. If the arrays are not aligned in this fashion, the cost of moving one or both of them is part of the cost of the operation. The choice of where to perform the operation then affects this cost. If an expression with several operands is to be evaluated, there may be many choices of where to perform the intermediate operations. An efficient algorithm is given to find the minimum-cost way to evaluate an expression, for several different data parallel architectures. This algorithm applies to any architecture in which the metric describing the cost of moving an array is robust. This encompasses most of the common data parallel communication architectures, including meshes of arbitrary dimension and hypercubes. Remarks are made on several variations of the problem, some of which are solved and some of which remain open.

Gilbert, John R.↗

A Radar Investigation of Precipitation Properties During Discrepancies Between GOES-16 GLM and LMA Observed Flash Rates in the Skyline Alabama Supercell of 22 April 2017

Conclusions: GLM (Geostationary Lightning Mapper) flash rates were 2 to 5 times lower than LMA (Lightning Mapping Array) in an Alabama supercell that was tracked using a combination of GLM flash initiation density and VIL (Vertically Integrated Liquid); Since most lightning was initiating at 8-9 kilometers (and not at low levels) according to LMA, flash height does not appear to be a primary factor in low GLM flash rates; When (LMA-GLM) flash rate differences were largest, the LMA observed flash areas were relatively small (and vice versa); Flash size may be a primary factor in low GLM flash rates due to detectability and/or flash clustering issues with small flashes within the coarse 8 kilometers by 8 kilometers resolution; High cloud liquid water droplet concentrations were inferred indirectly from riming necessary for large radar MESH (Maximum Expected Size of Hail), VIL and hail/graupel volumes. High cloud water droplet concentrations in supercells may decrease GLM detection efficiency due to optical extinction of near IR (near Infrared) emitted by lightning as it moves through cloud; Despite large flash rate differences, GLM & LMA lightning jumps during robust supercell generally agreed with each other and radar trends in HID (Radar Reflectivity and Hydrometeor Identification), MESH and VIL.However, more LMA jumps (than GLM) in developing supercell and more GLM jumps (than LMA) in weak to decaying supercell. Future work: improve GLM tracking.

GLM↗

A structured multi-block solution-adaptive mesh algorithm with mesh quality assessment

The dynamic solution adaptive grid algorithm, DSAGA3D, is extended to automatically adapt 2-D structured multi-block grids, including adaption of the block boundaries. The extension is general, requiring only input data concerning block structure, connectivity, and boundary conditions. Imbedded grid singular points are permitted, but must be prevented from moving in space. Solutions for workshop cases 1 and 2 are obtained on multi-block grids and illustrate both increased resolution of and alignment with the solution. A mesh quality assessment criteria is proposed to determine how well a given mesh resolves and aligns with the solution obtained upon it. The criteria is used to evaluate the grid quality for solutions of workshop case 6 obtained on both static and dynamically adapted grids. The results indicate that this criteria shows promise as a means of evaluating resolution.

Ingram, Clint L.↗

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.↗

Optimal expression evaluation for data parallel architectures

A data parallel machine represents an array or other composite data structure by allocating one processor per data item. A pointwise operation can be performed between two such arrays in unit time, provided their corresponding elements are allocated in the same processors. If the arrays are not aligned in this fashion, the cost of moving one or both of them is part of the cost of operation. The choice of where to perform the operation then affects this cost. If an expression with several operands is to be evaluated, there may be many choices of where to perform the intermediate operations. An efficient algorithm is given to find the minimum cost way to evaluate an expression, for several different data parallel architectures. The algorithm applies to any architecture in which the metric describing the cost of moving an array has a property called robustness. This encompasses most of the common data parallel communication architectures, including meshes of arbitrary dimension and hypercubes.

Gilbert, J. R.↗

Optimal expression evaluation for data parallel architectures

A data parallel machine represents an array or other composits data structure by allocating one processor per data item. A pointwise operation can be performed between two such arrays in unit time, provided their corresponding elements are allocated in the same processors. If the arrays are not aligned in this fashion, the cost of moving one or both of them is part of the cost of operation. The choice of where to perform the operation then affects this cost. If an expression with several operands is to be evaluated, there may be many choices of where to perform the intermediate operations. An efficient algorithm is given to find the minimum cost way to evaluate an expression, for several different data parallel architectures. The algorithm applies to any architecture in which the metric describing the cost of moving an array has a property called robustness. This encompasses most of the common data parallel communication architectures, including meshes of arbitrary dimension and hypercubes.

Gilbert, John R.↗

Finite Macro-Element Mesh Deformation in a Structured Multi-Block Navier-Stokes Code

A mesh deformation scheme is developed for a structured multi-block Navier-Stokes code consisting of two steps. The first step is a finite element solution of either user defined or automatically generated macro-elements. Macro-elements are hexagonal finite elements created from a subset of points from the full mesh. When assembled, the finite element system spans the complete flow domain. Macro-element moduli vary according to the distance to the nearest surface, resulting in extremely stiff elements near a moving surface and very pliable elements away from boundaries. Solution of the finite element system for the imposed boundary deflections generally produces smoothly varying nodal deflections. The manner in which distance to the nearest surface has been found to critically influence the quality of the element deformation. The second step is a transfinite interpolation which distributes the macro-element nodal deflections to the remaining fluid mesh points. The scheme is demonstrated for several two-dimensional applications.

Bartels, Robert E.↗