Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple time stepping”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Low-Thrust Trajectory Optimization with Simplified SQP Algorithm

The problem of low-thrust trajectory optimization in highly perturbed dynamics is a stressing case for many optimization tools. Highly nonlinear dynamics and continuous thrust are each, separately, non-trivial problems in the field of optimal control, and when combined, the problem is even more difficult. This paper de-scribes a fast, robust method to design a trajectory in the CRTBP (circular restricted three body problem), beginning with no or very little knowledge of the system. The approach is inspired by the SQP (sequential quadratic programming) algorithm, in which a general nonlinear programming problem is solved via a sequence of quadratic problems. A few key simplifications make the algorithm presented fast and robust to initial guess: a quadratic cost function, neglecting the line search step when the solution is known to be far away, judicious use of end-point constraints, and mesh refinement on multiple shooting with fixed-step integration.In comparison to the traditional approach of plugging the problem into a “black-box” NLP solver, the methods shown converge even when given no knowledge of the solution at all. It was found that the only piece of information that the user needs to provide is a rough guess for the time of flight, as the transfer time guess will dictate which set of local solutions the algorithm could converge on. This robustness to initial guess is a compelling feature, as three-body orbit transfers are challenging to design with intuition alone. Of course, if a high-quality initial guess is available, the methods shown are still valid.We have shown that endpoints can be efficiently constrained to lie on 3-body repeating orbits, and that time of flight can be optimized as well. When optimizing the endpoints, we must make a trade between converging quickly on sub-optimal endpoints or converging more slowly on end-points that are arbitrarily close to optimal. It is easy for the mission design engineer to adjust this trade based on the problem at hand.The biggest limitation to the algorithm at this point is that multi-revolution transfers (greater than 2 revolutions) do not work nearly as well. This restriction comes in because the relationship between node 1 and node N becomes increasingly nonlinear as the angular distance grows. Trans-fers with more than about 1.5 complete revolutions generally require the line search to improve convergence. Future work includes: Comparison of this algorithm with other established tools; improvements to how multiple-revolution transfers are handled; parallelization of the Jacobian computation; in-creased efficiency for the line search; and optimization of many more trajectories between a variety of 3-body orbits.

Parrish, Nathan L.↗
NASA NTRS2010

Human Operations Beyond LEO by the End of the Decade: An Affordable Near-Term "Stepping Stone"

For more than a decade, several teams have assessed designs for a long-duration free-space human habitat beyond low-Earth orbit (LEO), building upon years of hard-won experience with the International Space Station (ISS). These systems would enable multiple achievements for science and human space flight. Most were intended to be deployed using available or near-future capabilities within about a decade after funding begins and serve as the first major human "stepping stone" beyond LEO. Last year, Thronson and Talay summarized work up to that time on expandable or inflatable concepts for deployment at an Earth-Moon (E-M) L1 or L2 location. Here we summarize our team's more recent work both on a long-duration human habitat that could be deployed beyond LEO within a decade and on the priority goals that such a habitat might accomplish. Particulars of this and other concepts for human operations in cis-lunar space are posted on the web and will be presented at professional conferences, and detailed in future publications by our group.

Thronson, Harley↗
NASA NTRS1991

Characteristics of return stroke electric fields produced by lightning flashes at distances of 1 to 15 kilometers

Electric field derivative signals from single and multiple lightning strokes are presented. For about 25 pct. of all acquired waveforms, produced by return strokes, stepped leaders or intracloud discharges, type and distance of the signal source are known from the observations by an all sky video camera system. The analysis of the electric field derivative waveforms in the time domain shows a significant difference in the impulse width between return stroke signals and those of stepped leaders and intracloud discharges. In addition, the computed amplitude density spectrum of return stroke waveforms lies by a factor of 10 above that of stepped leaders and intracloud discharges in the frequency range from 50 to 500 kHz.

Hopf, CH.↗
NASA NTRS2014

The Combustion Experiment on the Sample Analysis at Mars (SAM) Instrument Suite on the Curiosity Rover

The combustion experiment on the Sample Analysis at Mars (SAM) suite on Curiosity will heat a sample of Mars regolith in the presence of oxygen and measure composition of the evolved gases using quadrupole mass spectrometry (QMS) and tunable laser spectrometry (TLS). QMS will enable detection of combustion products such as CO, CO2, NO, and other oxidized species, while TLS will enable precise measurements of the abundance and carbon isotopic composition (delta(sup 13)C) of the evolved CO2 and hydrogen isotopic composition (deltaD) of H2O. SAM will perform a two-step combustion to isolate combustible materials below approx.550 C and above approx.550 C. The combustion experiment on SAM, if properly designed and executed, has the potential to answer multiple questions regarding the origins of volatiles seen thus far in SAM evolved gas analysis (EGA) on Mars. Constraints imposed by SAM and MSL time and power resources, as well as SAM consumables (oxygen gas), will limit the number of SAM combustion experiments, so it is imperative to design an experiment targeting the most pressing science questions. Low temperature combustion experiments will primarily target the quantification of carbon (and nitrogen) contributed by SAM wet chemistry reagants MTBSTFA (N-Methyl-N-tert-butyldimethylsilyltrifluoroacetamide) and DMF (Dimethylformamide), which have been identified in the background of blank and sample runs and may adsorb to the sample while the cup is in the Sample Manipulation System (SMS). In addition, differences between the sample and "blank" may yield information regarding abundance and delta(sup 13)C of bulk (both organic and inorganic) martian carbon. High temperature combustion experiments primarily aim to detect refractory organic matter, if present in Cumberland fines, as well as address the question of quantification and deltaD value of water evolution associated with hydroxyl hydrogen in clay minerals.

Mars↗
NASA NTRS2013

Inter-Calibration and Concatenation of Climate Quality Infrared Cloudy Radiances from Multiple Instruments

A change in climate is not likely captured from any single instrument, since no single instrument can span decades of time. Therefore, to detect signals of global climate change, observations from many instruments on different platforms have to be concatenated. This requires careful and detailed consideration of instrumental differences such as footprint size, diurnal cycle of observations, and relative biases in the spectral brightness temperatures. Furthermore, a common basic assumption is that the data quality is independent of the observed scene and therefore can be determined using clear scene data. However, as will be demonstrated, this is not necessarily a valid assumption as the globe is mostly cloudy. In this study we highlight challenges in inter-calibration and concatenation of infrared radiances from multiple instruments by focusing on the analysis of deep convective or anvil clouds. TRMM/VIRS is potentially useful instrument to make correction for observational differences in the local time and foot print sizes, and thus could be applied retroactively to vintage instruments such as AIRS, IASI, IRIS, AVHRR, and HIRS. As the first step, in this study, we investigate and discuss to what extent AIRS and VIRS agree in capturing deep cloudy radiances at the same local time. The analysis also includes comparisons with one year observations from CrIS. It was found that the instruments show calibration differences of about 1K under deep cloudy scenes that can vary as a function of land type and local time of observation. The sensitivity of footprint size, view angle, and spectral band-pass differences cannot fully explain the observed differences. The observed discrepancies can be considered as a measure of the magnitude of issues which will arise in the comparison of legacy data with current data.

cloudy radiance↗
NASA NTRS1975

Prediction of health levels by remote sensing

Measures of the environment derived from remote sensing were compared to census population/housing measures in their ability to discriminate among health status areas in two urban communities. Three hypotheses were developed to explore the relationships between environmental and health data. Univariate and multiple step-wise linear regression analyses were performed on data from two sample areas in Houston and Galveston, Texas. Environmental data gathered by remote sensing were found to equal or surpass census data in predicting rates of health outcomes. Remote sensing offers the advantages of data collection for any chosen area or time interval, flexibilities not allowed by the decennial census.

Rush, M.↗
NASA NTRS1988

Trace element and isotope geochemistry of Cretaceous-Tertiary boundary sediments: identification of extra-terrestrial and volcanic components

Trace element and stable isotope analyses were performed on a series of sediment samples crossing the Cretaceous-Tertiary (K-T) boundary from critical sections at Aumaya and Sopelano, Spain. The aim is to possibly distinguish extraterrestrial vs. volcanic or authigenic concentration of platinum group and other elements in K-T boundary transitional sediments. These sediments also have been shown to contain evidence for step-wise extinction of several groups of marine invertebrates, associated with negative oxygen and carbon isotope excursions occurring during the last million years of the Cretaceous. These isotope excursions have been interpreted to indicate major changes in ocean thermal regime, circulation, and ecosystems that may be related to multiple events during latest Cretaceous time. Results to date on the petrographic and geochemical analyses of the Late Cretaceous and Early Paleocene sediments indicate that diagenesis has obviously affected the trace element geochemistry and stable isotope compositions at Zumaya. Mineralogical and geochemical analysis of K-T boundary sediments at Zumaya suggest that a substantial fraction of anomalous trace elements in the boundary marl are present in specific mineral phases. Platinum and nickel grains perhaps represent the first direct evidence of siderophile-rich minerals at the boundary. The presence of spinels and Ni-rich particles as inclusions in aluminosilicate spherules from Zumaya suggests an original, non-diagenetic origin for the spherules. Similar spherules from southern Spain (Caravaca), show a strong marine authigenic overprint. This research represents a new approach in trying to directly identify the sedimentary mineral components that are responsible for the trace element concentrations associated with the K-T boundary.

Margolis, S. V.↗
NASA NTRS2001

Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations

Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Brandt, Achi↗
NASA NTRS2002

A Closed-Loop Hardware Simulation of Decentralized Satellite Formation Control

In recent years, there has been significant interest in the use of formation flying spacecraft for a variety of earth and space science missions. Formation flying may provide smaller and cheaper satellites that, working together, have more capability than larger and more expensive satellites. Several decentralized architectures have been proposed for autonomous establishment and maintenance of satellite formations. In such architectures, each satellite cooperatively maintains the shape of the formation without a central supervisor, and processing only local measurement information. The Global Positioning System (GPS) sensors are ideally suited to provide such local position and velocity measurements to the individual satellites. An investigation of the feasibility of a decentralized approach to satellite formation flying was originally presented by Carpenter. He extended a decentralized linear-quadratic-Gaussian (LQG) framework proposed by Speyer in a fashion similar to an extended Kalman filter (EKE) which processed GPS position fix solutions. The new decentralized LQG architecture was demonstrated in a numerical simulation for a realistic scenario that is similar to missions that have been proposed by NASA and the U.S. Air Force. Another decentralized architecture was proposed by Park et al. using carrier differential-phase GPS (CDGPS). Recently, Busse et al demonstrated the decentralized CDGPS architecture in a hardware-in-the-loop simulation on the Formation Flying TestBed (FFTB) at Goddard Space Flight Center (GSFC), which features two Spirent Cox 16 channel GPS signal generator. Although representing a step forward by utilizing GPS signal simulators for a spacecraft formation flying simulation, only an open-loop performance, in which no maneuvers were executed based on the real-time state estimates, was considered. In this research, hardware experimentation has been extended to include closed-loop integrated guidance and navigation of multiple spacecraft formations using GPS receivers and real-time vehicle telemetry. A hardware closed-loop simulation has been performed using the decentralized LQG architecture proposed by Carpenter in the GPS test facility at the Center for Space Research (CSR). This is the first presentation using this type of hardware for demonstration of closed-loop spacecraft formation flying.

Ebimuma, Takuji↗
NASA NTRS2013

Evaluation of Aerosol Mixing State Classes in the GISS Modele-matrix Climate Model Using Single-particle Mass Spectrometry Measurements

Aerosol particles in the atmosphere are composed of multiple chemical species. The aerosol mixing state, which describes how chemical species are mixed at the single-particle level, provides critical information on microphysical characteristics that determine the interaction of aerosols with the climate system. The evaluation of mixing state has become the next challenge. This study uses aerosol time-of-flight mass spectrometry (ATOFMS) data and compares the results to those of the Goddard Institute for Space Studies modelE-MATRIX (Multiconfiguration Aerosol TRacker of mIXing state) model, a global climate model that includes a detailed aerosol microphysical scheme. We use data from field campaigns that examine a variety of air mass regimens (urban, rural, and maritime). At all locations, polluted areas in California (Riverside, La Jolla, and Long Beach), a remote location in the Sierra Nevada Mountains (Sugar Pine) and observations from Jeju (South Korea), the majority of aerosol species are internally mixed. Coarse aerosol particles, those above 1 micron, are typically aged, such as coated dust or reacted sea-salt particles. Particles below 1 micron contain large fractions of organic material, internally-mixed with sulfate and black carbon, and few external mixtures. We conclude that observations taken over multiple weeks characterize typical air mass types at a given location well; however, due to the instrumentation, we could not evaluate mass budgets. These results represent the first detailed comparison of single-particle mixing states in a global climate model with real-time single-particle mass spectrometry data, an important step in improving the representation of mixing state in global climate models.

climate models↗
NASA NTRS2015

Numerical Computation of a Continuous-thrust State Transition Matrix Incorporating Accurate Hardware and Ephemeris Models

A significant body of work exists showing that providing a nonlinear programming (NLP) solver with expressions for the problem constraint gradient substantially increases the speed of program execution and can also improve the robustness of convergence, especially for local optimizers. Calculation of these derivatives is often accomplished through the computation of spacecraft's state transition matrix (STM). If the two-body gravitational model is employed as is often done in the context of preliminary design, closed form expressions for these derivatives may be provided. If a high fidelity dynamics model, that might include perturbing forces such as the gravitational effect from multiple third bodies and solar radiation pressure is used then these STM's must be computed numerically. We present a method for the power hardward model and a full ephemeris model. An adaptive-step embedded eight order Dormand-Prince numerical integrator is discussed and a method for the computation of the time of flight derivatives in this framework is presented. The use of these numerically calculated derivatieves offer a substantial improvement over finite differencing in the context of a global optimizer. Specifically the inclusion of these STM's into the low thrust missiondesign tool chain in use at NASA Goddard Spaceflight Center allows for an increased preliminary mission design cadence.

Low Thrust↗
NASA NTRS1995

Combustion Module 1: Spacelab Racks Integrated at the Lewis Research Center for the First Time

The Combustion Module-1 (CM-1), NASA's largest (over 1800 lb) and one of the most sophisticated combustion experiments ever to fly on the Spacelab, will be carried on the first Microgravity Science Laboratory (MSL-1) mission aboard the space shuttle flight STS-84 in April 1997. The CM-1 project is a stepping stone to the space station era, because the hardware can support multiple investigators on the same mission and can be integrated at user locations and shipped to the launch site. CM-1 is being developed to accommodate microgravity combustion experiments that are designed to help explain and predict the behavior of combustion processes. Although the two principal investigators for CM-1 are both studying combustion processes, their investigations are quite different. Professor Paul D. Ronney of the University of Southern California will examine the Structures of Flame Balls at Low Lewis Numbers (SOFBALL), in which a variety of fuel-lean gaseous mixtures fill the combustion chamber and are ignited. Professor Gerard M. Faeth of the University of Michigan will investigate Laminar Soot Processes (LSP) by studying the key properties of burning gas jets of fuel, employing different fuels and nozzle sizes.

Source record↗
NASA NTRS2023

Harnessing the Digital Transformation for Development of Electrified Aircraft Propulsion Control Systems

Hybrid electric aircraft propulsion is an emerging technology that presents a variety of potential benefits along with technical integration challenges. Developing these new propulsion architectures with their complex control systems, and ultimately proving their benefit, is a multistep process. This process includes concept development and analysis, dynamic simulation, hardware-in-the-loop testing, full-scale testing, and so on. This effort is being revolutionized and indeed enabled by new digital tools that support increasing the technology readiness level throughout the maturation process. As part of this Digital Transformation, NASA has developed a suite of publicly available digital tools that facilitate the path from concept to implementation. This paper describes the NASA-developed tools and puts them in the context of control system development for hybrid electric aircraft propulsion. The three MATLAB®-based software packages are the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS), the Electrical Modeling and Thermal Analysis Toolbox (EMTAT), and the Thermal Systems Analysis Toolbox (TSAT). These tools are interactive, complementary, and compatible with each other. T-MATS is a modular thermodynamic modeling framework designed for creating custom component level models of jet engines. EMTAT is a modeling framework used to simulate a variety of power electronic devices, using both physics-based and power flow calculations. TSAT is a framework for modeling and analysis of dynamic heat transfer. These packages all consist of graphical, drag-and-drop, parameterizable building blocks representing various components of the system to be modeled, e.g., compressors, turbines, motors, energy storage devices, etc. They are designed to enable the user to model and simulate the end-to-end dynamic operation of a hybrid electric gas turbine engine powertrain at the timescale of the turbomachinery, capturing mechanical, electrical, and thermal interactions. This paper demonstrates through multiple examples how these tools have been used successfully in a variety of applications, including several of the early stages of hybrid electric gas turbine engine propulsion system development, from the initial system modeling to real-time interactive pilot-in-the-loop simulation to physical hardware-in-the-loop testing, each step bringing the technology closer to fruition.

hybrid electric propulsion↗
NASA NTRS2010

Diamond Heat-Spreader for Submillimeter-Wave Frequency Multipliers

The planar GaAs Shottky diode frequency multiplier is a critical technology for the local oscillator (LO) for submillimeter- wave heterodyne receivers due to low mass, tenability, long lifetime, and room-temperature operation. The use of a W-band (75-100 GHz) power amplifier followed by a frequency multiplier is the most common for submillimeter-wave sources. Its greatest challenge is to provide enough input power to the LO for instruments onboard future planetary missions. Recently, JPL produced 800 mW at 92.5 GHz by combining four MMICs in parallel in a balanced configuration. As more power at W-band is available to the multipliers, their power-handling capability be comes more important. High operating temperatures can lead to degradation of conversion efficiency or catastrophic failure. The goal of this innovation is to reduce the thermal resistance by attaching diamond film as a heat-spreader on the backside of multipliers to improve their power-handling capability. Polycrystalline diamond is deposited by hot-filament chemical vapor deposition (CVD). This diamond film acts as a heat-spreader to both the existing 250- and 300-GHz triplers, and has a high thermal conductivity (1,000-1,200 W/mK). It is approximately 2.5 times greater than copper (401 W/mK) and 20 times greater than GaAs (46 W/mK). It is an electrical insulator (resistivity approx. equals 10(exp 15) Ohms-cm), and has a low relative dielectric constant of 5.7. Diamond heat-spreaders reduce by at least 200 C at 250 mW of input power, compared to the tripler without diamond, according to thermal simulation. This superior thermal management provides a 100-percent increase in power-handling capability. For example, with this innovation, 40-mW output power has been achieved from a 250-GHz tripler at 350-mW input power, while the previous triplers, without diamond, suffered catastrophic failures. This breakthrough provides a stepping-stone for frequency multipliers-based LO up to 3 THz. The future work for this design is to apply the high output power from both the 250 and 300 GHz to multiple chains in order to generate milliwatts at 2.3 THz. Using the first generation of results for this innovation, 40 mW of output power were produced from a 240-GHz tripler at 350-mW input power, and 27- mW output power was produced from a 300-GHz tripler at 408-mW input power. This is two times higher than the current state-of-the-art output power capability. A finite-element thermal simulation also shows that 30-microns thick diamond dropped the temperature of the anodes by at least 200 C.

Lin, Robert H.↗
NASA NTRS1987

Grid generation and inviscid flow computation about cranked-winged airplane geometries

An algebraic grid generation procedure that defines a patched multiple-block grid system suitable for fighter-type aircraft geometries with fuselage and engine inlet, canard or horizontal tail, cranked delta wing and vertical fin has been developed. The grid generation is based on transfinite interpolation and requires little computational power. A finite-volume Euler solver using explicit Runge-Kutta time-stepping has been adapted to this grid system and implemented on the VPS-32 vector processor with a high degree of vectorization. Grids are presented for an experimental aircraft with fuselage, canard, 70-20-cranked wing, and vertical fin. Computed inviscid compressible flow solutions are presented for Mach 2 at 3.79, 7 and 10 deg angles of attack. Conmparisons of the 3.79 deg computed solutions are made with available full-potential flow and Euler flow solutions on the same configuration but with another grid system. The occurrence of an unsteady solution in the 10 deg angle of attack case is discussed.

Eriksson, L.-E.↗
NASA NTRS2011

Advanced Spacesuit Portable Life Support System Oxygen Regulator Development and Testing

The advanced spacesuit portable life support system (PLSS) oxygen regulators represent an evolutionary approach to regulator development. Several technology development prototypes have been produced that borrow much of the mechanical regulator design from the well proven Shuttle/ISS Extravehicular Mobility Unit (EMU) Secondary Oxygen Regulator, but incorporate a motor-settable pressure set-point feature that facilitates significantly greater operational flexibility. For example, this technology would enable EVA to begin at a higher suit pressure, which would reduce pre-breathe time, and then slowly step down to a lower pressure to increase suit mobility for the duration of the EVA. Comprehensive testing of the prototypes was performed on the component level as well as part of the PLSS 1.0 system level testing. Results from these tests characterize individual prototype performance and demonstrate successful operation during multiple nominal and contingency EVA modes

Campbell, Colin↗
NASA NTRS2015

Formation Flying for Satellites and Unmanned Aerial Vehicles

The shrinking size of satellites and unmanned aerial vehicles (UAVs) is enabling lower cost missions. As sensors and electronics continue to downsize, the next step is multiple vehicles providing different perspectives or variations for more precise measurements. While flying a single satellite or UAV autonomously is a challenge, flying multiple vehicles in a precise formation is even more challenging. The goal of this project is to develop a scalable mesh network between vehicles (satellites or UAVs) to share real-time position data and maintain formations autonomously. Newly available low-cost, commercial off-the-shelf credit card size computers will be used as the basis for this network. Mesh networking techniques will be used to provide redundant links and a flexible network. The Small Projects Rapid Integration and Test Environment Lab will be used to simulate formation flying of satellites. UAVs built by the Aero-M team will be used to demonstrate the formation flying in the West Test Area. The ability to test in flight on NASA-owned UAVs allows this technology to achieve a high Technology Readiness Level (TRL) (TRL-4 for satellites and TRL-7 for UAVs). The low cost of small UAVs and the availability of a large test range (West Test Area) dramatically reduces the expense of testing. The end goal is for this technology to be ready to use on any multiple satellite or UAV mission.

Merrill, Garrick↗
NASA NTRS2003

Numerical Models of Human Circulatory System under Altered Gravity: Brain Circulation

A computational fluid dynamics (CFD) approach is presented to model the blood flow through the human circulatory system under altered gravity conditions. Models required for CFD simulation relevant to major hemodynamic issues are introduced such as non-Newtonian flow models governed by red blood cells, a model for arterial wall motion due to fluid-wall interactions, a vascular bed model for outflow boundary conditions, and a model for auto-regulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models are solved iteratively using the pseudocompressibility method and dual time stepping. Moving wall boundary conditions from the first-order fluid-wall interaction model are used to study the influence of arterial wall distensibility on flow patterns and wall shear stresses during the heart pulse. A vascular bed modeling utilizing the analogy with electric circuits is coupled with an auto-regulation algorithm for multiple outflow boundaries. For the treatment of complex geometry, a chimera overset grid technique is adopted to obtain connectivity between arterial branches. For code validation, computed results are compared with experimental data for steady and unsteady non-Newtonian flows. Good agreement is obtained for both cases. In sin-type Gravity Benchmark Problems, gravity source terms are added to the Navier-Stokes equations to study the effect of gravitational variation on the human circulatory system. This computational approach is then applied to localized blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other model using an anatomical data set. A three- dimensional anatomical Circle of Willis configuration is reconstructed from human-specific magnetic resonance images using an image segmentation method. The blood flow through these Circle of Willis models is simulated to provide means for studying gravitational effects on the brain circulation under auto-regulation.

Kim, Chang Sung↗