Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discrete element modeling”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Equivalent-Continuum Modeling With Application to Carbon Nanotubes

A method has been proposed for developing structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with equivalent-continuum models. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As important examples with direct application to the development and characterization of single-walled carbon nanotubes and the design of nanotube-based devices, the modeling technique has been applied to determine the effective-continuum geometry and bending rigidity of a graphene sheet. A representative volume element of the chemical structure of graphene has been substituted with equivalent-truss and equivalent continuum models. As a result, an effective thickness of the continuum model has been determined. This effective thickness has been shown to be significantly larger than the interatomic spacing of graphite. The effective thickness has been shown to be significantly larger than the inter-planar spacing of graphite. The effective bending rigidity of the equivalent-continuum model of a graphene sheet was determined by equating the vibrational potential energy of the molecular model of a graphene sheet subjected to cylindrical bending with the strain energy of an equivalent continuum plate subjected to cylindrical bending.

Odegard, Gregory M.↗
NASA NTRSDate not supplied

Architectural Framework for Conceptualizing Exploration Class Medical Operations

Currently, discrete projects within NASA’s Human Research Program (HRP) take aim at providing future exploration missions with the capabilities necessary to operate within a progressively Earth-independent operative environment. Realizing an entire ecosystem that can accomplish these goals requires a clear visualization on how each project contributes to the overall objectives of Human Health and Performance (HH&P) in that context. The Exploration Medical Capability (ExMC) Element has developed an architectural framework for exploration class medical operations by utilizing a systems engineering approach using model-based systems engineering (MBSE). This model, a System of Systems Architecture (SoSA), captures a means of integrating various independent but related efforts into a comprehensive view associated with HRP’s goal of providing progressively Earth-Independent Medical Operations (EIMO). This level of visibility is realized by identifying interfaces between development efforts and existing programs that contribute to the overall mission of EIMO even if those extend past the traditional boundary of HRP. Potential interfaces that are considered important when constructing an architecture include program/project Needs, Goals, and Objectives (NGOs), capabilities, envisioned use cases/scenarios, and system/subsystem structure. This presentation will focus on the processes that were used to generate the SoSA along with a means to realize its utility amongst other MBSE products from ExMC. Evaluation of the best architecture moving forward can then be achieved that promotes cohesion between various functional areas of HH&P. Additionally, the SoSA is developed with a long-term vision for the dynamic integration of novel technological developments that may be implemented in the future. The SoSA is intended to enable Stakeholders the ability to view the broader interconnections that exist and begin to create avenues of communication between various efforts. The benefits of a SoSA provide an essential step towards addressing the unique challenges of attaining the goals of broader exploration class medical missions of the future.

Systems engineering↗
NASA NTRS2024

Architectural Framework for Conceptualizing Exploration Class Medical Operations

Currently, discrete projects within NASA’s Human Research Program (HRP) take aim at providing future exploration missions with the capabilities necessary to operate within a progressively Earth-independent operative environment. Realizing an entire ecosystem that can accomplish these goals requires a clear visualization on how each project contributes to the overall objectives of Human Health and Performance (HH&P) in that context. The Exploration Medical Capability (ExMC) Element has developed an architectural framework for exploration class medical operations by utilizing a systems engineering approach using model-based systems engineering (MBSE). This model, a System of Systems Architecture (SoSA), captures a means of integrating various independent but related efforts into a comprehensive view associated with HRP’s goal of providing progressively Earth-Independent Medical Operations (EIMO). This level of visibility is realized by identifying interfaces between development efforts and existing programs that contribute to the overall mission of EIMO even if those extend past the traditional boundary of HRP. Potential interfaces that are considered important when constructing an architecture include program/project Needs, Goals, and Objectives (NGOs), capabilities, envisioned use cases/scenarios, and system/subsystem structure. This presentation will focus on the processes that were used to generate the SoSA along with a means to realize its utility amongst other MBSE products from ExMC. Evaluation of the best architecture moving forward can then be achieved that promotes cohesion between various functional areas of HH&P. Additionally, the SoSA is developed with a long-term vision for the dynamic integration of novel technological developments that may be implemented in the future. The SoSA is intended to enable Stakeholders the ability to view the broader interconnections that exist and begin to create avenues of communication between various efforts. The benefits of a SoSA provide an essential step towards addressing the unique challenges of attaining the goals of broader exploration class medical missions of the future.

Systems engineering↗
NASA NTRS1984

Advanced numerics for multi-dimensional fluid flow calculations

In recent years, there has been a growing interest in the development and use of mathematical models for the simulation of fluid flow, heat transfer and combustion processes in engineering equipment. The equations representing the multi-dimensional transport of mass, momenta and species are numerically solved by finite-difference or finite-element techniques. However despite the multiude of differencing schemes and solution algorithms, and the advancement of computing power, the calculation of multi-dimensional flows, especially three-dimensional flows, remains a mammoth task. The following discussion is concerned with the author's recent work on the construction of accurate discretization schemes for the partial derivatives, and the efficient solution of the set of nonlinear algebraic equations resulting after discretization. The present work has been jointly supported by the Ramjet Engine Division of the Wright Patterson Air Force Base, Ohio, and the NASA Lewis Research Center.

Vanka, S. P.↗
NASA NTRSDate not supplied

A Lunar Ground Truth of Microbes that are Integral to Sustaining Biological Life Support Systems

A lunar outpost such as Gateway or a Mars transit vehicle will use a semi-closed to closed loop biological life support system (BLSS) adapted for microgravity and high levels of radiation. A moon habitat will likely exist under partial gravity conditions and sustained levels of high radiation as a semi-closed loop system able to get resources from Earth via Gateway. Additionally, the moon will act as a closed-loop testbed for Mars habitat operations and will employ various elements of in situ resource utilization (ISRU). Perhaps the longest running BLSS collective study to date is that of the European MELiSSA (Micro Ecological Life Support System Alternative) project, a circular life support system, established to gain knowledge on regenerative systems aimed at the highest degree of autonomy to produce food, water and oxygen from mission wastes. This setup has evaluated the use of discrete microbial compartments and a higher plant compartment to carry out the necessary life support functions to support crew. Such compartments utilize thermophilic anaerobes to break down human and inedible plant wastes, photoheterotrophic bacteria that can further metabolize volatile fatty acids, nitrifying bacteria that can convert ammonium to plant and microalgae available nitrates, and photoautotrophic bacteria and higher plants which will convert carbon dioxide to oxygen, purify water, and provide food for human consumption. This being the most advanced BLSS model example to date -aside from the Yuegong-1, Chinese Lunar Palace- with the most defined compartment composition, and the assumption that other BLSSs will build upon such models, it becomes relevant to study the effects of the lunar environment upon the microbes that are integral to a BLSS. In ramping up to advanced stage BLSS systems, relatively simple experiments can be conducted on Artemis missions using continuous culturing and sampling of the model microbes, subsequently subjected to DNA sequencing for mutational analysis monitoring and chemical analysis to assess the sustained ability to carry out their hallmark biochemical processes efficiently. A suite of microbes should be assessed, prioritizing those with the duel capacity to be utilized in an BLSS and the ability to biochemically facilitate ISRU goals, such as the ability to transform the biogeochemistry of Moon or Martian regolith into materials that can support crop growth or extract elements of industrial significance such as aluminum or iron.

Biological Life Support Systems↗
NASA NTRS1983

Wave envelope and infinite element schemes for fan noise radiation from turbofan inlets

Finite element models are presented for the calculation of radiated fan noise in the vicinity of turbofan inlets. The models involve the use of conventional axisymmetric finite elements in an inner region close to the inlet. The far field is represented by infinite elements or wave envelope elements. Theory and results are presented for the case with no mean flow. Comparisons of computed data with analytic solutions and measured values establish the utility of both the infinite element and wave envelope element schemes. The wave envelope scheme is shown to be effective also in the far field. Both schemes use meshes which are much sparser than would be required in conventional numerical discretizations and may be applied to real inlet configurations at modest computational cost.

Astley, R. J.↗
NASA NTRS1984

Interactive computer graphic surface modeling of three-dimensional solid domains for boundary element analysis

The establishment of the boundary element method (BEM) as a valid tool for solving problems in structural mechanics and in other fields of applied physics is discussed. The development of an integrated interactive computer graphic system for the application of the BEM to three dimensional problems in elastostatics is described. The integration of interactive computer graphic techniques and the BEM takes place at the preprocessing and postprocessing stages of the analysis process, when, respectively, the data base is generated and the results are interpreted. The interactive computer graphic modeling techniques used for generating and discretizing the boundary surfaces of a solid domain are outlined.

Perucchio, R.↗
NASA NTRS1997

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

Applications are described of high-performance computing methods to the numerical simulation of complete jet engines. The methodology focuses on the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field elements. New partitioned analysis procedures to treat this coupled three-component problem were developed. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. The NASA-sponsored ENG10 program was used for the global steady state analysis of the whole engine. This program uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 was developed as well as the capability for the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames.

Felippa, C. A.↗
NASA NTRS2007

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗
NASA NTRS1994

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

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

Chang, Sin-Chung↗
NASA NTRS1997

Ice velocity at the ice front of the Filchner-Ronne Ice Shelf, Antarctica, as observed with ERS interferometry

ERS images of the two ends of the ice front of the Filchner-Ronne ice shelf (Antarctica) were utilized interferometrically to study the deformation rate of the ice shelf in response to viscous creep. On the western bank of Berkner Island (BI), near Hemmen Ice Rise (HIR), a time series of ERS data acquired in February 1992, in both ascending and descending mode, and with a three-day time interval, were utilized to map the ice velocity in two dimensions. Finite-element ice-shelf flow simulations are compared with the ERS interferograms to interpret the ice motion in terms of the physical constraints on ice-shelf flow. The efforts to fit artificial interferograms generated with model velocity output suggest that the flow regime is strongly influenced by three processes. First, a void-creation process responsible for rifts at coastal margins tends to uncouple the ice shelf from the ice rise and neighboring coast of BI. Secondly, sea ice within the void space appears to act as a binding agent between discrete ice-shelf fragments, allowing rigid-body rotations of these fragments. Third, strain rates appear to be enhanced in a narrow zone adjacent to HIR, implying significant strain softening along the boundary. It is believed that synthetic aperture radar (SAR) interferogram/model intercomparison represents a powerful impetus towards the development of better, more physically realistic ice-shelf flow models.

Rignot, Eric↗
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRS1990

Scattering of impact wave by a crack in composite plate

The surface responses due to impact load on an infinite uniaxial graphite/epoxy plate with and without delamination cracks are investigated both in time and frequency domain by using a hybrid method combining the finite element discretization of the near-field with boundary integral representation of the field outside a contour enclosing completely the crack. This combined method leads to a set of linear unsymmetric complex matrix equations which are solved to obtain the response in the frequency domain by biconjugate gradient method. The time domain response is then obtained by using an FFT. In order to capture the time-domain characteristics accurately, high order finite elements have been used. Also, both the six node singular elements and eight node transition elements are used around the crack tips to model the crack-tip singularity. It is shown that from the numerical results for surface responses both depth and length of this crack can be identified.

Ju, T. H.↗
NASA NTRS2002

Numerical Modeling of Cavitating Venturi: A Flow Control Element of Propulsion System

In a propulsion system, the propellant flow and mixture ratio could be controlled either by variable area flow control valves or by passive flow control elements such as cavitating venturies. Cavitating venturies maintain constant propellant flowrate for fixed inlet conditions (pressure and temperature) and wide range of outlet pressures, thereby maintain constant, engine thrust and mixture ratio. The flowrate through the venturi reaches a constant value and becomes independent of outlet pressure when the pressure at throat becomes equal to vapor pressure. In order to develop a numerical model of propulsion system, it is necessary to model cavitating venturies in propellant feed systems. This paper presents a finite volume model of flow network of a cavitating venturi. The venturi was discretized into a number of control volumes and mass, momentum and energy conservation equations in each control volume are simultaneously solved to calculate one-dimensional pressure, density, and flowrate and temperature distribution. The numerical model predicts cavitations at the throat when outlet pressure was gradually reduced. Once cavitation starts, with further reduction of downstream pressure, no change in flowrate is found. The numerical predictions have been compared with test data and empirical equation based on Bernoulli's equation.

Majumdar, Alok↗
NASA NTRS2016

Introduction to COFFE: The Next-Generation HPCMP CREATE-AV CFD Solver

HPCMP CREATE-AV Conservative Field Finite Element (COFFE) is a modular, extensible, robust numerical solver for the Navier-Stokes equations that invokes modularity and extensibility from its first principles. COFFE implores a flexible, class-based hierarchy that provides a modular approach consisting of discretization, physics, parallelization, and linear algebra components. These components are developed with modern software engineering principles to ensure ease of uptake from a user's or developer's perspective. The Streamwise Upwind/Petrov-Galerkin (SU/PG) method is utilized to discretize the compressible Reynolds-Averaged Navier-Stokes (RANS) equations tightly coupled with a variety of turbulence models. The mathematics and the philosophy of the methodology that makes up COFFE are presented.

Glasby, Ryan S.↗
NASA NTRS1976

New convergence criteria for the vortex-lattice models of the leading-edge separation

The convergence criterion for the vortex-lattice technique which deals with delta wings exhibiting significant leading-edge separation was studied. It was shown that one can predict pressure distributions without irregularities which agree fairly well with experimental data (which show some irregularities of their own) by replacing the system of discrete vortex lines with a single concentrated core. This core has a circulation equal to the algebraic sum of the circulations around the discrete lines and is located at the centroid of these lines. Moreover, there is a requirement that the position and strength of the core must converge as the number of elements increases. Because the calculation of the position and strength of the core is much less involved than the calculation of the loads, this approach has the additional desirable feature of requiring less computational time.

Kandil, O. A.↗
NASA NTRSDate not supplied

Edge Based Viscous Method for Node-Centered Formulations

This paper presents a novel, efficient, conservative, edge-based method for evaluation of mean flow viscous fluxes and turbulence-model diffusion terms of the Reynolds-averaged Navier-Stokes equations on tetrahedral grids. The new method is implemented in a practical, node-centered, finite-volume computational fluid dynamics solver. The baseline finite-volume scheme that is equivalent to a second-order accurate finite-element Galerkin approximation of viscous stresses is reformulated. The order of operations to compute the cell-based Green-Gauss gradients is changed to combine the operations by edge, which leads to an equivalent formulation on tetrahedral grids, improves efficiency, and preserves the compact discretization stencil based on the nearest neighbors. The computational results presented in this paper verify the implementation of this edge-based method by comparing its accuracy and iterative convergence with those of the well verified and validated baseline formulation. Efficiency gains for residual and Jacobian evaluations result in significant reduction of time to solution. This novel edge-based formulation on tetrahedra can be seamlessly combined with the baseline formulation on cells of other types for computing solutions on mixed-element grids.

Edge Based↗