Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “adaptive timestepping”

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.

Lumped Parameter Models of the Central Nervous System for VIIP Research

INTRODUCTION: Current long-duration missions to the International Space Station and future exploration-class missions beyond low-Earth orbit, such as to Mars and asteroids, expose astronauts to increased risk of Visual Impairment and Intracranial Pressure (VIIP) syndrome [1]. It has been hypothesized that the headward shift of cerebral spinal fluid (CSF) and blood in microgravity may cause significant elevation of intracranial pressure (ICP), which in turn induces VIIP syndrome through biomechanical pathways [1, 2]. However, there is insufficient evidence to confirm this hypothesis. In this light, we are developing lumped-parameter models of fluid transport in the central nervous system (CNS) as a means to simulate the influence of microgravity on ICP. The CNS models will also be used in concert with the lumped parameter and finite element models of the eye described in the realted IWS abstracts submitted by Nelson et al., Feola et al. and Ethier et al. METHODS: We have developed a nine compartment CNS model (Figure 1) capable of both time-dependent and steady state fluid transport simulations, based on the works of Stevens et al. [3]. The breakdown of compartments within the model includes: vascular (3), CSF (2), brain (1) and extracranial (3). The boundary pressure in the Central Arteries [A] node is prescribed using an oscillating pressure function PA(t) simulating the carotid pulsatile pressure wave as developed by Linninger et al. [4]. For each time step, pressures are integrated through time using an adaptive-timestep 4th and 5th order Runga-Kutta solver. Once pressures are found, constitutive equations are used to solve for flowrates (Q) between each compartment. In addition to fluid flow between the different compartments, compliance (C) interactions between neighboring compartments are represented. We are also developing a second CNS model based on the works of Linninger et al. [4] which takes a more granular approach to represent the interactions of the intracranial and spinal compartments with the inclusion of arteries, arterioles, capillaries, venules, veins, venous sinus, and ventricles. The flow through the arteries, veins and CSF compartments are governed by continuity, momentum and distensibility balance equations. Furthermore, unlike the Stevens et al. approach, the Monro-Kellie doctrine of constant cranial volume and the bi-phasic nature of the brain parenchyma are implemented. These features appear to be more consistent with the physiologic and anatomical behavior of the CNS, and follow a modeling philosophy similar to the lumped parameter eye model that is intended to be integrated with the CNS model. However, Linninger’s approach has never been implemented to include hydrostatic gradient and microgravity simulation capabilities. Therefore, we aim at implement this modeling approach for spaceflight simulations and assess its overall applicability to VIIP research. OBJECTIVES: We will present verification and validation test results for both models, as well as head-to-head comparison to explore their strengths and limitations with respect to mathematical implementation and physiological significance for VIIP research. In doing so, we hope to provide some guidance to the HRP research community on how to appropriately leverage lumped parameter models for space biomedical research.

Vera, J.↗

Step-by-Step Simulation of Radiation Chemistry Using Green Functions for Diffusion-Influenced Reactions

Radiolytic species are formed approximately 1 ps after the passage of ionizing radiation through matter. After their formation, they diffuse and chemically react with other radiolytic species and neighboring biological molecules, leading to various oxidative damage. Therefore, the simulation of radiation chemistry is of considerable importance to understand how radiolytic species damage biological molecules [1]. The step-by-step simulation of chemical reactions is difficult, because the radiolytic species are distributed non-homogeneously in the medium. Consequently, computational approaches based on Green functions for diffusion-influenced reactions should be used [2]. Recently, Green functions for more complex type of reactions have been published [3-4]. We have developed exact random variate generators of these Green functions [5], which will allow us to use them in radiation chemistry codes. Moreover, simulating chemistry using the Green functions is which is computationally very demanding, because the probabilities of reactions between each pair of particles should be evaluated at each timestep [2]. This kind of problem is well adapted for General Purpose Graphic Processing Units (GPGPU), which can handle a large number of similar calculations simultaneously. These new developments will allow us to include more complex reactions in chemistry codes, and to improve the calculation time. This code should be of importance to link radiation track structure simulations and DNA damage models.

Plante, Ianik↗

Adaptive methods for nonlinear structural dynamics and crashworthiness analysis

The objective is to describe three research thrusts in crashworthiness analysis: adaptivity; mixed time integration, or subcycling, in which different timesteps are used for different parts of the mesh in explicit methods; and methods for contact-impact which are highly vectorizable. The techniques are being developed to improve the accuracy of calculations, ease-of-use of crashworthiness programs, and the speed of calculations. The latter is still of importance because crashworthiness calculations are often made with models of 20,000 to 50,000 elements using explicit time integration and require on the order of 20 to 100 hours on current supercomputers. The methodologies are briefly reviewed and then some example calculations employing these methods are described. The methods are also of value to other nonlinear transient computations.

Belytschko, Ted↗

Finite element methods for high speed flows

An explicit finite element based solution procedure for solving the equations of compressible viscous high speed flow is presented. The method uses domain splitting to advance the solution with different timesteps on different portions of the mesh. For steady inviscid flows, adaptive mesh refinement procedures are successfully employed to enhance the definition of discontinuities. Preliminary ideas on the application of adaptive mesh refinement to the solution of problems involving steady viscous flow are presented. Sample timings are given for the performance of the finite element code on modern supercomputers.

Loehner, R.↗

Generation and adaptation of 3-D unstructured grids for transient problems

Grid generation and adaptive refinement techniques suitable for the simulation of strongly unsteady flows past geometrically complex bodies in 3-D are described. The grids are generated using the advancing front technique. Emphasis is placed not to generate elements that are too small, as this would severely increase the cost of simulations with explicit flow solvers. The grids are adapted to an evolving flowfield using simple h-refinement. A grid change is performed every 5 to 10 timesteps, and only one level of refinement/coarsening is allowed per mesh change.

Loehner, Rainald↗

A Conservative, Scalable, Space-Time Blade Element Rotor Model for Multi-Rotor Vehicles

The development of a parallel blade-element rotor model and its implementation into an adaptive Cartesian method is described. The unsteady version of the rotor model applies a body force to all cells contained in the swept space-time volume at each timestep and special care is taken to maintain axisymmetry on the Cartesian grid. Mesh convergence of rotor thrust and torque is obtained with around 10000 cells in the disk for the steady model. Parallelization is accomplished using OpenMP and the rotor force computation is distributed across all available nodes. Simulations of an isolated XV-15 rotor in hover show good correlation with experimental data and predictions of multi-rotor thrust variation closely match previous high fidelity simulations. The final paper will also include results from the unsteady rotor model and parallel scaling tests.

Conservative↗

Adaptive explicit and implicit finite element methods for transient thermal analysis

The application of adaptive finite element methods to the solution of transient heat conduction problems in two dimensions is investigated. The computational domain is represented by an unstructured assembly of linear triangular elements and the mesh adaptation is achieved by local regeneration of the grid, using an error estimation procedure coupled to an automatic triangular mesh generator. Two alternative solution procedures are considered. In the first procedure, the solution is advanced by explicit timestepping, with domain decomposition being used to improve the computational efficiency of the method. In the second procedure, an algorithm for constructing continuous lines which pass only once through each node of the mesh is employed. The lines are used as the basis of a fully implicit method, in which the equation system is solved by line relaxation using a block tridiagonal equation solver. The numerical performance of the two procedures is compared for the analysis of a problem involving a moving heat source applied to a convectively cooled cylindrical leading edge.

Probert, E. J.↗

Unstructured grid methods for compressible flows

The implementation of the finite element method on unstructured triangular grids is described and the development of centered finite element schemes for the solution of the compressible Euler equation on general triangular and tetrahedral grids is discussed. Explicit and implicit Lax-Wendroff type methods and a method based upon the use of explicit multistep timestepping are considered. In the latter case, the convergence behavior of the method is accelerated by the incorporation of a fully unstructured multigrid procedure. The advancing front method for generating unstructured grids of triangles and tetrahedra is described and the application of adaptive mesh techniques to both steady and transient flow analysis is illustrated.

Morgan, K.↗

Virtual Petaflop Simulation: Parallel Potential Solvers and New Integrators for Gravitational Systems

"The orbit of any one planet depends on the combined motion of all the planets, not to mention the actions of all these on each other. To consider simultaneously all these causes of motion and to define these motions by exact laws allowing of convenient calculation exceeds, unless I am mistaken, the forces of the entire human intellect" -Isaac Newton 1687. Epochal surveys are throwing down the gauntlet for cosmological simulation. We describe three keys to meeting the challenge of N-body simulation: adaptive potential solvers, adaptive integrators and volume renormalization. With these techniques and a dedicated Teraflop facility, simulation can stay even with observation of the Universe. We also describe some problems in the formation and stability of planetary systems. Here, the challenge is to perform accurate integrations that retain Hamiltonian properties for 10(exp 13) timesteps.

Lake, George↗