Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “casing integrity”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Computer Simulation of the Forces Acting on a Submerged Polystyrene Probe as it Approaches the Succinonitrile Melt-Solid Interface

A Modeling approach to simulate both mesoscale and microscopic forces acting in a typical AFM experiment is presented. A mesoscale level interaction between the cantilever tip and the sample surface is primarily described by the balance of attractive Van der Waals and repulsive forces. Ultimately, the goal is to measure the forces between a particle and the crystal-melt interface. Two modes of AFM operation are considered in this paper - a stationary and a "tapping" one. The continuous mechanics approach to model tip-surface interaction is presented. At microscopic levels, tip contamination and details of tip-surface interaction are modeled using a molecular dynamics approach for the case of polystyrene - succinonitrile contact. Integration of the mesoscale model with a molecular dynamic model is discussed.

Bune, Andris V.↗
NASA NTRS2003

Neural Flight Control System

The Neural Flight Control System (NFCS) was developed to address the need for control systems that can be produced and tested at lower cost, easily adapted to prototype vehicles and for flight systems that can accommodate damaged control surfaces or changes to aircraft stability and control characteristics resulting from failures or accidents. NFCS utilizes on a neural network-based flight control algorithm which automatically compensates for a broad spectrum of unanticipated damage or failures of an aircraft in flight. Pilot stick and rudder pedal inputs are fed into a reference model which produces pitch, roll and yaw rate commands. The reference model frequencies and gains can be set to provide handling quality characteristics suitable for the aircraft of interest. The rate commands are used in conjunction with estimates of the aircraft s stability and control (S&C) derivatives by a simplified Dynamic Inverse controller to produce virtual elevator, aileron and rudder commands. These virtual surface deflection commands are optimally distributed across the aircraft s available control surfaces using linear programming theory. Sensor data is compared with the reference model rate commands to produce an error signal. A Proportional/Integral (PI) error controller "winds up" on the error signal and adds an augmented command to the reference model output with the effect of zeroing the error signal. In order to provide more consistent handling qualities for the pilot, neural networks learn the behavior of the error controller and add in the augmented command before the integrator winds up. In the case of damage sufficient to affect the handling qualities of the aircraft, an Adaptive Critic is utilized to reduce the reference model frequencies and gains to stay within a flyable envelope of the aircraft.

Gundy-Burlet, Karen↗
NASA NTRS2003

Optimal Feedback Control of Thermal Networks

An improved approach to the mathematical modeling of feedback control of thermal networks has been devised. Heretofore software for feedback control of thermal networks has been developed by time-consuming trial-and-error methods that depend on engineers expertise. In contrast, the present approach is a systematic means of developing algorithms for feedback control that is optimal in the sense that it combines performance with low cost of implementation. An additional advantage of the present approach is that a thermal engineer need not be expert in control theory. Thermal networks are lumped-parameter approximations used to represent complex thermal systems. Thermal networks are closely related to electrical networks commonly represented by lumped-parameter circuit diagrams. Like such electrical circuits, thermal networks are mathematically modeled by systems of differential-algebraic equations (DAEs) that is, ordinary differential equations subject to a set of algebraic constraints. In the present approach, emphasis is placed on applications in which thermal networks are subject to constant disturbances and, therefore, integral control action is necessary to obtain steady-state responses. The mathematical development of the present approach begins with the derivation of optimal integral-control laws via minimization of an appropriate cost functional that involves augmented state vectors. Subsequently, classical variational arguments provide optimality conditions in the form of the Hamiltonian equations for the standard linear-quadratic-regulator (LQR) problem. These equations are reduced to an algebraic Riccati equation (ARE) with respect to the augmented state vector. The solution of the ARE leads to the direct computation of the optimal proportional- and integral-feedback control gains. In cases of very complex networks, large numbers of state variables make it difficult to implement optimal controllers in the manner described in the preceding paragraph.

Papalexandris, Miltiadis↗
NASA NTRS2017

Vented Launch Vehicle Adaptor for a Manned Spacecraft with "Pusher" Launch Abort System

A system, method, and apparatus for a vented launch vehicle adaptor (LVA) for a manned spacecraft with a "pusher" launch abort system are disclosed. The disclosed LVA provides a structural interface between a commercial crew vehicle (CCV) crew module/service module (CM/SM) spacecraft and an expendable launch vehicle. The LVA provides structural attachment of the module to the launch vehicle. It also provides a means to control the exhaust plume from a pusher-type launch abort system that is integrated into the module. In case of an on-pad or ascent abort, which requires the module to jettison away from the launch vehicle, the launch abort system exhaust plume must be safely directed away from critical and dangerous portions of the launch vehicle in order to achieve a safe and successful jettison.

Vandervort, Robert E.↗
NASA NTRSDate not supplied

Characterizing the Impact of Equatorial Plasma Bubbles on GNSS Accuracy for Crop Yield Estimation in Brazil with NASA GOLD Data

In Brazil, equatorial plasma bubbles are the primary space weather phenomenon impacting the performance of Global Navigation Satellite Systems (GNSS), causing signal disruptions and positional inaccuracies in agricultural machinery used for mapping fields, guiding machinery, and ensuring accurate planting, spraying, and harvesting. As a result, farmers and stakeholders face challenges in maintaining accurate operational capabilities, often needing to shut down machinery to mitigate these effects. This project evaluated the feasibility of using NASA’s Heliophysics observations for a case study in March 2024. Data from NASA's Global-scale Observations of the Limb and Disk (GOLD) instrument were used to characterize equatorial plasma bubbles, while amplitude and phase indices from ground-based GNSS receivers confirmed their presence and impact on signal disruptions. Additionally, data from a local farm's rover receiver assessed the practical implications of these disruptions on Real-Time Kinematic (RTK) GNSS-based operations by calculating the 2-dimensional error. The project involved collaboration across multiple different organizations and sectors, including NASA’s Crustal Dynamics Data Information System (CDDIS) researchers, John Deere, a Brazilian farm for data acquisition, and Brazilian universities for data interpretation. The results aim to inform Brazilian farmers about position errors affecting their precision agriculture operations and provide insights into the needs of space weather data users and modeling. Lessons from this case study will guide future integration of space weather data and modeling with end users.

Melissa Marquez↗
DOE OSTI2023

Autonomous Energy Systems Functional Interoperability

Functional Interoperability is the capability of two or more networks, systems, devices, applications, or components to externally exchange and readily use information securely and effectively in order to achieve a stated objective. Our proposition is that NREL's ARIES capability is ideal to evaluate functional interoperability by defining use case (objective) and scenarios, by integrating systems, devices, and applications in a realistic laboratory environment, and by simulating scenarios.

ARIES↗
DOE OSTI2023

Performance and Cost Potential for Direct-Fired Supercritical CO2 Natural Gas Power Plants

Direct-fired supercritical CO2 (sCO2) power cycles are being explored as an attractive alternative to natural gas combined cycle (NGCC) plants with carbon capture and storage (CCS). Therefore, understanding their performance and cost potential is important for the commercialization of the technology. This study presents the techno-economic optimization results of natural gas-fired, utility-scale power plants based on the direct sCO2 power cycle, which are lacking in public literature. To identify the optimum plant configuration, the study considered multiple cases with varying levels of thermal integration with the plant air separation unit (ASU). Several design variables for each power cycle configuration were identified and optimized to minimize the levelized cost of electricity (LCOE) for each case. The optimization design variables include the sCO2 cooler outlet temperatures, recuperator approach temperatures, and pressure drops. High fidelity models for recuperators, coolers, and turbines were developed and used to capture the impact of design variables on plant efficiency and capital costs. The optimization was conducted using a combination of manual sensitivity analyses and automated derivative-free optimization algorithms available under NETL’s Framework for Optimization and Quantification of Uncertainty and Sensitivity platform. The optimized direct sCO2 power plants offered similar or slightly higher plant efficiencies than the reference NGCC plants based on the F-class gas turbine with CCS. The LCOE of the optimized direct sCO2 plants is 13 to 17% higher than the reference NGCC plants with CCS due to high capital costs associated with the ASU and sCO2 power block, though there is significant room for improvement due to the high uncertainty in component capital costs for these new plants. Recuperators make up over 50% of the sCO2 power block costs. Consequently, any research and development efforts to reduce the recuperator capital costs will benefit the technology’s commercialization. The study also presents preliminary results showing the impact of co-firing landfill gas and natural gas on plant efficiency, LCOE, and CO2 emissions.

Pidaparti, Sandeep↗
DOE OSTI2023

Performance and Cost Potential for Direct-Fired Supercritical CO2 Natural Gas Power Plants

Direct-fired supercritical CO2 (sCO2) power cycles are being explored as an attractive alternative to natural gas combined cycle (NGCC) plants with carbon capture and storage (CCS). Therefore, understanding their performance and cost potential is important for the commercialization of the technology. This study presents the techno-economic optimization results of natural gas-fired, utility-scale power plants based on the direct sCO2 power cycle, which are lacking in public literature. To identify the optimum plant configuration, the study considered multiple cases with varying levels of thermal integration with the plant air separation unit (ASU). Several design variables for each power cycle configuration were identified and optimized to minimize the levelized cost of electricity (LCOE) for each case. The optimization design variables include the sCO2 cooler outlet temperatures, recuperator approach temperatures, and pressure drops. High fidelity models for recuperators, coolers, and turbines were developed and used to capture the impact of design variables on plant efficiency and capital costs. The optimization was conducted using a combination of manual sensitivity analyses and automated derivative-free optimization algorithms available under NETL’s Framework for Optimization and Quantification of Uncertainty and Sensitivity platform. The optimized direct sCO2 power plants offered similar or slightly higher plant efficiencies than the reference NGCC plants based on the F-class gas turbine with CCS. The LCOE of the optimized direct sCO2 plants is 13 to 17% higher than the reference NGCC plants with CCS due to high capital costs associated with the ASU and sCO2 power block, though there is significant room for improvement due to the high uncertainty in component capital costs for these new plants. Recuperators make up over 50% of the sCO2 power block costs. Consequently, any research and development efforts to reduce the recuperator capital costs will benefit the technology’s commercialization. The study also presents preliminary results showing the impact of co-firing landfill gas and natural gas on plant efficiency, LCOE, and CO2 emissions.

Pidaparti, Sandeep↗
DOE OSTI2024

Improved Characterization of Soil Organic Matter by Integrating FT-ICR MS, Liquid Chromatography Tandem Mass Spectrometry, and Molecular Networking: A Case Study of Root Litter Decay under Drought Conditions

Understanding of how soil organic matter (SOM) chemistry is altered in a changing climate has advanced considerably; however, most SOM components remain unidentified, impeding the ability to characterize a major fraction of organic matter and predict what types of molecules, and from which sources, will persist in soil. Here we present a novel approach to better characterize SOM extracts by integrating information from three types of analyses, and we deploy this method to characterize decaying root-detritus soil microcosms subjected to either drought or normal conditions. To observe broad differences in composition, we employed direct infusion Fourier-transform ion cyclotron resonance mass spectrometry (DI-FT-ICR MS). We complemented this with liquid chromatography tandem mass spectrometry (LC-MS/MS) to identify components by library matching. Since libraries contain only a small fraction of SOM components, we also used fragment spectral cosine similarity scores to relate unknowns and library matches through molecular networks. This integrated approach allowed us to corroborate DI-FT-ICR MS molecular formulas using library matches, which included fungal metabolites and related polyphenolic compounds. We also inferred structures of unknowns from molecular networks and improved LC-MS/MS annotation rates from ~5 to 35% by considering DI-FT-ICR MS molecular formula assignments. Under drought conditions, we found greater relative amounts of lignin-like vs condensed aromatic polyphenol formulas and lower average nominal oxidation state of carbon, suggesting reduced decomposition of SOM and/or microbes under stress. Our integrated approach provides a framework for enhanced annotation of SOM components that is more comprehensive than performing individual data analyses in parallel.

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

A finite element: Boundary integral method for electromagnetic scattering

A method that combines the finite element and boundary integral techniques for the numerical solution of electromagnetic scattering problems is presented. The finite element method is well known for requiring a low order storage and for its capability to model inhomogeneous structures. Of particular emphasis in this work is the reduction of the storage requirement by terminating the finite element mesh on a boundary in a fashion which renders the boundary integrals in convolutional form. The fast Fourier transform is then used to evaluate these integrals in a conjugate gradient solver, without a need to generate the actual matrix. This method has a marked advantage over traditional integral equation approaches with respect to the storage requirement of highly inhomogeneous structures. Rectangular, circular, and ogival mesh termination boundaries are examined for two-dimensional scattering. In the case of axially symmetric structures, the boundary integral matrix storage is reduced by exploiting matrix symmetries and solving the resulting system via the conjugate gradient method. In each case several results are presented for various scatterers aimed at validating the method and providing an assessment of its capabilities. Important in methods incorporating boundary integral equations is the issue of internal resonance. A method is implemented for their removal, and is shown to be effective in the two-dimensional and three-dimensional applications.

Collins, J. D.↗
DOE OSTI2022

Towards physics-inspired data-driven weather forecasting: integrating data assimilation with a deep spatial-transformer-based U-NET in a case study with ERA5

Abstract. There is growing interest in data-driven weather prediction (DDWP), e.g., using convolutional neural networks such as U-NET that are trained on data from models or reanalysis. Here, we propose three components, inspired by physics, to integrate with commonly used DDWP models in order to improve their forecast accuracy. These components are (1) a deep spatial transformer added to the latent space of U-NET to capture rotation and scaling transformation in the latent space for spatiotemporal data, (2) a data-assimilation (DA) algorithm to ingest noisy observations and improve the initial conditions for next forecasts, and (3) a multi-time-step algorithm, which combines forecasts from DDWP models with different time steps through DA, improving the accuracy of forecasts at short intervals. To show the benefit and feasibility of each component, we use geopotential height at 500 hPa (Z500) from ERA5 reanalysis and examine the short-term forecast accuracy of specific setups of the DDWP framework. Results show that the spatial-transformer-based U-NET (U-STN) clearly outperforms the U-NET, e.g., improving the forecast skill by 45 %. Using a sigma-point ensemble Kalman (SPEnKF) algorithm for DA and U-STN as the forward model, we show that stable, accurate DA cycles are achieved even with high observation noise. This DDWP+DA framework substantially benefits from large (O(1000)) ensembles that are inexpensively generated with the data-driven forward model in each DA cycle. The multi-time-step DDWP+DA framework also shows promise; for example, it reduces the average error by factors of 2–3. These results show the benefits and feasibility of these three components, which are flexible and can be used in a variety of DDWP setups. Furthermore, while here we focus on weather forecasting, the three components can be readily adopted for other parts of the Earth system, such as ocean and land, for which there is a rapid growth of data and need for forecast and assimilation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Integrating Upper Class E Traffic Management (ETM) Operations into the National Airspace System: Use Cases and Research Questions

As new categories of vehicles are introduced in the National Airspace System, so too are novel concepts for a cooperative approach to traffic management environments. One of these new environments, Upper Class E Traffic Management (ETM), is expected to include a variety of high altitude, long endurance vehicles with a range of performance capabilities and mission profiles that operate in cooperative areas above 60,000 feet. In addition to developing the rules, architecture, and systems for operations within the ETM environment itself, it is also important to consider how ETM vehicles will integrate with traditional Air Traffic Management and interact with Air Traffic Control (ATC) as they traverse ATC-controlled airspace and transition in and out of cooperative ETM operating areas. As a first step toward future ETM demonstrations at the National Aeronautics and Space Administration (NASA) Ames Research Center’s Airspace Operations Laboratory, use cases with step-by-step procedures were developed to identify both nominal and off-nominal scenarios in which ETM operations will interact with ATC. As NASA prepares to develop a simulation platform to demonstrate ETM cooperative practices and ETM-ATC interactions, the procedures, ATC roles and responsibilities, data exchange requirements, and research questions that were identified as part of use case development will inform scenario and system architecture design. The upcoming simulation work will include initial prototype ETM-ATC coordination tools to support ATC controllers’ interactions with ETM operations. This paper will briefly discuss NASA’s upcoming ETM development work and then provide background on ETM-ATC interactions, describe each ETM-ATC interaction use case, and discuss open questions on concept, procedures, and assumptions.

ATC↗
NASA NTRSDate not supplied

Integrating Upper Class E Traffic Management (ETM) Operations into the National Airspace System: Use Cases and Research Questions

As new categories of vehicles are introduced in the National Airspace System, so too are novel concepts for a cooperative approach to traffic management environments. One of these new environments, Upper Class E Traffic Management (ETM), is expected to include a variety of high altitude, long endurance vehicles with a range of performance capabilities and mission profiles that operate in cooperative areas above 60,000 feet. In addition to developing the rules, architecture, and systems for operations within the ETM environment itself, it is also important to consider how ETM vehicles will integrate with traditional Air Traffic Management and interact with Air Traffic Control (ATC) as they traverse ATC-controlled airspace and transition in and out of cooperative ETM operating areas. As a first step toward future ETM demonstrations at the National Aeronautics and Space Administration (NASA) Ames Research Center’s Airspace Operations Laboratory, use cases with step-by-step procedures were developed to identify both nominal and off-nominal scenarios in which ETM operations will interact with ATC. As NASA prepares to develop a simulation platform to demonstrate ETM cooperative practices and ETM-ATC interactions, the procedures, ATC roles and responsibilities, data exchange requirements, and research questions that were identified as part of use case development will inform scenario and system architecture design. The upcoming simulation work will include initial prototype ETM-ATC coordination tools to support ATC controllers’ interactions with ETM operations. This paper will briefly discuss NASA’s upcoming ETM development work and then provide background on ETM-ATC interactions, describe each ETM-ATC interaction use case, and discuss open questions on concept, procedures, and assumptions.

ATC↗
NASA NTRS2005

Computation of Sound Propagation by Boundary Element Method

This report documents the development of a Boundary Element Method (BEM) code for the computation of sound propagation in uniform mean flows. The basic formulation and implementation follow the standard BEM methodology; the convective wave equation and the boundary conditions on the surfaces of the bodies in the flow are formulated into an integral equation and the method of collocation is used to discretize this equation into a matrix equation to be solved numerically. New features discussed here include the formulation of the additional terms due to the effects of the mean flow and the treatment of the numerical singularities in the implementation by the method of collocation. The effects of mean flows introduce terms in the integral equation that contain the gradients of the unknown, which is undesirable if the gradients are treated as additional unknowns, greatly increasing the sizes of the matrix equation, or if numerical differentiation is used to approximate the gradients, introducing numerical error in the computation. It is shown that these terms can be reformulated in terms of the unknown itself, making the integral equation very similar to the case without mean flows and simple for numerical implementation. To avoid asymptotic analysis in the treatment of numerical singularities in the method of collocation, as is conventionally done, we perform the surface integrations in the integral equation by using sub-triangles so that the field point never coincide with the evaluation points on the surfaces. This simplifies the formulation and greatly facilitates the implementation. To validate the method and the code, three canonic problems are studied. They are respectively the sound scattering by a sphere, the sound reflection by a plate in uniform mean flows and the sound propagation over a hump of irregular shape in uniform flows. The first two have analytical solutions and the third is solved by the method of Computational Aeroacoustics (CAA), all of which are used to compare the BEM solutions. The comparisons show very good agreements and validate the accuracy of the BEM approach implemented here.

Guo, Yueping↗
NASA NTRS2016

A Tale of Two Chambers: Iterative Approaches and Lessons Learned from Life Support Systems Testing in Altitude Chambers

The drive for the journey to Mars is in a higher gear than ever before. We are developing new spacecraft and life support systems to take humans to the Red Planet. The journey that development hardware takes before its final incarnation in a fully integrated spacecraft can take years, as is the case for the Orion environmental control and life support system (ECLSS). Through the Pressure Integrated Suit Test (PIST) series, NASA personnel at Johnson Space Center have been characterizing the behavior of a closed loop ECLSS in the event of cabin depressurization. This kind of testing - one of the most hazardous activities performed at JSC - requires an iterative approach, increasing in complexity and hazards). The PIST series, conducted in the Crew and Thermal Systems Division (CTSD) 11-ft Chamber, started with unmanned test precursors before moving to a human-in-the-loop phase, and continues to evolve with the eventual goal of a qualification test for the final system that will be installed on Orion. Meanwhile, the Human Exploration Spacecraft Testbed for Integration and Advancement (HESTIA) program is an effort to research and develop technologies that will work in concert to support habitation on Mars. September 2015 marked the first unmanned HESTIA test, with the goal of characterizing how ECLSS technologies work together in a closed environment. HESTIA will culminate in crewed testing, but it can benefit from the lessons learned from another test that is farther ahead in its development and life cycle. Discussing PIST and HESTIA, this paper illustrates how we approach testing, the kind of information that facility teams need to ensure efficient collaborations and successful testing, and how we can apply what we learn to execute future tests.

Callini, Gianluca↗
DOE OSTI2024

Technical and Economic Assessment and Gap Analysis of Advanced Nuclear Reactor Integration with a Reference Oil Refinery

Efforts to identify the most-economic methods to decarbonize several sectors of the U.S. economy are underway. Industrial processes such as crude-oil refining rely heavily on energy-dense and easily stored and transported fossil fuels for powering their operations. Refineries use large amounts of energy, primarily derived from fossil sources to separate crude-oil components, break down heavier hydrocarbons into lighter compounds, remove impurities, reform hydrocarbon molecules, and generate steam and electricity for pumps and compressors and other various auxiliary systems. Crude-oil refining operations such as distillation, cracking, desulfurization, reforming, utilities systems and some offsite facilities collectively account for most of the energy consumption. Other operations such as hydrocracking or hydrotreating also require hydrogen for developing hydrogenation reactions which involve substantial heating to keep the reactors at high-temperature and pressure levels. All heat and energy demands are typically provided by natural gas (NG), oil, or other fuels, which makes refinery industry one of the most-difficult sectors to decarbonize. Nuclear power is a viable and energy-dense source of clean electricity, heat, and hydrogen to provide the large, sustainable energy supply that the refining industry demands. The U.S. Department of Energy’s (DOE’s) Integrated Energy Systems (IES) program is working to perform research and development, design, economic siting, and risk analysis. This state-of-the-art work will enable the first on-site demonstrations and commercial deployments of advanced small modular nuclear reactors (SMNRs) integrated with industries such as chemical production, refining, iron and steel making, and more. IES seeks to demonstrate the ability of advanced nuclear reactors to meet the heat and power demands of these industries while reducing carbon emissions in a sustainable and cost-competitive way. The primary objective of this research effort is to analyze industrial-scale SMNR integration intended to decarbonize refining facilities. The foreseen outcome is the provision of reliable, cost-competitive, and sustainable clean energy, alongside a reduction of carbon emissions. Specifically, the focus of this work lies on meeting the reference facilities’ heat and electricity demands with nuclear power while also supplying clean hydrogen via integrated high-temperature steam electrolysis (HTSE). This report presents a comprehensive technical and economic assessment of the integration of advanced nuclear reactors into a reference refinery, leveraging financial incentives from the Inflation Reduction Act (IRA). The evaluation aims to explore the potential economic benefits and challenges associated with incorporating advanced nuclear reactors into refinery operations, particularly in terms of energy efficiency, economic implications and environmental impact. By examining both the technical feasibility and economic viability, this analysis seeks to identify existing gaps and propose solutions for successful nuclear integration implementation. The findings are intended to provide valuable insights for stakeholders considering the adoption of advanced nuclear reactors in the refining sector. A refinery reference-plant was developed, using an open-source refinery model, Petroleum Refinery Lifecycle Inventory Model (PRELIM) and expert assessment, as a base case for comparison with various nuclear integration options. The capacity of 100 kbd/day (KBD) of heavy crude-oil feed was selected to represent a general coking-type refinery with deep conversion capabilities (incorporating heavy-oil upgrading with FCC, coking, and associated hydrotreating process units), using a heavy crude-oil feed, which represents about 70% of U.S. refineries configurations. A summary of all cases considered in this study is shown in Table 1.

13 HYDRO ENERGY↗
DOE OSTI2020

Analytical integration of the tractions induced by non-singular dislocations on an arbitrary shaped triangular quadratic element

Here, an analytical model is proposed to evaluate the nodal force induced by a segment of dislocation upon an arbitrary shaped triangular element. This calculation is required in hybrid methods that associate dislocation dynamics to boundary or finite element to solve simultaneously the evolution of large ensembles of dislocations with complex boundary conditions. Nodal forces are defined as the triple integration of the unbalanced traction field induced by a straight dislocation upon the surface of the element. Following our previous approach (Queyreau et al 2014 Modelling Simul. Mater. Sci. Eng. 22 035004) on a simpler geometry and in the case of linear isotropic elasticity, triple integrals are solved by sequences of integration by parts that exhibit recurrence relations. The traction field is defined and finite everywhere even at the core of dislocations, thanks to the use of the non-singular stress expression formulated by Cai et al (2006 J. Mech. Phys. Solids 54 561–587). The nodal force expressions can be used when considering both a single convolution or double convolution of the Green's function with the core distribution. A solution is also proposed for the case of a semi-infinite segment through the study of the asymptotic behavior of the analytical expressions. The proposed approach is exact and very computationally efficient. The choice of an arbitrary shaped triangular element and quadratic shape functions allow the consideration of complex geometries and comply with automatic meshing procedures. These analytical expressions could also be employed to estimate dislocation interactions with interfaces.

42 ENGINEERING↗