SEARCH · Engineering Papers

Results for “Critical field”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2023

TROPHY: A Topologically Robust Physics-Informed Tracking Framework for Tropical Cyclones

Tropical cyclones (TCs) are among the most destructive weather systems. Realistically and efficiently detecting and tracking TCs are critical for assessing their impacts and risks. In particular, the eye is a signature feature of a mature TC. Therefore, knowing the eyes’ locations and movements is crucial for both operational weather forecasts and climate risk assessments. Recently, a multilevel robustness framework has been introduced to study the critical points of time-varying vector fields. The framework quantifies the robustness (i.e., structural stability) of critical points across varying neighborhoods. By relating the multilevel robustness with critical point tracking, the framework has demonstrated its potential in cyclone tracking. An advantage is that it identifies cyclonic features using only 2D wind vector fields, which is encouraging as most tracking algorithms require multiple dynamic and thermodynamic variables at different altitudes. A disadvantage is that the framework does not scale well computationally for datasets containing a large number of cyclones. Herein this paper introduces a topologically robust physics-informed tracking framework (TROPHY) for TC tracking. The main idea is to integrate physical knowledge of TC to drastically improve the computational efficiency of multilevel robustness framework for large-scale climate datasets. First, during preprocessing, we propose a physics-informed feature selection strategy to filter 90% of critical points that are short-lived and have low stability, thus preserving good candidates for TC tracking. Second, during in-processing, we impose constraints during the multilevel robustness computation to focus only on physics-informed neighborhoods of TCs. We apply TROPHY to 30 years of 2D wind fields from reanalysis data in ERA5 and generate a number of TC tracks. In comparison with the observed tracks, we demonstrate that TROPHY can capture TC characteristics (e.g., frequency, intensity, duration, latitudes with maximum intensity, and genesis) that are comparable to and sometimes even better than a well-validated TC tracking algorithm that requires multiple dynamic and thermodynamic scalar fields.

97 MATHEMATICS AND COMPUTING↗

NASA NTRSDate not supplied

Surface Hardness Testing of Mobile Launcher 1 (ML-1) Vehicle Support Posts (VSP) Following Artemis I Launch Exposure

Following the successful launch of Artemis I on November 16th, 2022 (Fig 1), a surface hardness assessment was performed on the eight primary VSPs (Figs. 2, 3) located on the Zero Deck of the ML. This assessment was performed as part of an instrument investigation for strain gauges installed on VSP interior surfaces that had a “strange reaction”. This report does not address the VSP strain gauge investigation, but rather VSP structural integrity for future service (Artemis II, III, …). Surface hardness measurements were taken to assess the strength of the ASTM A148 steel castings to confirm thermal exposure from launch had not reannealed the VSPs. However unlikely, it is possible high temperature may have heated and tempered (reannealed) the VSPs, resulting in a change of built-in stresses. This loss in built-in stress may explain the change in off-set measured by the strain gauges. Initial hardness measurements were made at two interior surfaces of VSPs near the strain gauges and the corresponding exterior surfaces (engine hole-facing). Field hardness measurements on the exterior, engine hole surface-facing surfaces for all eight VSPs ranged from 175 to 226 HB/P, well below the minimum hardness of 302 Brinell (HB). Opposing interior surface hardness measurements near strain gauges ranged from 264 to 358 HB/P. In all cases for flight VSPs, a drop in interior-to-exterior hardness ranging from 70 to 132 HB/P was measured. Surface hardness measurements were then made on the two spare VSPs not exposed to launch temperatures for baseline comparison. These two spare VSPs had a drop in Brinell hardness from interior to exterior surfaces ranging from 3 to 17 HB/P. The large interior-to-exterior delta hardness for eight flight VSPs but not for the two spare VSPs suggests launch exposure temperatures softened the exterior surfaces, raising two questions: 1) Are the VSPs compromised for future use, and 2) How deep is the exterior surface softening effect? To address concerns of low surface hardness values at non-critical exterior surfaces, additional field surface hardness measurements were taken at high-stress locations on VSP3: Both door jams, back radius, and internal left and right radius near the base. These values came in well above the required minimum hardness (317 to 365 HB/P). For a more representative pre-launch versus post-launch surface hardness comparison, VSP Acceptance Data Packages (ADPs) were reviewed, and hardness measurements were made at the same eight original locations on VSP3 where hardness measurements were made at the foundry. To determine the degree of softening, 1/8” of material was machined away at three exterior surfaces of engine hole-facing side (Fig. 4). A significant increase in surface hardness was measured ranging from 50 to 86 HB/P. Values at these three locations were still slightly below the required minimum hardness, so an additional 1/8” material was machined away, and hardness values were retaken. At the ¼” depth, hardness values increased to 321 to 324 HB/P. These values fall above the required minimum, concluding surface hardness testing. Rationale for VSP continued use was reached with no restrictions during a Chief Engineering Review based on hardness values exceeding minimum at critical locations and subsurface locations for VSP #3.

Artemis I↗

DOE OSTI2025

TARP Identifying Critical Thresholds for Acute Response of Plants and Ecosystems to Water Stress at Walker Branch Watershed, 2002-2005

This dataset contains data from a manipulative field study aimed at identifying critical thresholds for acute response of plants and ecosystems to water stress (TARP) that took place at Walker Branch Watershed in Oak Ridge, Tennessee from 2002-2005 (2002-06-20 to 2005-12-16). The study used understory tents for the removal of 100% of the growing-season throughfall and stem flow, to provide data on the impact of acute drought on mechanisms responsible for growth and mortality of deciduous forest canopy trees representative of common plant functional types (Liriodendron and Quercus). Through three years of manipulation (2003, 2004 and 2005; pretreatment measurements in 2002) various measures of tree response to surface moisture deficits were recorded including root and leaf traits, plant nonstructural carbohydrates status, hourly sapflow, basal area, and periodic observations of foliar photosynthesis and conductance. Additionally, environmental data such as air and soil temperature, soil water content, and soil matric potential were recorded. This dataset contains data 16 files in comma-separate (*.csv) format. Additional metadata are provided: 16 data dictionaries and a file-level metadata file in comma separate (*.csv) format and a user guide in PDF (*.pdf) format.

basal area↗

DOE OSTI2021

Effect of heat treatments on superconducting properties and connectivity in K-doped BaFe2As2

Abstract Fe-based superconductors and in particular K-doped BaFe 2 As 2 (K-Ba122) are materials of interest for possible future high-field applications. However the critical current density ( J c ) in polycrystalline Ba122 is still quite low and connectivity issues are suspected to be responsible. In this work we investigated the properties of high-purity, carefully processed, K-Ba122 samples synthesized with two separate heat treatments at various temperatures between 600 and 825 °C. We performed specific heat characterization and T c -distribution analysis up to 16 T and we compared them with magnetic T c and J c characterizations, and transmission-electron-microscopy (TEM) microstructures. We found no direct correlation between the magnetic T c and J c , whereas the specific heat T c -distributions did provide valuable insights. In fact the best J c -performing sample, heat treated first at 750 °C and then at 600 °C, has the peak of the T c -distributions at the highest temperatures and the least field sensitivity, thus maximizing H c2 . We also observed that the magnetic T c onset was always significantly lower than the specific heat T c : although we partially ascribe the lower magnetization T c to the small grain size (< λ , the penetration depth) of the K-Ba122 phase, this behaviour also implies the presence of some grain-boundary barriers to current flow. Comparing the T c -distribution with J c , our systematic synthesis study reveals that increasing the first heat treatment above 750 °C or the second one above 600 °C significantly compromises the connectivity and suppresses the vortex pinning properties. We conclude that high-purity precursors and clean processing are not yet enough to overcome all J c limitations. However, our study suggests that a higher temperature T c -distribution, a larger H c2 and a better connectivity could be achieved by lowering the second heat treatment temperature below 600 °C thus enhancing, as a consequence, J c .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Measurement-induced criticality in random quantum circuits

We investigate the critical behavior of the entanglement transition induced by projective measurements in (Haar) random unitary quantum circuits. Using a replica approach, we map the calculation of the entanglement entropies in such circuits onto a two-dimensional statistical-mechanics model. In this language, the area- to volume-law entanglement transition can be interpreted as an ordering transition in the statistical-mechanics model. We derive the general scaling properties of the entanglement entropies and mutual information near the transition using conformal invariance. We analyze in detail the limit of infinite on-site Hilbert space dimension in which the statistical-mechanics model maps onto percolation. In particular, we compute the exact value of the universal coefficient of the logarithm of subsystem size in the n th Rényi entropies for n ≥ 1 in this limit using relatively recent results for the conformal field theory describing the critical theory of two-dimensional (2D) percolation, and we discuss how to access the generic transition at finite on-site Hilbert space dimension from this limit, which is in a universality class different from 2D percolation. Here, we also comment on the relation to the entanglement transition in random tensor networks, studied previously in Vasseur et al. [Phys. Rev. B 100, 134203 (2019)].

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2010

Australian Soil Moisture Field Experiments in Support of Soil Moisture Satellite Observations

Large-scale field campaigns provide the critical fink between our understanding retrieval algorithms developed at the point scale, and algorithms suitable for satellite applications at vastly larger pixel scales. Retrievals of land parameters must deal with the substantial sub-pixel heterogeneity that is present in most regions. This is particularly the case for soil moisture remote sensing, because of the long microwave wavelengths (L-band) that are optimal. Yet, airborne L-band imagers have generally been large, heavy, and required heavy-lift aircraft resources that are expensive and difficult to schedule. Indeed, US soil moisture campaigns, have been constrained by these factors, and European campaigns have used non-imagers due to instrument and aircraft size constraints. Despite these factors, these campaigns established that large-scale soil moisture remote sensing was possible, laying the groundwork for satellite missions. Starting in 2005, a series of airborne field campaigns have been conducted in Australia: to improve our understanding of soil moisture remote sensing at large scales over heterogeneous areas. These field data have been used to test and refine retrieval algorithms for soil moisture satellite missions, and most recently with the launch of the European Space Agency's Soil Moisture Ocean Salinity (SMOS) mission, to provide validation measurements over a multi-pixel area. The campaigns to date have included a preparatory campaign in 2005, two National Airborne Field Experiments (NAFE), (2005 and 2006), two campaigns to the Simpson Desert (2008 and 2009), and one Australian Airborne Cal/val Experiment for SMOS (AACES), just concluded in the austral spring of 2010. The primary airborne sensor for each campaign has been the Polarimetric L-band Microwave Radiometer (PLMR), a 6-beam pushbroom imager that is small enough to be compatible with light aircraft, greatly facilitating the execution of the series of campaigns, and a key to their success. An L-band imaging radar is being added to the complement to provide simultaneous active-passive L-band observations, for algorithm development activities in support of NASA's upcoming Soil Moisture Active Passive (.S"M) mission. This paper will describe the campaigns, their objectives, their datasets, and some of the unique advantages of working with small/light sensors and aircraft. We will also review the main scientific findings, including improvements to the SMOS retrieval algorithm enabled by NAFE observations and the evaluation of the Simpson Desert as a calibration target for L-band satellite missions. Plans for upcoming campaigns will also be discussed.

Kim, Edward↗

DOE OSTI2026

Deriving effective electrode–ion interactions from free-energy profiles at electrochemical interfaces

Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.

Roncoroni, Fabrice [Lawrence Berkeley National Lab↗

NASA NTRS1999

The Effect of a Rotating Magnetic Field on Rayleigh-Benard Convection

A destabilizing vertical temperature gradient and a rotating magnetic field have been applied to a cylindrical column of liquid gallium. The convective flows which arise as a function of these driving parameters are identified. For small magnetic field strengths, a region of stationary flow is observed. This region is bounded by critical values of the Rayleigh and Hartmann numbers. As the rotating magnetic field is increased, the critical Rayleigh number can increase by more than a factor of 10. The rotating magnetic field itself induces an instability at a critical value of the Hartmann number independent of the Rayleigh number. The nature of the bifurcations (whether subcritical or supercritical) and the convective flows occurring at the critical Rayleigh numbers are dependent upon the Hartmann number. In the stationary flow region, the rotating magnetic field can induce fluid motion with velocities on the order of cm/s. Thus, a rotating magnetic field can be used to achieve the benefits of forced convection during crystal growth processes without triggering deleterious instabilities.

Volz, Martin P.↗

NASA NTRS2008

Proceedings of the Sixth NASA Langley Formal Methods (LFM) Workshop

Today's verification techniques are hard-pressed to scale with the ever-increasing complexity of safety critical systems. Within the field of aeronautics alone, we find the need for verification of algorithms for separation assurance, air traffic control, auto-pilot, Unmanned Aerial Vehicles (UAVs), adaptive avionics, automated decision authority, and much more. Recent advances in formal methods have made verifying more of these problems realistic. Thus we need to continually re-assess what we can solve now and identify the next barriers to overcome. Only through an exchange of ideas between theoreticians and practitioners from academia to industry can we extend formal methods for the verification of ever more challenging problem domains. This volume contains the extended abstracts of the talks presented at LFM 2008: The Sixth NASA Langley Formal Methods Workshop held on April 30 - May 2, 2008 in Newport News, Virginia, USA. The topics of interest that were listed in the call for abstracts were: advances in formal verification techniques; formal models of distributed computing; planning and scheduling; automated air traffic management; fault tolerance; hybrid systems/hybrid automata; embedded systems; safety critical applications; safety cases; accident/safety analysis.

Rozier, Kristin Yvonne↗

DOE OSTI2025

Singly occupied 4⁢𝑓 antiferromagnetic insulators: CePO 4 and CeVO 4

Rare-earth containing wide band gap oxides, which provide spin-photon interface and narrow linewidth optical emission, are getting significant attention as the most promising candidate materials in advancing quantum transduction and memories. Here, in this work, from ab initio calculations, we identify antiferromagnetic ground states in structurally preferred monoclinic CePO 4 and tetragonal CeVO 4 exhibiting localized occupied and unoccupied Ce 4⁢𝑓 states with 4⁢𝑓−4⁢𝑓 transition characteristics. Interestingly, in CePO 4 , O 2⁢𝑝 and P 3⁢𝑝 states hybridize negligibly with Ce 4⁢𝑓 states, while in CeVO 4 , V 3⁢𝑑 and O 2⁢𝑝 states hybridize and appear as extended states in between the occupied and unoccupied Ce 4⁢𝑓 states. Here, phonon calculations and analysis identify and differentiate Raman active phonon modes along with the spin phonon coupling of Ce in both CePO 4 and CeVO 4 that ultimately lead to different 4⁢𝑓 ground-state crystal-field multiplets, which are critical to accurately describe electronic transitions for foundational quantum transduction and memories. Further, the identified 𝐶 1 site symmetry of Ce, lacking inversion symmetry in CePO 4 , is relevant for quantum memories and 𝐷 2⁢𝑑 site symmetry of Ce exhibiting inversion symmetry in CeVO 4 is relevant for quantum transduction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Impact of CH 4 addition on the electron properties and electric field dynamics in a Ar nanosecond-pulsed dielectric barrier discharge

Non-equilibrium plasmas derive their low temperature reactivity from producing and driving energetic electrons and active species under large electric fields. Therefore, the impact of reactants on the plasma properties including electron number density, electric field, and electron temperature is critical for applications such as plasma CH 4 reforming. Due to experimental complexity, electron properties and the electric field are rarely measured together in the same discharge. Here, we combine time-resolved Thomson scattering and electric field induced second harmonic generation (EFISH) to probe electron temperature, electron density, and electric field strength in a 60 Torr CH 4 /Ar nanosecond-pulsed dielectric barrier discharge (ns-DBD) while varying the CH 4 mole fraction from 0 to 8%. These measurements are compared to a 1-D numerical model to benchmark its predictions and identify areas of uncertainty. Nonlinear coupling between CH 4 addition, electron temperature, electron density, and the electric field was directly observed. Contrary to previous measurements in He, the electron temperature increased with CH 4 mole fraction. This rise in electron temperature is identified as electron heating by residual electric fields that increased with larger CH 4 mole fraction. Moreover, the electron number density has been found to decrease rapidly with the increase of methane mole fraction. Comparison of these measurements with the model yielded better agreement at higher CH 4 mole fractions and with the usage of ab initio calculated Ar electron-impact cross-sections from the B-spline R-matrix (BSR) database. Furthermore, the calculated plasma properties are shown to be sensitive to the residual surface charge implanted on the quartz dielectric surfaces. Without considering surface charge in the simulations, the calculated electric field profiles agreed well with the measurements, but the electron properties were underpredicted by more than a factor of three. Therefore, measurements of either the electric field or electron properties measurements alone are insufficient to fully validate modelling predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1989

The topology and vorticity dynamics of a three-dimensional plane compressible wake

The three-dimensional aspects of transition in a low Mach number plane compressible wake are studied numerically. Comparisons are made between the topology of the velocity field and the vorticity dynamics of the flow based on results from direct numerical simulations of the full compressible Navier-Stokes equations. The velocity field is integrated to obtain instantaneous streamlines at different stages in the evolution. A generalized three-dimensional critical point theory is applied to classify the critical points of the velocity field.

Chen, Jacqueline H.↗

DOE OSTI2020

Unusual zwitterionic catalytic site of SARS–CoV-2 main protease revealed by neutron crystallography

The main protease (3CL M pro ) from SARS–CoV-2, the etiological agent of COVID-19, is an essential enzyme for viral replication. 3CL M pro possesses an unusual catalytic dyad composed of Cys 145 and His 41 residues. A critical question in the field has been what the protonation states of the ionizable residues in the substrate-binding active-site cavity are; resolving this point would help understand the catalytic details of the enzyme and inform rational drug development against this pernicious virus. Here, we present the room-temperature neutron structure of 3CL M pro , which allowed direct determination of hydrogen atom positions and, hence, protonation states in the protease. We observe that the catalytic site natively adopts a zwitterionic reactive form in which Cys 145 is in the negatively charged thiolate state and His 41 is doubly protonated and positively charged, instead of the neutral unreactive state usually envisaged. The neutron structure also identified the protonation states, and thus electrical charges, of all other amino acid residues and revealed intricate hydrogen-bonding networks in the active-site cavity and at the dimer interface. The fine atomic details present in this structure were made possible by the unique scattering properties of the neutron, which is an ideal probe for locating hydrogen positions and experimentally determining protonation states at near-physiological temperature. Our observations provide critical information for structure-assisted and computational drug design, allowing precise tailoring of inhibitors to the enzyme's electrostatic environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Development of radiopharmaceuticals for targeted alpha therapy: Where do we stand?

Targeted alpha therapy is an oncological treatment, where cytotoxic doses of alpha radiation are locally delivered to tumor cells, while the surrounding healthy tissue is minimally affected. This therapeutic strategy relies on radiopharmaceuticals made of medically relevant radionuclides chelated by ligands, and conjugated to targeting vectors, which promote the drug accumulation in tumor sites. This review discusses the state-of-the-art in the development of radiopharmaceuticals for targeted alpha therapy, breaking down their key structural components, such as radioisotope, targeting vector, and delivery formulation, and analyzing their pros and cons. Moreover, we discuss current drawbacks that are holding back targeted alpha therapy in the clinic, and identify ongoing strategies in field to overcome those issues, including radioisotope encapsulation in nanoformulations to prevent the release of the daughters. Lastly, we critically discuss potential opportunities the field holds, which may contribute to targeted alpha therapy becoming a gold standard treatment in oncology in the future.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2021

Transfer-Learnt Energy Models for Predicting Electricity Consumption in Buildings with Limited and Sparse Field Data

Modeling energy consumption is critical for energy-efficient utilization of the electric appliances in a building, smart grid programs (like demand-response), and many other smart home applications. State-of-the-art energy modeling techniques either rely on theoretical models, or extensive instrumentation of the building envelope to gather ``big" data to train a deep neural network. While theoretical models are often limited by their estimation accuracy, it is not always feasible to gather a significant amount of field data. In this paper, we explore transfer learning-based strategies to train much more accurate model for energy estimation when using a sparse field data. We transferred knowledge, in the form of data and parameters, from the simulation framework to the field data. We evaluated the efficacy of our approach on field data collected from six commercial buildings and our results indicate that transfer learning-based models trained over one month data can perform comparative (and in some cases better) than the state-of-the-art machine learning and deep learning solutions.

Jain, Milan↗

NASA NTRS2023

Effects of Optimisation Parameters on Data-Driven Magnetofrictional Modelling of Active Regions

Context . The solar magnetic field plays an essential role in the formation, evolution, and dynamics of large-scale eruptive structures in the corona. The estimation of the coronal magnetic field, the ultimate driver of space weather, particularly in the ‘low’ and ‘middle’ corona, is presently limited due to practical difficulties. Data-driven time-dependent magnetofrictional modelling (TMFM) of active region magnetic fields has been proven to be a useful tool to study the corona. The input to the model is the photospheric electric field that is inverted from a time series of the photospheric magnetic field. Constraining the complete electric field, that is, including the non-inductive component, is critical for capturing the eruption dynamics. We present a detailed study of the effects of optimisation of the non-inductive electric field on the TMFM of AR 12473. Aims . We aim to study the effects of varying the non-inductive electric field on the data-driven coronal simulations, for two alternative parametrisations. By varying parameters controlling the strength of the non-inductive electric field, we wish to explore the changes in flux rope formation and their early evolution and other parameters, for instance, axial flux and magnetic field magnitude. Methods . We used the high temporal and spatial resolution cadence vector magnetograms from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). The non-inductive electric field component in the photosphere is critical for energising and introducing twist to the coronal magnetic field, thereby allowing unstable configurations to be formed. We estimated this component using an approach based on optimising the injection of magnetic energy. Results . Our data show that flux ropes are formed in all of the simulations except for those with the lower values of these optimised parameters. However, the flux rope formation, evolution and eruption time varies depending on the values of the optimisation parameters. The flux rope is formed and has overall similar evolution and properties with a large range of non-inductive electric fields needed to determine the non-inductive electric field component that is critical for energising and introducing twist to the coronal magnetic field. Conclusions . This study shows that irrespective of non-inductive electric field values, flux ropes are formed and erupted, which indicates that data-driven TMFM can be used to estimate flux rope properties early in their evolution without needing to employ a lengthy optimisation process.

A. Kumari↗

DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗

DOE OSTI2023

Large deformation and brittle failure calculated using the dual-domain material point method

The dual domain material point (DDMP) method is explored as a candidate to be implemented in a general purpose code to perform simulations of materials with complex geometry that undergo large history-dependent deformation and failure. To test its candidacy, we study its mesh convergence, its sensitivity to mesh orientation, and its ability to handle softening and failure of a material. Simulations of large deformation and simulations of mechanical failure are performed using both DDMP and the material point method (MPM). When cell-crossing of material points is not an issue and when there are a sufficient number of material points in each computation cell, the numerical error decreases with the square of the cell size as expected for both MPM and DDMP. DDMP has reduced error compared with MPM when there are many instances of material points crossing cell boundaries due to the continuous nature of the modified gradient of the shape functions. Simulations of a specimen under tension are also performed where the background mesh is aligned and misaligned with the tension direction. MPM displays a significant mesh-dependent stress field, DDMP shows negligible mesh dependency. Despite a mesh orientation-dependent stress field from MPM, the critical tension and failure mode from both MPM and DDMP calculations have negligible mesh dependency when using a non-local failure model. If only the failure mode is important (i.e., local stresses are unimportant), MPM with a non-local failure model is a suitable method for modeling failure with small deformations. However, if local stresses are also important or if there are large deformations with many cell-crossings before failure, DDMP should be the method that is used. A needed improvement for DDMP is identified from our numerical simulations.

36 MATERIALS SCIENCE↗