Tag Archives: Multiphysics

Ten research papers utilizing Energy2D published in the past two years

Screenshots from recent papers that use Energy2D
Energy2D simulation of fire
Energy2D is a multiphysics simulation program that was created from scratch and is still under development (though its progress has slowed down significantly because my priority has been given to its Energy3D cousin). The software was originally intended to be a teaching and learning tool for high school students who are interested in studying engineering. Over the past two years, however, we have seen 10 research papers published in various journals and conferences that involved significant applications of Energy2D as a scientific research tool for modeling natural phenomena and engineering systems. The problems that these researchers simulated range from solar energy, industrial processes, geophysics, and building science. The authors come from universities from all over the world, including top-notch institutions in US, Europe, and China.

Energy2D simulation of thermal bridge
Among them, researchers from Delft University of Technology, Technical University of Darmstadt, and Eindhoven University of Technology wrote in their recent paper about the validity of Energy2D: "The software program Energy2D is used to solve the dynamic Fourier heat transfer equations for the Convective Concrete case. Energy2D is a relatively new program (Xie, 2012) and is not yet widely used as a building performance simulation tool. To gain more confidence in the predictions with Energy2D, an analytical validation study was therefore carried out first, inspired by the approach described in Hensen and Nakhi (1994). Those analytical solutions and the simulation results of the dynamic response to a 20°C temperature step change on the surface of a concrete construction with the following properties were compared for this research." They concluded that "the simulation results never divert from the exact solution more than 0.45°C and it is therefore considered acceptable to further use this model."

The publication of these papers and very positive user feedback suggest that Energy2D seems to have found itself an interesting niche market. Many scientists and engineers are unable to invest a lot of time and money on its complicated commercial counterparts. But they nonetheless need a handy simulation tool that is much more flexible, intuitive, and capable than formulas in books to deal with realistic geometry -- at least in 2D. This is where Energy2D comes into play.

Reaching this milestone is critically important to the free and open-source Energy2D software, whose future will be reliant on community support. Its modest popularity among scientists is a valid demonstration of the broader impact expected by the National Science Foundation that funded its development. One can only imagine that there are many more users who used the software in their workplace but didn't publish. Now that good words about it have spread, we expect the usage to continue and even accelerate. To better support our users, we have added a community forum recently. We also plan to work with Professor Bob Hanson to port the Java code to JavaScript through his SwingJS translator so that the program can run on more devices.

The list of these papers is as follows:
  1. Mahfoud Abderrezek & Mohamed Fathi, Experimental Study of the Dust Effect on Photovoltaic Panels' Energy Yield, Solar Energy, Volume 142, pp 308-320, 2017
  2. Dennis de Witte, Marie L. de Klijn-Chevalerias, Roel C.G.M. Loonen, Jan L.M. Hensen, Ulrich Knaack, & Gregor Zimmermann, Convective Concrete: Additive Manufacturing to Facilitate Activation of Thermal Mass, Journal of Facade Design and Engineering, Volume 5, No. 1, 2017
  3. Javier G. Monroy & Javier Gonzalez-Jimenez, Gas Classification in Motion: An Experimental Analysis, Sensors and Actuators B: Chemical, Volume 240, pp 1205-1215, 2017
  4. Tom Rainforth, Tuan Anh Le, Jan-Willem van de Meent, Michael A. Osborne, & Frank Wood, Bayesian Optimization for Probabilistic Programs, 30th Conference on Neural Information Processing Systems, Barcelona, Spain, 2016
  5. E. Rozos, I. Tsoukalas, & C. Makropoulos, Turning Black into Green: Ecosystem Services from Treated Wastewater, 13th IWA Specialized Conference on Small Water and Wastewater Systems, Athens, Greece, 2016
  6. W. Taylor Shoulders, Richard Locke, & Romain M. Gaume, Elastic Airtight Container for the Compaction of Air-Sensitive Materials, Review of Scientific Instruments, Volume 87, 063908, 2016
  7. Zachary R. Adam, Temperature Oscillations near Natural Nuclear Reactor Cores and the Potential for Prebiotic Oligomer Synthesis, Origins of Life and Evolution of Biospheres, Volume 46, Issue 2, pp 171-187, 2016
  8. Jiarui Chen, Shuyu Qin, Xinglong Wu, & Paul K Chu, Morphology and Pattern Control of Diphenylalanine Self-Assembly via Evaporative Dewetting, ACS Nano, Volume 10, No. 1, pp 832-838, 2016
  9. Atanas Vasilev, Geothermal Evolution of Gas Hydrate Deposits: Bulgarian Exclusive Economic Zone in the Black Sea, Comptes rendus de l‘Académie bulgare des Sciences, Volume 68, No. 9, pp 1135-1144, 2015
  10. Pedro A. Hernández, et al., Magma Emission Rates from Shallow Submarine Eruptions Using Airborne Thermal Imaging, Remote Sensing of Environment, Volume 154, pp 219-225, November 2014

Scanning radiation flux with moving sensors in Energy2D

Figure 1: Moving sensors facing a rectangular radiator.
The heat flux sensor in Energy2D can be used to measure radiative heat flux, as well as conductive and convective heat fluxes. Radiative heat flux depends on not only the temperature of the object the sensor measures but also the angle at which it faces the object. The latter is known as the view factor.

In radiative heat transfer, a view factor between two surfaces A and B is the proportion of the radiation which leaves surface A that strikes surface B. If the two surfaces face each other directly, the view factor is greater than the case in which they do not. If the two surfaces are closer, the view factor is greater.

Figure 2: Rotating sensors inside and outside a ring radiator.
To conveniently visualize the effect of a view factor, Energy2D allows you to attach a heat flux sensor to a moving or rotating particle, with a settable linear or angular velocity. In this way, we can set up sensors to automatically "scan" the field of radiation heat flux like a radar.

Figure 1 shows a moving sensor and a rotating sensor, as well as the data they record. A third sensor is also placed to the right of an object that is being heated by the radiator. This object has an emissivity of one so it also radiates. Its radiation flux is recorded by the third sensor whose data shows a slowly increasing heat flux as the object slowly warms up.

As an interesting test case, Figure 2 shows two rotating sensors, one placed precisely at the center of a ring radiator and the other outside. The almost steady line recorded by the first sensor suggests that the view factor at the center does not change, which makes sense. The small sawtooth shape is due to the limitation of discretization in our numerical simulation.

Multiphysics simulations of inelastic collisions with Energy2D

Figure 1. Mechano-thermal simulation of inelastic collision.
Many existing simulations of inelastic collisions show the changes of speeds and energy of the colliding objects without showing what happens to the lost energy, which is often converted into thermal energy that spreads out through heat transfer. With the new multiphysics modeling capabilities, the Energy2D software can show the complete picture of energy transfer from the mechanical form to the thermal form in a single simulation.

Figure 2. Thermal marks left by collisions.
Figure 1 shows the collisions of three identical balls (mass = 10 kg, speed = 1 m/s) with three fixed objects that have different elasticities (0, 0.5, and 1). The results show that, in the case of the completely inelastic collision, all the kinetic energy of the ball (5 J) is converted into thermal energy of the rectangular hit object (at this point, the particles in Energy2D do not hold thermal energy, but this will be changed in a future version), whereas in the case of completely elastic collision, the ball B1 does not lose any kinetic energy to the hit object. In the cases of inelastic collisions, you can see the thermal marks created by the collisions. The thermometers placed in the objects also register a rise of temperatures. This view resembles infrared images of floors taken immediately after being hit by tennis balls.

Figure 3. Collisions in Energy2D.
Energy2D supports particle collisions with all the 2D shapes that it provides: rectangles, ellipses, polygons, and blobs. Figure 2 shows the thermal marks on two blobs created by a few bouncing particles. And Figure 3 shows another simulation of collision dynamics with a lot of particles bouncing off complex shapes (boy, it took me quite a while in this July 4 weekend to hunt down most of the bugs in the collision code).

The multiphysics functionality of Energy2D is an exciting new feature as it allows more realistic modeling of natural phenomena. Even in science classrooms, realism of simulations is not just something that is nice to have. If computer simulations are to rival real experiments, it must produce not only the expected effects but also the unexpected side effects. Capable of achieving just that, a multiphysics simulation can create a deep and wide learning space just like real experiments. For engineering design, this depth and breadth are not options -- there is no open-endedness without this depth and breadth and there is no engineering without open-endedness.

Towards a multiphysics Energy2D

Figure 1: Particle motions driven by convective flow.
Up to yesterday, our Energy2D software has been a program for simulating, mostly, fluid and heat flows. But there are also objects in the world that are not fluids. To simulate that part of the world, we have to incorporate some other physics. A simple addition is to couple particles with fluids. This technique is commonly known as discrete phase modeling in the CFD community. It is used to model things such as suspension particles in fluids.

Figure 2: Heat traces of fireballs.
The latest version of Energy2D has a particle solver and a particle editor. Particles in Energy2D observe collision dynamics among themselves and interact with fluid and heat flows: particles can not only be moved by the fluid but also exert reaction force and transfer heat to the fluid. Figure 1 shows the motion of two types of particles driven by a convective flow. Depending on its density (relative to the fluid density), a particle may be buoyant enough to flow with the fluid or so heavy that it must sink to the bottom. This is shown in Figure 1: The black particles are the heavy ones and the white ones are the light ones; the convective force is not strong enough to move the black ones.

Particles can also transfer physical properties such as energy and momentum to the fluid while they are moving. Figure 2 shows the heat traces left by fireballs of different sizes.

Figure 3: Thermophoresis (Soret's effect)
With this new capacity, we can simulate phenomena such as thermophoresis, in which the different particle types in a mixture respond to a temperature gradient differently and thereby can be separated by just heating them up.

If you are enticed enough to want to see these simulations at work, click the links below the figures.

These new features represent an overdue step towards making Energy2D a versatile multiphysics simulation system. For engineering simulations, multiphysics is essential as real-world problems are often complicated by more than one mechanisms, each driven by its own physics.

The particle dynamics shown here is very simple (just a weekend's work). In the long run, I expect that a generic contact dynamics engine such as that of Box2D will be implemented in Energy2D. Coupling the Eulerian and Lagrangian reference frames, this integration will make Energy2D more interesting and useful. That would be a critical step towards our goal for Energy2D to simulate as many energy-related natural phenomena as possible.