Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Mixed finite element models for free vibrations of thin-walled beams

Simple, mixed finite element models are developed for the free vibration analysis of curved thin-walled beams with arbitrary open cross section. The analytical formulation is based on a Vlasov's type thin-walled beam theory with the effects of flexural-torsional coupling, transverse shear deformation and rotary inertia included. The fundamental unknowns consist of seven internal forces and seven generalized displacements of the beam. The element characteristic arrays are obtained by using a perturbed Lagrangian-mixed variational principle. Only C(sup o) continuity is required for the generalized displacements. The internal forces and the Lagrange multiplier are allowed to be discontinuous at interelement boundaries. Numerical results are presented to demonstrate the high accuracy and effectiveness of the elements developed. The standard of comparison is taken to be the solutions obtained by using 2-D plate/shell models for the beams.

Noor, Ahmed K.↗
NASA NTRS1997

Soft-Decision Decoding of Binary Linear Block Codes Based on an Iterative Search Algorithm

This correspondence presents a suboptimum soft-decision decoding scheme for binary linear block codes based on an iterative search algorithm. The scheme uses an algebraic decoder to iteratively generate a sequence of candidate codewords one at a time using a set of test error patterns that are constructed based on the reliability information of the received symbols. When a candidate codeword is generated, it is tested based on an optimality condition. If it satisfies the optimality condition, then it is the most likely (ML) codeword and the decoding stops. If it fails the optimality test, a search for the ML codeword is conducted in a region which contains the ML codeword. The search region is determined by the current candidate codeword and the reliability of the received symbols. The search is conducted through a purged trellis diagram for the given code using the Viterbi algorithm. If the search fails to find the ML codeword, a new candidate is generated using a new test error pattern, and the optimality test and search are renewed. The process of testing and search continues until either the MEL codeword is found or all the test error patterns are exhausted and the decoding process is terminated. Numerical results show that the proposed decoding scheme achieves either practically optimal performance or a performance only a fraction of a decibel away from the optimal maximum-likelihood decoding with a significant reduction in decoding complexity compared with the Viterbi decoding based on the full trellis diagram of the codes.

Lin, Shu↗
NASA NTRS2006

Airframe Noise Prediction by Acoustic Analogy: Revisited

The present work follows a recent survey of airframe noise prediction methodologies. In that survey, Lighthill s acoustic analogy was identified as the most prominent analytical basis for current approaches to airframe noise research. Within this approach, a problem is typically modeled with the Ffowcs Williams and Hawkings (FW-H) equation, for which a geometry-independent solution is obtained by means of the use of the free-space Green function (FSGF). Nonetheless, the aeroacoustic literature would suggest some interest in the use of tailored or exact Green s function (EGF) for aerodynamic noise problems involving solid boundaries, in particular, for trailing edge (TE) noise. A study of possible applications of EGF for prediction of broadband noise from turbulent flow over an airfoil surface and the TE is, therefore, the primary topic of the present work. Typically, the applications of EGF in the literature have been limited to TE noise prediction at low Mach numbers assuming that the normal derivative of the pressure vanishes on the airfoil surface. To extend the application of EGF to higher Mach numbers, the uniqueness of the solution of the wave equation when either the Dirichlet or the Neumann boundary condition (BC) is specified on a deformable surface in motion. The solution of Lighthill s equation with either the Dirichlet or the Neumann BC is given for such a surface using EGFs. These solutions involve both surface and volume integrals just like the solution of FW-H equation using FSGF. Insight drawn from this analysis is evoked to discuss the potential application of EGF to broadband noise prediction. It appears that the use of a EGF offers distinct advantages for predicting TE noise of an airfoil when the normal pressure gradient vanishes on the airfoil surface. It is argued that such an approach may also apply to an airfoil in motion. However, for the prediction of broadband noise not directly associated with a trailing edge, the use of EGF does not appear to offer any advantages over the use of FSGF at the present stage of development. It is suggested here that the applications of EGF for airframe noise analysis be continued. As an example pertinent to airframe noise prediction, the Fast Scattering Code of NASA Langley is utilized to obtain the EGF numerically on the surface of a three dimensional wing with a flap and leading edge slat in uniform rectilinear motion. The interpretation and use of these numerical Green functions are then discussed.

Farassat, F.↗
NASA NTRS2011

180-GHz Interferometric Imager

A 180-GHz interferometric imager uses compact receiver modules, combined high- and low-gain antennas, and ASIC (application specific integrated circuit) correlator technology, enabling continuous, all-weather observations of water vapor with 25-km resolution and 0.3-K noise in 15 minutes of observation for numerical weather forecasting and tropical storm prediction. The GeoSTAR-II prototype instrument is broken down into four major subsystems: the compact, low-noise receivers; sub-array modules; IF signal distribution; and the digitizer/correlator. Instead of the single row of antennas adopted in GeoSTAR, this version has four rows of antennas on a coarser grid. This dramatically improves the sensitivity in the desired field of view. The GeoSTAR-II instrument is a 48-element, synthetic, thinned aperture radiometer operating at 165-183 GHz. The instrument has compact receivers integrated into tiles of 16 elements in a 4x4 arrangement. These tiles become the building block of larger arrays. The tiles contain signal distribution for bias controls, IF signal, and local oscillator signals. The IF signals are digitized and correlated using an ASIC correlator to minimize power consumption. Previous synthetic aperture imagers have used comparatively large multichip modules, whereas this approach uses chip-scale modules mounted on circuit boards, which are in turn mounted on the distribution manifolds. This minimizes the number of connectors and reduces system mass. The use of ASIC technology in the digitizers and correlators leads to a power reduction close to an order of magnitude.

Kangaslahti, Pekka P.↗
NASA NTRS2016

Atmospheric Diabatic Heating in Different Weather States and the General Circulation

Analysis of multiple global satellite products identifies distinctive weather states of the atmosphere from the mesoscale pattern of cloud properties and quantifies the associated diabatic heating/cooling by radiative flux divergence, precipitation, and surface sensible heat flux. The results show that the forcing for the atmospheric general circulation is a very dynamic process, varying strongly at weather space-time scales, comprising relatively infrequent, strong heating events by ''stormy'' weather and more nearly continuous, weak cooling by ''fair'' weather. Such behavior undercuts the value of analyses of time-averaged energy exchanges in observations or numerical models. It is proposed that an analysis of the joint time-related variations of the global weather states and the general circulation on weather space-time scales might be used to establish useful ''feedback like'' relationships between cloud processes and the large-scale circulation.

Cooling↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
DOE OSTI2024

CEBAF Injector for K Long Beam Conditions

The Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab concurrently operates four experimental Halls with distinct bunch charge specifications and repetition rates. Numerous critical beam parameters within CEBAF are configured in the injector, some remaining unchanged throughout the accelerator. Consequently, the injector plays a crucial role in determining final beam characteristics, including bunch structure, beam sizes, bunch lengths, energy spread, and beam transmission. The Jefferson Lab KL experiment is scheduled to take place at CEBAF in Hall D, featuring a much lower bunch repetition rate of 7.80 MHz or 15.59 MHz, below the nominal values of 249.5 MHz or 499 MHz. Although the proposed average current of 5 ?A or 10 ?A is low compared to the maximum CEBAF cur- rent of approximately 180 ?A, the corresponding bunch charge is unusually high for CEBAF injector operation. This study focuses on the behavior of low-repetition-rate, high-bunch- charge (0.32 to 0.64 pC) beams in the CEBAF injector. We investigated the evolution and transmission of low-charge beams to space-charge dominated high-charge beams in the front end of the CEBAF injector for two configurations: the pre-existing CEBAF Phase 1 injector upgrade, operated at 130 kV, and the existing CEBAF Phase 2 injector upgrade, operated at 140 kV, 180 kV, and to be operated at 200 kV. The electron beam through the CEBAF injector is characterized using beam dynamics simulations and comparisons with the available measurements performed at 130 kV. Multi-objective genetic optimizations of the CEBAF injector were performed to determine the operating magnetic elements and RF settings for the evolution and transmission of low, moderate, and high charge beams in the CEBAF injector at 180 kV and 200 kV DC gun voltages. Subsequently, simulations at the same voltages were conducted to obtain the beam characteristics at the front end of the CEBAF injector. The laser spot size and laser pulse length at the cathode were varied to observe their effects on beam transmission in the injector at different voltages (130 kV, 180 kV, and 200 kV). Experimental studies at 130 kV, 140 kV, and 180 kV validate the simulations. Beam study measurements are carried out using EPICS tools, while optimizations and simulations are facilitated by General Particle Tracer. Based on the findings, optimal parameters for the upcoming Jefferson Lab KL experiment are proposed, utilizing a lower repetition rate and higher bunch charge

Pokharel, Sunil↗
DOE OSTI2024

Model-Based Approaches to Generate Knowledge from Data in a Plant Reliability Context

One challenge that nuclear power plant system engineers are facing is continuous generation of an extremely large amount of equipment reliability (ER) data. These data elements come in textual (e.g., condition reports) and numeric (e.g., generated by monitoring systems) forms. They provide system engineers with valuable insights and information by discovering anomalous behaviors or degradation trends, identifying possible causes behind such behaviors and trends, and predicting their direct consequences. This paper directly targets the knowledge generation from ER data by putting “data into context.” We employ model-based system engineering (MBSE) of systems and assets to represent and capture their architecture and functional (i.e., cause-effect) relations. ER data elements are processed by first identifying which of the developed MBSE elements they are referring to. This task is harder for textual data since the information contained in issue or maintenance reports needs to be “understood” by a computational tool. We called this process “knowledge extraction” since our methods extract knowledge from textual data. Last, once numeric and textual ER data elements have been processed and “understood,” we discover possible cause-effect relations among them. This is performed by observing whether a logical connection through the MBSE models exists, and if there is a temporal relationship among them. The logic and temporal are the two main ingredients to perform “machine reasoning” from ER data.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS1982

The CMC:3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 1: Theoretical

An order-of-magnitude analysis of the subsonic three dimensional steady time averaged Navier-Stokes equations, for semibounded aerodynamic juncture geometries, yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement of first order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. Numerical results document accuracy, convergence, and utility of the developed finite element algorithm, and the CMC:3DPNS computer code applied to an idealized wing-body juncture region. Additional results document accuracy aspects of the algorithm turbulence closure model.

Baker, A. J.↗
NASA NTRSDate not supplied

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗
NASA NTRS2023

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗
NASA NTRS1973

A field investigation and numerical simulation of coastal fog

A field investigation of the microphysical and micrometeorological features of fogs occurring near Los Angeles and Vandenberg, California was conducted. Observations of wind speed and direction, temperature, dew point, vertical wind velocity, dew deposition, drop-size distribution, liquid water content, and haze and cloud nucleus concentration were obtained. These observations were initiated in late evening prior to fog formation and continued until the time of dissipation in both advection and radiation fogs. Data were also acquired in one valley fog and several dense haze situations. The behavior of these parameters prior to and during fog are discussed in detail. A two-dimensional numerical model was developed to investigate the formation and dissipation of advection fogs under the influence of horizontal variations in surface temperature. The model predicts the evolution of potential temperature, water vapor content, and liquid water content in a vertical plane as determined by vertical turbulent transfer and horizontal advection. Results are discussed from preliminary numerical experiments on the formation of warm-air advection fog and dissipation by natural and artificial heating from the surface.

Mack, E. J.↗
NASA NTRS2011

Recent Changes in the Greenland Ice Sheet as Seen from Space

Many changes in the Greenland Ice Sheet have been reported in the recent scientific literature and have been attributed to various responses of the ice sheet due to regional (and global) warming. Because melting of the ice sheet would contribute approximately 7 m to sea-level rise, the lives and habitat of hundreds of millions of people worldwide would be directly and indirectly affected if continued ice-sheet melting occurs. As mean-annual global temperatures have increased, there has been an increasing focus on studying the Greenland Ice Sheet using available satellite data, and numerous expeditions have been undertaken. Regional "clear-sky" surface temperature increases since the early 1980s in the Arctic, measured using Advanced Very High Resolution Radiometer (AVHRR) infrared data, range from 0.57+/-0.02 C to 0.72+/-0.10 C per decade. Arctic warming has important implications for ice-sheet mass balance because much of the periphery of the Greenland Ice Sheet is already near O C during the melt season, and is thus vulnerable to more extensive melting if temperatures continue to increase. An increase in melting of the ice sheet would accelerate sea-level rise, an issue of increasing concern to billions of people worldwide. The surface temperature of the ice sheet has been studied in even greater detail using Moderate-Resolution Imaging Spectroradiometer (MODIS) data in the six individual drainage basins as well as for the ice sheet as a whole. Surface temperature trends in the decade of the 2000s have not been strong, according to the MODIS measurements. In addition to surface-temperature increases over the last few decades as measured by AVHRR, other changes have been observed such as accelerated movement of many of Greenland's outlet glaciers and sudden draining of supraglacial lakes. Decreasing mass of the ice sheet since (at least) 2002 has been measured using Gravity Recovery and Climate Experiment (GRACE) data, along with an build-up of ice at the higher elevations and a decrease of ice at the lower elevations as measured using airborne Lidar and Ice, Cloud and Land Elevation Satellite (ICESat) data. The seminar will address the above issues using a variety of NASA satellite data and ground observations.

Hall, Dorothy K.↗
DOE OSTI2023

Report on FY 2023 Research and Development on Specially Designed Creep-fatigue Experiments on Alloy 617 in Support of Improving Creep-fatigue Evaluation Approaches

Experimental and numerical studies in support of developing the integrated Elastic–Perfectly Plastic (EPP) plus Simplified Model Test (SMT) design methodology, referred to as the EPP+SMT method, continued in FY 2023. This report focuses on the methods for extrapolating the EPP+SMT creep-fatigue life curves at long hold times and low strain ranges at elevated temperatures. In this work, the available uniaxial creep-fatigue failure data on Alloy 617 at temperatures of 950°C and 850°C were analyzed to provide a guidance on the development of the extrapolation method and the creep-fatigue failure criteria. A viscoplastic constitutive model for Alloy 617 was adopted to extrapolate the mechanical responses to low strain ranges and long hold times. A set of design curves of Alloy 617 at temperatures of 950°C, 850°C, and 800°C with tensile hold times of 1 hr, 100 hr, and 1,000 hr are developed. Furthermore, creep-fatigue testing on two notch specimen geometries, shallow-notch and sharp V-notch, on Alloy 617 was performed 950°C to understand the multiaxial stress relaxation behavior. The experimental and numerical results on the notch specimens were compared with those on the standard uniaxial smooth bar specimens. The effect of multi-axial stress state combined with elastic follow-up on the stress relaxation behavior was investigated in this report.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Numerical Modeling of CO2 Sequestration within a Five-Spot Well Pattern in the Morrow B Sandstone of the Farnsworth Hydrocarbon Field: Comparison of the TOUGHREACT, STOMP-EOR, and GEM Simulators

The objectives of this study were (1) to assess the fate and impact of CO2 injected into the Morrow B Sandstone in the Farnsworth Unit (FWU) through numerical non-isothermal reactive transport modeling, and (2) to compare the performance of three major reactive solute transport simulators, TOUGHREACT, STOMP-EOR, and GEM, under the same input conditions. The models were based on a quarter of a five-spot well pattern where CO2 was injected on a water-alternating-gas schedule for the first 25 years of the 1000 year simulation. The reservoir pore fluid consisted of water with or without petroleum. The results of the models have numerous broad similarities, such as the pattern of reservoir cooling caused by the injected fluids, a large initial pH drop followed by gradual pH neutralization, the long-term persistence of an immiscible CO2 gas phase, the continuous dissolution of calcite, very small decreases in porosity, and the increasing importance over time of carbonate mineral CO2 sequestration. The models differed in their predicted fluid pressure evolutions; amounts of mineral precipitation and dissolution; and distribution of CO2 among immiscible gas, petroleum, formation water, and carbonate minerals. The results of the study show the usefulness of numerical simulations in identifying broad patterns of behavior associated with CO2 injection, but also point to significant uncertainties in the numerical values of many model output parameters.

STOMP↗
NASA NTRS2000

Numerical Calculation of the Drag Force Acting on an Insoluble Particle Moving in Front of a Solidifying Interface

The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. Whether a particle will be engulfed or continuously pushed by SLI is essentially determined by the balance of forces acting on the particle. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that not only the steady-state approach is inappropriate to model the interaction process but also that even at steady-state the theoretical solution underestimates the drag force. A regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gave a relationship of the form F(sub d) (exp num) = (square root of 3) pi eta V(sub p)((R(sub p) / d) (exp 10 gamma/3)) as compared to the theoretical solution F(sub d) (exp theor) = 3 X (square root of 2) pi eta V(sub p)((R(sub p) / d) (exp 3/2)) where F(sub d) is the drag force, eta is the dynamic viscosity of the fluid, V(sub p) is the particle velocity, R(sub p) is the particle radius, gamma is Euler's constant (gamma is approximately equal to 0.577), and d is the width of the gap between the particle and the SLI. Model validation as well as the influence of the interface shape on the value of F(sub d) will also be discussed.

Catalina, Adrian V.↗
NASA NTRS1989

Review of FD-TD numerical modeling of electromagnetic wave scattering and radar cross section

Applications of the finite-difference time-domain (FD-TD) method for numerical modeling of electromagnetic wave interactions with structures are reviewed, concentrating on scattering and radar cross section (RCS). A number of two- and three-dimensional examples of FD-TD modeling of scattering and penetration are provided. The objects modeled range in nature from simple geometric shapes to extremely complex aerospace and biological systems. Rigorous analytical or experimental validatons are provided for the canonical shapes, and it is shown that FD-TD predictive data for near fields and RCS are in excellent agreement with the benchmark data. It is concluded that with continuing advances in FD-TD modeling theory for target features relevant to the RCS problems and in vector and concurrent supercomputer technology, it is likely that FD-TD numerical modeling will occupy an important place in RCS technology in the 1990s and beyond.

Taflove, Allen↗
DOE OSTI2020

Linear thermomagnetic energy harvester for low-grade thermal energy harvesting

Low-grade thermal energy, either from waste heat or from natural resources, constitutes an enormous energy reserve that remains to be fully harvested. Harvesting low-grade heat is challenging because of the low Carnot efficiency. Among various thermal energy harvesting mechanisms available for capturing low-grade heat (temperature less than 100 degrees C), the thermomagnetic effect has been found to be quite promising. In this study, we demonstrate a scalable thermomagnetic energy harvester architecture that exhibits 140% higher power density compared to the previously published spring-mass designs. The alternating force required to oscillate the thermomagnetic mass is generated through the interaction between two magnetic forces in opposite directions. We employed numerical modeling to illustrate the behavior of a thermomagnetic device under different operating conditions and to obtain the optimal hot-side and cold-side temperatures for continuous mode operations. A miniaturized thermomagnetic harvester was fabricated and experiments were conducted to systematically evaluate the performance. The prototype was found to exhibit an oscillation frequency of 0.33 Hz, a work output of 0.6 J/kg/cycle, and a power density of 0.2 W/kg of gadolinium under the temperature difference of 60 K.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗