Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complex Loading”

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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Automated Work Package: Initial Wireless Communication Platform Design, Development, and Evaluation

The Department of Energy's Light Water Reactor Sustainability Program is developing the scientific basis to ensure long-term reliability, productivity, safety, and security of the nuclear power industry in the United States. The Instrumentation, Information, and Control (II and C) pathway of the program aims to increase the role of advanced II and C technologies to achieve this objective. One of the pathway efforts at Idaho National Laboratory (INL) is to improve the work packages execution process by replacing the expensive, inefficient, bulky, complex, and error-prone paper-based work orders with automated work packages (AWPs). An AWP is an automated and dynamic presentation of the work package designed to guide the user through the work process. It is loaded on a mobile device, such as a tablet, and is capable of communicating with plant equipment and systems to acquire plant and procedure states. The AWP replaces those functions where a computer is more efficient and reliable than a human. To enable the automatic acquisition of plant data, it is necessary to design and develop a prototype platform for data exchange between the field instruments and the AWP mobile devices. The development of the platform aims to reveal issues and solutions generalizable to large-scale implementation of a similar system. Topics such as bandwidth, robustness, response time, interference, and security are usually associated with wireless communication. These concerns, along with other requirements, are listed in an earlier INL report. Specifically, the targeted issues and performance aspects in this work are relevant to the communication infrastructure from the perspective of promptness, robustness, expandability, and interoperability with different technologies. The development of a prompt system entailed significant performance optimization, including the use of parallel processing, parallel communication, modularity, noise filtering, and frequency tuning. A robust system was achieved by enforcing a persistent continuous data exchange that automatically detects and conveys a failure of any link or component in any part of the system, and that automatically restores the communication once the link or component is restored. An expandable topology was enabled by modularity and load sharing. The number of servers, mobile devices, and field instruments is not limited from a system structure point of view. The system interoperability was achieved by enforcing the independence of the field communication application from the mobile device application. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

SULI End of Summer Presentation

This presentation made on Power Point is a brief summary of my work enhancing the Supervisory Control and Data Acquisition (SCADA) architecture for Idaho National Lab’s Critical Infrastructure Test Range Complex (CITRC). This work was part of an on-going project to create a digital twin of CITRC that will be implemented for resiliency research. Resiliency is growing more important as extreme weather, aging infrastructure, and increasing load all exert stress on the United States power grid. This presentation highlights how I documented the existing devices and data flow on site and then updated specific device settings to better support the needs of the project. It also examines future use of this data in the digital twin data warehouse (DeepLynx), and future Artificial Intelligence and Advanced Distribution Management System applications. It will be presented to the B751 Water and Energy Systems Analysis group at INL during my last week.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Origins of Clustering of Metalate–Extractant Complexes in Liquid–Liquid Extraction

Effective and energy-efficient separation of precious and rare metals is very important for a variety of advanced technologies. Liquid-liquid extraction (LLE) is a relatively less energy intensive separation technique, widely used in separation of lanthanides, actinides, and platinum group metals (PGMs). In LLE, the distribution of an ion between an aqueous phase and an organic phase is determined by enthalpic (coordination interactions) and entropic (fluid reorganization) contributions. The molecular scale details of these contributions are not well understood. Preferential extraction of an ion from the aqueous phase is usually correlated with the resulting fluid organization in the organic phase, as the longer-range organization increases with metal loading. However, it is difficult to determine the extent to which organic phase fluid organization causes, or is caused by, metal loading. In this study, we demonstrate that two systems with the same metal loading may impart very different organic phase organizations and investigate the underlying molecular scale mechanism. Small-angle X-ray scattering shows that the structure of a quaternary ammonium extractant solution in toluene is affected differently by the extraction of two metalates (octahedral PtCl 6 2- and square-planar PtCl 6 2- ), although both are completely transferred into the organic phase. The aggregates formed by the metalate-extractant complexes (approximated as reverse micelles) exhibit a more long-range order (clustering) with PtCl 6 2- compared to that with PtCl 6 2- . Vibrational sum frequency generation spectroscopy and complementary atomistic molecular dynamics simulations on model Langmuir monolayers indicate that the two metalates affect the interfacial hydration structures differently. Furthermore, the interfacial hydration is correlated with water extraction into the organic phase. Overall, these results support a strong relationship between the organic phase organizational structure and the different local hydration present within the aggregates of metalate-extractant complexes, which is independent of metalate concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Scatter in Carbon/Silicon Carbide (C/SiC) Composites Quantified

Carbon-fiber-reinforced silicon carbide matrix (C/SiC) composites processed by chemical vapor infiltration are candidate materials for aerospace thermal structures. Carbon fibers can retain properties at very high temperatures, but they are known to have poor oxidation resistance in adverse, high-temperature environments. Nevertheless, the combination of CVI-SiC matrix with higher stiffness and oxidation resistance, the interfacial coating, and additional surface-seal coating provides the necessary protection to the carbon fibers, and makes the material viable for high-temperature space applications operating under harsh environments. Furthermore, C/SiC composites, like other ceramic matrix composites (CMCs), exhibit graceful non-catastrophic failure because of various inherent energy dissipating mechanisms. The material exhibits nonlinearity in deformation even at very low stress levels. This is the result of the severe matrix microcracking present in the as processed composite because of large differences between the coefficients of thermal expansion of the fiber and the matrix. Utilization of these advanced composites in next generation space vehicles will require innovative structural configurations, updated materials, and refined analyses. Structural safety issues for these vehicles are in direct competition with performance and cost. One would have to quantify the uncertainties associated with the design using formal probabilistic methods. Specifically four fundamental aspects on which analyses are based-- (1) loading conditions, (2) material behavior, (3) geometrical configurations, and (4) structural connections between the composite components and baseline structure--are stochastic in nature. A direct way to formally account for uncertainties is to develop probabilistic structural analysis methods where all participating variables are described by appropriate probability density functions. The present work, however, focuses on analyzing the stochastic material behavior of these advanced composites using formal probabilistic analysis methods. Often, some of the desirable property characteristics that allow composites to offer advantages over conventional structural materials (like tailoring of composite properties) and the complexity are in fact responsible for their greater statistical variability and the requirements for more characterization tests. Composite properties are anisotropic as well, having different properties in different directions. This means that characterization of a property such as stiffness--which will vary greatly depending on the orientation of the fiber relative to the direction of the testing--must be repeated for several different directions and loading conditions. The fabrication process for composites also introduces statistical variations in properties and geometry. A composite part is produced in a number of steps, each of which introduces statistical variability. The matrix is usually produced from a combination of raw materials; and the fiber, which has its own set of properties, is often coated or surface treated, introducing yet another source of variability.

Murthy, Pappu L. N.↗
DOE OSTI2026

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗
DOE OSTI2024

Evaluating SWAT + model uncertainties for human and natural outcomes: Application in a Great Lakes agricultural watershed

Nutrient exports from agricultural lands in the Great Lakes Region pose significant threats to water quality and ecological health through eutrophication, hypoxia, and harmful algal blooms. Climate change and agricultural adaptation practices complicate future nutrient loading due to intensified hydrologic cycles and land use decisions. Our research focuses on evaluating the Soil and Water Assessment Tool (SWAT) plus model parametric uncertainties for human and natural outcomes across different scales. These factors are integral to ensuring a balance between productive agricultural practices and maintaining the health of watershed hydrology. However, uncertainties in modeling such complex interactions pose significant challenges, limiting our ability to precisely determine critical factors that influence crop yield and soil moisture. Our analysis employs Sobol global sensitivity analysis to evaluate first-order, second order, and total-order indices for SWAT crop growth parameter, ensuring comprehensive assessment of individual and interactive effects on model outputs. The objective is to identify the parameters that significantly affect model outputs for crop yield and soil moisture and improve our understanding of their interactions at the basin and hydrological response unit (HRU) scale. Our case study, the Portage River Watershed, which drains into Lake Erie, is chosen to better capture finer scale interactions crucial for predicting nutrient loading under future climate scenarios. This foundational work is aimed at setting the stage for the future development of an agent-based model (ABM). The ABM model would incorporate SWAT outputs to dynamically simulate decision-making processes.

Bunyon, Enock↗
NASA NTRS2009

Influence of Combined Whole-Body Vibration Plus G-Loading on Visual Performance

Recent engineering analyses of the integrated Ares-Orion stack show that vibration levels for Orion crews have the potential to be much higher than those experienced in Gemini, Apollo, and Shuttle vehicles. Of particular concern to the Constellation Program (CxP) is the 12 Hz thrust oscillation (TO) that the Ares-I rocket develops during the final ~20 seconds preceding first-stage separation, at maximum G-loading. While the structural-dynamic mitigations being considered can assure that vibration due to TO is reduced to below the CxP crew health limit, it remains to be determined how far below this limit vibration must be reduced to enable effective crew performance during launch. Moreover, this "performance" vibration limit will inform the operations concepts (and crew-system interface designs) for this critical phase of flight. While Gemini and Apollo studies provide preliminary guidance, the data supporting the historical limits were obtained using less advanced interface technologies and very different operations concepts. In this study, supported by the Exploration Systems Mission Directorate (ESMD) Human Research Program, we investigated display readability-a fundamental prerequisite for any interaction with electronic crew-vehicle interfaces-while observers were subjected to 12 Hz vibration superimposed on the 3.8 G loading expected for the TO period of ascent. Two age-matched groups of participants (16 general population and 13 Crew Office) performed a numerical display reading task while undergoing sustained 3.8 G loading and whole-body vibration at 0, 0.15, 0.3, 0.5, and 0.7 g in the eyeballs in/out (x-axis) direction. The time-constrained reading task used an Orion-like display with 10- and 14-pt non-proportional sans-serif fonts, and was designed to emulate the visual acquisition and processing essential for crew system monitoring. Compared to the no-vibration baseline, we found no significant effect of vibration at 0.15 and 0.3 g on task error rates (ER) or response times (RT). Significant degradations in both ER and RT, however, were observed at 0.5 and 0.7 g for 10-pt, and at 0.7 g for 14-pt font displays. These objective performance measures were mirrored by participants' subjective ratings. Interestingly, we found that the impact of vibration on ER increased with distance from the center of the display, but only for vertical displacements. Furthermore, no significant ER or RT aftereffects were detected immediately following vibration, regardless of amplitude. Lastly, given that our reading task required no specialized spaceflight expertise, our finding that effects were not statistically distinct between our two groups is not surprising. The results from this empirical study provide initial guidance for evaluating the display readability trade-space between text-font size and vibration amplitude. However, the outcome of this work should be considered preliminary in nature for a number of reasons: 1. The single 12 Hz x-axis vibration employed was based on earlier load-cycle models of the induced TO environment at the end of Ares-I first stage flight. Recent analyses of TO mitigation designs suggest that significant concurrent off-axis vibration may also occur. 2. The shirtsleeve environment in which we tested fails to capture the full kinematic and dynamic complexity of the physical interface between crewmember and the still-to-bematured helmet-suit-seat designs, and the impact these will have for vibration transmission and consequent performance. 3. By examining performance in this reading and number processing task, we are only assessing readability, a first and necessary step that in itself does not directly address the performance of more sophisticated operational tasks such as vehicle-health monitoring or manual control of the vehicle.

readability↗
DOE OSTI2022

New large-strain FFT-based formulation and its application to model strain localization in nano-metallic laminates and other strongly anisotropic crystalline materials

This paper presents a new robust large-strain (LS) elasto-viscoplastic (EVP) formulation based on Fast Fourier Transforms (FFTs) for the prediction of the micro-mechanical response and microstructure evolution of polycrystalline and multiphase materials, with emphasis on the effect of strong crystallographic and/or morphologic anisotropy on localization of plastic deformation. In this work, the novel LS-EVPFFT formulation allows treatment of complex initial geometries and large deformations considering three grids of material points: a regular grid in the reference configuration, where FFTs can be performed; an irregular grid in the initial configuration, created by applying a stress-free displacement field to the reference regular grid; and an irregular grid in the current configuration, undergoing large strains and rotations as the material is loaded. Further numerical stability of the new formulation also required the use of a novel expression for the discrete modified Green’s operator, which reduces spurious field oscillations. After presenting and validating the new formulation by comparison with preexisting implementations and analytical solutions, LS-EVPFFT is applied to the prediction of slip and kink bands formation in polycrystalline columnar ice, and kink bands in single crystal zinc wires, showing good agreement with classic experiments. Finally, the model is used to study kink band formation during compression of Cu–Nb nano-metallic laminates (NMLs), in which accurate treatment of the complex geometry associated with the tortuosity of interfaces and large deformations become critical, showing consistency with corresponding micropillar experiments.

36 MATERIALS SCIENCE↗
NASA NTRS1989

The measurement of boundary layers on a compressor blade in cascade. Volume 1: Experimental technique, analysis and results

Measurements were made of the boundary layers and wakes about a highly loaded, double-circular-arc compressor blade in cascade. These laser Doppler velocimetry measurements have yielded a very detailed and precise data base with which to test the application of viscous computational codes to turbomachinery. In order to test the computational codes at off-design conditions, the data were acquired at a chord Reynolds number of 500,000 and at three incidence angles. Moreover, these measurements have supplied some physical insight into these very complex flows. Although some natural transition is evident, laminar boundary layers usually detach and subsequently reattach as either fully or intermittently turbulent boundary layers. These transitional separation bubbles play an important role in the development of most of the boundary layers and wakes measured in this cascade and the modeling or computing of these bubbles should prove to be the key aspect in computing the entire cascade flow field. In addition, the nonequilibrium turbulent boundary layers on these highly loaded blades always have some region of separation near the trailing edge of the suction surface. These separated flows, as well as the subsequent near wakes, show no similarity and should prove to be a challenging test for the viscous computational codes.

Zierke, William C.↗
DOE OSTI2014

Introduction of π-Complexation into Porous Aromatic Framework for Highly Selective Adsorption of Ethylene over Ethane

We report herein a strategy of incorporating air stable Ag(I) ions into water stable, high surface area porous organic polymer (POP) affording significant increase in ethylene uptake capacity and extremely high Qst for ethylene (over 100 kJ/mol at low ethylene load-ing) as illustrated in the context of Ag(I) ion functionalized PAF-1, PAF-1-SO3Ag. IAST calculations using single-component-isotherm data and equimolar ethylene/ethane ratio at 296 K reveal PAF-1-SO3Ag shows exceptionally high ethylene/ethane adsorption selectivi-ty (Sads: 27 to 125), far surpassing benchmark zeolite and any other MOF reported in literature. This alongside excellent water/air stability, high ethylene uptake capacity, and mild regeneration requirements make PAF-1-SO3Ag hold promise for adsorption-based eth-ylene/ethane separations, paving a way to develop Ag(I) ion function-alized POPs as a new platform for highly selective adsorption of eth-ylene over ethane.

Li, Baiyan↗
NASA NTRS2018

Efficacy of FTIR Analysis in Determining CO2 Loading on Diglycolamine

In support of advanced air revitalization technologies to enable human spaceflight beyond low earth orbit, performance studies have been conducted using a liquid amine, Diglycolamine (DGA) between teams at NASA?s Johnson Spaceflight Center (JSC) and Ames Research Center (ARC). Liquid amines have been used in regenerable earth-based systems to remove CO2 from industrial systems as well as for closed-environment air revitalization because they can be regenerated at lower temperatures than solid sorbent systems. As an additional advantage to solid sorbent based systems, liquid sorbents can be cycled between an adsorbing contactor and degassing chamber, thereby reducing system complexity by operation in a continuous loop. In an effort to inform a regeneration system design for micro-gravity applications, ARC has performed a number of tests to characterize the degas mechanics of DGA. In order to accurately measure the amount of CO2 captured or released by the amine, methods such as gravimetric weighing and chemical decomposition are reasonable, however the initial design for a scaled down degas system required analysis on small sample sizes. In order to accomplish CO2 loading measurements with a sample size of approximately 500µL, Fourier-transform infrared spectroscopy (FTIR) was used. Calibration against acid decomposition showed relatively good correlation and test data showed reasonable adherence to expected trends, however more extensive testing should be conducted to fully validate the usage of FTIR to determine CO2 loading on DGA.

Diglycolamine↗
DOE OSTI2014

Statistics of reversible transitions in two-state trajectories in force-ramp spectroscopy

A possible way to extract information about the reversible dissociation of a molecular adhesion bond from force fluctuations observed in force ramp experiments is discussed. For small loading rates the system undergoes a limited number of unbinding and rebinding transitions observable in the so-called force versus extension (FE) curves. The statistics of these transient fluctuations can be utilized to estimate the parameters for the rebinding rate. This is relevant in the experimentally important situation where the direct observation of the reversed FE-curves is hampered, e.g., due to the presence of soft linkers. I generalize the stochastic theory of the kinetics in two-state models to the case of time-dependent kinetic rates and compute the relevant distributions of characteristic forces. While for irreversible systems there is an intrinsic relation between the rupture force distribution and the population of the free-energy well of the bound state, the situation is slightly more complex if reversible systems are considered. For a two-state model, a “stationary” rupture force distribution that is proportional to the population can be defined and allows to consistently discuss quantities averaged over the transient fluctuations. While irreversible systems are best analyzed in the soft spring limit of small pulling device stiffness and large loading rates, here I argue to use the stiffness of the pulling device as a control parameter in addition to the loading rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Sandia fracture challenge 2: Sandia California's modeling approach

The second Sandia Fracture Challenge illustrates that predicting the ductile fracture of Ti-6Al-4V subjected to moderate and elevated rates of loading requires thermomechanical coupling, elasto-thermo-poro-viscoplastic constitutive models with the physics of anisotropy and regularized numerical methods for crack initiation and propagation. We detail our initial approach with an emphasis on iterative calibration and systematically increasing complexity to accommodate anisotropy in the context of an isotropic material model. Blind predictions illustrate strengths and weaknesses of our initial approach. We then revisit our findings to illustrate the importance of including anisotropy in the failure process. Furthermore, mesh-independent solutions of continuum damage models having both isotropic and anisotropic yields surfaces are obtained through nonlocality and localization elements.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Thermal, Rheological, Mechanical and Morphological Behavior of High Density Polyethylene and Carboxymethyl Cellulose Blend

In this study water soluble sodium carboxymethyl cellulose (CMC) was blended with high density polyethylene (HDPE) by peroxide-initiated melt compounding technique. The compatibility of the blended polymers were carried out by silane crosslinking agent. A series of blends were prepared by varying the CMC contents up to a maximum of 50 phr. The physical properties of non-crosslinked and crosslinked blends were investigated in detail. FTIR analysis of crosslinked blend confirmed the presence of Si–O–Si and Si–O–C absorption peaks at 1050 and 1159 cm{sup −1}. Thermal stability of crosslinked blends improved as compared to its non-crosslinked congener. Rheological study of crosslinked blends illustrated high complex viscosity and dynamic shear storage modulus. The tensile strength of virgin polyethylene was 8.1 MPa whereas the maximum tensile strength of 19.6 MPa was observed in crosslinked blend. Similarly lower deformation was observed in crosslinked blends under static load. Scanning electron microscopy of crosslinked formulations also showed strong adhesion between the polymers interface. The compatibility of HDPE and CMC is attributed to both free radical and condensation reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Influence of build direction on the fracture mechanism of 3D printed octet lattices

Here, we investigate the effects of 3D printing build directions on the fracture properties of octet lattice metamaterials made of polylactic acid (PLA) and how these effects vary with the relative density of the lattices. Single-edge notch bend samples are 3D printed in two orthogonal build directions at various relative densities. Our results show that the work of fracture for octet lattices with build directions parallel to the crack plane is significantly higher than those with build directions perpendicular to the crack plane. We also observed that the ratio of specific work of fracture between the two build directions remains nearly constant across different relative densities. In contrast, the ratio of peak load between the two build directions decreases as relative density increases. Furthermore, the build direction dictates the fracture mechanism. While the perpendicular build direction predominantly results in a brittle (mode I) fracture, the parallel build direction leads to a complex, delamination-dominated (mode II type) fracture. This phenomenon is largely governed by the weak interfaces formed between the printed layers and their interaction with the lattice geometry. These results reveal that the build direction governs the fracture mechanism and work of fracture of these lattice metamaterials and is therefore an important design consideration.

3D printed PLA↗
DOE OSTI2025

Direct ink writing of aqueous-based Gadolinium (III) oxide slurries

Gadolinium (III) oxide (gadolinia, Gd 2 O 3 ) has recently been identified as an intriguing material for applications in the medical, solid oxide fuel cell, and nuclear industries. This interest drives the need for developing and understanding manufacturing techniques that can produce dense Gd 2 O 3 structures. Direct ink writing (DIW), an extrusion-based additive manufacturing method, has also garnered interest because of its capability to produce dense ceramic parts with increased complexity in an economical manner. In this study, DIW was explored as a manufacturing technique for Gd 2 O 3 . Experiments were performed to develop Gd 2 O 3 bearing inks capable of being processed via DIW. Ink solids loading and sintering temperatures were varied to assess their impact on the final density and microstructure. Optimum sintering conditions are proposed and were experimentally verified at a dwell temperature of 1500°C. Gd 2 O 3 samples were successfully manufactured using DIW, achieving densities greater than 96 % of the theoretical density.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bayesian blacksmithing: discovering thermomechanical properties and deformation mechanisms in high-entropy refractory alloys

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Dynamic load balance scheme for the DSMC algorithm

The direct simulation Monte Carlo (DSMC) algorithm, devised by Bird, has been used over a wide range of various rarified flow problems in the past 40 years. While the DSMC is suitable for the parallel implementation on powerful multi-processor architecture, it also introduces a large load imbalance across the processor array, even for small examples. The load imposed on a processor by a DSMC calculation is determined to a large extent by the total of simulator particles upon it. Since most flows are impulsively started with initial distribution of particles which is surely quite different from the steady state, the total of simulator particles will change dramatically. The load balance based upon an initial distribution of particles will break down as the steady state of flow is reached. The load imbalance and huge computational cost of DSMC has limited its application to rarefied or simple transitional flows. In this paper, by taking advantage of METIS, a software for partitioning unstructured graphs, and taking the total of simulator particles in each cell as a weight information, the repartitioning based upon the principle that each processor handles approximately the equal total of simulator particles has been achieved. The computation must pause several times to renew the total of simulator particles in each processor and repartition the whole domain again. Thus the load balance across the processors array holds in the duration of computation. The parallel efficiency can be improved effectively. The benchmark solution of a cylinder submerged in hypersonic flow has been simulated numerically. Besides, hypersonic flow past around a complex wing-body configuration has also been simulated. The results have displayed that, for both of cases, the computational time can be reduced by about 50%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗