Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modeling software”

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 433 records · Page 24

Modifications of Hinge Mechanisms for the Mobile Launcher

The further development and modifications made towards the integration of the upper and lower hinge assemblies for the Exploration Upper Stage umbilical are presented. Investigative work is included to show the process of applying updated NASA Standards within component and assembly drawings for selected manufacturers. Component modifications with the addition of drawings are created to precisely display part geometries and geometric tolerances, along with proper methods of fabrication. Comparison of newly updated components with original Apollo era components is essential to correctly model the part characteristics and parameters, i.e. mass properties, material selection, weldments, and tolerances. 3-Dimensional modeling software is used to demonstrate the necessary improvements. In order to share and corroborate these changes, a document management system is used to store the various components and associated drawings. These efforts will contribute towards the Mobile Launcher for Exploration Mission 2 to provide proper rotation of the Exploration Upper Stage umbilical, necessary for providing cryogenic fill and drain capabilities.

Exploration Upper Stage Umbilical↗

Effect of Near-Interface Compensation of CdSeTe Absorber Layers on Solar Cell Performance

Arsenic has been shown to be an effective p-type dopant of CdTe. However, challenges remain in the fabrication of high efficiency CdTe solar cells using As. As-doped CdTe is prone to self-compensation and observed accumulation of dopant atoms in CdTe near the interface with MgZnO (MZO) suggests that this could be occurring. In this study, we use SCAPS 1-D modeling software to investigate the effect of near-interface compensation. We consider three cases: shallow donors, deep recombination centers, and a thin layer of excess positive charge accumulation. All of these cases are shown to have significant effects on the current-voltage characteristics, while the thin charge layer also affects capacitance-voltage measurements.

arsenic↗

Three-Dimensional Heat Flux and Thermal Analysis of Angled Tungsten Samples on DIII-D

ITER-grade tungsten and dispersoid-strengthened tungsten samples with the top surface angled at ~15° towards the incident plasma flux were exposed to 9 H-mode discharges with edge-localized modes (ELMs) in the lower divertor of DIII-D tokamak using the Divertor Material Evaluation System (DiMES). Surface damage included cracking and flaking of material on the two samples farthest away from the plasma strike point, and significant melting of the two samples closest to the strike point. Heat flux and thermal analysis tools new to DIII-D have been applied to better understand this material response and to help optimize the exposure conditions for future experiments. SMITER field-line tracing simulations based on IRTV data and EFIT equilibria estimate an average inter-ELM perpendicular heat flux, 𝑞⊥,𝑖nter−𝐸LM , on the angled surfaces of 10.1 – 19.6 MW/m² for a majority of the 9 discharges, increasing to 15.6 – 24.5 MW/m² for the single, higher-power shot where samples melted. Fast camera data showed shallow intra-ELM melting and re-solidification, which transitioned to bulk inter-ELM melting with melt motion in the 𝐽⃗ 𝑥 𝐵⃗ direction. About 50% of the protruding volume of the most affected sample was displaced via melt-motion. SIERRA thermal modeling software was able to reproduce an onset time of melting consistent with fast camera data and final sample conditions, within < 200 ms. Maximum surface temperatures of 3122 K and 2787 K are estimated for the samples farthest away from the strike point, while the closest samples achieve melting at 4067 ms and 4750 ms into the ~5000 ms plasma exposure. A +10% increase in both the SMITER 𝑞⊥,𝑖nter−𝐸LM calculations and the estimated ELM heat loads 𝑞⊥, 𝐸LM was required to achieve this result, which is within the uncertainty of the diagnostic data but likely accounts for non-ideal geometry effects plus other physics uncertainties not included in this first iteration of modeling. This work provided valuable estimates of the 3D temperature evolution to help better understand the observed surface morphology and internal recrystallization of samples, which are discussed in detail in a complementary manuscript [1]. Benchmarking efforts with more diagnosed DIII-D experiments are underway to further refine the SMITER and SIERRA models for DiMES. Future use of these tools will enable researchers to precisely target heat flux exposure conditions in DIII-D to test, but not exceed, the thermomechanical limitations of novel plasma-facing materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations

The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested. In this study, we compare three classes of methods for antibody/antigen (Ab/Ag) binding affinity calculations: (i) a method that relies on the physical separation of the Ab/Ag complex in equilibrium molecular dynamics (MD) simulations, (ii) a collection of 18 scoring functions that act on an ensemble of structures created using homology modeling software, and (iii) methods based on the molecular mechanics-generalized Born surface area (MM-GBSA) energy decomposition, in which the individual contributions of the energy terms are scaled to optimize agreement with the experiment. When applied to a set of 49 antibody mutations in two Ab/HIV gp120 complexes, all of the methods are found to have modest accuracy, with the highest Pearson correlations reaching about 0.6. In particular, the most computationally intensive method, i.e., MD simulation, did not outperform several scoring functions. The optimized energy decomposition methods provided marginally higher accuracy, but at the expense of requiring experimental data for parametrization. Within each method class, we examined the effect of the number of independent computational replicates, i.e., modeled structures or reinitialized MD simulations, on the prediction accuracy. We suggest using about ten modeled structures for scoring methods, and about five simulation replicates for MD simulations as a rule of thumb for obtaining reasonable convergence. We anticipate that our study will be a useful resource for practitioners working to incorporate binding affinity calculations within their protein design and optimization process.

59 BASIC BIOLOGICAL SCIENCES↗

Statechart Analysis with Symbolic PathFinder

We report here on our on-going work that addresses the automated analysis and test case generation for software systems modeled using multiple Statechart formalisms. The work is motivated by large programs such as NASA Exploration, that involve multiple systems that interact via safety-critical protocols and are designed with different Statechart variants. To verify these safety-critical systems, we have developed Polyglot, a framework for modeling and analysis of model-based software written using different Statechart formalisms. Polyglot uses a common intermediate representation with customizable Statechart semantics and leverages the analysis and test generation capabilities of the Symbolic PathFinder tool. Polyglot is used as follows: First, the structure of the Statechart model (expressed in Matlab Stateflow or Rational Rhapsody) is translated into a common intermediate representation (IR). The IR is then translated into Java code that represents the structure of the model. The semantics are provided as "pluggable" modules.

Pasareanu, Corina S.↗

Parametric Cost Modeling of a Mid-Lift-to-Drag Ratio Vehicle for Human Mars Entry, Descent, and Landing

Cost is one of the biggest obstacles to sending humans to Mars. However, spacecraft costs are typically estimated after the preliminary vehicle and mission concepts have been designed. By automating the cost estimation process, the effect of any change in vehicle or mission design on the mission cost can be determined more rapidly. This paper uses the tool SAPE-C (Systems Analysis for Planetary Entry, Descent, and Landing-Cost). SAPE-C integrates the cost modeling software SEER-H (System Estimation and Evaluation of Resources-Hardware) with a number of systems analysis tools and enables large design tradespace exploration early in the design process. SAPE-C is used to analyze several tradespaces of the Mid-Lift-to-Drag ratio rigid entry vehicle for human Mars missions. Key findings include quantifying the impacts of payload, vehicle length, and main engine specific impulse on the overall cost of the vehicle, and how the payload per lander and number of landers affects the cost of a campaign to Mars.

Mid L/D↗

Storage Futures Study: Grid Operational Impacts of Widespread Storage Deployment

This report, the fifth in the Storage Futures Study series, uses cost-driven scenarios from the ReEDS model as a starting point to examine the operational impacts of grid-scale storage deployment and relationships between this deployment and the contribution of variable renewable energy. We use commercial production cost modeling software to evaluate hourly operation of five scenarios that reach between 210 gigawatts (GW) and 930 GW of installed storage by 2050. We find that storage plays an important role in these power systems between now and 2050 - by storing the lowest-marginal cost generation (often, overgeneration from solar or wind plants) and generating energy during the highest net load periods of the day and year. Storage helps with the integration of variable renewable energy and by providing an important resource to provide continued reliable power.

25 ENERGY STORAGE↗

A Model for Assessing the Liability of Seemingly Correct Software

Current research on software reliability does not lend itself to quantitatively assessing the risk posed by a piece of life-critical software. Black-box software reliability models are too general and make too many assumptions to be applied confidently to assessing the risk of life-critical software. We present a model for assessing the risk caused by a piece of software; this model combines software testing results and Hamlet's probable correctness model. We show how this model can assess software risk for those who insure against a loss that can occur if life-critical software fails.

Voas, Jeffrey M.↗

Aviation Simulators for the Desktop: Panel and Demonstrations

Panel Members are: Christine M. Mitchell (Georgia Tech), Michael T. Palmer (NASA Langley), Greg Pisani (NASA Ames), and Amy R. Pritchett (MIT). The Panel members are affiliated with aviation human factors groups from NASA Ames, NASA Langley, MITCHELL Department of Aerospace and Aeronautical Engineering, and Georgia Technics Center for Human-Machine Systems Research. Panelists will describe the simulator(s) used in their respective institutions including a description of the FMS aircraft models, software, hardware, and displays. Panelists will summarize previous, on-going, and planned empirical studies conducted with the simulators. Greg Pisanich will describe two NASA Ames simulation systems: the Stone Soup Simulator (SSS), and the Airspace Operations Human Factors Simulation Laboratory. The the Stone Soup Simulator is a desktop-based, research flight simulator that includes mode control, flight management, and datalink functionality. It has been developed as a non-proprietary simulator that can be easily distributed to academic and industry researchers who are collaborating on NASA research projects. It will be used and extended by research groups represented by at least two panelists (Mitchell and Palmer). The Airspace Operations Simulator supports the study of air traffic control in conjunction with the flight deck. This simulator will be used provide an environment in which many AATT and free flight concepts can be demonstrated and evaluated. Mike Palmer will describe two NASA Langley efforts: The Langley Simulator and MD-11 extensions to the NASA Amesbury simulator. The first simulator is publicly available and combines a B-737 model with a high fidelity flight management system. The second simulator enhances the S3 simulator with MD-11 electronic flight displays together with modifications to the flight and FMS models to emulate MD-11 dynamics and operations. Chris Mitchell will describe GT-EFIRT (Georgia Tech-Electronic Flight Instrument Research Tool) and B-757 enhancements to the NASA Ames S3. GT-EFIRT is a medium fidelity simulator used to conduct preliminary studies of the CATS (crew activity tracking system). Like the Langley efforts with S3, the Georgia Tech enhancements will allow it to emulate the dynamics and operations of a widely used glass cockpit. Amy Pritchett will describe the MIT simulator(s) that have been used in a range of research investigating cockpit displays, warning devices, and flight deck-ATC interaction.

Pisanich, Greg↗

Quantifying Errors in 3D CME Parameters Derived from Synthetic Data Using White-Light Reconstruction Techniques

Current efforts in space weather forecasting of CMEs have been focused on predicting their arrival time and magnetic structure. To make these predictions, methods have been developed to derive the true CME speed, size, position, and mass, among others. Difficulties in determining the input parameters for CME forecasting models arise from the lack of direct measurements of the coronal magnetic fields and uncertainties in estimating the CME 3D geometric and kinematic parameters after eruption. White-light coronagraph images are usually employed by a variety of CME reconstruction techniques that assume more or less complex geometries. This is the first study from our International Space Science Institute (ISSI) team “Understanding Our Capabilities in Observing and Modeling Coronal Mass Ejections”, in which we explore how subjectivity affects the 3D CME parameters that are obtained from the Graduated Cylindrical Shell (GCS) reconstruction technique, which is widely used in CME research. To be able to quantify such uncertainties, the “true” values that are being fitted should be known, which are impossible to derive from observational data. We have designed two different synthetic scenarios where the “true” geometric parameters are known in order to quantify such uncertainties for the first time. We explore this by using two sets of synthetic data: 1) Using the ray-tracing option from the GCS model software itself, and 2) Using 3D magnetohydrodynamic (MHD) simulation data from the Magnetohydrodynamic Algorithm outside a Sphere code. Our experiment includes different viewing configurations using single and multiple viewpoints. CME reconstructions using a single viewpoint had the largest errors and error ranges overall for both synthetic GCS and simulated MHD white-light data. As the number of viewpoints increased from one to two, the errors decreased by approximately 4° in latitude, 22° in longitude, 14° in tilt, and 10° in half-angle. Our results quantitatively show the critical need for at least two viewpoints to be able to reduce the uncertainty in deriving CME parameters. We did not find a significant decrease in errors when going from two to three viewpoints for our specific hypothetical three spacecraft scenario using synthetic GCS white-light data. As we expected, considering all configurations and numbers of viewpoints, the mean absolute errors in the measured CME parameters are generally significantly higher in the case of the simulated MHD white-light data compared to those from the synthetic white-light images generated by the GCS model. We found the following CME parameter error bars as a starting point for quantifying the minimum error in CME parameters from white-light reconstructions: Δθ (latitude)=6° +2° -3° , Δϕ (longitude)=11° +18° -6° , Δγ (tilt)=25° +8° -7° , Δx (half-angle)=10° +12° -6° , Δh (height)=0.6 +1.2 -0.4 R ⨀ , and Δκ (ratio)=0.1 +0.03 -0.02 .

Coronal mass ejections↗

Quantifying Errors in 3D CME Parameters Derived From Synthetic Data Using White-Light Reconstruction Techniques

Current efforts in space weather forecasting of CMEs have been focused on predicting their arrival time and magnetic structure. To make these predictions, methods have been developed to derive the true CME speed, size, position, and mass, among others. Difficulties in determining the input parameters for CME forecasting models arise from the lack of direct measurements of the coronal magnetic fields and uncertainties in estimating the CME 3D geometric and kinematic parameters after eruption. White-light coronagraph images are usually employed by a variety of CME reconstruction techniques that assume more or less complex geometries. This is the first study from our International Space Science Institute (ISSI) team “Understanding Our Capabilities in Observing and Modeling Coronal Mass Ejections”, in which we explore how subjectivity affects the 3D CME parameters that are obtained from the Graduated Cylindrical Shell (GCS) reconstruction technique, which is widely used in CME research. To be able to quantify such uncertainties, the “true” values that are being fitted should be known, which are impossible to derive from observational data. We have designed two different synthetic scenarios where the “true” geometric parameters are known in order to quantify such uncertainties for the first time. We explore this by using two sets of synthetic data: 1) Using the ray-tracing option from the GCS model software itself, and 2) Using 3D magnetohydrodynamic (MHD) simulation data from the Magnetohydrodynamic Algorithm outside a Sphere code. Our experiment includes different viewing configurations using single and multiple viewpoints. CME reconstructions using a single viewpoint had the largest errors and error ranges overall for both synthetic GCS and simulated MHD white-light data. As the number of viewpoints increased from one to two, the errors decreased by approximately 4° in latitude, 22° in longitude, 14° in tilt, and 10° in half-angle. Our results quantitatively show the critical need for at least two viewpoints to be able to reduce the uncertainty in deriving CME parameters. We did not find a significant decrease in errors when going from two to three viewpoints for our specific hypothetical three spacecraft scenario using synthetic GCS white-light data. As we expected, considering all configurations and numbers of viewpoints, the mean absolute errors in the measured CME parameters are generally significantly higher in the case of the simulated MHD white-light data compared to those from the synthetic white-light images generated by the GCS model. We found the following CME parameter error bars as a starting point for quantifying the minimum error in CME parameters from white-light reconstructions: Δθ (latitude)=6° +2° -3° , Δϕ (longitude)=11° +18° -6° , Δγ (tilt)=25° +8° -7° , Δx (half-angle)=10° +12° -6° , Δh (height)=0.6 +1.2 -0.4 R ⨀ , and Δκ (ratio)=0.1 +0.03 -0.02 .

Coronal mass ejections↗

Investigations into mirror fabrication metrology analysis

This final report describes the work performed under this delivery order from June 1993 through August 1994. The scope of work included three distinct tasks in support of the AXAF-I program. The objective of the first task was to perform investigations of the grinding and polishing characteristics of the zerodur material by fabricating several samples. The second task was to continue the development of the integrated optical performance modeling software for AXAF-I. The purpose of third and final task was to develop and update the database of AXAF technical documents for an easy and rapid access. The MSFC optical and metrology shops were relocated from the B-wing of Building 4487 to Room BC 144 of Building 4466 in the beginning of this contract. This included dismantling, packing, and moving the equipment from its old location, and then reassembling it at the new location. A total of 65 zerodur samples, measuring 1 inch x 2 inches x 6 inches were ground and polished to a surface figure of lambda/10 p-v, and a surface finish of 5A rms were fabricated for coating tests. A number of special purpose tools and metal mirrors were also fabricated to support various AXAF-I development activities. In the metrology area, the ZYGO Mark 4 interferometer was relocated and also upgraded with a faster and more powerful processor. Surface metrology work was also performed on the coating samples and other optics using ZYGO interferometer and WYKO profilometer. A number of new features have been added to the GRAZTRACE program to enhance its analysis and modeling capabilities. A number of new commands have been added to the command mode GRAZTRACE program to provide a better control to the user on the program execution and data manipulation. Some commands and parameter entries have been reorganized for a uniform format. The command mode version of the convolution program CONVOLVE has been developed. An on-line help system and a user's manual have also been developed for the benefit of the users. The database of AXAF technical documents continues to progress. The titles, company name, date, and location of over 390 documents have been entered in this database. This database provides both a data search and retrieval function, and a data adding function. These functions allow a user to quickly search the data files for documents or add new information. A detailed user's guide has also been prepared. This user guide includes a document classification guide, a list of abbreviations, and a list of acronyms, which have been used in compiling this database of AXAF-I technical documents.

Dimmock, John O.↗

Lander Lighting Solution to Reduce Pilot & Autonomous Approach Errors

The south pole lighting environment will have harsh low inclination sunlight, making overhead judgement of surface features difficult. Autonomous solutions are great, however, the need for visual monitoring and independent go/no-go decisions remain. Our project proposes that lunar landing systems will be better served by including a powerful landing light system that improves visibility of surfaces from overhead by illuminating the ground at critical distances for the crew to make last minute decisions regarding an approach. The project utilized computer-based optical modeling software to predict requirements for a potential landing light system. The analysis based the lamp prediction from commercially available LED chip sets and lamp optics. The goal was to illustrate a method to raise the surface contrast of a landing site within an acceptable contrast threshold for most camera systems and human observers to recognize hazards that would not be noticed with low inclination sunlight alone. The Apollo lunar landings benefitted from overhead sun or dark conditions. The surface lighting at the Lunar South Pole is a harsh environment where surfaces are lit from a low inclination angle by the sun (from the side). This change in lighting condition precipitates a need for updated lunar landing systems that facilitate improved recognition of landing sites, and thereby increase pilot awareness of landing hazards. The reliance on LIDAR and other autonomous mechanisms alone is risky given the known usage of visual monitoring for operator concurrence on current spacecraft programs and present-day autonomous land-based vehicles. Visual monitoring via cameras or windows requires the surface contrast to be within 3 orders of magnitude for reasonable recognition of objects. Artificial overhead illumination, when sufficiently sized, provides a means to even out contrast problems created by low inclination sunlight, potentially reducing piloting errors. Current vehicle requirements do not specify this type of guidance for the purpose of increasing mission success. An optical ray-trace simulation model was developed in Zemax Optics Studio to predict the best combination of LED power, LED optics, lamp quantity, and lamp location to raise the surface contrast to within 2 orders of magnitude from 3 orders required to further increased visibility and reduce risk. The project considered the following design constraints: potential base diameter of lander, approach distance(s) for a go-no-go decision point (200 meter), solar inclination angle (2-7), lunar surface reflectance (10%), LED chip sets, LED focusing optics, LED power, lamp quantity, lamp locations, and illumination diameter of lunar surface landing zone. The results can be used to establish minimum design constraints for vehicle landing light systems. With a solar inclination angle ranging from 2-7 degrees, the horizontal illumination of the lunar surface is attenuated by about 10% when compared to overhead illumination from the Sun. This modifies the sun's maximum of 130,000 lux to 13,000 lux horizontal illuminance. The artificial lighting system was designed to provide an 18-meter-wide illumination zone, to create viewing clearances around a 6-meter-wide lander. The system provides an average illuminance of 300 lux, meeting the 2 orders of magnitude criteria. The solution utilized modern Chip On Board LEDs, that each utilized 17 watts, with focusing Total Internal Reflection (TIR) lenses. A lighting system of 300 LEDs was arrayed along the "bottom" of a “lander”. With 17 watts per LED, the system is estimated to require 5100 watts. This is a large amount of power, but it would only be needed during critical phases during the landing. LED lighting systems can be dimmed, and it is assumed that as the lander arrives closer to the landing site, the lighting system power can be adjusted as needed to produce the necessary surface illuminance. The designed reduction of contrast improves reliability of safety assessments using real time visible light camera systems and out the window viewing by the crew.

T A Clark↗

Effects of Absorber Near-Interface Compensation on Cd(Se,Te) Solar Cell Performance

Arsenic (As) has been shown to be an effective p-type dopant for CdTe, although high performance in As-doped devices remains difficult to achieve. Arsenic is prone to self-compensation in CdTe, as evidenced by the accumulation of dopant atoms in CdTe/Cd(Se,Te) near the interface with MgxZn1-xO (MZO). In this study, we use SCAPS 1D modeling software to investigate the effect of near-interface compensation, helping elucidate loss pathways in present-day As-doped devices and informing future growth directions. We consider three possible results of As accumulation: shallow donors, deep recombination centers, and a thin layer of excess acceptor accumulation. The reduction in near-interface carrier concentration caused by shallow donors is shown to improve open-circuit voltages (VOC), whereas deep levels are detrimental to all performance parameters. The thin charge layer affects capacitance-voltage (CV) measurements by reducing the depletion width while maintaining the same carrier concentration, replicating CV behavior that has been observed in actual devices. These results illustrate the importance of monitoring dopant accumulation within 100 nm of the interface, and suggest that reducing or eliminating the As concentration in this region would be beneficial. An undoped Cd(Se,Te) layer at the interface is suggested as a possible device structure to boost performance.

absorber↗

Optimal Sizing and Dispatch of Solar Power with Storage

Designers of utility-scale solar plants with storage, seeking to maximize some aspect of plant performance, face multiple challenges. In many geographic locations, there is significant penetration of photovoltaic generation, which depresses energy prices during the hours of solar availability. An energy storage system affords the opportunity to dispatch during higher-priced time periods, but complicates plant design and dispatch decisions. Solar resource variability compounds these challenges, because determining optimal system sizes requires simultaneously considering how the plant will be operated under the imposed market and weather conditions. We develop an approach to analyze the economic performance of hybrid and single-technology solar power plants, which incorporates optimal dispatch, and considers the expected electricity market and weather conditions. We utilize the System Advisor Model software package to simulate the operation of multiple renewable generation and energy storage technologies, in conjunction with hourly-fidelity generation decisions determined by a revenue-maximizing, mixed-integer linear program. We show that, under our assumed market and weather conditions, the lifetime benefit-to-cost ratio can be improved by 6 to 19 percent, relative to a baseline design without optimizing, and that a concentrating solar power with thermal energy storage design produces significantly more energy per year, but is less profitable under our cost assumptions.

black-box optimization↗

Scientific computing plan for the ECCE detector at the Electron Ion Collider

The Electron Ion Collider (EIC) is the next generation of precision QCD facility to be built at Brookhaven National Laboratory in conjunction with Thomas Jefferson National Laboratory. There are a significant number of software and computing challenges that need to be overcome at the EIC. During the EIC detector proposal development period, the ECCE consortium began identifying and addressing these challenges in the process of producing a complete detector proposal based upon detailed detector and physics simulations. Here, in this document, the software and computing efforts to produce this proposal are discussed; furthermore, the computing and software model and resources required for the future of ECCE are described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Pyomo: Accidentally outrunning the bear

Pyomo is an open-source optimization modeling software that has undergone significant evolution since its inception in 2008. Pyomo has evolved to enhance flexibility, solver integration, and community engagement. Modern collaborative tools for open-source software have facilitated the development of new Pyomo functionality and improved our development process through automated testing and performance-tracking pipelines. However, Pyomo faces challenges typical of research software, including resource limitations and knowledge retention. The Pyomo team’s commitment to better development practices and community engagement reflects a proactive approach to these issues. We describe Pyomo’s development journey, highlighting both successes and failures, in the hopes that other open-source research software packages may benefit from our experiences.

automation↗

Micromagnetic Modeling of a Magnetically Unstable Zone and Its Geological Significances

Abstract Recent micromagnetic simulations have found that particles in the transitional zone between the single domain (SD) and single vortex (SV) zone are prone to thermal and external magnetic field instabilities that could adversely affect the accuracy of interpretations of paleomagnetic recordings. In this study, we attempt to evaluate the internal magnetization characteristics of these magnetically unstable (MU) particles and the influence on paleomagnetic observations by simulating the magnetic behavior of 68–104 nm truncated octahedral magnetite particles using the MERRILL modeling software. We found that: (a) The size region of the “MU zone” for grains of truncated octahedron shape is different from cubic octahedrons and spheres, indicating that the zone may be controlled by the geometry and shape of particles; (b) The MU zone has a range of 79–97 nm region, which is dominated by a hard‐axis aligned single vortex (HSV); and (c) MU particles are unstable as a function of temperature and have low coercive fields. Finally, the numerical fitting of hysteresis parameters for experimental data suggests that the influence of such MU particles in samples should not be ignored, especially for samples with fine‐grained magnetic minerals as the primary magnetic recording carriers. This research has extended our understanding of the behavior of the “MU zone” and its significance on paleomagnetic records.

Wang, Yuqin↗