Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh data”

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,315 records · Page 73

The computation of steady nozzle flow by a time-dependent method

The steady flow in two-dimensional and axisymmetric nozzles was computed using a time-dependent method. In this method the interior mesh points were computed using the MacCormack finite-difference scheme, while a characteristic scheme was used to calculate the boundary mesh points. No explicit artificial viscosity term was included. The fluid was assumed to be a perfect gas. This method was used to compute the flow in a 45 deg - 15 deg conical, converging-diverging nozzle, a 15 deg conical, converging nozzle, and a 10 deg conical, plug nozzle. Good agreement between the numerical solution and experimental data was found. In contrast to previous time-dependent methods, the computational times were less than one minute on a CDC 6600 computer.

Cline, M. C.↗

Transonic Navier-Stokes wing solution using a zonal approach. Part 1: Solution methodology and code validation

A fast diagonalized Beam-Warming algorithm is coupled with a zonal approach to solve the three-dimensional Euler/Navier-Stokes equations. The computer code, called Transonic Navier-Stokes (TNS), uses a total of four zones for wing configurations (or can be extended to complete aircraft configurations by adding zones). In the inner blocks near the wing surface, the thin-layer Navier-Stokes equations are solved, while in the outer two blocks the Euler equations are solved. The diagonal algorithm yields a speedup of as much as a factor of 40 over the original algorithm/zonal method code. The TNS code, in addition, has the capability to model wind tunnel walls. Transonic viscous solutions are obtained on a 150,000-point mesh for a NACA 0012 wing. A three-order-of-magnitude drop in the L2-norm of the residual requires approximately 500 iterations, which takes about 45 min of CPU time on a Cray-XMP processor. Simulations are also conducted for a different geometrical wing called WING C. All cases show good agreement with experimental data.

Flores, J.↗

Transonic Navier-Stokes wing solution using a zonal approach. Part 1: Solution methodology and code validation

A fast diagonalized Beam-Warming algorithm is coupled with a zonal approach to solve the three dimensional Euler/Navier-Stokes equations. The computer code, called Transonic Navier-Stokes (TNS), uses a total of four zones for wing configurations (or can be extended to complete aircraft configurations by adding zones). In the inner blocks near the wing surface, the thin-layer Navier-Stokes equations are solved, while in the outer two blocks the Euler equations are solved. The diagonal algorithm yields a speedup of as much as a factor of 40 over the original algorithm/zonal method code. The TNS code, in addition, has the capability to model wind tunnel walls. Transonic viscous solutions are obtained on a 150,000-point mesh for a NACA 0012 wing. A three-order-of-magnitude drop in the L2-norm of the residual requires approximately 500 iterations, which takes about 45 min of CPU time on a Cray-XMP processor. Simulations are also conducted for a different geometrical wing called WING C. All cases show good agreement with experimental data.

Flores, J.↗

Wing flutter boundary prediction using an unsteady Euler aerodynamic method

Modifications to an existing three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1) for the aeroelastic analysis of wings are described. These modifications, which were previously added to CFL3D Version 1.0, include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time-integration with the government flow equations. The paper gives a brief description of these modifications and presents unsteady calculations which check the modifications to the code. Euler flutter results for an isolated 45 degree swept-back wing are compared with experimental data for seven freestream Mach numbers which define the flutter boundary over a range of Mach number from 0.499 to 1.14. These comparisons show good agreement in flutter characteristics for freestream Mach numbers below unity. For freestream Mach numbers above unity, the computed aeroelastic results predict a premature rise in the flutter boundary as compared with the experimental boundary. Steady and unsteady contours of surface Mach number and pressure are included to illustrate the basic flow characteristics of the time-marching flutter calculations and to aid in identifying possible causes for the premature rise in the computational flutter boundary.

Lee-Rausch, Elizabeth M.↗

Wing flutter boundary prediction using unsteady Euler aerodynamic method

Modifications to an existing 3D implicit upwind Euler/Navier-Stokes code for the aeroelastic analysis of wings are described. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time-integration with the governing flow equations. The paper gives a brief description of these modifications and presents unsteady calculations which check the modifications to the code. Euler flutter results for an isolated 45 deg swept-back wing are compared with experimental data for seven freestream Mach numbers which define the flutter boundary over a range of Mach number from 0.499 to 1.14. These comparisons show good agreement in flutter characteristics for freestream Mach numbers below unity. For freestream Mach numbers above unity, the computed aeroelastic results predict a premature rise in the flutter boundary as compared with the experimental boundary. Steady and unsteady contours of surface Mach number and pressure are included to illustrate the basic flow characteristics of the time-marching flutter calculations and to aid in identifying possible causes for the premature rise in the computational flutter boundary.

Lee-Rausch, Elizabeth M.↗

Remote Sensing-Driven Hydrodynamic Modeling in Data-Scarce Regions: Integrating ICESat-2, Sentinel-2, SWOT and Re-analysis Models for Coastal Monitoring

Hydrodynamic models in coastal and estuarine systems are typically constrained by sparse bathymetry, boundary, and validation data, especially in regions where field campaigns are costly or impractical. Here we develop and test a fully satellite-driven framework for hydrodynamic modeling in South Africa’s Langebaan Lagoon without using any local in situ measurements. Bathymetry is derived by training multispectral Sentinel-2 reflectance against ICESat-2 ATL24 photon-derived depths using an XGBoost model optimized with Bayesian search. The final satellite-derived bathymetry reproduces independent ATL24 points with RMSE = 0.45 m and R 2 = 0.97. This bathymetry was used in a depth-averaged Delft3D Flexible Mesh model driven at the open boundary by TPXO tidal harmonics and by ERA5 winds. We validate modeled water surface elevation against 16 SWOT low-rate (250 m, unsmoothed) passes in 2023. SWOT–model comparisons yield an overall RMSE of 0.11 m and R 2 = 0.61, with typical point differences <0.10 m (∼7% of the 1.5 m tidal range), and showed consistent spatial gradients in water level from the offshore boundary, through Saldanha Bay, and into the lagoon. At the offshore boundary, TPXO and SWOT sea surface heights agree closely (R 2 = 0.86). A simple phase adjustment of ∼26,min between TPXO and SWOT lowers the RMSE from 0.18,m to 0.11,m, showing that phase offset accounts for some of the discrepancy, with additional errors likely linked to non-tidal signals. Our results demonstrate that combining passive optical, photon-counting LiDAR, radar interferometry, and global tidal/atmospheric models enables robust, transferrable hydrodynamic modeling in data-scarce coastal systems, offering a cost-effective pathway for monitoring.

ICESat-2↗

Enhancements to BISON U-Zr Metallic Fuel X447 Example Problem

As development of a metallic fuel modeling capability in BISON has progressed, the need for an example problem used as a comparison basis was observed. Collaborative work between researchers at Los Alamos National Laboratory (LANL) and Idaho National Laboratory (INL) then proceeded to determine a viable rod to use as the basis and create a BISON input deck utilizing as many metallic fuel models as feasible. The basis chosen was what would be considered a generic rod from subassembly X447, an assembly irradiated in EBR-II towards the end of its operating life, heavily based on reported data for fuel pin DP11. Thus, the approach was adopted to use flow characteristics from subassembly X447 as a basis for the convective heat transfer solution, power history and axial power profiles that are representative of rod DP11 from subassembly X447. The rod simulated is a U-10Zr wt% (U-22.5Zr at%) composition. A 2D-RZ mesh would be used to capture axial thermal hydraulic effects, axial swelling and stress-strain calculations over the full length of the rod. After initial work was invested, a refinement of the various models and input parameters was conducted to ensure consistency between operator-declared conditions, model input requirements and those represented in the example problem. This report serves as a synopsis of the enhancements and refinements to the example problem conducted throughout the 2016 fiscal year.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗

Magnetic-field measurements for the Lewis Research Center cyclotron

The magnetic field of the Lewis Center cyclotron was mapped by using a Hall-effect magnetic-field transducer. Main-field Fourier coefficients were determined on a polar mesh of 40 radii for each of seven levels of main-field coil current. Incremental fields for eight sets of trim coils and two sets of harmonic coils were also determined at four of these main-field levels. A stored-program, digital computer was used to perform the measurements. The process was entirely automatic; all data-taking and data-reduction activities were specified by the computer programs. A new method for temperature compensation of a Hall element was used. This method required no temperature control of the element. Measurements of the Hall voltage and Hall-element resistance were sufficient to correct for temperature effects.

Fessler, T. E.↗

Large eddy simulation of homogeneous isotropic turbulence

This paper presents results from a comprehensive investigation of large-eddy simulations of homogeneous isotropic turbulence. Calculations have been made using grid meshes, a number of flowfield filters, and two simple subgrid scale turbulence models. Particular attention has been paid to the degree of isotropy and the ability of the approach to predict higher-order statistical quantities (skewness and flatness). Directions for needed improvement in the simulation approach are indicated, and a plea is made to experimenters to process their data in a way which will facilitate comparison with large eddy simulations.

Ferziger, J. H.↗

Free-edge delamination: Laminate width and loading conditions effects

The width and loading conditions effects on free-edge stress fields in composite laminates are investigated using a three-dimensional finite element analysis. This analysis includes a special free-edge region refinement or superelement with progrssive substructuring (mesh refinement) and finite thickness interply layers. The different loading conditions include in-plane and out-of-plane bending, combined axial tension and in-plane shear, twisting, uniform temperature and uniform moisture. Results obtained indicate that: axial tension causes the smallest magnitude of interlaminar free edge stress compared to other loading conditions; free-edge delamination data obtained from laboratory specimens cannot be scaled to structural components; and composite structural components are not likely to delaminate.

Murthy, P. L. N.↗

Free-edge delamination - Laminate width and loading conditions effects

The width and loading conditions effects on free-edge stress fields in composite laminates are investigated using a three-dimensional finite element analysis. This analysis includes a special free-edge region refinement or superelement with progressive substructuring (mesh refinement) and finite thickness interply layers. The different loading conditions include in-plane and out-of-plane bending, combined axial tension and in-plane shear, twisting, uniform temperature and uniform moisture. Results obtained indicate that: axial tension causes the smallest magnitude of interlaminar free edge stress compared to other loading conditions; free-edge delamination data obtained from laboratory specimens cannot be scaled to structural components; and composite structural components are not likely to delaminate.

Murthy, Pappu L. N.↗

Personal communications via ACTS satellite HBR transponders

The concept of a fully meshed network of briefcase-sized terminals is presented for personal communications over Ka-band satellite transponders. In this concept, undesirable double-hop delays are avoided for voice communications. The bandwidth and power resources of the transponder are efficiently shared by users in a simple demand-assigned manner via code-division multiple access (CDMA). Voice, data, and facsimile are statistically multiplexed at each terminal. In order to minimize terminal costs, frequency-precorrected, and level-preadjusted continuous-wave tones are sent from the central network control station in each beam so that the terminals in each down-link beam can use these pilots as references for antenna acquisition and tracking, as reliable frequency sources, and as indicators of signal fade for up-link power control (ULPC). The potential CDMA 'near-far' problem due to up-link fades is mitigated by using ULPC. Quasi-burst mode transmission is employed to minimize the potential clock and pseudorandom number code synchronization.

Fang, Russell J. F.↗

Measurements of the Influence of Integral Length Scale on Stagnation Region Heat Transfer

The purpose was twofold: first, to determine if a length scale existed that would cause the greatest augmentation in stagnation region heat transfer for a given turbulence intensity and second, to develop a prediction tool for stagnation heat transfer in the presence of free stream turbulence. Toward this end, a model with a circular leading edge was fabricated with heat transfer gages in the stagnation region. The model was qualified in a low turbulence wind tunnel by comparing measurements with Frossling's solution for stagnation region heat transfer in a laminar free stream. Five turbulence generating grids were fabricated; four were square mesh, biplane grids made from square bars. Each had identical mesh to bar width ratio but different bar widths. The fifth grid was an array of fine parallel wires that were perpendicular to the axis of the cylindrical leading edge. Turbulence intensity and integral length scale were measured as a function of distance from the grids. Stagnation region heat transfer was measured at various distances downstream of each grid. Data were taken at cylinder Reynolds numbers ranging from 42,000 to 193,000. Turbulence intensities were in the range 1.1 to 15.9 percent while the ratio of integral length scale to cylinder diameter ranged from 0.05 to 0.30. Stagnation region heat transfer augmentation increased with decreasing length scale. An optimum scale was not found. A correlation was developed that fit heat transfer data for the square bar grids to within +4 percent. The data from the array of wires were not predicted by the correlation; augmentation was higher for this case indicating that the degree of isotropy in the turbulent flow field has a large effect on stagnation heat transfer. The data of other researchers are also compared with the correlation.

Vanfossen, G. James↗

The NASA Soil Moisture Active Passive (SMAP) Mission: Overview

The Soil Moisture Active Passive (SMAP) mission is one of the first Earth observation satellites being developed by NASA in response to the National Research Council?s Decadal Survey [1]. Its mission design consists of L-band radiometer and radar instruments sharing a rotating 6-m mesh reflector antenna to provide high-resolution and high-accuracy global maps of soil moisture and freeze/thaw state every 2-3 days. The combined active/passive microwave soil moisture product will have a spatial resolution of 10 km and a mean latency of 24 hours. In addition, the SMAP surface observations will be combined with advanced modeling and data assimilation to provide deeper root zone soil moisture and net ecosystem exchange of carbon. SMAP is expected to launch in the late 2014 - early 2015 time frame.

O'Neill, Peggy↗

CFD Based Design of a Filming Injector for N+3 Combustors

An effort was undertaken to perform CFD analysis of fluid flow in Lean-Direct Injection (LDI) combustors with axial swirl-venturi elements for next-generation LDI-3 combustor design. The National Combustion Code (NCC) was used to perform non-reacting and two-phase reacting flow computations for a newly-designed pre-filming type fuel injector LDI-3 injector, in a single-injector and a five-injector array configuration. All computations were performed with a consistent approach of mesh-optimization, spray-modeling, ignition and kinetics-modeling. Computational predictions of the aerodynamics of the single-injector were used to arrive at an optimized main-injector design that meets effective area and fuel-air mixing criteria. Emissions (EINOx) characteristics were predicted for a medium-power engine cycle condition, and will be compared with data when it is made available from experimental measurements. The use of a PDF-like turbulence-chemistry interaction model with NCC's Time-Filtered Navier-Stokes (TFNS) solver is shown to produce a significant impact on the CFD results, when compared with a laminar-chemistry TFNS approach for the five-injector computations.

combustion chambers↗

Preliminary Scale-Resolving Simulations of Laminar-To-Turbulent Transition in Swept-Wing Flow

Swept wings and control surfaces are common elements of modern aircraft, and it has been shown both experimentally and theoretically that laminar-to-turbulent transition of the three-dimensional boundary layer that develops over them is highly sensitive to surface roughness. Numerous studies have been conducted on the effect of discrete roughness elements or distributed roughness elements on swept flow transition, however so far limited computational effort has been dedicated to the study of transition over swept wings with randomly distributed micron-sized roughness. In the present work, we set up to reproduce the extensive experimental data base generated by Dagenhart et al for the infinite swept wing NLF(2)-0415. To this purpose, we perform scale-resolving simulations of flow transition over smooth and rough surfaces using a high-order space-time spectral-element Discontinuous-Galerkin solver. Different types of surface roughnesses are implemented by elastically deforming the original mesh. The study shows that the experimental results cannot be accounted for by a perfectly smooth wing and reveals a strong sensitivity of the transition process to the representation of the surface roughness. The crossflow patterns and transition location approach those measured for some of the surface profiles, however a correlation between the wavenumber spectrum of the surface, grid resolution and boundary layer stability is yet to be established.

Denison, Marie↗

High-fidelity Multidisciplinary Sensitivity Analysis and Design Optimization for Rotorcraft Applications

A multidisciplinary sensitivity analysis of rotorcraft simulations involving tightly coupled high-fidelity computational fluid dynamics and comprehensive analysis solvers is presented and evaluated. A sensitivity-enabled fluid dynamics solver and a nonlinear flexible multibody dynamics solver are coupled to predict aerodynamic loads and structural responses of helicopter rotor blades. A discretely consistent adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows and unstructured dynamic overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system is validated by conducting simulations for a benchmark rotorcraft model and comparing solutions with established analyses and experimental data. Sensitivities of lift computed by the multidisciplinary sensitivity analysis are verified by comparison with the sensitivities obtained by complex-variable simulations. Finally the multidisciplinary sensitivity analysis is applied to a constrained gradient-based design optimization for a HART-II rotorcraft configuration.

Wang, Li↗