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 829 records · Page 46

Evaluation of Hybrid FPOG Applications in Regulated and Deregulated Markets Using HERON

Recent changes in the U.S. energy market, such as low natural gas prices and increased electricity production for variable renewable energy (VRE) sources, have led to an economic crisis for existing light-water reactor (LWR) nuclear power plants (NPP). Many owners and operators of LWRs have elected to decommission these plants rather than continue using them as consistent sources of clean baseload power. This has led to exploration of various possibilities to increase the economic viability of these units, including market restructuring to monetize benefits LWRs already provide to the grid through ancillary markets, load following and economic dispatch, and possible integration of secondary systems directly to the NPP for production of additional products through technologies such as hydrogen electrolysis or water desalination. Previous studies have considered the technologies associated with these Integrated Energy Systems (IES) activities, and the analysis of markets for these secondary products. To analyze the economic viability of various system configurations including IES, especially given the uncertainty surrounding load demand, electricity prices, and the availability of VRE resources, the stochastic technoeconomic analysis package HERON (Heuristic Energy Resource Optimization Network) was released earlier this year as an extension of the risk analysis framework RAVEN (Risk Analysis Virtual Environment). HERON focuses foremost on making the complex uncertainty quantification analysis tools approachable for energy systems analysts, also providing general dispatch optimization algorithms for those workflows. HERON continues to be improved and tested as a significant part of the IES viability analyses performed in this work. HERON is not a capacity expansion model. To consider market and grid energy system development in a variety of scenarios, HERON is best used in coupling with modelling tools such as US-REGEN, which sacrifice some of the uncertainty analysis and resolution of HERON's modelling for the ability to efficiently predict the change in the grid energy system's profile due to economic drivers over decades. HERON can then use this information to explore the economic viability of introducing changes to the predicted outcomes, such as the introduction of an IES. In this work, experts at EPRI using US-REGEN provide six projection scenarios for use in HERON stochastic technoeconomic analysis (STEA) in considering the options available for increasing LWR economic viability through introduction of a hydrogen-centric IES using a high-temperature steam electrolysis plant (HTSE), hydrogen storage, and a constant-rate contracted hydrogen consumer. The results obtained are differential in nature; they do not report expected profits for any configuration, but rather report on the possible increase in the NPV of a configuration with respect to a baseline no-IES configuration. Due to the uncertainty captured in the variable net load of the systems, there is likewise uncertainty in the mean values reported. We consider this viability both in terms of a regulated market, where the energy producers and IES are owned and operated by single entity, as well as a deregulated market, where the IES chooses its bid for electricity generation and is then dispatched by the grid system operator. Results indicate that for deregulated markets, the inclusion of the IES is often statistically beneficial. This is especially true in policies that are not favorable towards nuclear, as nuclear is less often dispatched and is forced to deal with frequent idle capacity. In the nominal case as well as the case of carbon tax policies, inclusion of the IES clearly benefited the economic performance of the NPP. In the regulated case, however, there was a trend towards minimizing the IES, likely due to the optimal sizing performed by US-REGEN of the NPP within the system as well as the lack of penalty for idle capacity at the NPP in the regulated market analyses.

99 GENERAL AND MISCELLANEOUS↗

Scalable Pattern Matching in Metadata Graphs via Constraint Checking

Pattern matching is a fundamental tool for answering complex graph queries. Unfortunately, existing solutions have limited capabilities: They do not scale to process large graphs and/or support only a restricted set of search templates or usage scenarios. Moreover, the algorithms at the core of the existing techniques are not suitable for today’s graph processing infrastructures relying on horizontal scalability and shared-nothing clusters, as most of these algorithms are inherently sequential and difficult to parallelize. In this article we present an algorithmic pipeline that bases pattern matching on constraint checking. The key intuition is that each vertex and edge participating in a match has to meet a set of constraints implicitly specified by the search template. These constraints can be verified independently and typically are less expensive to compute than searching the full template. The pipeline we propose generates these constraints and iterates over them to eliminate all the vertices and edges that do not participate in any match, thus reducing the background graph to a subgraph that is the union of all template matches—the complete set of all vertices and edges that participate in at least one match. Additional analysis can be performed on this annotated, reduced graph, such as full match enumeration, match counting, or computing vertex/edge centrality. Furthermore, a vertex-centric formulation for constraint checking algorithms exists, and this makes it possible to harness existing high-performance, vertex-centric graph processing frameworks. This technique (i) enables highly scalable pattern matching in metadata (labeled) graphs; (ii) supports arbitrary patterns with 100% precision; (iii) enables tradeoffs between precision and time-to-solution, while always selects all vertices and edges that participate in matches, thus offering 100% recall; and (iv) supports a set of popular data analytics scenarios. We implement our approach on top of HavoqGT, an open-source asynchronous graph processing framework, and demonstrate its advantages through strong and weak scaling experiments on massive scale real-world (up to 257 billion edges) and synthetic (up to 4.4 trillion edges) labeled graphs, respectively, and at scales (1,024 nodes / 36,864 cores), orders of magnitude larger than used in the past for similar problems. This article serves two purposes: First, it synthesises the knowledge accumulated during a long-term project. Second, it presents new system features, usage scenarios, optimizations, and comparisons with related work that strengthen the confidence that pattern matching based on iterative pruning via constraint checking is an effective and scalable approach in practice. The new contributions include the following: (i) We demonstrate the ability of the constraint checking approach to efficiently support two additional search scenarios that often emerge in practice, interactive incremental search and exploratory search. (ii) We empirically compare our solution with two additional state-of-the-art systems, Arabsque and TriAD. (iii) We show the ability of our solution to accommodate a more diverse range of datasets with varying properties, e.g., scale, skewness, label distribution, and match frequency. (iv) We introduce or extend a number of system features (e.g., work aggregation, load balancing, and the ability to cap the generated traffic) and design optimizations and demonstrate their advantages with respect to improving performance and scalability. (v) We present bottleneck analysis and insights into artifacts that influence performance. (vi) We present a theoretical complexity argument that motivates the performance gains we observe.

97 MATHEMATICS AND COMPUTING↗

Synthesis of palm-based polyurethane-LiClO{sub 4} via prepolymerization

Palm-based polyurethane (pPU) with varying lithium salt (LiClO{sub 4}) content was synthesized. Higher loading percentage of LiClO{sub 4} in the pPU led to the inhibition of prepolymerization process from taking place. Hydrogen bonded C=O was detected in the FTIR spectrum indicating the hydrogen bonding between the urethane bonds. Ordered complexed C=O was observed in the FTIR spectrum confirming the complex formation between urethane bond and Li{sup +} ion. DSC thermogram showed the increase in the LiClO{sub 4} content could increase the glass transition temperature. SEM micrographs exhibited that more bubbles were formed when the LiClO{sub 4} increased from 10 to 30wt% indicating the reaction between free isocyanate groups with moisture presence in the salt due to the hygroscopic properties of LiClO{sub 4}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Iron Availability within the Leaf Vasculature Determines the Magnitude of Iron Deficiency Responses in Source and Sink Tissues in Arabidopsis

Abstract Iron (Fe) uptake and translocation in plants are fine-tuned by complex mechanisms that are not yet fully understood. In Arabidopsis thaliana, local regulation of Fe homeostasis at the root level has been extensively studied and is better understood than the systemic shoot-to-root regulation. While the root system is solely a sink tissue that depends on photosynthates translocated from source tissues, the shoot system is a more complex tissue, where sink and source tissues occur synchronously. In this study, and to gain better insight into the Fe deficiency responses in leaves, we overexpressed Zinc/Iron-regulated transporter-like Protein (ZIP5), an Fe/Zn transporter, in phloem-loading cells (proSUC2::AtZIP5) and determined the timing of Fe deficiency responses in sink (young leaves and roots) and source tissues (leaves). Transgenic lines overexpressing ZIP5 in companion cells displayed increased sensitivity to Fe deficiency in root growth assays. Moreover, young leaves and roots (sink tissues) displayed either delayed or dampened transcriptional responses to Fe deficiency compared to wild-type (WT) plants. We also took advantage of the Arabidopsis mutant nas4x-1 to explore Fe transcriptional responses in the opposite scenario, where Fe is retained in the vasculature but in an unavailable and precipitated form. In contrast to proSUC2::AtZIP5 plants, nas4x-1 young leaves and roots displayed a robust and constitutive Fe deficiency response, while mature leaves showed a delayed and dampened Fe deficiency response compared to WT plants. Altogether, our data provide evidence suggesting that Fe sensing within leaves can also occur locally in a leaf-specific manner.

Cell Biology↗

Literature Review on the Theory and Measurement of Dry Sliding Friction of Metals

Understanding the friction behavior between two sliding bodies can inform the design of machines, processing of materials, and simulation of dynamic processes. Kinetic friction is a complex phenomenon that depends on a multitude of factors such as sliding velocity, normal force, contact area, surface roughness, material properties, lubrication conditions, and thermal effects. This literature review covers the major known effects of sliding velocity, normal load, and surface roughness on the measured kinetic friction coefficient in the context of microscopic friction phenomena. Classic macroscale friction models are reviewed to illustrate approaches for simulating friction behavior. Prominent experimental friction setups within the literature are discussed with respect to achievable velocity and pressure regimes. The background information gathered here will be used to inform experimental procedures and modeling strategies of the exploratory research (ER) project titled “Measurement of Dynamic Friction via Kolsky Bar” (20200418ER).

36 MATERIALS SCIENCE↗

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&C) pathway of the program aims to increase the role of advanced II&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.

Automated Work Packages↗

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 ↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗