Simple 1D Model Research Articles

Numerical approach and experimental investigation of a stratified hot water storage tank for domestic hot water production

Thermal stratification within storage tanks plays a critical role in the efficiency of thermal energy storage systems. In this paper, one-dimensional numerical model of a stratified storage tank is developed to simulate the thermal stratification process within a storage tank during the charging and the discharging processes, capturing the temperature distribution inside the storage tank over time. This work investigates the evolution and transient behavior of stratified flows during the heating and cooling process within sensible-heat energy storage systems under ReD=167 and ReD=616 for heating process. For that purpose, an experimental campaign was performed in an instrumented water storage tank to measure the stratification process when a hot water round jet was injected into a slender cylindrical tank filled with an initially uniform hydrostatic water field stabilized at a lower temperature. This information was used to develop a simple 1D model able to be applied during long-term charging/discharging transients of these kind of systems. The 1D model developed was validated against the experiments performed following prototypical charging/discharging tests according to the standard EN16147 and can be used to describe the energy performance of this kind of systems during long-term charging/discharging sequences. The present model provides a reduction of about 5 % of the error of the time of the water cooling process with respect to other models in the literature.

Room temperature electron beam sensitive viscoplastic response of ultra-ductile amorphous olivine films

The mechanical properties of amorphous olivine (a-olivine) deformed at room temperature are investigated in situ in a TEM under uniaxial tension using a Push-to-Pull (PTP) device. Thin films of a-olivine were produced by pulsed laser deposition (PLD). With or without electron irradiation, a-olivine films deform plastically, with a gradual transition that makes impossible the determination of a precise threshold. The strength attains values up to 2.5 GPa. The increasing strain-rate in load control results in an apparent softening with stress drop. The fracture strain reaches values close to 30 % without e-beam irradiation. Under electron illumination at 200 kV, the strength is lower, around 1.7 GPa, while higher elongations close to 36 % are obtained. Alternating beam-off and beam-on sequences lead to exceptionally large fracture strains equal to 68 % at 200 kV and 139 % at 80 kV. EELS measurements were performed to characterize the interaction between the electron beam and a-olivine. At a voltage of 80 kV, radiolysis accompanied by oxygen release dominates whereas at high voltage (300 kV) the interaction is dominated by knock-on type defects. Radiolysis is also the main interaction mechanism at 200 kV with low exposition which corresponds to most of our in situ TEM deformation experiments. To interpret the mechanical data, a simple 1D model has been developed to rationalize the load transfer between the PTP device and the specimen. The strain-rate sensitivity is 6 to 10 times higher when a-olivine is deformed under electron irradiation.

2D PIC modeling of the helical scrape-off layer current driven by hybrid divertor biased targets in tokamak plasmas

Abstract The heat flux control of the divertor plate via strike-point splitting generated by biased targets was proposed in the HL-2A tokamak (Cui et al 2021 Fusion Eng. Des. 173 112963). To understand the helical scrape-off layer (SOL) currents driven by hybrid biasing, two SOL current models (model A and B) are employed. Model A is a simplified 2D model that focuses on investigating the effect of biasing on the sheath and elucidating the fundamental physical mechanism of bias-driven SOL current paths. The potential, charge density, electric field and current densities are calculated. Model B takes into account the actual tokamak geometry and computes the resonant magnetic perturbations (RMPs) generated by bias-driven linear decay currents. Additionally, strike-point splitting is observed in the HL-2A tokamak, indicating that the SOL currents generated by hybrid biasing are capable of generating strong RMPs and consequently influence the magnetic topology. These results confirm the potential of heat/particle flux control by hybrid divertor biased targets.

How Landmass Distribution Influences the Atmospheric Dynamics of Tidally Locked Terrestrial Exoplanets

Interpretation of the ongoing efforts to simulate the atmospheres of potentially habitable terrestrial exoplanets requires that we understand the underlying dynamics and chemistry of such objects to a much greater degree than 1D or even simple 3D models enable. Here, for the tidally locked habitable-zone planet TRAPPIST-1e, we explore one effect which can shape the dynamics and chemistry of terrestrial planets: the inclusion of an Earth-like land–ocean distribution with orography. To do this we use the Earth-system model WACCM6/CESM2 to run a pair of TRAPPIST-1e models with N2–O2 atmospheres and with the substellar point fixed over either land or ocean. The presence of orography shapes atmospheric transport, and in the case of Earth-like orography, breaks the symmetry between the Northern and Southern Hemispheres which was previously found in slab ocean models. For example, peak zonal jet speeds in the Southern Hemisphere are 50%–100% faster than similar jets in the Northern Hemisphere. This also affects the meridional circulation, transporting equatorial material toward the south pole. As a result we also find significant changes in the atmospheric chemistry, including the accumulation of potentially lethal quantities of ozone at both the south pole and the surface. Future studies which investigate the effects of landmass distribution on the dynamics of exoplanetary atmospheres should pay close attention to both the dayside land fraction as well as the orography of the land. Simply modeling a flat landmass will not give a complete picture of its dynamical impact.

Improving a 1D Hydraulic Model to Include Bridges as Internal Boundary Conditions

The paper describes the implementation of internal boundary conditions in the 1D ORSADEM hydraulic model to simulate the effect of a hydraulic in-line structure. The proposed model introduces a simplified representation of the bridge geometry by imposing an equivalent narrowing, computed according to the opening size and characteristics, combined with the mass and energy balance at the structure. The model is then applied to a series of experimental tests concerning the propagation of shock waves through different types of bridges, representing different flow conditions, from free surface flow to overflow. The tests are also simulated with the original 1D ORSADEM model, including the standard head losses and the cross-section narrowing due to the presence of a structure. The comparison with the experimental measurements shows that the proposed model can simulate the shock wave flow through the bridges with a higher accuracy than the standard formulation. These findings highlight the possibility of properly evaluating the backwater effect at bridges even with a simple 1D model if the physical narrowing of the cross-section is modeled.

Tetrahedral grids in Monte Carlo radiative transfer

Context. To understand the structures of complex astrophysical objects, 3D numerical simulations of radiative transfer processes are invaluable. For Monte Carlo radiative transfer, the most common radiative transfer method in 3D, the design of a spatial grid is important and non-trivial. Common choices include hierarchical octree and unstructured Voronoi grids, each of which has advantages and limitations. Tetrahedral grids, commonly used in ray-tracing computer graphics, can be an interesting alternative option. Aims. We aim to investigate the possibilities, advantages, and limitations of tetrahedral grids in the context of Monte Carlo radiative transfer. In particular, we want to compare the performance of tetrahedral grids to other commonly used grid structures. Methods. We implemented a tetrahedral grid structure, based on the open-source library TetGen, in the generic Monte Carlo radiative transfer code SKIRT. Tetrahedral grids can be imported from external applications or they can be constructed and adaptively refined within SKIRT. We implemented an efficient grid traversal method based on Plücker coordinates and Plücker products. Results. The correct implementation of the tetrahedral grid construction and the grid traversal algorithm in SKIRT were validated using 2D radiative transfer benchmark problems. Using a simple 3D model, we compared the performance of tetrahedral, octree, and Voronoi grids. With a constant cell count, the octree grid outperforms the tetrahedral and Voronoi grids in terms of traversal speed, whereas the tetrahedral grid is poorer than the other grids in terms of grid quality. All told, we find that the performance of tetrahedral grids is relatively poor compared to octree and Voronoi grids. Conclusions. Although the adaptively constructed tetrahedral grids might not be favourable in most media representative of astrophysical simulation models, they still form an interesting unstructured alternative to Voronoi grids for specific applications. In particular, they might prove useful for radiative transfer post-processing of hydrodynamical simulations run on tetrahedral or unstructured grids.

Massive fabrication of functional hepatic cancer spheroids by micropatterned GelMA hydrogel chip for drug screening

Since hepatic cancer incidence and mortality continue to grow worldwide, it is necessary to develop the biomimetic tumor models for drug development and tumor therapeutics. Cellular spheroids as an excellent simple 3D model can bridge the gap between 2D cell culture and live tissue. In this study, we proposed a biological methacrylated gelatin (GelMA) hydrogel-based microplatform for the massive generation of hepatocellular spheroids and downstream investigation of drug resistance. Micropatterned GelMA hydrogel microwell chip (GHM-chip) with tunable array was easily achieved in standard 24-culture well plates through the micro-molding fabrication strategy. The fabricated GHM-chip induced multicellular self-assembly behavior within the defined topography and further formed spheroidal structure. By regulating cell seeding density and designing microwell size, uniform hepatic cancer spheroids with tunable diameters were obtained in a simplicity, stability and controllable manner. In addition, the screening chemotherapy study of anti-cancer drug was completed through non-destructive recovery of spheroids from the GHM-chip. Beyond that, the recovered functional spheroids have potential application value in various biomedical fields such as tumor biology, pharmacology, and tissue microengineering. Finally, the proposed GHM-chip incorporated into standard cell culture plates with easy to manufacture and operate properties, may be an efficient culture microplatform for cancer research applications.

COMPARING SCIENTIFIC COMPUTING ENVIRONMENTS FOR SIMULATING 2D NON-BUOYANT FLUID PARCEL TRAJECTORY UNDER INERTIAL OSCILLATION: A PRELIMINARY EDUCATIONAL STUDY

This study presents a preliminary numerical investigation of the two-dimensional trajectory of a non-buoyant fluid parcel subjected to inertial oscillations and abrupt external forcing events. The simulations were implemented using Python, GNU Octave, R, Julia, and Fortran open-source scientific computing environments. By running 1,000 iterations in each environment, we evaluated the computational performance of these languages in tackling this idealized problem. The results, visualized through static plots and animation, validate the numerical model's ability to represent the fundamental physics governing fluid motion. Statistical analysis using the Kruskal-Wallis test and Dunn's post-hoc test with Bonferroni correction revealed that Fortran exhibits significantly faster execution times than other environments. However, the choice of programming language should also consider factors such as coding expertise, library availability, and scalability requirements. This study focuses on the performance of scientific computing environments within each language rather than the languages themselves. The observed execution times should be interpreted in the context of the specific environments used, as they often leverage optimized libraries written in lower-level languages. Despite the limitations of this work, such as the simplified 2D model and the use of a single hardware configuration, this study provides valuable insights into selecting appropriate computational tools. It contributes to educational resources for teaching idealized fluid dynamics models. Future studies could explore more complex scenarios, a more comprehensive range of programming environments, and the impact of different numerical schemes and physical parameterizations.

Towards strain gauge 2.0: Substituting the electric resistance routinely deposited on polyimide film by the optimal pattern for full‐field strain measurement

AbstractThe checkerboard constitutes the best pattern for full‐field strain measurement because it maximizes image gradient. In the experimental mechanics community, employing this pattern is currently strongly limited because depositing it on the surface of specimens raises practical difficulties. A recent study shows, however, that it is technically possible by using a laser engraver. The present paper aims to push this solution forward by printing a checkerboard pattern on a thin polymeric film and then gluing the resulting laser‐engraved film on the specimen surface. The underlying idea is to separate the manufacturing process of this optical strain gauge on the one hand and its use on the other hand to help spread this strain measuring tool in the experimental mechanics community. The polymeric film employed here is the same as that used in the manufacturing process of classic electrical gauges, so one can rely on the know‐how of classic strain gauge bonding to glue this optical strain gauge on the specimen surface. The main difference between the proposed tool and classic electrical gauges is that the strain field beneath the polymeric support is measured instead of localized strain values. The paper is a proof of concept for this strain field measuring tool. The manufacturing and bonding processes are described in the paper. The localized spectrum analysis, a spectral technique developed for processing images of periodic patterns, is used to retrieve the strain fields from checkerboard images. Through two complementary examples, we show the ability of this new type of strain gauge to detect and quantify local details in the strain field beneath. A simplified 1D model is also proposed to assess the minimum width of the strain peak that can reliably be measured with this technique.

Numerical simulations of exocomet transits: Insights from β Pic and KIC 3542116

In recent years, the topic of existence and exploration of exocomets has been gaining increasing attention. The asymmetrical decrease in the star’s brightness due to the passage of a comet-like object in front of the star was successfully predicted. It was subsequently confirmed on the basis of the light curves of stars observed by Kepler and TESS orbital telescopes. Since then, there have been successful attempts to fit the asymmetrical dips observed in the stars’ light curves utilizing a simple 1D model of an exponentially decaying optically thin dust tail. In this work, we propose fitting the photometric profiles of some known exocomet transits based on a Monte Carlo approach to build up the distribution of dust particles in a cometary tail. As the exocomet prototypes, we used the physical properties of certain Solar System comets belonging to the different dynamical groups and moving at heliocentric distances of 0.6 au, 1.0 au, 5.0 au, and 5.5 au. We obtained a good agreement between the observed and modeled transit light curves. We also show that the physical characteristics of dust particles, such as the particle size range, the power index of dust size distribution, the particle terminal velocity, and distance to the host star affect the shape of the transit light curve, while the dust productivity of the comet nucleus and the impact parameter influence its depth and duration. The estimated dust production rates of the transiting exocomets are at the level of the most active Solar System comets.

Research on the Influence of Cone Angle on the Sealing Performance of Pipeline Joint

In order to evaluate the sealing performance of the pipe connectors, we developed a simplified 3D model by adjusting the cone angle of the nozzle and analyzing the contact pressure of the sealing structure, the width of the contact surface, and the equivalent stress as the response quantities. Based on a finite element analysis approach, we found that increasing the cone angle of the nozzle improves the sealing performance, but may negatively affect structural reliability and maintenance. As such, the optimal range for the nozzle cone angle is relatively reasonable between 55° and 60°. Furthermore, we use a sealed interface leakage model to calculate the joint leakage rate under theoretical conditions.

3D modeling and measurement of HTS tape stacks in linear superconducting magnetic bearings

Superconducting magnetic bearings (SMBs) are among the possible new technologies to be incorporated in maglev vehicles. Stacks of high-temperature superconductor (HTS) tapes can be used as an alternative to bulks, because stacks offer better mechanical properties, a better thermal conductivity and a simpler production process. Numerical modeling has been employed as a cost-effective, fast and reliable tool for improving the performance of SMBs. Several scenarios can be simulated with fast and relatively simple 2D models; however, in some cases using 3D models is inevitable. In this study, we use a full 3D model to solve the problem of magnetization of the tape stacks and obtaining the hysteresis force loop between a permanent magnet and the tape stacks. For this purpose, we employ an energy minimization-based method called minimum electromagnetic entropy production in 3D, combined with a homogenization technique and the Jc(B,θ) dependence of the HTS tape as input. The modeling results agree very well with the experiment both in the zero-field cooled and field-cooled conditions. The presented approach offers significant computational advantages, delivering faster and more efficient results compared to previously proposed 3D methods.

Calculation and prediction of divertor detachment via impurity seeding by using one-dimensional model

Abstract Achieving the divertor detachment helps to alleviate excessive heat load and sputtering problems on the target plates, thereby prolonging the lifetime of divertor components for fusion devices. In order to provide a rapid but relatively reliable prediction of plasma parameters along the flux tube for future devices design, a one-dimensional (1D) modeling code for the impurity seeded detached divertor operational point has been developed based on the Python language, which is a fluid model based on the previous work.[1] The experimental observation of the onset of divertor detachment by neon (Ne) and argon (Ar) seeding in EAST is well reproduced by using the 1D modeling code. The comparison of the 1D modeling and 2D simulation by SOLPS-ITER code for CFETR detachment operation with Ne and Ar seeding also shows good consistency. The predictions of the radiative power loss and requirement of the corresponding impurity concentration for achieving divertor detachment via different impurity seeding for the high heating power scenarios on EAST and CFETR Phase II by the 1D model are also presented. Based on the predictions, the optimized parameter space for divertor detachment operation on EAST and CFETR has also been determined. Such a simple but reliable 1D model can provide a reasonable parameter inputs for a detailed and accurate analysis by 2D or 3D modeling tools through rapid parameter scanning.

Numerical Study of a Latent Heat Storage System’s Performance as a Function of the Phase Change Material’s Thermal Conductivity

The thermal conductivities of most commonly used phase change materials (PCMs) are typically fairly low (in the range of 0.2 to 0.4 W/m·K) and are an important consideration when designing latent heat energy storage systems (LHESSs). Because of that, material scientists have been asking the following question: “by how much does the thermal conductivity of a PCM needs to be increased to positively impact the design and performance of a LHESS?” The answer to this question is not straightforward as the performance of a LHESS depends on the PCM’s thermal conductivity, other PCM thermophysical properties, the type of heat exchange system geometry used, the mode of operation, and the targeted power/energy storage of the LHESS. This paper presents work related to this question through a numerical study based on a simplified 2D model of an experimental setup studied previously in the authors’ laboratory. A model created in COMSOL Multiphysics, based on conduction and accounting for the solid-liquid phase change process, was initially validated against experimental results and then used to study the impact of the PCM’s thermal conductivity (dodecanoic acid) on the discharging power of the LHESS. The results show that even increasing the thermal conductivity of the PCM by a factor of 50 only leads to a maximum instantaneous power increase by a factor of 2 or 3 depending on the LHESS configurations.

Numerical modeling of magnetic cores with combined electrical steel grades

PurposeThis paper aims to investigate the impact of combining grain-oriented electrical steel (GOES) grades on specific iron losses and the flux density distribution within a single-phase magnetic core.Design/methodology/approachThis paper presents the results of finite-element method (FEM) simulations investigating the impact of mixing two different GOES grades on losses of a single-phase magnetic core. The authors used different models: a 3D model with a highly detailed geometry including both saturation and anisotropy, as well as a simplified 2D model to save computation time. The behavior of the flux distribution in the mixed magnetic core is analyzed. Finally, the results from the numerical simulations are compared with experimental results.FindingsThe specific iron losses of a mixed magnetic core exhibit a nonlinear decrease with respect to the GOES grade with the lowest losses. Analyzing the magnetic core behavior using 2D and 3D FEM shows that the rolling direction of the GOES grades plays a critical role on the nonlinearity variation of the specific losses.Originality/valueThe novelty of this research lies in achieving an optimum trade-off between the manufacturing cost and the core efficiency by combining conventional and high-performance GOES grade in a single-phase magnetic core.

Topology optimization for flow machine rotor design considering resonance and low mass density flows

When designing structures that take part in rotating fluid flow machines, one important effect that may be taken into account is resonance, which is evaluated through the modal analysis. Under rotating conditions, the formulation and behavior of the structure changes, requiring a new formulation. Thus, this work formulates the fluid flow topology optimization for rotating resonance, while also considering a low mass density (compressible or incompressible) fluid flow. Therefore, three different effects are considered in the fluid flow-based design: from the modal analysis, from the solid structure, and from the low mass density fluid flow. The topology optimization problem is formulated for a simplified 2D model of a rotating fluid flow device, and a nodal design variable. The topology optimization is performed through an integer linear programming algorithm. Numerical examples are presented for the developed topology optimization formulation.

Development of Thermal Performance Prediction for Large Planar Military Antenna with Multi-Cooling Channels

Large planar military antenna boasts a range of electrical components, including TRA(Transmit-Receive Assembly), signal processors, etc. which engage in computations and calculations. These processes generate a significant amount of heat, leading to unforeseen consequences for the equipment. To mitigate these adverse effects, it's imperative to implement a cooling system that can effectively reduce heat-related issues. Given the antenna's intricate nature and the multitude of components it houses, a two-step estimation process is necessary. The first step involves a comprehensive model calculation to determine the total flow characteristics, while the second step entails a thermal analysis of individual TRA set. In this study, we depicted an antenna set using simplified 3D models of its components, considering their material and thermal properties. The sequential analysis process facilitated the calculation of branched flow rates, providing insights into the individual TRA. This approach also allowed us to design a cooling system for the TRA set, assessing its thermal stability in high-temperature environments. To ensure the optimal performance of TRA, breaking down the analysis into stages based on the cooling system's structure can assist operators in predicting numerical results more effectively.

Heating and Evaporation of Sessile Droplets: Simple and Advanced Models.

New advanced and simple two-dimensional (2D) models of sessile droplet heating and cooling and evaporation are suggested. In contrast to the earlier developed one-dimensional (1D) model, based on the assumption that heat supplied from the supporting surface is homogeneously and instantaneously spread throughout the droplet, both new 2D models consider the spatial distribution of this heat. The advanced 2D model is based on the numerical solution to the equations of conservation of mass, momentum, vapor mass fraction, and energy with standard boundary and initial conditions, using COMSOL Multiphysics code. Simple 2D and 1D models assume that droplets retain their truncated spherical shapes during the evaporation process. In the 1D model, the analytical solution to the 1D heat conduction equation inside the droplet is implemented into a numerical code. In the simple 2D model, the 2D version of this equation is solved numerically using COMSOL Multiphysics code. Droplet deformation, temperature gradients along the droplet surface, and the Marangoni effect are not considered in this model. The predictions of all three models are validated using in-house experimental data obtained from studies of sessile droplets of distilled water with initial volumes of 5.2, 3.2, and 2.2 μL, at an ambient temperature of 298.15 K, and at atmospheric pressure. The observed values of normalized droplet radii squared are shown to be close to those predicted by all three models. This allows us to recommend the application of the simplest 1D model for predicting this parameter. The time dependences of the droplet average surface temperature predicted by the advanced 2D model are shown to be close to those observed experimentally. The simple 2D and 1D models can correctly predict the initial rapid decrease in droplet average surface temperature followed by its gradual increase, in agreement with experimental data.

Efficiency of timber-concrete composite floors with discrete perfobond connections: A numerical study

To evaluate the behaviour of a single timber-concrete composite (TCC) T-beam, a simplified 2D finite element (FE) model in Abaqus is developed and validated using experimental tests and the Gamma method given in Annex B of EN 1995–1-1. The shear connection used in the beam consist of perfobond connectors, which are adhesively bonded to the timber beam. In the developed model, beam and shell finite elements have been used to simulate/reproduce the timber beam and the concrete slab, respectively. These have then been connected using connector elements which incorporate the mechanical behaviour of the tested shear connections, obtained in a previous research. After the successful validation of the FE model, this is used to further analyse different parameters such as: i) different number of connections (varying from 8 to 20), ii) different types of loading (point load, distributed load and line load), iii) different concrete slab widths and iv) connection models with elastic perfectly plastic behaviour applied to TCC T-beam and a 5 m by 5 m floor system. The results of the conducted study provide an in-depth understanding of the effect and impact of the investigated parameters which can support a future optimisation of the structural solution. Moreover, the results show that the use of a simplified 2D model can support designers to overcome the limitations of analytical approaches, such as the Gamma method, to analyse large spacing and to take benefit of bidirectional effects (load sharing) if analysed as a whole instead of a single T-beam.

Implementing a simple 2D constitutive model for rocks into finite element method

Our research group has developed a 2D constitutive model for rocks that incorporates a strain-dependent elastic modulus. This model efficiently represents nonlinear stress–strain curves using just three equations and four parameters. In this study, we have implemented the simple 2D model into finite element codes to tackle engineering problems. We have calculated and compared stresses and displacements around a pressurized thick-walled hollow cylinder, a lined non-circular tunnel, and a rock slope using the elastic, simple 2D, and elasto-plastic models. The results confirmed that the simple 2D model closely matches the stress distribution obtained from the elastic solution at low-stress levels. Additionally, the simple 2D model produces a plastic zone with a larger inward displacement than the elasto-plastic model at higher stress levels. Only the simple 2D model captured sidewall convergence, roof sag, and floor heave in non-circular tunnels, and the toppling failure for 90° slope and toe translational slide for 60° slope in rock slopes with smooth critical slip surface. We did not encounter any convergence difficulties while solving the simple 2D model.