Search NASASearch

SEARCH · Search NASA

Results for “Finite State Machine”

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 109 records · Page 6

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials

Pore-scale visualization of natural hydrate-bearing sediments

Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.

Liu, Mengwei [Oak Ridge Institute for Science and

Design and Evaluation of Composite Fuselage Panels Subjected to Combined Loading Conditions

Methodologies used in industry for designing transport aircraft composite fuselage structures are discussed. Several aspects of the design methodologies are based on assumptions from metallic fuselage technology which requires that full-scale structures be tested with the actual loading conditions to validate the designs. Composite panels which represent crown and side regions of a fuselage structure are designed using this approach and tested in biaxial tension. Descriptions of the state-of-the-art test facilities used for this structural evaluation are presented. These facilities include a pressure-box test machine and a D-box test fixture in a combined loads test machine which are part of a Combined Loads Test System (COLTS). Nonlinear analysis results for a reference shell and a stiffened composite panel tested in the pressure-box test machine with and without damage are presented. The analytical and test results are compared to assess the ability of the pressure-box test machine to simulate a shell stress state with and without damage. A combined loads test machine for testing aircraft primary structures is described. This test machine includes a D-box test fixture to accommodate curved stiffened panels and the design features of this test fixture are presented. Finite element analysis results for a curved panel to be tested in the D-box test fixture are also discussed.

Ambur, Damodar R.

Verification of Numerical Algorithms

The following strategy is suggested for specification and proof: (1) Defer the construction of a formal program specification with respect to I/O assertions unit the correctness of the program with respect to an abstract mathematical model of program intent is demonstrated. (2) Prove that an abstract machine (using infinite precision arithmetic) would compute that object exactly. (3) Prove that the computational sequences of arithmetic operations that occur in the abstract machine must be precisely the same at every step as those occurring on an actual machine (with finite precision arithmetic), executing the same program. (4) Use a Verification Conditions VC-generator that knows about the semantics of arithmetic operations to annotate the program with assertions that bound (or in some circumstances estimate) the difference between the actual machine state variables and the corresponding ones of the abstract machine. Construct the formal program specification by combining the verification conditions into theorems about computational error that can be proved with mechanical assistance.

Source record

Solving high-dimensional partial integral differential equations: The finite expression method

Partial integro-differential equations (PIDEs) have broad applications in the sciences, from electro-magnetism to options pricing. Here, in this paper, we introduce a new finite expression method (FEX) to solve PIDEs. This approach builds upon the original FEX and its inherent advantages with new advances: 1) A novel method of parameter grouping is proposed to reduce the number of coefficients in high-dimensional function approximation; 2) A Taylor series approximation method is implemented to significantly improve the computational efficiency and accuracy of the evaluation of the integral terms of PIDEs. The new FEX based method, denoted FEX-PG to indicate the addition of the parameter grouping (PG) step to the algorithm, provides both high accuracy and interpretable numerical solutions, with the outcome being an explicit equation that facilitates intuitive understanding of the underlying solution structures. These features are often absent in traditional methods, such as finite element methods (FEM) and finite difference methods, as well as in deep learning-based approaches. To benchmark our method against recent advances, we apply the new FEX-PG to solve benchmark PIDEs in the literature. In high-dimensional settings, FEX-PG exhibits strong and robust performance, achieving relative errors on the order of single precision machine epsilon, significantly outperforming existing approaches based on neural networks.

Combinatorial optimization

Constrained or unconstrained? Neural-network-based equation discovery from data

Throughout many fields, practitioners often rely on differential equations to model systems. Yet, for many applications, the theoretical derivation of such equations and/or the accurate resolution of their solutions may be intractable. Instead, recently developed methods, including those based on parameter estimation, operator subset selection, and neural networks, allow for the data-driven discovery of both ordinary and partial differential equations (PDEs), on a spectrum of interpretability. The success of these strategies is often contingent upon the correct identification of representative equations from noisy observations of state variables and, as importantly and intertwined with that, the mathematical strategies utilized to enforce those equations. Specifically, the latter has been commonly addressed via unconstrained optimization strategies. Representing the PDE as a neural network, we propose to discover the PDE (or the associated operator) by solving a constrained optimization problem and using an intermediate state representation similar to a physics-informed neural network (PINN). The objective function of this constrained optimization problem promotes matching the data, while the constraints require that the discovered PDE is satisfied at a number of spatial collocation points. We present a penalty method and a widely used trust-region barrier method to solve this constrained optimization problem, and we compare these methods on numerical examples. Our results on several example problems demonstrate that the latter constrained method outperforms the penalty method, particularly for higher noise levels or fewer collocation points. This work motivates further exploration into using sophisticated constrained optimization methods in scientific machine learning, as opposed to their commonly used, penalty-method or unconstrained counterparts. For both of these methods, we solve these discovered neural network PDEs with classical methods, such as finite difference methods, as opposed to PINNs-type methods relying on automatic differentiation. Here, we briefly highlight how simultaneously fitting the data while discovering the PDE improves the robustness to noise and other small, yet crucial, implementation details.

Data-driven discovery

Active magnetic bearing control loop modeling for a finite element rotordynamics code

A mathematical model of an active electromagnetic bearing which includes the actuator, the sensor and the control system is developed and implemented in a specialized finite element code for rotordynamic analysis. The element formulation and its incorporation in the model of the machine are described in detail. A solution procedure, based on a modal approach in which the number of retained modes is controlled by the user, is then shown together with other procedures for computing the steady-state response to both static and unbalance forces. An example of application shows the numerical results obtained on a model of an electric motor suspended on a five active-axis magnetic suspension. The comparison of some of these results with the experimental characteristics of the actual system shows the ability of the present model to predict its performance.

Genta, Giancarlo

Local Analysis-Test Correlation Of Tow-Steered Composite Shells With Small Cutouts

The prebuckling and postbuckling behavior of two composite tow-steered shells with small cutouts is assessed using nonlinear dynamic finite element analyses and compared in detail with experimental measurements. The cylindrical shells were manufactured without cutouts using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the axially stiff crown and keel, to ±45 degrees on the shear-stiff sides. The first shell with overlaps has laminate thickness variations on the crown and keel that result from application of all 24 tows during each pass of the fiber placement system. The second shell without overlaps uses the fiber placement system’s tow drop/add capability to achieve a more uniform shell wall thickness. An unreinforced cutout representing a passenger door on a commercial aircraft fuselage is machined into the side of each of the two shells. These shells with cutouts were tested in axial compression and buckled elastically in previous work. Detailed nonlinear finite element analysis results are compared with their corresponding measured local load-displacement and load-strain responses in prebuckling, at global buckling, and into a stable postbuckled state. Test data from displacement transducers, strain gages, and digital image correlation are extracted at the centers of the crown and keel, and at the middles of the top and left edges of the cutout. The overall agreement between these measured and analytical responses is excellent in prebuckling through global buckling, and very good from global buckling through postbuckling. As such, the excellent correlation observed here increases confidence in applying tow-steered composites in operational vehicles.

Composites

Local Analysis-Test Correlation Of Tow-Steered Composite Shells With Small Cutouts

The prebuckling and postbuckling behavior of two composite tow-steered shells with small cutouts is assessed using nonlinear dynamic finite element analyses and compared in detail with experimental measurements. The cylindrical shells were manufactured without cutouts using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the axially stiff crown and keel, to ±45 degrees on the shear-stiff sides. The first shell with overlaps has laminate thickness variations on the crown and keel that result from application of all 24 tows during each pass of the fiber placement system. The second shell without overlaps uses the fiber placement system’s tow drop/add capability to achieve a more uniform shell wall thickness. An unreinforced cutout representing a passenger door on a commercial aircraft fuselage is machined into the side of each of the two shells. These shells with cutouts were tested in axial compression and buckled elastically in previous work. Detailed nonlinear finite element analysis results are compared with their corresponding measured local load-displacement and load-strain responses in prebuckling, at global buckling, and into a stable postbuckled state. Test data from displacement transducers, strain gages, and digital image correlation are extracted at the centers of the crown and keel, and at the middles of the top and left edges of the cutout. The overall agreement between these measured and analytical responses is excellent in prebuckling through global buckling, and very good from global buckling through postbuckling. As such, the excellent correlation observed here increases confidence in applying tow-steered composites in operational vehicles.

Composites

LeWRON: Learning ElectroWeak phase tRansitiON with agentic architecture

An agent to analyze electroweak phase transition. LeWRON turns a BSM model description or a reproduction target into a structured run: symbolic setup, effective-potential artifacts, finite-temperature machinery, generated model code, a scientific report, and an interactive exploration session. It keeps both machine-readable artifacts and human-readable notes, so a run can be resumed, audited, revised, and shared.

Wang, Isaac [Fermi National Accelerator Laboratory

Three dimensional magnetic fields in extra high speed modified Lundell alternators computed by a combined vector-scalar magnetic potential finite element method

A 3D finite element (FE) approach was developed and implemented for computation of global magnetic fields in a 14.3 kVA modified Lundell alternator. The essence of the new method is the combined use of magnetic vector and scalar potential formulations in 3D FEs. This approach makes it practical, using state of the art supercomputer resources, to globally analyze magnetic fields and operating performances of rotating machines which have truly 3D magnetic flux patterns. The 3D FE-computed fields and machine inductances as well as various machine performance simulations of the 14.3 kVA machine are presented in this paper and its two companion papers.

Demerdash, N. A.

Analysis of composite ablators using massively parallel computation

In this work, the feasibility of using massively parallel computation to study the response of ablative materials is investigated. Explicit and implicit finite difference methods are used on a massively parallel computer, the Thinking Machines CM-5. The governing equations are a set of nonlinear partial differential equations. The governing equations are developed for three sample problems: (1) transpiration cooling, (2) ablative composite plate, and (3) restrained thermal growth testing. The transpiration cooling problem is solved using a solution scheme based solely on the explicit finite difference method. The results are compared with available analytical steady-state through-thickness temperature and pressure distributions and good agreement between the numerical and analytical solutions is found. It is also found that a solution scheme based on the explicit finite difference method has the following advantages: incorporates complex physics easily, results in a simple algorithm, and is easily parallelizable. However, a solution scheme of this kind needs very small time steps to maintain stability. A solution scheme based on the implicit finite difference method has the advantage that it does not require very small times steps to maintain stability. However, this kind of solution scheme has the disadvantages that complex physics cannot be easily incorporated into the algorithm and that the solution scheme is difficult to parallelize. A hybrid solution scheme is then developed to combine the strengths of the explicit and implicit finite difference methods and minimize their weaknesses. This is achieved by identifying the critical time scale associated with the governing equations and applying the appropriate finite difference method according to this critical time scale. The hybrid solution scheme is then applied to the ablative composite plate and restrained thermal growth problems. The gas storage term is included in the explicit pressure calculation of both problems. Results from ablative composite plate problems are compared with previous numerical results which did not include the gas storage term. It is found that the through-thickness temperature distribution is not affected much by the gas storage term. However, the through-thickness pressure and stress distributions, and the extent of chemical reactions are different from the previous numerical results. Two types of chemical reaction models are used in the restrained thermal growth testing problem: (1) pressure-independent Arrhenius type rate equations and (2) pressure-dependent Arrhenius type rate equations. The numerical results are compared to experimental results and the pressure-dependent model is able to capture the trend better than the pressure-independent one. Finally, a performance study is done on the hybrid algorithm using the ablative composite plate problem. It is found that there is a good speedup of performance on the CM-5. For 32 CPU's, the speedup of performance is 20. The efficiency of the algorithm is found to be a function of the size and execution time of a given problem and the effective parallelization of the algorithm. It also seems that there is an optimum number of CPU's to use for a given problem.

Shia, David

Absolute stability analysis of attitude control systems for large boosters.

A method for performing absolute stability analyses of attitude control systems for large launch vehicles is presented. Absolute stability of these systems is shown in a finite region of the state space. The regions are computed by using the Lur'e-Postnikov Liapunov function. This function is chosen to provide additional information about the exponential property of absolute stability. Significant advantages of the method proposed in this paper are: it is independent of the order of the system; algebraic operations involved in the computations are relatively simple and convenient for machine implementation; and the obtained results are valid not only for a particular nonlinearity but also for an entire class of nonlinear characteristics that satisfy certain general conditions. A system model representing the Saturn V launch vehicle is used to illustrate the method.

Seltzer, S. M.

Progress Toward the Critical Design of the Superconducting Rotor for NASA's 1.4 MW High-Efficiency Electric Machine

NASA Glenn Research Center is developing a 1.4 MW high-efficiency electric machine for future electrified aircraft to reduce energy consumption, emissions, and noise. This wound-field, synchronous machine employs a self-cooled, superconducting rotor to achieve excellent specific power and efficiency. This paper presents the progress made over the past year toward a critical design of the superconducting rotor, including refined analysis and risk reduction testing on prototype superconducting coils. A 3D finite element model of the rotor is developed to simulate the combined loading imposed by the initial cool down of the rotor from room temperature to cryogenic temperature and the subsequent rotation of the rotor up to the design speed. The model is evaluated to better predict the stress state in the critical components and the reduction in physical gap between the stator and rotor due to the resulting radial deflections. This model includes an enhanced model of the superconducting coil that considers the coil as an assembly of four individual coil layers. An improved coil fabrication process is described in detail. Two sub-scale test coils are fabricated and then tested under repeated thermal cycling between room temperature and liquid nitrogen temperature (77 K). It is demonstrated that the proposed coil fabrication procedure can produce no-insulation, high temperature superconducting coils that can survive several thermal cycles without apparent degradation.

electric machines

Towards real-time simulation of large space structures: Stabilization of fluid/thermal/structure interactions and implementation on high performance supercomputers

Within the Center for Space Construction, the SIMSTRUC project's objectives center around the development of simulation tools for the realistic analysis of large space structures. The word 'tools' is the broad sense; it designates mathematical models, finite element/finite difference formulations, computational algorithms, implementations on advanced computer architectures, and visualization capabilities. The results of our activities during the first year within the SIMSTRUC project are reported. On the modeling side, an alternative approach to fluid/thermal/structure interaction analysis that is a departure from the 'loosely coupled' and 'unified' approaches that are being currently practiced are described. The advantages of our approach both in terms of accuracy and computational efficiency were demonstrated. On the computational side, a software architecture for parallel/vector and massively parallel supercomputers that speeds up finite element and finite difference computations by several orders of magnitude is presented. As an example, the simulation of the deployment of a space structure that used to require over six hours of a workstation using a conventional finite element software, now runs on a multiprocessor using a parallel computation strategy in less than three seconds. In order to promote the physical understanding of the simulation behavior, a real-time visualization capability on the Connection Machine, which allows the analyst to watch the graphical animation of the results at the same time these are generated, was also developed. It is believed that by combining efficient analytical formulations with the state-of-the-art high performance computer implementations and superfast visualization capabilities, SIMSTRUC is moving fast towards the real-time simulation of large space structures. The designers as well as the researchers will certainly benefit from this technology.

Farhat, C.

Delamination-informed lifecycle decisions: A dielectric and machine learning framework for composite sorting and recycling

Composite materials are widely used in aerospace, marine, and automotive sectors due to their high strength-to-weight ratio and durability. However, their long-term reliability can be compromised by damage accumulation. Specifically, delamination initiation serves as a precursor to structural failure, which is often difficult to detect during damage inspection. Identifying and sorting delamination initiation in samples not only increases operational safety while providing critical information for end-of-life decisions, which influences both the service life extension value and the efficiency of fiber extraction during recycling. This research addresses two challenges: (1) developing a nondestructive, ex-situ framework to sort composite materials based on damage severity, particularly delamination, and (2) understanding how damage in composites influences resin removal during pyrolysis. Both experimental work and finite element analysis were performed to predict critical stress levels that are associated with delamination onset. Based on these results, three loading levels 50 %, 75 %, and 90 % of maximum stress, were selected for controlled experiments, generating composite samples with varying extents of damage for machine learning model training. Microscopic imaging of these samples confirmed the damage progression from matrix cracking to delamination, validating the computational predictions. We explored supervised machine learning using dielectric measurements to classify damage states. Preliminary results show an artificial neural network can identify early delamination which is a potential precursor to failure, with 94.44 % accuracy on our dataset. A parallel investigation into the effect of damage severity on pyrolysis recycling showed that heavily delaminated samples required significantly less energy for comparable matrix removal than undamaged samples.

dielectric variables

Application of two-point difference schemes to the conservative Euler equations for one-dimensional flows

An implicit finite difference method is presented for obtaining steady state solutions to the time dependent, conservative Euler equations for flows containing shocks. The method used the two-point differencing approach of Keller with dissipation added at supersonic points via the retarded density concept. Application of the method to the one-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions shows the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. It is shown that the scheme offers certain advantages over the more widely used three-point schemes, especially in regard to application of boundary conditions.

Wornom, S. F.

An implementation of neural simulation-based inference for parameter estimation in ATLAS

Neural simulation-based inference (NSBI) is a powerful class of machine-learning-based methods for statistical inference that naturally handles high-dimensional parameter estimation without the need to bin data into low-dimensional summary histograms. Such methods are promising for a range of measurements, including at the Large Hadron Collider, where no single observable may be optimal to scan over the entire theoretical phase space under consideration, or where binning data into histograms could result in a loss of sensitivity. This work develops a NSBI framework for statistical inference, using neural networks to estimate probability density ratios, which enables the application to a full-scale analysis. It incorporates a large number of systematic uncertainties, quantifies the uncertainty due to the finite number of events in training samples, develops a method to construct confidence intervals, and demonstrates a series of intermediate diagnostic checks that can be performed to validate the robustness of the method. As an example, the power and feasibility of the method are assessed on simulated data for a simplified version of an off-shell Higgs boson couplings measurement in the four-lepton final states. This approach represents an extension to the standard statistical methodology used by the experiments at the Large Hadron Collider, and can benefit many physics analyses.

frequentist statistics