SEARCH · Engineering Papers

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Global Impacts of Marine Methanethiol Emissions and Chemistry in the Atmosphere

Oceanic emissions of dimethyl sulfide (DMS) have long been known to influence aerosol particle composition, cloud condensation nuclei (CCN) concentration, and Earth’s radiative budget. However, the impact of oceanic emissions of methanethiol (MeSH), a sulfur compound produced by the same oceanic precursor as DMS, has been relatively less explored. The gas-phase oxidation of MeSH has a higher effective yield of SO 2 and a shorter oxidative lifetime compared to DMS, highlighting the relevance of this pathway for the modeled representation of particle formation, growth, and CCN abundance in the marine atmosphere. Here, we use the global chemical transport model GEOS-Chem to explore possible scenarios representative of specific environmental conditions and MeSH emission schemes based on previous experimental studies. We further implement and test previously reported chemical mechanisms for MeSH oxidation, along with additional improvements, highlighting key uncertainties and sensitivities for regional and global sulfur budgets. We place our results in the context of recent modeling updates to DMS chemistry and cloud processing, which further impact SO 2 production in the marine atmosphere in parallel with MeSH oxidation. Within the overall marine sulfur budget, our findings highlight that MeSH plays a significant role in SO 2 production in the marine atmosphere, contributing to regional surface layer concentration increases of up to 40–60%. These results point to the importance of MeSH for efforts aimed at improving the modeled representation of sulfur spatiotemporal patterns relevant to air quality predictions and climate impact assessments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Fracture Capabilities in Grizzly with the extended Finite Element Method (X-FEM)

Efforts are underway to develop fracture mechanics capabilities in the Grizzly code to enable it to be used to perform deterministic fracture assessments of degraded reactor pressure vessels (RPVs). A capability was previously developed to calculate three-dimensional interaction- integrals to extract mixed-mode stress-intensity factors. This capability requires the use of a finite element mesh that conforms to the crack geometry. The eXtended Finite Element Method (X-FEM) provides a means to represent a crack geometry without explicitly fitting the finite element mesh to it. This is effected by enhancing the element kinematics to represent jump discontinuities at arbitrary locations inside of the element, as well as the incorporation of asymptotic near-tip fields to better capture crack singularities. In this work, use of only the discontinuous enrichment functions was examined to see how accurate stress intensity factors could still be calculated. This report documents the following work to enhance Grizzly’s engineering fracture capabilities by introducing arbitrary jump discontinuities for prescribed crack geometries; X-FEM Mesh Cutting in 3D: to enhance the kinematics of elements that are intersected by arbitrary crack geometries, a mesh cutting algorithm was implemented in Grizzly. The algorithm introduces new virtual nodes and creates partial elements, and then creates a new mesh connectivity; Interaction Integral Modifications: the existing code for evaluating the interaction integral in Grizzly was based on the assumption of a mesh that was fitted to the crack geometry. Modifications were made to allow for the possibility of a crack front that passes arbitrarily through the mesh; and Benchmarking for 3D Fracture: the new capabilities were benchmarked against mixed-mode three-dimensional fracture problems with known analytical solutions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Large deformation and brittle failure calculated using the dual-domain material point method

The dual domain material point (DDMP) method is explored as a candidate to be implemented in a general purpose code to perform simulations of materials with complex geometry that undergo large history-dependent deformation and failure. To test its candidacy, we study its mesh convergence, its sensitivity to mesh orientation, and its ability to handle softening and failure of a material. Simulations of large deformation and simulations of mechanical failure are performed using both DDMP and the material point method (MPM). When cell-crossing of material points is not an issue and when there are a sufficient number of material points in each computation cell, the numerical error decreases with the square of the cell size as expected for both MPM and DDMP. DDMP has reduced error compared with MPM when there are many instances of material points crossing cell boundaries due to the continuous nature of the modified gradient of the shape functions. Simulations of a specimen under tension are also performed where the background mesh is aligned and misaligned with the tension direction. MPM displays a significant mesh-dependent stress field, DDMP shows negligible mesh dependency. Despite a mesh orientation-dependent stress field from MPM, the critical tension and failure mode from both MPM and DDMP calculations have negligible mesh dependency when using a non-local failure model. If only the failure mode is important (i.e., local stresses are unimportant), MPM with a non-local failure model is a suitable method for modeling failure with small deformations. However, if local stresses are also important or if there are large deformations with many cell-crossings before failure, DDMP should be the method that is used. A needed improvement for DDMP is identified from our numerical simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2018

XMOF2D

XMOF2D is a library implementing the X-MOF (eXtended Moments of Fluid) interface reconstruction method on unstructured polygonal meshes with convex cells. Given a computational mesh with multi-material cells, for which volume fractions and centroids of contained materials are provided, piecewise linear material interfaces are constructed inside such cells using the Moment-of-Fluid (MOF) method. Additionally, a hierarchical mesh is created, where every multi-material cell of the computational mesh contains a local mesh with fully established local topology, as well as the parentage information to relate mesh entities from the local and the computational meshes.

Kikinzon, Eugene↗

DOE OSTI2017

Report on UQ Assessments to support SESAME wire-wrapped bundle experiment

This work assesses the influence of assumptions made when generating a mesh of a wire-wrapped geometry. The contact region between a wire and its adjacent pin is commonly modeled by either embedding the wire to the adjacent pin or trimming the wire so that a gap separates the wire from its adjacent pin. These models are referred to as close-gap and open-gap approaches herein and are applied to two geometries. The first geometry consists of a single pin wire-wrapped subchannel. A polyhedral mesh and a hexahedral mesh are generated. The second and third geometry are a 7-pin and a 19-pinwire-wrapped bundles meshed with polyhedral elements only. Pressure drops are obtained with the STAR-CCM+computational fluid dynamic package. Sensitivity analyses of the mesh density, the mesh type, and the turbulent models are performed. Numerical results show that the best match to the experimental data and to the Cheng-Todreas correlation is obtained with the combination of a hexahedral mesh, the shear stress transport (SST) turbulent model, and the open-gap approach. In the case of the 7-pin geometry, the best results are obtained with the open-gap approach and the SST turbulent model. The 19-pin geometry yields contradictory results to the 7-pin geometry results, and thus will require further investigations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Development of a Novel Accelerator for Neutron Transport Solution Using the Galerkin Spectral Element Methods (Final Report)

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to develop and implement advanced acceleration schemes for the DOE NEAMS neutronics code PROTEUSSN. The project team included the University of Massachusetts Lowell, The Ohio State University, University of Michigan, and Argonne National Laboratory. The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transport method (NODAL). The DSA scheme has been implemented in PROTEUS to speed up SN transport calculations. The existing DSA scheme employs a consistent finite element formulation on the same fine mesh structure as the SN solution While the use of consistent discretizations makes the DSA method effective, it does not necessarily make it efficient because the numerical solution of a discretized elliptic diffusion problem itself can be costly, particularly when the problem size becomes large. We developed a new discontinuous Galerkin (DG) discretization of the diffusion equation, called DG-DSA, which can effectively and efficiently accelerate the SN transport iterations. As compared with the previous work, the novelty of our method is that the diffusion equation is solved on a coarse-mesh grid using the DG methods, and the DG diffusion discretization incorporates local hp adaptation, i.e., local adaptation of mesh size and/or polynomial degree, based on local total cross section (or optical thickness). Therefore, the resulting number of degrees of freedom (DOF) of the DG discretization is much less than the conventional consistent DSA discretizations, and thus DG-DSA can achieve significant improvement in computational efficiency. We implemented this scheme in PROTEUS-SN. In addition, we developed a new stabilization scheme named linear prolongation CMFD (lpCMFD). A novel feature of this scheme is that the conventional flat flux ratio–based scaling approach is replaced with a linear interpolation of the scalar flux differences at the coarse-mesh cell edges between the neutron transport and CMFD calculations. Fourier convergence analysis and numerical results show that lpCMFD is unconditionally stable and more effective than other CMFD based on acceleration schemes such as pCMFD, odCMFD. We also developed a new nonlinear diffusion acceleration scheme for solving neutron transport equations. This scheme, called LR-NDA, employs a local refinement approach on the framework of CMFD by solving a local boundary value problem of the scalar flux on the coarse-mesh structure to replace the piecewise constant scalar flux obtained by CMFD. The refined flux is then used to update the scalar flux in the neutron transport source iteration. It has been demonstrated that that LR-NDA is much more effective and stable than CMFD for a wide range of optical thicknesses. LR-NDA is a local adaptive method, which means LR-NDA does not necessarily require local refinement for all the coarse-mesh cells on the problem domain, i.e., it can be used only for relatively optically thick regions where the standard CMFD scheme would encounter the convergence problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Mapping of orthogonal 2D flux coordinates for two nearby magnetic X-points to logically rectangular domains

The behavior of plasmas in configurations with magnetic null points is of potential interest for a number of systems. Such nulls, sometimes referred to as X-points, a name strictly valid only for wellseparated nulls where magnetic separatrices cross at a right angle, may occur in only one component of the magnetic field, such as the poloidal field in a tokamak. Numerical simulations are often used to study and predict turbulence and transport in such systems. For cases with strong magnetic field such that the ion and electron gyro-frequencies are much larger that particle collision frequencies, there is a strong anisotropy in the plasma transport along the magnetic field compared to across the field. Because of the strong anisotropy, it is often advantageous for simulations to use a coordinate system based on magnetic flux surfaces together with a coordinate orthogonal to the flux surfaces. This report presents the topological considerations needed to construct the mapping on an orthogonal flux-surface mesh onto a logically rectangular domain often used in simulation codes of the divertor region of tokamaks having two nearby X-points, thus extending the mapping now common single-null or well-separated double nullpoint cases with one null located at the top and bottom of the device. This extension is especially relevant to study the capability of snowflake divertors in tokamaks for reducing the peak heat flux of exhaust plasma on the target plates [1-3]. The analysis presented here follows the analytic description of the poloidal magnetic field, B pol , near two closely-spaced X-points given by Ryutov [1-3]. This analysis is applicable to experimental devices where the two X-points are generated by distributed poloidal field coils and the plasma core current. This analytic description is used to identify topological features of an orthogonal mesh. A key parameter that specifies the various topology changes is the poloidal angle between the primary X-point and the second X-point, and the changes are associated with a set of critical angles. Understanding these changes leads to the final step in the mesh construction by identifying the locations of various regions or patches in configuration space and in a corresponding logically rectangular index space. UEDGE [4] and other 2-D fluid plasma simulation codes [5-7] describe the edge plasma of toroidally symmetric tokamak configurations with one or more X-points in a divertor region. The 2-D spatial mesh is usually based on a magnetic configuration produced by an MHD equilibrium code such as EFIT and written out as one or more “eqdsk” files. As illustrated in Fig. 1, the domain of the simulated plasma consists of several flux tubes bounded by poloidal magnetic flux surfaces. The number and spacing of the flux surfaces defines the “radial” resolution of the mesh. The “radial” direction is locally orthogonal to flux surfaces; the “poloidal” direction is locally parallel to flux surfaces. The “poloidal” resolution of the mesh is defined by the number of cells in the poloidal direction between adjacent flux surfaces. For boundary plasma transport codes used for tokamak configurations such as UEDGE it is convenient, though not necessary, to assume that all flux tubes have the same number of poloidal cells so the mesh is logically rectangular in index space, but does possess internal boundaries linking different sub-rectangles in the computational domain. The following sections of this report describe indexing schemes that map tokamak magnetic configurations with one or more x-points in the divertor region(s) onto a rectangular index space. Simulation results from the UEDGE code illustrate the schemes for a few configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Nalu-Wind and OpenFAST: A high-fidelity modeling and simulation environment for wind energy. Milestone ECP-Q2-FY19

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. This milestone represents the culmination of several parallel development activities towards the goal of establishing a full-physics simulation capability for modeling wind turbines operating in turbulent atmospheric inflow conditions. The demonstration simulation performed in this milestone is the first step towards the "ground truth" simulation and includes the following components: neutral atmospheric boundary layer inflow conditions generated using a precursor simulation, a hybrid RANS/LES simulation of the wall-resolved turbine geometry, hybridization of the turbulence equations using a blending function approach to transition from the atmospheric scales to the blade boundary layer scales near the turbine, fluid-structure interaction (FSI) that accounts for the complete set of blade deformations (bending, twisting and pitch motion, yaw and tower displacements) by coupling to a comprehensive turbine dynamics code (OpenFAST). The use of overset mesh methodology for the simulations in this milestone presents a significant deviation from the previous efforts where a sliding mesh approach was employed to model the rotation of the turbine blades. The choice of overset meshes was motivated by the need to handle arbitrarily large deformations of the blade and to allow for blade pitching in the presence of a controller and the ease of mesh generation compared to the sliding mesh approach. FSI and the new timestep algorithm used in the simulations were developed in partnership with the A2e High-Fidelity Modeling project. The individual physics components were verified and validated (V%V) through extensive code-to-code comparisons and with experiments where possible. The detailed V&V efforts provide confidence in the final simulation where these physics models were combined together even though no detailed experimental data is available to perform validation of the final configuration. Taken together, this milestone successfully demonstrated the most advanced simulation to date that has been performed with Nalu-Wind.

17 WIND ENERGY↗

DOE OSTI2020

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Replace Human Intelligence with Fast and Smart Geometric Reasoning and Graph Neural Network to Accelerate Next Gen ModSim Workflows

We present an agent-guided approach to CAD geometry decomposition that automates hex/hybrid meshing with graph neural networks (GNNs) to accelerate next-generation ModSim workflows. Our end-to-end pipeline (i) reduces 3D boundary-representation (B-Rep) models to a 2D chordal axis skeleton (CAT) and then to a 1D bipartite graph of surface and curve nodes, (ii) assigns per node labels as Cubit® WebCut actions, (iii) trains a multi-action GNN under supervised learning, and (iv) predicts five surface-node and three curve-node actions on out-of-distribution test geometries. Each graph node carries geometric, topological, and meshing attributes drawn from the B-Rep “skin” and CAT “skeleton,” with two-way mappings across 3D↔2D↔1D representations to maintain traceability back to 3D CAD. The supervised learning model exhibits stable convergence of the binary cross-entropy loss and achieves 98.7% accuracy on unseen lattice models. To operationalize decision-making, we rank predicted commands by geometric significance and prototyped the agent-guided workflow through the Cubit® Meshing PowerTool GUI. As a stretch goal, we explore reinforcement learning (RL) to reduce or remove label requirements and to learn policies for action sequences that maximize total reward (e.g., size of hex-meshable regions and resulting hex mesh quality). When all-hex meshing is not feasible, the agent assists in producing hybrid meshes—prioritizing hex in critical regions and transitioning to tetrahedral elements (tets) elsewhere—maintaining fidelity while ensuring robustness. The overarching objective is to replace manual, heuristics-based decomposition with data-driven, reproducible automation, cutting meshing turnaround time by orders of magnitude. We anticipate direct impact on simulation workflows through intelligent, scalable decomposition of complex CAD models into hex-meshable subdomains.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Closing the Loop between In Situ Stress Complexity and EGS Fracture Complexity

42 ENGINEERING↗

DOE OSTI2022

Effects of Microbial-Mineral Interactions on Organic Carbon Stabilization in a Ponderosa Pine Root Zone: A Micro-Scale Approach

Soil microbial communities affect the formation of micro-scale mineral-associated organic matter (MAOM) where complex processes, including adhesion, aggregate formation, microbial mineral weathering and soil organic matter stabilization occur in a narrow zone of large biogeochemical gradients. Here we designed a field study to examine carbon stabilization mechanisms by using in-growth mesh bags containing biotite that were placed in a ponderosa pine root zone for 6 months and compared to the surrounding bulk soil. We sought to determine the composition of the microbial community in the mesh bags compared to the surrounding soils, analyze the direct interactions between microbes and biotite, and finally identify the nature of the newly formed MAOM within the mesh-bags. Our results revealed that minerals in the mesh bags were colonized by a microbial community that produced organic matter in situ. The 16S rRNA gene sequencing and ITS2 region characterization showed phylogenetic similarity between the mesh bag and bulk soil archaea/bacteria and fungi microbiomes, with significant differences in alpha- and beta-diversity and species abundances. Organic matter pools in the mesh bags, analyzed by Fourier transform ion cyclotron resonance mass spectrometry, contained protein- (peptides) and lipid-like compounds while the bulk soil OM was comprised of lignin-like and carboxyl-rich alicyclic molecules. These results support that the newly formed biotite associated organic compounds have a microbial signature in the mesh bags. High-resolution electron microscopy documented strongly adhered organic compounds to biotite surfaces, formation of microaggregates, elemental uptake at the microbe (organic matter)-mineral interface, and distortion of biotite layers. Overall, this study shows the direct and indirect involvement of soil microbial communities from the root zone of ponderosa pine in the formation of MAOM, soil organic carbon stabilization, microaggregation, and mineral weathering at micro- and nano-scales.

58 GEOSCIENCES↗

DOE OSTI2016

SANTA BARBARA CLUSTER COMPARISON TEST WITH DISPH

The Santa Barbara cluster comparison project revealed that there is a systematic difference between entropy profiles of clusters of galaxies obtained by Eulerian mesh and Lagrangian smoothed particle hydrodynamics (SPH) codes: mesh codes gave a core with a constant entropy, whereas SPH codes did not. One possible reason for this difference is that mesh codes are not Galilean invariant. Another possible reason is the problem of the SPH method, which might give too much “protection” to cold clumps because of the unphysical surface tension induced at contact discontinuities. In this paper, we apply the density-independent formulation of SPH (DISPH), which can handle contact discontinuities accurately, to simulations of a cluster of galaxies and compare the results with those with the standard SPH. We obtained the entropy core when we adopt DISPH. The size of the core is, however, significantly smaller than those obtained with mesh simulations and is comparable to those obtained with quasi-Lagrangian schemes such as “moving mesh” and “mesh free” schemes. We conclude that both the standard SPH without artificial conductivity and Eulerian mesh codes have serious problems even with such an idealized simulation, while DISPH, SPH with artificial conductivity, and quasi-Lagrangian schemes have sufficient capability to deal with it.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1991

Generation and use of unstructured grids for turbomachinery calculations

A wavefront mesh generator for two dimensional triangular meshes as well as a brief description of the solution method used with these meshes are presented. The interest is in creating meshes for solving the equations of fluid mechanics in complex turbomachinery problems, although the mesh generator and flow solver may be used for a larger variety of applications. The focus is on the flexibility and power of the mesh generation method for triangulating extremely complex geometries and in changing the geometry to create a new mesh. Two turbomachinery applications are presented which take advantage of this method: the analysis of pylon/strut and pylon/OGV interaction in the bypass of a turbofan.

Lindquist, Dana R.↗

NASA NTRS1992

AXAF VETA-I mirror encircled energy measurements and data reduction

The AXAF VETA-I mirror encircled energy was measured with a series of apertures and two flow gas proportional counters at five X-ray energies ranging from 0.28 to 2.3 keV. The proportional counter has a thin plastic window with an opaque wire mesh supporting grid. Depending on the counter position, this mesh can cause the X-ray transmission to vary as much as +/-9 percent, which directly translates into an error in the encircled energy. In order to correct this wire mesh effect, window scan measurements were made, in which the counter was scanned in both horizontal (Y) and vertical (Z) directions with the aperture fixed. Post VETA measurement of the VXDS setup were made to determine the exact geometry and position of the mesh grid. Computer models of the window mesh were developed to simulate the X-ray transmission based on this measurement. The window scan data were fitted to such mesh models and corrections were made. After this study, the mesh effect was well understood and the final results of the encircled energy were obtained with an uncertainty of less than 0.8 percent.

Zhao, Ping↗

NASA NTRS1993

Development of an unstructured solution adaptive method for the quasi-three-dimensional Euler and Navier-Stokes equations

A general solution adaptive scheme based on a remeshing technique is developed for solving the two-dimensional and quasi-three-dimensional Euler and Favre-averaged Navier-Stokes equations. The numerical scheme is formulated on an unstructured triangular mesh utilizing an edge-based pointer system which defines the edge connectivity of the mesh structure. Jameson's four-stage hybrid Runge-Kutta scheme is used to march the solution in time. The convergence rate is enhanced through the use of local time stepping and implicit residual averaging. As the solution evolves, the mesh is regenerated adaptively using flow field information. Mesh adaptation parameters are evaluated such that an estimated local numerical error is equally distributed over the whole domain. For inviscid flows, the present approach generates a complete unstructured triangular mesh using the advancing front method. For turbulent flows, the approach combines a local highly stretched structured triangular mesh in the boundary layer region with an unstructured mesh in the remaining regions to efficiently resolve the important flow features. One-equation and two-equation turbulence models are incorporated into the present unstructured approach. Results are presented for a wide range of flow problems including two-dimensional multi-element airfoils, two-dimensional cascades, and quasi-three-dimensional cascades. This approach is shown to gain flow resolution in the refined regions while achieving a great reduction in the computational effort and storage requirements since solution points are not wasted in regions where they are not required.

Jiang, Yi-Tsann↗

NASA NTRS1993

Research in unsteady aerodynamics and computational aeroelasticity at the NASA Langley Research Center

This paper presents recent results in the unsteady aerodynamics and computational aeroelasticity research programs at the NASA Langley Research Center. These programs include development of two types of computational methods: methods that use structured computational meshes and those that use unstructured meshes. Results show that an aeroelastic analysis method that uses unsteady transonic small disturbance (TSD) potential aerodynamics and structured, Cartesian meshes is capable of accurate analysis of complex aircraft configurations. The paper describes recent enhancements to the TSD method that allow analysis of vehicles with swept, flexible vertical surfaces and flexible fuselages and presents selected results that verify the accuracy of the new capabilities. Modifications to a structured-mesh Euler/Navier-Stokes method to allow aeroelastic analysis are described, and a wing flutter analysis using the resulting method is presented. Advantages of using unstructured meshes for the analysis of complex configurations are discussed. The paper presents development of unstructured-mesh Euler/Navier-Stokes methods for unsteady aerodynamics and aeroelastic analysis. Spatial and temporal adaption methods on unstructured meshes are described, and selected results are presented.

Whitlow, Woodrow, Jr.↗

NASA NTRS1993

AXAF VETA-I mirror encircled energy measurements and data reduction

The AXAF VETA-I mirror encircled energy was measured with a series of apertures and two flow gas proportional counters at five X-ray energies ranging from 0.28 to 2.3 keV. The proportional counter has a thin plastic window with an opaque wire mesh supporting grid. Depending on the counter position, this mesh can cause the X-ray transmission to vary as much as +/- 9 percent, which directly translates into an error in the encircled energy. In order to correct this wire mesh effect, window scan measurements were made, in which the counter was scanned in both horizontal (Y) and vertical (Z) directions with the aperture fixed. Post VETA measurement of the VXDS setup were made to determine the exact geometry and position of the mesh grid. Computer models of the window mesh were developed to simulate the X-ray transmission based on this measurement. The window scan data were fitted to such mesh models and corrections were made. After this study, the mesh effect was well understood and the final results of the encircled energy were obtained with an uncertainty of less than 0.8 percent.

Zhao, Ping↗