SEARCH · Search NASA

Results for “embedded methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Fluid & Thermal Analysis of a Manifold Microchannel Heat Sink

A three-dimensional (3-D) numerical model has been developed to study the fluid flow and heat transfer through a manifold-microchannel (MMC) heat sink, a microfluidic cooling system designed as a heat dissipation solution for compact, high-powered electronic systems. An MMC consists of a 3-D manifold structure that distributes working fluid through alternating inlet and outlet manifold channels that guide the fluid to and from the microchannels, shortening the flow length and thereby, reducing the pressure drop across the system. Wide bandgap semiconductors such as silicon carbide (SiC) are increasingly utilized due to their ability to enable smaller electronic devices and components to operate more efficiently at higher frequencies. These semiconductors can greatly benefit from the MMC, which can extract all the heat produced by compact, high-powered electronics, enabling them to achieve optimal performance. The MMC also utilizes an embedded cooling approach which allows for fabrication directly within SiC-based electronics, resulting in significant mass and weight savings and lower thermal resistance than conventional heat sink designs. This paper outlines the methods used to conduct fluid and thermal analysis of different MMC configurations for flow parameters such as flow rate and inlet flow conditions. This study involves both steady state and transient thermal analysis as a gradually increasing heat flux is applied to the system over time, which is representative of the heat dissipation from the electronics over a given period. Results from the analysis show that the MMC is capable of effectively dissipating heat flux of ~100 W/cm2 while maintaining considerably low pressure drop (< 1 kPa) and thermal resistance (< 1 K/W) for certain single-phase working fluids.

manifold microchannel↗

NASA NTRSDate not supplied

Evaluating Lunar Water Processing System Model Configurations for Small Scale Oxygen and Hydrogen Production Within JAXA’s ISRU Technology Demonstration Plant

In-Situ Resource Utilization (ISRU) refers to novel methods of extracting and processing local resources for use in life support and propulsion systems, reducing or eliminating the required consumables to be transferred from Earth. Current estimates of water-ice availability embedded in regolith within the Moon’s permanently shadowed regions (PSR’s)range between 1-5% by weight. However, the composition and characteristics of the “wet” regolith is unknown. Alternate ISRU excavation techniques and Concept of Operations (ConOps) must be explored to optimize surface system operations based on these factors. To assess the feasibility of different ISRU subsystem technologies and compare system architecture configurations, an interchangeable system model was generated to incorporate technologies spanning excavation of raw materials to storage of products and determine optimal arrangement of total system processing needs. Total Mass, Volume, and Power (M/V/P) requirements were computed for 168 design iterations of this water processing plant.

ISRU↗

NASA NTRSDate not supplied

Evaluating Lunar Water Processing System Model Configurations for Small Scale Oxygen and Hydrogen Production Within JAXA'S ISRU Technology

Introduction: In-Situ Resource Utilization (ISRU) refers to novel methods of extracting and processing local resources for use in life support and propulsion systems, reducing or eliminating the required consumables to be transferred from Earth. Current estimates of water-ice availability embedded in regolith within the Moon’s permanently shadowed regions (PSR’s) range between 1-5% by weight. However, the composition and characteristics of the “wet” regolith is unknown. Alternate ISRU excavation techniques and Concept of Operations (ConOps) must be explored to optimize surface system operations based on these factors. To assess the feasibility of different ISRU subsystem technologies and compare system architecture configurations, an interchangeable system model was generated to incorporate technologies spanning excavation of raw materials to storage of products and determine optimal arrangement of total system processing needs. Total Mass, Volume, and Power (M/V/P) requirements were computed for 168 design iterations of this water processing plant. System Model: In FY24, the System Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) team developed a lunar water processing system model using the Mission Analysis and Integration Tool (MAIT) to estimate the M/V/P for ISRU subsystems operating under a wide range of Hydrogen (H2) and Oxygen (O2) production targets for the Space Technology Mission Directorate (STMD) [1]. Based on Japan Aerospace Exploration Agency’s (JAXA) surface operational requirements, this system architecture was modified to include the ability to excavate consolidated icy regolith (versus granular ice excavation using Kennedy Space Center’s (KSC) ISRU Pilot Excavator, IPEx) and explore the feasibility of processing the lunar water both inside and outside of the PSR. For the consolidated icy regolith case study, excavation was performed via a mobility transport chassis outfitted with The Regolith Ice Drill for Exploring New Terrain (TRIDENT) for drilling [2] and the Cold Operable Lunar Deployable Arm (COLDArm) [3] for regolith transfer. The system model determines the required rover and payload. M/V/P to handle the required regolith processing rates. The regolith is then sorted and heated to sublimate the ice (via an auger dryer). The exiting high temperature, low pressure vapor is cleaned of volatiles (via cold trap) and electrolyzed to produce H2 and O2. These products are then dried, liquified with 20 K and 90 K cryocoolers (for H2 and O2, respectively), and stored in cylindrical tanks. Study Goals: Due to the different ConOps options of regolith transport to the ridge for processing versus processing it directly inside the PSR, as well as the unknown regolith/water-ice composition, new excavation techniques and their power configurations are being evaluated within a ISRU system architecture for production targets less than NASA’s pilot plant (1 mT). This analysis investigates the feasibility of numerous excavation techniques, power architectures, and logistical operations and determines an optimal system configuration with regards to M/V/P. It aims to investigate which parameters, both locally and globally, have the greatest effect on each subsystem within the plant. This can be used to identify the most critical components of the plant, and guide future decisions on allocating funding for research and development. The results from this study may provide subsystem developers with appropriate interfaces with excavation subsystems and downstream processes, and assessing the overall feasibility of each excavation technique, power architecture, and logistical timeframe. References: [1] Carlson, A. et. al. (2024) ICES. [2] Zacny, K., et. al. (2024) “ASCE Earth and Space”. [3] McCormick, R., et. Al. (2024) IEEE Xplore.

ISRU↗

DOE OSTI2024

Modeling bicarbonate formation in an alkaline solution with multi-level quantum mechanics/molecular dynamics simulations

Understanding carbonate speciation and how it may be modulated is essential for the advancement of carbon dioxide (CO 2 ) capture and storage technologies, which often rely on the transformation of CO 2 into carbonate, e.g. via the formation of carbonate minerals. To date, few atomic-level, quantum-mechanics-based simulations have been carried out to characterize how carbonic acid (H 2 CO 3 ) and bicarbonate ($HCO^{-}_{3}$) form in aqueous solution, and how pH affects this process. Recently, Martirez and Carter utilized rare-event sampling density functional theory molecular dynamics simulations in combination with multi-level embedded correlated wavefunction theory, thus accounting for both solvent dynamics and electron correlation accurately, to elucidate the mechanism of H 2 CO 3 formation in neutral solution (J. Am. Chem. Soc., 145, 12561, 2023). Here, we perform a complementary simulation using the same method to map out the energetics of $HCO^{-}_{3}$ formation from dissolved CO 2 in basic solution. We find that, as in H 2 CO 3 formation, including water dynamics is important to obtain an accurate prediction of the energetics for the aforementioned reaction. Furthermore, only with MD did we identify the correct pathway for the reaction, in which water – not hydroxide – acts as the initial nucleophile and only at the transition state does it lose a proton.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1995

Method and apparatus for measuring temperatures in fabrics and flexible thermal insulations

A temperature sensor uses a type R thermocouple wire element in a ceramic sheath to sense temperatures up to 3,200 deg F., and is particularly suitable for flexible insulations. The sensor includes a thermocouple wire embedded in a sheath having two sections disposed at right angles to each other. The junction of the thermocouple is located at one end of one of the sections and the lead wires extend from the other section. The section which includes the junction is secured to a flexible surface with ceramic cement.

Kourtides, Demetrius A.↗

NASA NTRS1980

Method for producing a biological reagent

Polymeric functional microspheres containing metal or metal compounds are formed by addition polymerization of a covalently bondable olefinic monomer such as hydroxyethylmethacrylate in the presence of finely divided metal or metal oxide particles, such as iron, gold, platinum or magnetite, which are embedded in the resulting microspheres. The microspheres can be covalently bonded to chemotherapeutic agents, antibodies, or other proteins providing a means for labeling or separating labeled cells. Labeled cells or microspheres can be concentrated at a specific body location such as in the vicinity of a malignant tumor by applying a magnetic field to the location and then introducing the magnetically attractable microspheres or cells into the circulatory system of the subject. Labeled cells can be separated from a cell mixture by applying a predetermined magnetic field to a tube in which the mixture is flowing. After collection of the labeled cells, the magnetic field is discontinued and the labeled sub-cell population recovered.

Yen, Shiao-Ping S.↗

NASA NTRS2014

Automated Testcase Generation for Numerical Support Functions in Embedded Systems

We present a tool for the automatic generation of test stimuli for small numerical support functions, e.g., code for trigonometric functions, quaternions, filters, or table lookup. Our tool is based on KLEE to produce a set of test stimuli for full path coverage. We use a method of iterative deepening over abstractions to deal with floating-point values. During actual testing the stimuli exercise the code against a reference implementation. We illustrate our approach with results of experiments with low-level trigonometric functions, interpolation routines, and mathematical support functions from an open source UAS autopilot.

Metrics↗

NASA NTRS1973

Transonic flow past lifting wings.

Work conducted by Norstrud (1971) has been extended to lifting flows with the inclusion of embedded continuous supercritical regions. The approach taken follows some fundamental steps proposed for two-dimensional flows by Oswatitsch (1950). The governing integral equation is replaced by a system of nonlinear algebraic equations. The method of parametric differentiation is applied to the solution of this system of equations. The analytical analysis is discussed together with a numerical analysis in the case of a wing configuration with arbitrary thickness distribution.

Norstrud, H.↗

NASA NTRS1982

Relaxation solution of the full Euler equations

A numerical procedure for the relaxation solution of the full steady Euler equations is described. By embedding the Euler system in a second order surrogate system, central differencing may be used in subsonic regions while retaining matrix forms well suited to iterative solution procedures and convergence acceleration techniques. Hence, this method allows the development of stable, fully conservative differencing schemes for the solution of quite general inviscid flow problems. Results are presented for both subcritical and shocked supercritical internal flows. Comparisons are made with a standard time dependent solution algorithm.

Johnson, G. M.↗

NASA NTRS1982

Relaxation solution of the full Euler equations

Johnson, G. M.↗

NASA NTRS2001

Using Neural Networks in Decision Making for a Reconfigurable Electro Mechanical Actuator (EMA)

The objectives of this project were to demonstrate applicability and advantages of a neural network approach for evaluating the performance of an electro-mechanical actuator (EMA). The EMA in question was intended for the X-37 Advanced Technology Vehicle. It will have redundant components for safety and reliability. The neural networks for this application are to monitor the operation of the redundant electronics that control the actuator in real time and decide on the operating configuration. The system we proposed consists of the actuator, sensors, control circuitry and dedicated (embedded) processors. The main purpose of the study was to develop suitable hardware and neural network capable of allowing real time reconfiguration decisions to be made. This approach was to be compared to other methods such as fuzzy logic and knowledge based systems considered for the same application. Over the course of the project a more general objective was the identification of the other neural network applications and the education of interested NASA personnel on the topic of Neural Networks.

Latino, Carl D.↗

NASA NTRS1997

Mechanical signals in plant development: a new method for single cell studies

Cell division, which is critical to plant development and morphology, requires the orchestration of hundreds of intracellular processes. In the end, however, cells must make critical decisions, based on a discrete set of mechanical signals such as stress, strain, and shear, to divide in such a way that they will survive the mechanical loads generated by turgor pressure and cell enlargement within the growing tissues. Here we report on a method whereby tobacco protoplasts swirled into a 1.5% agarose entrapment medium will survive and divide. The application of a controlled mechanical load to agarose blocks containing protoplasts orients the primary division plane of the embedded cells. Photoelastic analysis of the agarose entrapment medium can identify the lines of principal stress within the agarose, confirming the hypothesis that cells divide either parallel or perpendicular to the principal stress tensors. The coincidence between the orientation of the new division wall and the orientation of the principal stress tensors suggests that the perception of mechanical stress is a characteristic of individual plant cells. The ability of a cell to determine a shear-free orientation for a new partition wall may be related to the applied load through the deformation of the matrix material. In an isotropic matrix a uniaxial load will produce a rotationally symmetric strain field, which will define a shear-free plane. Where high stress intensities combine with the loading geometry to produce multiaxial loads there will be no axis of rotational symmetry and hence no shear free plane. This suggests that two mechanisms may be orienting the division plane, one a mechanism that works in rotationally symmetrical fields, yielding divisions perpendicular to the compressive tensor, parallel to the long axis of the cell, and one in asymmetric fields, yielding divisions parallel to the short axis of the cell and the compressive tensor.

NASA Discipline Number 40-50↗

NASA NTRS2004

New Effective Material Couple--Oxide Ceramic and Carbon Nanotube-- Developed for Aerospace Microsystem and Micromachine Technologies

The prime driving force for using microsystem and micromachine technologies in transport vehicles, such as spacecraft, aircraft, and automobiles, is to reduce the weight, power consumption, and volume of components and systems to lower costs and increase affordability and reliability. However, a number of specific issues need to be addressed with respect to using microsystems and micromachines in aerospace applications--such as the lack of understanding of material characteristics; methods for producing and testing the materials in small batches; the limited proven durability and lifetime of current microcomponents, packaging, and interconnections; a cultural change with respect to system designs; and the use of embedded software, which will require new product assurance guidelines. In regards to material characteristics, there are significant adhesion, friction, and wear issues in using microdevices. Because these issues are directly related to surface phenomena, they cannot be scaled down linearly and they become increasingly important as the devices become smaller. When microsystems have contacting surfaces in relative motion, the adhesion and friction affect performance, energy consumption, wear damage, maintenance, lifetime and catastrophic failure, and reliability. Ceramics, for the most part, do not have inherently good friction and wear properties. For example, coefficients of friction in excess of 0.7 have been reported for ceramics and ceramic composite materials. Under Alternate Fuels Foundation Technologies funding, two-phase oxide ceramics developed for superior high-temperature wear resistance in NASA's High Operating Temperature Propulsion Components (HOTPC) project and new two-layered carbon nanotube (CNT) coatings (CNT topcoat/iron bondcoat/quartz substrate) developed in NASA's Revolutionary Aeropropulsion Concepts (RAC) project have been chosen as a materials couple for aerospace applications, including micromachines, in the nanotechnology lubrication task because of their potential for superior friction and wearf properties in air and in an ultrahigh vacuum, spacelike environment. At the NASA Glenn Research Center, two-phase oxide ceramic eutectics, Al2O3/ZrO2(Y2O3), were directionally solidified using the laser-float-zone process, and carbon nanotubes were synthesized within a high-temperature tube furnace at 800 C. Physical vapor deposition was used to coat all quartz substrates with 5-nm-thick iron as catalyst and bondcoat, which formed iron islands resembling droplets and serving as catalyst particles on the quartz. A series of scanning electron micrographs showing multiwalled carbon nanotubes directionally grown as aligned "nanograss" on quartz is presented. Unidirectional sliding friction eperiments were conducted at Glenn with the two-layered CNT coatings in contact with the two-phase Al2O3/ZrO2(Y2O3) eutectics in air and in ultrachigh vacuum. The main criteria for judging the performance of the materials couple for solid lubrication and antistick applications in a space environment were the coefficient of friction and the wear resistance (reciprocal of wear rate), which had to be less than 0.2 and greater than 10(exp 5) N(raised dot)/cubic millimetes, respectively, in ultrahigh vacuum. In air, the coefficient of friction for the CNT coatings in contact with Al2O3/ZrO2 (Y2O3) eutectics was 0.04, one-fourth of that for quartz. In an ultrahigh vacuum, the coefficient of friction for CNT coatings in contact with Al2O3/ZrO2 (Y2O3) was one-third of that for quartz. The two-phase Al2O3/ZrO2 (Y2O3) eutectic coupled with the two-layered CNT coating met the coefficient of friction and wear resistance criteria both in air and in an ultrahigh vacuum, spacelike environment. This material's couple can dramatically improve the stiction (or adhesion), friction, and wear resistance of the contacting surfaces, which are major issues for microdevices and micromachines.

Miyoshi, Kazuhisa↗

NASA NTRS2004

Resistivity of Carbon-Carbon Composites Halved

Carbon-carbon composites have become the material of choice for applications requiring strength and stiffness at very high temperatures (above 2000 C). These composites comprise carbon or graphite fibers embedded in a carbonized or graphitized matrix. In some applications, such as shielding sensitive electronics in very high temperature environments, the performance of these materials would be improved by lowering their electrical resistivity. One method to lower the resistivity of the composites is to lower the resistivity of the graphite fibers, and a proven method to accomplish that is intercalation. Intercalation is the insertion of guest atoms or molecules into a host lattice. In this study the host fibers were highly graphitic pitch-based graphite fibers, or vapor-grown carbon fibers (VGCF), and the intercalate was bromine. Intercalation compounds of graphite are generally thought of as being only metastable, but it has been shown that the residual bromine graphite fiber intercalation compound is remarkably stable, resisting decomposition even at temperatures at least as high as 1000 C. The focus of this work was to fabricate composite preforms, determine whether the fibers they were made from were still intercalated with bromine after processing, and determine the effect on composite resistivity. It was not expected that the resistivity would be lowered as dramatically as with graphite polymer composites because the matrix itself would be much more conductive, but it was hoped that the gains would be substantial enough to warrant its use in high-performance applications. In a collaborative effort supporting a Space Act Agreement between the NASA Glenn Research Center and Applied Sciences, Inc. (Cedarville, OH), laminar preforms were fabricated with pristine and bromine-intercalated pitch-based fibers (P100 and P100-Br) and VGCF (Pyro I and Pyro I-Br). The green preforms were carbonized at 1000 C and then heat treated to 3000 C. To determine whether the fibers in the samples were still intercalated after composite fabrication, they were subjected to X-ray diffraction. The composites containing intercalated graphite fibers showed much higher background scatter than that of pristine fibers, indicating the presence of bromine in the samples. More importantly, faint features indicative of intercalation were visible in the diffraction pattern, showing that the fibers were still intercalated.

Gaier, James R.↗

NASA NTRS2005

Reduced-Order Modeling of a Heaving Airfoil

A reduced-order model of a flapping airfoil is developed using Proper Orthogonal Decomposition (POD). The proper basis functions, developed from snapshots of full Navier-Stokes simulations, are used for a Galerkin projection of the governing equations. The resulting coupled, nonlinear ordinary di.erential equations have a low dimension because the first few basis members capture most of the energy of the flow. The reduced-order model is used to simulate heaving motions that are both similar to and different from the motion(s) used to generate the basis functions, and the errors in the model are quantified. Several methods are used to generate mode sets that can be used over a range of heaving parameters, including snapshots from one, two, and multiple Navier-Stokes simulations. As snapshots from additional simulations are added to the decomposition, the mode sets become richer and can simulate a wider range of parameter space, at some computational cost. Whereas the POD method is fully applicable in three dimensions, the simulation technique based on a body-fixed and body-fitted grid suffers large overhead when extended to three dimensions. To reduce the overhead, an embedding technique is discussed which embeds the solid wing into a fixed Cartesian grid. The wing, which can now have multiple pieces and also be flexible, is represented by a distribution of body forces. This distribution is determined to give exactly the flow around a flapping wing.

Haj-Hariri, H.↗

NASA NTRS2009

Evolutionary Computing Methods for Spectral Retrieval

A methodology for processing spectral images to retrieve information on underlying physical, chemical, and/or biological phenomena is based on evolutionary and related computational methods implemented in software. In a typical case, the solution (the information that one seeks to retrieve) consists of parameters of a mathematical model that represents one or more of the phenomena of interest. The methodology was developed for the initial purpose of retrieving the desired information from spectral image data acquired by remote-sensing instruments aimed at planets (including the Earth). Examples of information desired in such applications include trace gas concentrations, temperature profiles, surface types, day/night fractions, cloud/aerosol fractions, seasons, and viewing angles. The methodology is also potentially useful for retrieving information on chemical and/or biological hazards in terrestrial settings. In this methodology, one utilizes an iterative process that minimizes a fitness function indicative of the degree of dissimilarity between observed and synthetic spectral and angular data. The evolutionary computing methods that lie at the heart of this process yield a population of solutions (sets of the desired parameters) within an accuracy represented by a fitness-function value specified by the user. The evolutionary computing methods (ECM) used in this methodology are Genetic Algorithms and Simulated Annealing, both of which are well-established optimization techniques and have also been described in previous NASA Tech Briefs articles. These are embedded in a conceptual framework, represented in the architecture of the implementing software, that enables automatic retrieval of spectral and angular data and analysis of the retrieved solutions for uniqueness.

Terrile, Richard↗

NASA NTRS2010

Proceedings of the Second NASA Formal Methods Symposium

This publication contains the proceedings of the Second NASA Formal Methods Symposium sponsored by the National Aeronautics and Space Administration and held in Washington D.C. April 13-15, 2010. Topics covered include: Decision Engines for Software Analysis using Satisfiability Modulo Theories Solvers; Verification and Validation of Flight-Critical Systems; Formal Methods at Intel -- An Overview; Automatic Review of Abstract State Machines by Meta Property Verification; Hardware-independent Proofs of Numerical Programs; Slice-based Formal Specification Measures -- Mapping Coupling and Cohesion Measures to Formal Z; How Formal Methods Impels Discovery: A Short History of an Air Traffic Management Project; A Machine-Checked Proof of A State-Space Construction Algorithm; Automated Assume-Guarantee Reasoning for Omega-Regular Systems and Specifications; Modeling Regular Replacement for String Constraint Solving; Using Integer Clocks to Verify the Timing-Sync Sensor Network Protocol; Can Regulatory Bodies Expect Efficient Help from Formal Methods?; Synthesis of Greedy Algorithms Using Dominance Relations; A New Method for Incremental Testing of Finite State Machines; Verification of Faulty Message Passing Systems with Continuous State Space in PVS; Phase Two Feasibility Study for Software Safety Requirements Analysis Using Model Checking; A Prototype Embedding of Bluespec System Verilog in the PVS Theorem Prover; SimCheck: An Expressive Type System for Simulink; Coverage Metrics for Requirements-Based Testing: Evaluation of Effectiveness; Software Model Checking of ARINC-653 Flight Code with MCP; Evaluation of a Guideline by Formal Modelling of Cruise Control System in Event-B; Formal Verification of Large Software Systems; Symbolic Computation of Strongly Connected Components Using Saturation; Towards the Formal Verification of a Distributed Real-Time Automotive System; Slicing AADL Specifications for Model Checking; Model Checking with Edge-valued Decision Diagrams; and Data-flow based Model Analysis.

Munoz, Cesar↗

DOE OSTI2026

Latent heat thermal energy storage performance maps enabling fast & accurate building energy simulations

Thermal energy storage (TES) using phase change materials (PCMs) has gained attention as an effective approach to manage energy demand fluctuations and shift peak building loads. PCM embedded heat exchangers (PCM-HXs) offer high energy storage density and low temperature variation during phase change, being suitable for load-shifting applications. However, this component is typically evaluated using computationally expensive methods, which present significant challenges when the ultimate goal is to assess the performance of PCM-HX integrated thermal energy storage systems in the full building context. In this paper, we present a methodology to generate highly accurate and computationally efficient PCM-HX performance maps which can be easily integrated into building energy simulation tools to analyze the feasibility of space conditioning systems with latent heat PCM-based TES. The performance maps are generated using a computationally efficient PCM-HX simulation tool based on a Generalized Resistance-Capacitance Model (GRCM) which can simulate arbitrary PCM-HXs with high accuracy and significantly less computational effort compared to full CFD simulations. The methodology was verified for a case study considering a 5-ton (~17.5 kW) air-to-water heat pump-thermal energy storage system (HP-TES), which was co-simulated in Modelica for a DOE prototype small-office building in Vienna, Austria, using Spawn of EnergyPlus™. The TES performance maps provided accurate predictions of PCM-HX behavior when used as Modelica component, with deviations within 2-4% while also achieving at least 103 computational time reduction. Leveraging this faster prediction capability, four PCMs with different melting temperatures for cooling (12°C, 16°C) and heating (31°C, 36°C) were assessed to investigate their impact on system performance. This work highlights the importance of robust PCM-HX models for efficient and high-fidelity building-level simulations, presenting new opportunities for advanced control strategy development and parametric analysis of TES configurations in a computationally efficient manner

Modelica Building Simulations↗