SEARCH · Engineering Papers

Results for “coarse 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 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Mesh generation/refinement using fractal concepts and iterated function systems

A novel method of mesh generation is proposed which is based on the use of fractal concepts to derive contractive, affine transformations. The transformations are constructed in such a manner that the attractors of the resulting maps are a union of the points, lines and surfaces in the domain. In particular, the mesh nodes may be generated recursively as a sequence of points which are obtained by applying the transformations to a coarse background mesh constructed from the given boundary data. A Delaunay triangulation or similar edge connection approach can then be performed on the resulting set of nodes in order to generate the mesh. Local refinement of an existing mesh can also be performed using the procedure. The method is easily extended to three dimensions, in which case the Delaunay triangulation is replaced by an analogous 3D tesselation.

Bova, S. W.↗

NASA NTRS2004

A Computational Icing Effects Study for a Three-Dimensional Wing: Comparison with Experimental Data and Investigation of Spanwise Variation

Steady-state solutions of the Reynolds-averaged Navier-Stokes (RANS) equations were computed using the Colbalt flow solver for a constant-section, rectangular wing based on an extruded two-dimensional glaze ice shape. The one equation Spalart-Allmaras turbulence model was used. The results were compared with data obtained from a recent wind tunnel test. Computed results indicate that the steady RANS solutions do not accurately capture the recirculating region downstream of the ice accretion, even after a mesh refinement. The resulting predicted reattachment is farther downstream than indicated by the experimental data. Additionally, the solutions computed on a relatively coarse baseline mesh had detailed flow characteristics that were different from those computed on the refined mesh or the experimental data. Steady RANS solutions were also computed to investigate the effects of spanwise variation in the ice shape. The spanwise variation was obtained via a bleeding function that merged the ice shape with the clean wing using a sinusoidal spanwise variation. For these configurations, the results predicted for the extruded shape provided conservative estimates for the performance degradation of the wing. Additionally, the spanwise variation in the ice shape and the resulting differences in the flow fields did not significantly change the location of the primary reattachment.

Thompson, D.↗

NASA NTRS1975

High-order numerical solutions using cubic splines

The cubic spline collocation procedure for the numerical solution of partial differential equations was reformulated so that the accuracy of the second-derivative approximation is improved and parallels that previously obtained for lower derivative terms. The final result is a numerical procedure having overall third-order accuracy for a nonuniform mesh and overall fourth-order accuracy for a uniform mesh. Application of the technique was made to the Burger's equation, to the flow around a linear corner, to the potential flow over a circular cylinder, and to boundary layer problems. The results confirmed the higher-order accuracy of the spline method and suggest that accurate solutions for more practical flow problems can be obtained with relatively coarse nonuniform meshes.

Rubin, S. G.↗

NASA NTRS1997

Load Balancing Sequences of Unstructured Adaptive Grids

Mesh adaption is a powerful tool for efficient unstructured grid computations but causes load imbalance on multiprocessor systems. To address this problem, we have developed PLUM, an automatic portable framework for performing adaptive large-scale numerical computations in a message-passing environment. This paper makes several important additions to our previous work. First, a new remapping cost model is presented and empirically validated on an SP2. Next, our load balancing strategy is applied to sequences of dynamically adapted unstructured grids. Results indicate that our framework is effective on many processors for both steady and unsteady problems with several levels of adaption. Additionally, we demonstrate that a coarse starting mesh produces high quality load balancing, at a fraction of the cost required for a fine initial mesh. Finally, we show that the data remapping overhead can be significantly reduced by applying our heuristic processor reassignment algorithm.

Biswas, Rupak↗

DOE OSTI2024

Peridynamic Models for Random Media Found by Coarse Graining

Using coarse graining, the upscaled mechanical properties of a solid with small scale heterogeneities are derived. The method maps internal forces at the small scale onto peridynamic bond forces in the coarse grained mesh. These upscaled bond forces are used to calibrate a peridynamic material model with position-dependent parameters. These parameters incorporate mesoscale variations in the statistics of the small scale system. The upscaled peridynamic model can have a much coarser discretization than the original small scale model, allowing larger scale simulations to be performed efficiently. The convergence properties of the method are investigated for representative random microstructures. In conclusion, a bond breakage criterion for the upscaled peridynamic material model is also demonstrated.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Transient simulation of biomass combustion in a circulating fluidized bed riser

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has sky-rocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. Herin this paper summarizes the modeling effort on a 50 kW th CFB combustor designed, built, and operated at CanmetENERGY in Ottawa, Canada. The numerical model employs the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture both the physics and chemistry of a CFB combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

09 BIOMASS FUELS↗

DOE OSTI2022

Frequency Transform Method for Transient Analysis of Nuclear Reactors in Monte Carlo

Not Available

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Development of fast-running eighth core nodal solution

A fast-running, eighth-core solution technique has been developed for the quick core calculations in ANC, the Westinghouse core analysis code based on the nodal expansion method, and the BEACON{sup TM} Core Monitoring System (BEACON). This technique, referred to as 'Fast Eighth' (Fast8), enables users to obtain satisfactory macroscopic core calculation results such as reactivity, axial offset and assembly power much more quickly than calculations performed in the standard full-core or quarter-core geometry. Fast8 uses three techniques: folding the radial geometry into eighth-core geometry, reducing the axial geometry into a coarse axial mesh geometry, and smearing the spacer grids into nodes on-the-fly after the first calculation in the standard geometry. These geometry adjustments in Fast8 contribute to the reduction of the number of neutronic and material nodes while maintaining the reaction rates in each node. The axial mesh condition dependency in Fast8 is investigated in a xenon transient calculation in single assembly geometry, and the Fast8 maximum axial mesh size is determined based on these results. Fast8 is then applied to the load swing calculation, and the results are compared with the reference, full-core geometry, solution. The critical boron concentration and axial offset of Fast8 show good agreement with the reference solutions, and that Fast8 is a good technique for obtaining macroscopic core calculation results in significantly shorter time. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Biomass Combustion in a Circulating Fluidized Bed Combustor

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has skyrocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing resulting in higher combustion efficiency and the ability to use lower quality fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions. SO 2 emissions are conveniently controlled through the addition of Ca-based sulphur sorbents within the combustor. This report summarizes the current modeling effort on a 50 kW th CFB combustor with a diameter of 10 cm and a height of 5 m; designed, built, and operated at CanmetENERGY in Ottawa, Canada employing the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using 9 kg of olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture the physics and chemistry of a circulating fluidized bed combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

A hybrid finite volume method and smoothed particle hydrodynamics approach for efficient and accurate blast simulations

Modeling strong shock waves in fluids remains a persistent challenge in computational physics. Essential to research efforts in industry and defense, numerous methods have been devised to improve the accuracy and efficiency of shock simulations. A novel, hybrid Finite Volume Method (FVM)-Smoothed Particle Hydrodynamics (SPH) approach is capable of further improving efficiency and retaining accuracy by exploiting the favorable characteristics of each respective method. This hybrid approach is presented for shock capturing in compressible fluids. The Python framework Pyro2 is employed to simulate a coarse FVM mesh, while the Python framework PySPH is utilized to model the fluid in regions with high gradients through SPH particles. The performance of the hybrid FVM-SPH scheme, compared to the individual FVM and SPH methods, is assessed in 1 kt and 10 kt blast simulations. Our results indicate that the hybrid approach offers higher computational efficiency than SPH while preserving its accuracy and characteristics. The hybrid approach had a relative speedup of 11.3x and 22.3x over the FVM and SPH approaches for the 1 kt simulation and a relative speedup of 14.7x and 20.9x over the FVM and SPH approaches for the 10 kt simulation. The hybrid SPH algorithm enables future compressible fluid simulations with more extensive capabilities than grid-based methods alone, presenting potential applications in modeling fluid-structure interactions and solid deformation and fracturing in blast simulations.

Myers, Conner↗

NASA NTRS1973

Convection in the earth's mantle

The computer is used to solve for thermal convection within the earth's mantle. A review of the knowledge of surface displacements and of the present understanding of the mantle and its relevant physical and chemical properties is contained in the paper. Applicable equations assume a Newtonian fluid layer heated from below and within, with gravity acting downward. The numerical method employs finite differences and was constructed with a view toward the faithful simulation of coupling mechanisms. It enables surveying the effect of a parameter using a relatively coarse computing mesh. Some of the results obtained are presented.

Turcotte, D. L.↗

NASA NTRS1994

Development of an adaptive hp-version finite element method for computational optimal control

In this research effort, the usefulness of hp-version finite elements and adaptive solution-refinement techniques in generating numerical solutions to optimal control problems has been investigated. Under NAG-939, a general FORTRAN code was developed which approximated solutions to optimal control problems with control constraints and state constraints. Within that methodology, to get high-order accuracy in solutions, the finite element mesh would have to be refined repeatedly through bisection of the entire mesh in a given phase. In the current research effort, the order of the shape functions in each element has been made a variable, giving more flexibility in error reduction and smoothing. Similarly, individual elements can each be subdivided into many pieces, depending on the local error indicator, while other parts of the mesh remain coarsely discretized. The problem remains to reduce and smooth the error while still keeping computational effort reasonable enough to calculate time histories in a short enough time for on-board applications.

Hodges, Dewey H.↗

NASA NTRS1996

Subsonic Analysis of 0.04-Scale F-16XL Models Using an Unstructured Euler Code

The subsonic flow field about an F-16XL airplane model configuration was investigated with an inviscid unstructured grid technique. The computed surface pressures were compared to wind-tunnel test results at Mach 0.148 for a range of angles of attack from 0 deg to 20 deg. To evaluate the effect of grid dependency on the solution, a grid study was performed in which fine, medium, and coarse grid meshes were generated. The off-surface vortical flow field was locally adapted and showed improved correlation to the wind-tunnel data when compared to the nonadapted flow field. Computational results are also compared to experimental five-hole pressure probe data. A detailed analysis of the off-body computed pressure contours, velocity vectors, and particle traces are presented and discussed.

Lessard, Wendy B.↗

NASA NTRS2012

Implementation of a Smeared Crack Band Model in a Micromechanics Framework

The smeared crack band theory is implemented within the generalized method of cells and high-fidelity generalized method of cells micromechanics models to capture progressive failure within the constituents of a composite material while retaining objectivity with respect to the size of the discretization elements used in the model. An repeating unit cell containing 13 randomly arranged fibers is modeled and subjected to a combination of transverse tension/compression and transverse shear loading. The implementation is verified against experimental data (where available), and an equivalent finite element model utilizing the same implementation of the crack band theory. To evaluate the performance of the crack band theory within a repeating unit cell that is more amenable to a multiscale implementation, a single fiber is modeled with generalized method of cells and high-fidelity generalized method of cells using a relatively coarse subcell mesh which is subjected to the same loading scenarios as the multiple fiber repeating unit cell. The generalized method of cells and high-fidelity generalized method of cells models are validated against a very refined finite element model.

Pineda, Evan J.↗

NASA NTRS2012

Progressive Failure of a Unidirectional Fiber-Reinforced Composite Using the Method of Cells: Discretization Objective Computational Results

Pineda, Evan J.↗

DOE OSTI2022

Linear source method of characteristics in HELIOS2

The Studsvik Scandpower code HELIOS2 is a fast lattice transport code that has gained wide use for reactor fuel cycle analysis in commercial, regulatory, and research environments. It contains two transport solvers - Collision Probabilities (CP) and Method of Characteristics (MoC) - with a general geometry ability and generalized depletion for a variety of fuel types, experimental facilities, and isotope production applications. The HELIOS2 lattice code has been extended with the development of a linear source MoC option. The linear source equations are derived within the scope of the HELIOS2 general geometry and non-cyclic ray tracing. The implementation was tested for fine and coarse spatial mesh models of the 2D C5G7 MOX computational benchmark, which provides a pure and extreme test of a transport solver. The linear source yields at least comparable accuracy with reduced computational resource requirements relative to the flat-source approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Choosing Transport Events for Initiating Splitting and Rouletting

A study was performed to determine which transport events should be used to initiate a weight window lookup to achieve the best variance reduction performance. A weight window lookup potentially triggers particle splitting (in important regions of phase space) or rouletting (in unimportant regions), thereby optimizing computational effort. Potential initiating transport events include collisions (both pre- and post-collision), geometry surface crossings, traversing a mean-free path, and streaming across a weight window boundary. Permutations of these initiating events were tested on an urban model with background radiation sources and a spent fuel cask with a neutron dose mesh tally. Generally, all methods perform better with finer weight window meshes. Tracking on weight windows performs well for coarse weight window meshes, while a combination of splitting each mean-free path, geometric surface crossing, and before collisions performs well for fine weight window meshes.

42 ENGINEERING↗

NASA NTRSDate not supplied

Cartesian Mesh Simulations for the Third AIAA Sonic Boom Prediction Workshop

Simulation results are presented for all cases from the Third AIAA Sonic Boom Prediction Workshop. An inviscid, embedded-boundary Cartesian-mesh flow solver is used in conjunction with adjoint-based mesh adaptation to compute nearfield pressure signatures. Specialized techniques are applied to maximize accuracy and minimize cost on Cartesian meshes. The Richardson-based error estimate highlights regions of the signatures most sensitive to mesh refinement. Timing results and coarse, medium, and fine mesh sizes for nearfield cases demonstrate that the parallel decomposition approach is efficient in both computational time and wall-clock. Pressure signals are propagated to the ground using an augmented Burgers’ equation solver to predict boom carpets. Ground signatures and loudness metrics are presented for a standard atmosphere as well as more realistic atmospheric profiles, which affect overall noise levels and can significantly widen the boom carpet. Mesh convergence studies show that high sampling frequencies, around 500 kHz, are required for propagation, and the sampling frequency increases at large off-track angles with longer acoustic ray paths and propagation times. The numerical methods yield accurate results for predicting low sonic boom signatures while being among the least computationally expensive of the workshop.

ARMD↗