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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Statistics of internal stress fluctuations in dislocated crystals and relevance to density-based dislocation dynamics models

A statistical analysis of internal stress fluctuations, defined as the difference between the local mean stress and stress on dislocations, is presented for deforming crystals with 3D discrete dislocation systems. Dislocation realizations are generated using dislocation dynamics simulations and the associated stress field is computed as a superposition of a regularized stress field of dislocation lines within the domain of the solution and a complementary stress field computed via a finite-element boundary value problem. The internal stress fluctuations of interest are defined by an ensemble of the difference between the stress on dislocation lines and the local mean field stress in the crystal. The latter is established in a piecewise fashion over small voxels in the crystal thus allowing the difference between the local average stress and stress on segments to be easily estimated. The results show that the Schmid stress (resolved shear stress) and Escaig stress fluctuations on various slip systems sampled over a random set of points follow a Cauchy (Lorentz) distribution at all strain levels, with the amplitude and width of the distribution being dependent on the strain. Finally, the implications of the Schmid and Escaig internal stress fluctuations are discussed from the points of view of dislocation cross-slip and the dislocation motion in continuum dislocation dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Validation of the SCEPTRE Boltzmann-CSD Solver

A new Boltzmann-CSD solver has been developed within the SCEPTRE radiation-transport code, based on the 1st-order form of the transport equation, using discontinuous finite elements in space and energy and discrete ordinates in angle. The Boltzmann-CSD solver has been validated against experimental data for electron energy deposition distributions and for electron emission spectra. Comparison of the calculated results with experimental data shows excellent agreement for many of the test configurations and reasonable agreement for other test configurations. The tests have also been modeled with the ITS Monte Carlo code, which also shows excellent to reasonable agreement with the SCEPTRE results and experimental data. The SCEPTRE Boltzmann-CSD solver relies on electron cross sections generated by the legacy CEPXS code, which currently is limited to electron-only Boltzmann-CSD cross sections. Performing full electron-photon radiation transport with the Boltzmann-CSD solver will require further development in the cross section generating code. For the energy- deposition calculations, neglecting photon transport results in at most about 5% overprediction of the energy deposition for high-energy electrons on high-Z targets, and relatively insignificant difference for the other test configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

Preliminary design of an intermittent smoke flow visualization system

A prototype intermittent smoke flow visualization system for studying the flowfield of an aircraft in flight has been constructed and demonstrated. It provides discrete pulses of dense white smoke suitable for video imaging to determine the unsteady vortex core trajectory of fluid elements in a high angle-of-attack flowfield. Two methods of pulsing the smoke were initially investigated: (1) periodically diverting the smoke between two exit ports and (2) completely blocking the smoke flow for short times. System dynamics have been modeled mathematically, data have been collected in a wind tunnel with blockage times up to 80 milliseconds, and the prototype is currently being flown on a general aviation airplane to collect three-dimensional video data. Three different plenum chamber sizes are available. Data collected so far are consistent and repeatable, though care must be taken to provide adequate contrast levels for accurate video resolution. Camera frame rates of at least 180 frames/second and wide angle lenses for the video cameras are needed to acquire meaningful vortex core velocities and accelerations for the general aviation test aircraft installation.

Ward, D. T.↗
NASA NTRS1993

Integrated analysis and design of thick composite structures for optimal passive damping characteristics

The development of novel composite mechanics for the analysis of damping in composite laminates and structures and the more significant results of this effort are summarized. Laminate mechanics based on piecewise continuous in-plane displacement fields are described that can represent both intralaminar stresses and interlaminar shear stresses and the associated effects on the stiffness and damping characteristics of a composite laminate. Among other features, the mechanics can accurately model the static and damped dynamic response of either thin or thick composite laminates, as well as, specialty laminates with embedded compliant damping layers. The discrete laminate damping theory is further incorporated into structural analysis methods. In this context, an exact semi-analytical method for the simulation of the damped dynamic response of composite plates was developed. A finite element based method and a specialty four-node plate element were also developed for the analysis of composite structures of variable shape and boundary conditions. Numerous evaluations and applications demonstrate the quality and superiority of the mechanics in predicting the damped dynamic characteristics of composite structures. Finally, additional development was focused on the development of optimal tailoring methods for the design of thick composite structures based on the developed analytical capability. Applications on composite plates illustrated the influence of composite mechanics in the optimal design of composites and the potential for significant deviations in the resultant designs when more simplified (classical) laminate theories are used.

Saravanos, D. A.↗
NASA NTRS1998

Aeroservoelastic and Structural Dynamics Research on Smart Structures Conducted at NASA Langley Research Center

An overview of smart structures research currently underway at the NASA Langley Research Center in the areas of aeroservoelasticity and structural dynamics is presented. Analytical and experimental results, plans, potential technology pay-offs, and challenges are discussed. The goal of this research is to develop the enabling technologies to actively and passively control aircraft and rotorcraft vibration and loads using smart devices. These enabling technologies and related research efforts include developing experimentally-validated finite element and aeroservoelastic modeling techniques; conducting bench experimental tests to assess feasibility and understand system trade-offs; and conducting large-scale wind- tunnel tests to demonstrate system performance. The key aeroservoelastic applications of this research include: active twist control of rotor blades using interdigitated electrode piezoelectric composites and active control of flutter, and gust and buffeting responses using discrete piezoelectric patches. In addition, NASA Langley is an active participant in the DARPA/ Air Force Research Laboratory/ NASA/ Northrop Grumman Smart Wing program which is assessing aerodynamic performance benefits using smart materials. Keywords: aeroelasticity, smart structures, piezoelectric actuators, active fiber composites, rotorcraft, buffet load alleviation, individual blade control, aeroservoelasticity, shape memory alloys, damping augmentation, piezoelectric power consumption

McGowan, Anna-Maria Rivas↗
NASA NTRS1997

Aeroservoelastic and Structural Dynamics Research on Smart Structures Conducted at NASA Langley Research Center

An overview of smart structures research currently underway at the NASA Langley Research Center in the areas of aeroservoelasticity and structural dynamics is presented. Analytical and experimental results, plans, potential technology pay-offs, and challenges are discussed. The goal of this research is to develop the enabling technologies to actively and passively control aircraft and rotorcraft vibration and loads using smart devices. These enabling technologies and related research efforts include developing experimentally-validated finite element and aeroservoelastic modeling techniques; conducting bench experimental tests to assess feasibility and understand system trade-offs; and conducting large-scale wind tunnel tests to demonstrate system performance. The key aeroservoelastic applications of this research include: active twist control of rotor blades using interdigitated electrode piezoelectric composites and active control of flutter, and gust and buffeting responses using discrete piezoelectric patches. In addition, NASA Langley is an active participant in the DARPA/Air Force Research Laboratory/NASA/Northrop Grumman Smart Wing program which is assessing aerodynamic performance benefits using smart materials.

McGowan, Anna-Maria Rivas↗
DOE OSTI2020

CFD-DEM and PR-DNS studies of low-temperature densely packed beds

Over the past few decades, granular media is gaining attention as a viable option for heat transfer fluids (HTFs). Several research efforts are studying the use of particle-based heat transfer fluids in a wide variety of applications. With this motivation, the current work focusses on analyzing the different heat transfer mechanisms in low-temperature mono-sized densely packed granular media. To study the heat transfer behavior of granular media at different scales, the current work employs a two-way coupled computational strategy. The motion of particles is solved using the Discrete Element Method (DEM) and the interstitial air is solved using a Finite-Volume (CFD) approach. The Open-Source library CFDEM Coupling® is used in the current study to join the Finite Volume PISO solver of OpenFOAM® and the DEM solver of LIGGGHTS®. Typically, particle-particle contact conduction and particle-air convection are the most popular closure models. But recent research identified a different heat transfer phenomenon in packed beds that cannot be identified by conduction or convection models. While closure models were developed to implement this on a CFD-DEM framework, they did not capture the effect of intra-particulate thermal gradients on this phenomenon. Hence the current work also employs Particle-Resolved Direct Numerical Simulations (PR-DNS) to gain valuable insights allowing for the modification of existing models. A new closure model is then proposed here and is implemented in the CFD-DEM framework. This model provides key insights into the different heat transfer mechanism of packed beds.

42 ENGINEERING↗
NASA NTRS2014

Aeroelastic Optimization Study Based on X-56A Model

A design process which incorporates the object-oriented multidisciplinary design, analysis, and optimization (MDAO) tool and the aeroelastic effects of high fidelity finite element models to characterize the design space was successfully developed and established. Two multidisciplinary design optimization studies using an object-oriented MDAO tool developed at NASA Armstrong Flight Research Center were presented. The first study demonstrates the use of aeroelastic tailoring concepts to minimize the structural weight while meeting the design requirements including strength, buckling, and flutter. A hybrid and discretization optimization approach was implemented to improve accuracy and computational efficiency of a global optimization algorithm. The second study presents a flutter mass balancing optimization study. The results provide guidance to modify the fabricated flexible wing design and move the design flutter speeds back into the flight envelope so that the original objective of X-56A flight test can be accomplished.

flutter constraints↗
DOE OSTI2025

Is tokenization needed for masked particle modeling?

In this work, we significantly enhance masked particle modeling (MPM), a self-supervised learning scheme for constructing highly expressive representations of unordered sets relevant to developing foundation models for high-energy physics. In MPM, a model is trained to recover the missing elements of a set, a learning objective that requires no labels and can be applied directly to experimental data. We achieve significant performance improvements over previous work on MPM by addressing inefficiencies in the implementation and incorporating a more powerful decoder. We compare several pre-training tasks and introduce new reconstruction methods that utilize conditional generative models without data tokenization or discretization. We show that these new methods outperform the tokenized learning objective from the original MPM on a new test bed for foundation models for jets, which includes using a wide variety of downstream tasks relevant to jet physics, such as classification, secondary vertex finding, and track identification.

conditional generative models↗
NASA NTRS1991

Numerical aspects of computing high Reynolds number flows on unstructured meshes

An edge-data structure describing a mesh edge-wise given the vertices of each edge and neighboring cell information is used developing algorithms for the Navier-Stokes equations on triangular meshes. Edge formulas for the Galerkin and finite-element discretization of gradient, divergence, Hessian, and Laplacian operators are derived. A simple edge formula is derived for the discretization of the Laplacian operator, where precise theoretical conditions for a discrete maximum principle can be ascertained. Practical issues associated with solving the Navier-Stokes equations on unstructured meshes are addressed, along with issues concerning the generation of highly stretched triangular meshes and the modeling of turbulence on unstructured meshes. A turbulence modeling strategy is proposed, and numerical results for a high-Reynolds-number flow about single- and multielement airfoils are discussed.

Barth, Timothy J.↗
NASA NTRS2026

Development and Experimental Validation of a Path-Dependent Spin Forming Finite Element Model

Spin forming is an advanced manufacturing process widely used in the aerospace and defense sectors to produce lightweight, high-strength cylindrical components with tight dimensional tolerances. This study explores the applicability of the path-dependent Mechanical Threshold Stress (MTS) constitutive model by simulating the evolution of geometry, machining forces, and plastic deformation during the spin forming of a 10-mm thick 6061-O aluminum cylinder. While numerical modeling of spin forming has advanced substantially over the past decade, systematic verification and experimental validation of material models remain limited, particularly in predicting through-thickness process evolution. The MTS model, incorporating a Voce hardening rule, is employed for its ability to represent cyclic loading, rapidly varying temperature fields, and strain rates characteristic of spin forming. Numerical convergence analysis indicates discretization uncertainties between 0.3% and 9.2% for key quantities of interest. Experimental validation demonstrates that the MTS model, when implemented with a verified mesh, accurately reproduces both elastic and plastic behavior of 6061-O aluminum, predicting peak roller loads within 11–18% of measurements, geometric tolerances within 3%, and plastic strain distributions within 10% of experimental values. Collectively, these results establish a validated computational framework for predictive spin-forming simulations with quantified confidence, providing a foundation for extension to other alloys, geometries, and forming conditions.

Spin forming↗
NASA NTRS2004

The Influence of Mesh Density on the Impact Response of a Shuttle Leading-Edge Panel Finite Element Simulation

A study was performed to examine the influence of varying mesh density on an LS-DYNA simulation of a rectangular-shaped foam projectile impacting the space shuttle leading edge Panel 6. The shuttle leading-edge panels are fabricated of reinforced carbon-carbon (RCC) material. During the study, nine cases were executed with all possible combinations of coarse, baseline, and fine meshes of the foam and panel. For each simulation, the same material properties and impact conditions were specified and only the mesh density was varied. In the baseline model, the shell elements representing the RCC panel are approximately 0.2-in. on edge, whereas the foam elements are about 0.5-in. on edge. The element nominal edge-length for the baseline panel was halved to create a fine panel (0.1-in. edge length) mesh and doubled to create a coarse panel (0.4-in. edge length) mesh. In addition, the element nominal edge-length of the baseline foam projectile was halved (0.25-in. edge length) to create a fine foam mesh and doubled (1.0- in. edge length) to create a coarse foam mesh. The initial impact velocity of the foam was 775 ft/s. The simulations were executed in LS-DYNA version 960 for 6 ms of simulation time. Contour plots of resultant panel displacement and effective stress in the foam were compared at five discrete time intervals. Also, time-history responses of internal and kinetic energy of the panel, kinetic and hourglass energy of the foam, and resultant contact force were plotted to determine the influence of mesh density. As a final comparison, the model with a fine panel and fine foam mesh was executed with slightly different material properties for the RCC. For this model, the average degraded properties of the RCC were replaced with the maximum degraded properties. Similar comparisons of panel and foam responses were made for the average and maximum degraded models.

Jackson, Karen E.↗
NASA NTRS1989

Propagation and scattering of elastic waves in fiber-reinforced plates

The objectives of this program are as follows: modelling of guided waves in fiber-reinforced plates in terms of different modes; and analysis of scattering by transverse cracks using modal representation. A hybrid numerical method combining the finite element representation of a region around the crack with the modal representation in the exterior region will be used in this program. Modes will be obtained using the through-the-thickness discretization of the displacement field.

Bratton, Robert L.↗
NASA NTRSDate not supplied

High-Lift Common Research Model: RANS, HRLES and WMLES Perspectives for CLmax Prediction Using LAVA

A unified assessment of three turbulence treatments: Reynolds Averaged Navier-Stokes(RANS), Hybrid RANS/LES (HRLES) and Equilibrium Wall-Modelled Large Eddy Simulation(WMLES) is presented for the High-Lift Common Research Model (CRM-HL). For the free-air configuration, steady-state RANS simulations show very accurate drag polar predictions in the low-𝛼linear regime. However, strong grid sensitivity is reported near the maximum lift-state(𝐶𝐿max), with finer-grids showing larger errors and predicting erroneous flow topologies on the wing. Our RANS simulations show that several corrections for the Spalart-Allmaras (SA)turbulence model widely used in the community lead to more erroneous results compared to the baseline closure, without exception. Both scale-resolving methods (HRLES and WMLES)address these drawbacks and predict an outboard separation pattern on the main element that is in good agreement with the oil flow photographs taken from the QinetiQ wind tunnel experiments, when LES-appropriate grids and numerical discretizations are used. While RANS simulations with the baseline SA closure do not show any wing-root separation post𝐶𝐿max, both HRLES and WMLES show onset of corner flow separation with varying degrees of progression, along with a weak pitch break in the wing-contribution of the overall pitching moment. This post-𝐶𝐿max pitch break seen in the free-air simulations is weaker than the break observed in experiments, with a weaker break reported in WMLES for each iteration of grid-refinement. In-tunnel simulations using both SA-baseline RANS and WMLES show a much stronger post-𝐶𝐿max break with the WMLES predictions showing excellent agreement with the experiment in terms of both the flow-topology observed and the pressure-coefficients at various spanwise stations. Sensitivity to the tunnel wall boundary layer is characterized via comparisons between viscous and inviscid treatments for the tunnel walls. WMLES predictions show moderate sensitivity at the predicted inboard flow-state at 𝐶𝐿max along with the progression towards a post-𝐶𝐿max stall; however, this stalled state at 𝛼≈20◦(inside the tunnel) obtained with both tunnel wall treatments appears to be largely identical.

TTT↗
NASA NTRSDate not supplied

FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.

CFD↗
NASA NTRS1986

Remote sensing of Earth Terrain

The Mueller matrix and polarization covariance matrix for polarimetric radar systems were studied. The clutter is modelled by a layer of random permittivity, described by a three dimensional correlation function, with variance, and horizontal and vertical correlation lengths. This model is applied, using the wave theory with Born approximations carried to the second order, to find the backscattering elements of the polarimetric matrices. It is found that 8 out of 16 elements of the Mueller matrix are identically zero, corresponding to a covariance matrix with four zero elements. Theoretical predictions are matched with experimental data for vegetation fields. A general mixing formula is derived for discrete scatterers immersed in a host medium. The inclusion particles are assumed to be ellipsoidal. The electric field inside the scatterers is determined by quasistatic analysis, assuming the diameter of the inclusion particles to be much smaller than the wavelength. The results are applicable to general multiphase mixtures.

Kong, J. A.↗
NASA NTRS1987

Finite octree meshing through topologically driven geometric operators

The octree technique is developed into the finite octree, and an overview is given. Modeler requirements are given. The octree discretization is discussed along with geometric communication operators. Geometric communication operators returning topological associativity and geometric communication operators returning spatial data are also discussed and illustrated. The advantages are given of the boundary representation and of geometric communication operators. The implementation plays an important role in the integration with a variety of geometric modelers. The capabilities of closed loop processes within a complete finite element system are presented.

Grice, Kurt R.↗
NASA NTRS1993

Stiffness, thermal expansion, and thermal bending formulation of stiffened, fiber-reinforced composite panels

A method is presented for formulating stiffness terms and thermal coefficients of stiffened, fiber-reinforced composite panels. The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. Nonlinear, temperature and load dependent constitutive material data of each laminate are used to 'build-up' the stiffened panel membrane, bending, and membrane-bending coupling stiffness terms and thermal coefficients. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN analysis program (FEA) is revealed that provides accurate solutions of entire airframes or engines with coarsely meshed models. An example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. The erroneous results and the correct ones produced from this method are compared to analysis from discretely meshed three-dimensional FEA.

Collier, Craig S.↗