Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pointing Performance Analysis”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Impact of Dopants on Suppressing Gas Evolution from High-Nickel Layered Oxide Cathodes

Gas release from high-Ni layered-oxide cathodes (LiNi x Mn 1-x-y-z Co y Al z O 2 ; x > 0.8) can jeopardize the overall performance and safety characteristics of the cell. A comprehensive assessment of rational cathode design with common dopants, such as Ni, Co, Al, and Mn, to suppress gas evolution is crucial for battery safety, yet it remains to be conducted. Furthermore, we present an in-situ gas analysis on nine high-Ni cathode materials with online electrochemical mass spectrometry (OEMS). We show that regardless of the dopant, reactive oxygen release from the cathode lattice remains a critical process for gas evolution. A series of comparisons reveals that the intensity and onset point of gas release is strongly dependent on the cathode composition. Notably, Al and Mn are the most effective dopants to suppress gas evolution from the cathode at 4.4 V. We further highlight lattice stability limits, across these nine compositions, between 85 % to 93 % state of charge.

25 ENERGY STORAGE↗
DOE OSTI2020

Integrated and Advanced Lighting Controls: The Next Frontier in Energy Savings?

Lighting energy use has significantly declined as a portion of overall building energy usage due largely to the uptake of more efficient light-emitting diode technology. Although the per unit cost of lighting controls has declined, the installation and use of lighting controls are low and the installation and use of advanced controls are abysmally low. Increased use of advanced and connected lighting controls, and the ability to integrate lighting with other building systems (e.g., heat and air, plug loads, miscellaneous loads) that make the building “smarter,” offer the potential for additional energy savings. Studies into the performance of new, advanced lighting control technologies in commercial buildings have generated new findings. Although each field evaluation has interesting results individually, an analysis of their collective findings across a spectrum of commercial building types will help establish the current state of the technology and implications related to market readiness. This paper will analyze results from available field evaluation studies (recently completed as well as ongoing) of retrofit sites and share findings and insights including: • The amount of energy savings that may be expected from these systems. • Cost effectiveness, monetization, and cost points of these systems. • Ideal characteristics of these systems the most viable. • Lessons learned from the installation, commissioning, and occupant response of these systems. The results will support Department of Energy research investments and help inform efficiency programs of the potential of these emerging technologies and how they may fit within the clean energy future.

Myer, Michael↗
DOE OSTI2021

Leveraging In-Network Computing and Programmable Switches for Streaming Analysis of Scientific Data

With the emergence of programmable network devices that match the performance of fixed function devices, several recent projects have explored in-network computing, where the processing that is traditionally done outside the network is offloaded to the network devices. In-network computing has typically been applied to network functions (e.g., load balancing, NAT, and DNS), caching, data reduction/aggregation, and coordination/consensus functions. In some cases it has been used to accelerate stream-processing tasks that involve small payloads and simple operations. In this work we focus on leveraging in-network computing for stream processing of scientific datasets with large payloads that require complex operations such as floating-point computations and logarithmic functions. We demonstrate in-network computing for a real-world scientific application performing streaming normalization of a 2-D image from a light source experiment. We discuss the challenges we encountered and potential approaches to address them.

Sankaran, Ganesh↗
DOE OSTI2025

Probing Operando Electrochemical Strain Generation in α-NaFeO 2 Composite Cathodes during Cycling of Na-Ion Batteries

The transition metal oxide (TMO) cathodes in Na-ion batteries suffer from low-capacity retention. Chemo-mechanical instabilities lead to the deterioration of the electrochemical performance of TMO cathodes in Li-ion batteries. However, there is not much known about the chemo-mechanical instabilities in the TMO cathodes for Na-ion batteries. Understanding the governing forces behind the interplay between the electrochemical performance and mechanical stability in TMO cathodes is critical for the development of Na-ion batteries. Here, we synchronize the digital image correlation (DIC) technique with electrochemical analysis to capture the real-time deformation behavior of the α-NaFeO 2 cathodes during cycling. When the charge cutoff voltage is 3.6 V, the cathode experiences reversible deformations (except for the first cycle). There is negative strain (shrinkage) generation during Na extraction and positive strain (expansion) generation during the subsequent Na insertion. A detailed analysis of the potential-dependent strain rate evolution points out complicated phase transformations and nonequilibrium conditions in the α-NaFeO 2 cathodes during cycling. When the charge cutoff voltage was increased to 4.2 V, there was a rapid capacity loss and large plastic deformations in the α-NaFeO 2 cathodes. We provide an in-depth discussion about the possible mechanisms behind the chemo-mechanical instabilities in the α-NaFeO 2 . In conclusion, the correlation is critical to develop material-based strategies to mitigate instability mechanisms in TMO cathodes for Na-ion batteries.

Wable, Minal [University of Maryland Baltimore Cou↗
DOE OSTI2021

Alignment of magnetic particles in anisotropic Nd–Fe–B bonded magnets

High-performance anisotropic bonded magnets need a high volumetric filling fraction of magnetic powder in a polymer binder and good alignment while possessing adequate mechanical properties. The correlation between the degree of alignment (DOA) and the volumetric filling fraction of magnetic powder, the type of binder and the processing parameters have been studied using Nd–Fe–B anisotropic bonded magnets. Here, thermomagnetic analysis and differential scanning calorimetry results show that the magnetic powder can be well aligned in a bonded magnet under an external magnetic field at a temperature higher than the binders' melting point by 20–60 K, depending on the type of magnetic powder and binder. The DOA results mainly from the interplay between the Zeeman energy (E app ) and the inter-particle static energy (E stat ) in bonded magnets at the chosen alignment temperature and the magnetic fields. A good alignment can be achieved by an alignment magnetic field which is about twice its coercivity in an anisotropic Nd–Fe–B bonded magnet, which is confirmed by our experimental and modeling results.

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

Successive magnetic orderings in the Ising spin chain magnet DyNi 5 Ge 3

In this report, we investigated a new rare-earth-based one-dimensional Ising spin chain magnet DyNi 5 Ge 3 by means of magnetization, specific heat, and powder neutron diffraction measurements. Due to the crystalline electrical field splitting, the magnetic Dy ions share an Ising-like ground doublet state. Owning to the local point symmetry, these Ising moments form into two canted magnetic sublattices, which were further confirmed by the angle-dependent magnetization measurement. In zero fields, two successive antiferromagnetic phase transitions were found at temperatures T N1 =6K and T N2 =5K, respectively. Only part of the moments are statically ordered in this intermediate state between T N1 and T N2 . Powder neutron diffraction experiments at different temperatures were performed as well. Further, an incommensurate magnetic propagation vector of k m =(0.5,0.4,0.5) was identified. The refined spin configurations through the irreducible representation analysis confirmed that these Ising spins are canted in the crystal ab plane.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spatiotemporal Analyses of Groundwater and Shoreline Cr(VI) Concentrations in the 100 Areas at Hanford

Cleanup efforts have been ongoing since the late 1990s to remediate contaminated waste sites and groundwater in the 100 Areas at the U.S. Department of Energy (DOE) Hanford Site. One of the primary contaminants of concern is hexavalent chromium (Cr(VI)), which was used as a corrosion inhibitor in cooling water for nuclear reactors that formerly operated along the shoreline of the Columbia River. Cleanup efforts have included 1) removal, treatment (as needed), and disposal of contaminated sediments; 2) in situ redox manipulation as a permeable reactive barrier; 3) pump-and-treat; 4) soil flushing; and 5) monitored natural attenuation. DOE’s annual groundwater monitoring reports document the significant reductions in Cr(VI) plume areas that have occurred over the past 10 years or more as a result of these cleanup efforts. The Record of Decision for the 100-HR-3 operable unit specified a cleanup level (CUL) for Cr(VI) in groundwater of 48 µg/L to protect human receptors, and a surface water CUL of 10 µg/L to protect aquatic organisms in the Columbia River. The Record of Decision did not specify point-of-compliance locations for the surface water CUL. Data for 2019 from the six groundwater operable units (OUs) in the 100 Areas indicate that the 48 μg/L groundwater CUL has been achieved in 100% of the wells in the 100-BC and 100-NR OUs, and in 89- 97% of the wells in the other OUs (100-KR, 100-HR-D, 100-HR-H, 100-FR). Data for 2019 indicate that 100% of the aquifer tubes monitored for Cr(VI) in the 100 Areas have concentrations below the 48 μg/L groundwater CUL. However, the 10 μg/L standard has not yet been consistently achieved for both inland groundwater monitoring wells and shoreline aquifer tubes. This report describes a series of data analyses performed to identify consistent relationships, if any, between inland well and shoreline Cr(VI) concentrations within the 100 Areas. To this end, select monitoring data for Cr(VI) measured in groundwater and aquifer tubes at the 100 Areas were analyzed for a 10-year period—2010 to 2019. Relationships between inland groundwater plumes and surface-water points of discharge in and along the Columbia River were examined through several analyses that included inland well and aquifer tube concentrations as a function of distance from the shoreline, evaluation of cumulative probability plots, trend analysis, correlation analysis, cluster analysis, and identification of plume trajectories for each of the 100 Areas. The analyses did not identify consistent relationships between inland groundwater Cr(VI) concentrations and shoreline concentrations within the 100 Areas due to several confounding factors influencing groundwater flow directions and Cr(VI) concentrations. The proximity of groundwater Cr(VI) plumes to the river, and the highly dynamic nature of the river, influence the transport behavior of the plumes and create challenges for quantifying attenuation of Cr(VI) between the inland monitoring wells and shoreline concentrations. Other factors contributing to temporal and spatial Cr(VI) concentrations, as supported by some of the data analyses, include the presence of vadose zone sources, variable sorption behavior, and complexities associated with Cr(VI) mass transfer between the upper and lower aquifers and their interactions with the river. Hence, monitoring to assess compliance with target CULs will need to be determined for each area individually since several factors influence Cr(VI) concentrations in the 100 Areas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Aero-Thermal Characterization of Accelerating and Diffusing Passages Downstream of Rotating Detonation Combustors

Cycle benefits of rotating detonation engines show up to five percentage points of efficiency gain for low-pressure ratio engines. An optimal integration between the combustor and the turbine needs to be guaranteed to realize this potential gain. The rotating detonation combustor (RDC) exhausts transonic flow with shocks rotating at frequencies ranging from a few to tens of kilohertz depending on the number of present waves. Hence, the turbine design requires precise knowledge of the fluctuations and losses downstream of the combustor. This paper focuses on the quantification of fluctuations and losses for accelerating and diffusing passages. The analysis of the combustor is performed via reactive unsteady Reynolds Averaged Navier-Stokes (URANS) simulations. The unsteady RANS equations are solved via CFD++ from Metacomp with a one-step reaction mechanism for an H2-air mixture. The resolving of the boundary layer is achieved with a structured mesh of around 36 million cells. Inlet pressure of 10 bar and two different back pressures are applied to the combustor to model the interconnection with downstream turbines. Finally, we present and assess a methodology to reduce the computational time to model these passages ten times.

Braun, James↗
DOE OSTI2025

Techno-economic comparison of sCO 2 cycles for particle-based CSP at design-point conditions

In this work, we compare the techno-economic performance of supercritical carbon dioxide power cycles integrated in a particle CSP system. We model four core cycle configurations: simple (with optional bypass), recompression (with optional bypass), partial cooling, and turbine split flow, which each demonstrate different benefits in a CSP system, such as high efficiency, low cost, or large HTF temperature differences. We parametrically sweep cycle design variables for each configuration. The set of power cycle performance results are then combined with a design point particle CSP system model which calculates the system specific cost. The simple cycle and turbine split flow cycles have the best performance in the baseline results, with system specific costs of 5,912 and 5,899 $\$$/kWe respectively. In addition to the baseline set of results, we also vary key parameters and costs in a sensitivity study. The cycle designs with the best system performance limit their efficiency to ~45 %, despite demonstrating higher maximum efficiencies, due to the rapid increase in cost of recuperation as efficiencies rise. The simple cycle has strong performance in the analysis and is on average only 1.4 % worse than the optimal configuration. Lowered turbine inlet temperatures from the sensitivity study improve performance by reducing the PHX and turbine cost. Decreasing the inlet temperature from 700 to 625°C results in an >8 % decrease in system specific cost. Future work should expand sensitivity analyses to colder turbine inlet temperatures and calculate system performance by simulating annual performance with off-design solar and cycle component models.

14 SOLAR ENERGY↗
DOE OSTI2026

Hard X-ray imaging with sub-10 nm resolution by scanning bonded 2D multilayer Laue lenses

Multilayer Laue lenses (MLLs) offer significant advantages over traditional diffractive focusing optics, such as zone plates (ZPs), in terms of efficiency and durability when focusing hard X-rays to 10 nm and below. Typically, a pair of 1D MLLs is aligned independently to achieve point focusing. This approach only permits scanning microscopy by moving the sample stage, which presents a limitation when it comes to scanning large or bulky objects at high speed. In this study, we present an experimental demonstration of scanning monolithically bonded 2D MLLs while keeping the sample stationary, effectively addressing the limitation of the sample’s size and weight. We characterized the system’s positional and angular stabilities during scanning, which are critical parameters to the nanofocusing performance of MLLs. In addition, we provided the X-ray measurement results and corresponding image resolution analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of undoped substrates on high performance silicon solar cells

Today's highest-efficiency silicon solar cells typically operate near the threshold between low-level and high-level injection. It is not well understood if pushing further into a regime in which the cell operating point is solidly in high-level injection at all times of the day has further benefits for the solar cell performance. From a reliability perspective, cells fabricated on lower doped silicon have a larger breakdown voltage. This advantage can affect the design of modules allowing higher voltages and a relaxation of the number of bypass diodes needed. In this project, we present a comprehensive assessment, both experimental and using simulation, of how bulk resistivity, light intensity, and operation temperature impact the performance of the silicon solar cell. This work incorporates a comprehensive device physics analysis assisted by numerical simulation. The source code is now available under General Public License (GPL3), and to further leverage the findings of this project and outreach people outside of the scientific community, we are working with www.pveducation.org (which receives >1 million visitors a year) to create interactive content using our simulation results. The simulation results indicate that high bulk resistivity wafers (>>10 Ωcm) require bulk Shockley-Read-Hall (SRH) lifetimes in the millisecond range to outperform wafers with standard bulk resistivities (<10 Ωcm). Additionally, above bulk resistivities of 10 Ωcm (the exact value depends on the bulk characteristics of the wafer), the cell efficiency is weakly dependent on the bulk resistivity. As a result, ingot manufactures may have an opportunity to further reduce wafer cost by growing higher resistivity ingots that are more tolerant to resistivity variations. This project is particularly relevant today, as solar cell architectures with improved surface passivation and milliseconds lifetimes wafers are commercially available, leveraging potential benefits of using higher bulk resistivities. Outside of interdigitated back contact (IBC) cell, reported studies on high resistivity silicon (>100 Ωcm) are very limited. To the best of our knowledge, this project provides for the first-time experimental insight on solar cells fabricated on wafers with bulk resistivities up to several thousand Ωcm, delivering a comprehensive vision of their performance under real-world temperature and light intensity operation conditions. We manufactured and characterized solar cells with bulk resistivities in the range of 1 Ωcm to >15k Ωcm. Under standard testing conditions (STC), we measured solar cells efficiencies over 20% over the entire range of bulk resistivities, using our baseline cell processing. To evaluate the cell performance in real-world operation conditions, the solar cells were measured at different temperatures (25-80°C) and at different light intensities (0.1-1 suns). The measurements show that the bulk resistivity does not impact the solar cell response to temperature and light intensity. Similar thermal coefficients (TC) were measured for standard and high bulk resistivities, and they are comparable with the TC values reported in the literature for standard bulk resistivities <10 Ωcm. After light soaking, the solar cell didn’t show signs of light-induced degradation (LID). This result was expected since n-type float zone (FZ) wafers were used in this work, i.e. low traces of boron and low concentration of oxygen (oxygen is typically found in the seed end of Czochralski (CZ) ingots). We measured for high bulk resistivities (>10 Ωcm) extremely high breakdown voltages (>1000V). In conclusion, the insight provided by this project can positively impact the levelized cost of energy (LCOE) of the photovoltaic systems through its effect on cell and ingot manufacturing yield, silicon cell power output, and module reliability.

14 SOLAR ENERGY↗
DOE OSTI2020

Impedance Analysis and PHIL Demonstration of Reactive Power Oscillations in a Wind Power Plant Using a 4-MW Wind Turbine

This paper presents impedance-based analysis, mitigation, and power-hardware-in-the-loop (PHIL) demonstration of reactive power oscillations in a wind power plant using a 4-MW Type III wind turbine drivetrain. Because such low-frequency oscillations result from interactions among slower control loops of wind turbines regulating phasor quantities—active and reactive power output of the wind turbine and the magnitude of voltages at the point of interconnection (POI)—a new type of admittance is defined in terms of phasor quantities for their analysis. The so-called power-domain admittance of a wind turbine is defined as the transfer function from the frequency and magnitude of voltages at the POI to the active and reactive power output of the turbine. The power-domain admittance responses of the 4-MW wind turbine are measured using a 7-MVA grid simulator to identify the source of the reactive power oscillations. Power-domain impedance analysis and PHIL experiments are performed to explain how a resonant mode manifests as turbine-to-turbine and plant-to-grid reactive power oscillations. It is discovered that weaker grids exhibiting high inductive impedance mitigate oscillations in the reactive power output of wind power plants; however, a higher grid impedance does not help in damping turbine-to-turbine reactive power oscillations. This paper presents a simple droop-based solution to eliminate both turbine-to-turbine and plant-to-grid reactive power oscillations.

17 WIND ENERGY↗
DOE OSTI2025

A 28 nm multiply-accumulate ASIC architecture for on-chip data compression in MHz frame rate X-ray and electron pixel detectors

Modern X-ray detector systems urgently require compact, efficient, and fast data compression schemes to handle the transmission of big data from pixel arrays, enabling frame rates in the MHz regime. Here, in this work, a data compression ASIC that implements a streaming fixed-length lossy compression scheme is introduced and analyzed, proving the feasibility and benefits of on-chip compression. The compression scheme utilizes a vector matrix product logic, which performs a number of floating-point multiplications, additions, and accumulations. The logic is verified, synthesized, and shown to fit in the area resource available for the X-ray detector under study, which comprises 192 × 168 pixels each of 12-bit width, and having a total area of 20 mm× 20 mm, about 2 mm× 20 mm of which are available for the digital logic. Several system architectures, precisions, and compression ratios ranging from 100 to 250 were analyzed to pave the way for on-chip fixed-length compression (e.g., principal component analysis, singular value decomposition) and data reduction (e.g., azimuthal integration) for X-ray and electron detectors.

Data compression↗
DOE OSTI2021

The Tianlai dish pathfinder array: design, operation, and performance of a prototype transit radio interferometer

ABSTRACT The Tianlai Dish Pathfinder Array is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization universe as a means for measuring large-scale cosmic structure. It performs drift scans of the sky at constant declination. We describe the design, calibration, noise level, and stability of this instrument based on the analysis of about 5% of 6200 h of on-sky observations through 2019 October. Beam pattern determinations using drones and the transit of bright sources are in good agreement, and compatible with electromagnetic simulations. Combining all the baselines, we make maps around bright sources and show that the array behaves as expected. A few hundred hours of observations at different declinations have been used to study the array geometry and pointing imperfections, as well as the instrument noise behaviour. We show that the system temperature is below 80 K for most feed antennas and that noise fluctuations decrease as expected with integration time, at least up to a few hundred seconds. Analysis of long integrations, from 10 nights of observations of the North Celestial Pole (NCP), yielded visibilities with amplitudes of 20–30 mK, consistent with the expected signal from the NCP radio sky with ${\lt}10\,$ mK precision for 1 MHz × 1 min binning. Hi-pass filtering the spectra to remove smooth spectrum signal yields a residual consistent with zero signal at the $0.5\,$ mK level.

(cosmology:) large-scale structure of the Universe↗
DOE OSTI2023

Closed Loop Geothermal Working Group: GeoCLUSTER App, Subsurface Simulation Results, and Publications

To better understand the heat production, electricity generation performance, and economic viability of closed-loop geothermal systems in hot-dry rock, the Closed-Loop Geothermal Working Group -- a consortium of several national labs and academic institutions has tabulated time-dependent numerical solutions and levelized cost results of two popular closed-loop heat exchanger designs (u-tube and co-axial). The heat exchanger designs were evaluated for two working fluids (water and supercritical CO2) while varying seven continuous independent parameters of interest (mass flow rate, vertical depth, horizontal extent, borehole diameter, formation gradient, formation conductivity, and injection temperature). The corresponding numerical solutions (approximately 1.2 million per heat exchanger design) are stored as multi-dimensional HDF5 datasets and can be queried at off-grid points using multi-dimensional linear interpolation. A Python script was developed to query this database and estimate time-dependent electricity generation using an organic Rankine cycle (for water) or direct turbine expansion cycle (for CO2) and perform a cost assessment. This document aims to give an overview of the HDF5 database file and highlights how to read, visualize, and query quantities of interest (e.g., levelized cost of electricity, levelized cost of heat) using the accompanying Python scripts. Details regarding the capital, operation, and maintenance and levelized cost calculation using the techno-economic analysis script are provided. This data submission will contain results from the Closed Loop Geothermal Working Group study that are within the public domain, including publications, simulation results, databases, and computer codes. GeoCLUSTER is a Python-based web application created using Dash, an open-source framework built on top of Flask that streamlines the building of data dashboards. GeoCLUSTER provides users with a collection of interactive methods for streamlining the exploration and visualization of an HDF5 dataset. The GeoCluster app and database are contained in the compressed file geocluster_vx.zip, where the "x" refers to the version number. For example, geocluster_v1.zip is Version 1 of the app. This zip file also contains installation instructions. **To use the GeoCLUSTER app in the cloud, click the link to "GeoCLUSTER on AWS" in the Resources section below. To use the GeoCLUSTER app locally, download the geocluster_vx.zip to your computer and uncompress this file. When uncompressed this file comprises two directories and the geocluster_installation.pdf file. The geo-data app contains the HDF5 database in condensed format, and the GeoCLUSTER directory contains the GeoCLUSTER app in the subdirectory dash_app, as app.py. The geocluster_installation.pdf file provides instructions on installing Python, the needed Python modules, and then executing the app.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

University of Missouri Research Reactor (MURR) Design Demonstration Element End Fitting Structural Rigidity Analysis

The primary objective of this work is to assess the extent to which the stiffness (measured by means of maximum displacement) of the end fittings in the DDE contributes to the stiffness of the entire element, and how it compares to the equivalent stiffness of the end fittings in the LEU element. Structural analysis of both the LEU element and the DDE were performed using COMSOL 5.3a finite element software. Supporting combs are used on the leading and trailing edges of fuel plates for both the LEU element and the DDE. Therefore, simulations with and without combs are performed as two bounding boundary conditions on the leading edge of the fuel plates. Three types of loads are analyzed in this work: the hydraulic load due to the channel flow disparity-induced pressure differential, the thermal load due to the thermal expansion of the fuel plates, and a point load equal in magnitude to the LEU element’s weight.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗