Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

An Energy Conserving Finite Difference Scheme for Simulation of Collisions

Nonlinear phenomena play an essential role in the sound production process of many musical instruments. A common source of these effects is object collision, the numerical simulation of which is known to give rise to stability
issues. This paper presents a method to construct numerical schemes that conserve the total energy in simulations of one-mass systems involving collisions, with no conditions imposed on any of the physical or numerical parameters.
This facilitates the adaptation of numerical models to experimental data, and allows a more free parameter adjustment in sound synthesis explorations. The energy preservedness of the proposed method is tested and demonstrated though several examples, including a bouncing ball and a non-linear oscillator, and implications regarding the wider applicability are discussed.

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

Fourier´s law from a chain of coupled anharmonic oscillators under energy conserving shot noise.

We analyze the transport of heat along a chain of particles interacting through anharmonic potentials consisting of quartic terms in addition to harmonic quadratic terms and subject to heat reservoirs at its ends. Each particle is also subject to an impulsive shot noise with exponentially distributed waiting times whose effect is to change the sign of its velocity, thus conserving the energy of the chain. We show that the introduction of this energy conserving stochastic noise leads to Fourier's law. That is for large system size L the heat current J behaves as J ‘approximately’ 1/L, which amounts to say that the conductivity k is constant. The conductivity is related to the current by J = kΔT/L, where ΔT is the difference in the temperatures of the reservoirs. The behavior of heat conductivity k for small intensities¸ of the shot noise and large system sizes L are obtained by assuming a scaling behavior of the type k = ‘L POT a Psi’(L’lambda POT a/b’) where a and b are scaling exponents. For the pure harmonic case a = b = 1, characterizing a ballistic conduction of heat when the shot noise is absent. For the anharmonic case we found values for the exponents a and b smaller then 1 and thus consistent with a superdiffusive conduction of heat without the shot noise. We also show that the heat conductivity is not constant but is an increasing function of temperature.

Heat conduction in a chain of harmonic and anharmonic oscillators under the presence of an energy conserving noise

Reproducing Fourier's law of heat conduction from a microscopic stochastic model is a long standing challenge in statistical physics. As was shown by Rieder, Lebowitz and Lieb many years ago, a chain of harmonically coupled oscillators connected to two heat baths at different temperatures does not reproduce the diffusive behaviour of Fourier's law, but instead a ballistic one with an infinite thermal conductivity. Since then, there has been a substantial effort from the scientific community in identifying the key mechanism necessary to reproduce such diffusivity, which usually revolved around anharmonicity and the effect of impurities. Recently, it was shown by Dhar, Venkateshan and Lebowitz that Fourier's law can be recovered by introducing an energy conserving noise, whose role is to simulate the elastic collisions between the atoms and other microscopic degrees of freedom, which one would expect to be present in a real solid. For a one-dimensional chain this is accomplished numerically by randomly flipping - under the framework of a Poisson process with a variable “rate of collisions" - the sign of the velocity of an oscillator. In this poster we present Langevin simulations of a one-dimensional chain of oscillators coupled to two heat baths at different temperatures. We consider both harmonic and anharmonic (quartic) interactions, which are studied with and without the energy conserving noise. With these results we are able to map in detail how the heat conductivity k is influenced by both anharmonicity and the energy conserving noise. We also present a detailed analysis of the behaviour of k as a function of the size of the system and the rate of collisions, which includes a finite-size scaling method that enables us to extract the relevant critical exponents. Finally, we show that for harmonic chains, k is independent of temperature, both with and without the noise. Conversely, for anharmonic chains we find that k increases roughly linearly with the temperature of a given reservoir, while keeping the temperature difference fixed.

Energy budget procedures and performance criteria for energy conserving building illumination systems /

"April 1980, issued May 1980."

Performance of passive solar and energy conserving houses in California : final subcontract report /

November 1983.

Kinetic energy conserving integrators for Gaussian thermostatted SLLOD

A new integration scheme is developed for nonequilibrium molecular dynamics simulations where the temperature is constrained by a Gaussian thermostat. The utility of the scheme is demonstrated by its application to the SLLOD algorithm which is the standard nonequilibrium molecular dynamics algorithm for studying shear flow. Unlike conventional integrators, the new integrators are constructed using operator-splitting techniques to ensure stability and that little or no drift in the kinetic energy occurs. Moreover, they require minimum computer memory and are straightforward to program. Numerical experiments show that the efficiency and stability of the new integrators compare favorably with conventional integrators such as the Runge-Kutta and Gear predictor-corrector methods. (C) 1999 American Institute of Physics. [S0021-9606(99)50125-6].

Changes in resting and walking energy expenditure and walking speed during pregnancy in obese women

Background Energy conserving processes reported in undernourished women during pregnancy are a recognised strategy to provide energy required to support fetal development. Women who are obese before conceiving arguably have sufficient fat stores to support the energy demands of pregnancy without the need to provoke energy conserving mechanisms. Objective We tested the hypothesis that obese women would demonstrate behavioural adaptation (i.e. decrease in self-selected walking (SSW) speed) but not metabolic compensation (i.e. decrease in resting metabolic rate (RMR) or metabolic cost of walking) during gestation. Design RMR, SSW speed, metabolic cost of walking, and anthropometry were measured in 23 women (BMI: 33.6 ± 2.5 kg/m2; 31 ± 4 years) at approximately weeks 15 (wk 15) and 30 (wk 30) of gestation. RMR was also measured in two cohorts of non-pregnant controls matched for age, weight and height of the pregnant cohort at wk 15 (N=23) and wk 30 (N=23). Results GWG varied widely (11.3 ± 5.4 kg) and 52% of women gained more weight than is recommended. RMR increased significantly by an average 177 ± 176 kcal/d (11±12%; P<0.0001); however the within-group variability was large. Both the metabolic cost of walking and SSW speed decreased significantly (P<0.01). While RMR increased in >80% of the cohort, the net oxygen cost of walking decreased in the same proportion of women. Conclusions While the increase in RMR was greater than was explained by weight gain, there was evidence of both behavioural and biological compensation in the metabolic cost of walking in obese women during gestation.

Information capacity of energy harvesting sensor nodes

Sensor nodes with energy harvesting sources are gaining popularity due to their ability to improve the network life time and are becoming a preferred choice supporting `green communication'. We study such a sensor node with an energy harvesting source and compare various architectures by which the harvested energy is used. We find its Shannon capacity when it is transmitting its observations over an AWGN channel and show that the capacity achieving energy management policies are related to the throughput optimal policies. We also obtain the capacity when energy conserving sleep-wake modes are supported and an achievable rate for the system with inefficiencies in energy storage.