The main goal of this work was to develop numerical models for studying the behaviour of fibres in an accelerated flow. This is of special interest for e.g. papermaking. The early stage of the paper manufacturing process determines most of the final properties of a paper sheet. The complexity of studying the flow of fibre suspensions both experimentally and numerically emphasises a need for new ideas and developments. By means of solving the evolution of a convective-dispersion equation, i.e. the Fokker-Planck equation, a fully 3D approach with respect to the position and the two fibre angles, polar and azimuthal angles, following a streamline is presented. As an input to the fibre orientation model the turbulent flow field is solved by Computational Fluid Dynamics (CFD) with second-order closure in the turbulence model. In this work two new hypotheses have been presented for the variation of the non-dimensional rotational diffusivity with non-dimensional fibre length, Lf /η and the Reynolds number based on the Taylor micro-scale of the turbulence, Reλ Parameters for the two new hy- potheses and earlier models are determined with the aim of achieving a general relation and a value of the rotational dispersion coeffcient of stiff fibres in an anisotropic turbulent fluid flow. Earlier modelling work has been focused on solving the planar approach, i.e. assuming all fibres to be in one plane. This planar approach is discussed and compared with the fully 3D approach and its validity is evaluated. The optimization of parameters for the different hypotheses correlated on a central streamline, showed a good agreement with an independent experimental result in the undisturbed region. Moreover, it is particularly interesting that the boundary layer region and the wake region are predicted fairly well and the phenomena are well described, which has not been the case earlier. It seems that the new hypothesis based on the variation of the non-dimensional fibre length, Lf /η gives the best correlation in these shear-layer regions. Further- more it was established that the planar approach fails to predict shear layers, i.e. the boundary layer and the wake regions. As emphasized in the theory section, the planar formulation is strictly valid only if all fibres are oriented in one plane, which is not the case in the shear layers. In the undisturbed region, the 3D and the planar approaches, agree in their results. This leads to the conclusion that both approaches are suitable when shear layers are not studied. / QC 20100812

QC 20120613

Process industries are in general depending, in one way or the other, on fluid mechanics.Specifically, paper manufacturing, which probably is the dominant processindustry in Sweden, is depending on the flow of cellulose fibres suspended in water.As a part of the process the suspension, consisting of fibres in water, is spread out onor between two moving permeable weaves, i.e. wires. The speed of this is usually 10–30 m/s and the suspension is spread out by a plane jet issuing from a headbox nozzle.It has been show that the conditions in the headbox and jet have a large influence onthe quality of the final paper sheet. Primarily, streaks in the paper sheet are believedto be the result of streamwise streaks in the headbox jet.The thesis is aimed at the flow phenomena which occur in the headbox jet. Theinvestigations have been made with numerical calculations, stability theory and modelexperiments using water, as well as experiments with a real paper machine headboxand fibre suspension. In the thesis an introduction to the hydrodynamics of planeliquid jets is presented together with a description of the paper forming process andthe fluid mechanics of headbox flow.The basic flow and stability of a two-dimensional plane liquid jet has been investigatedby numerical calculations, stability theory and experiments. The calculationsof the laminar basic flow is successfully compared to pitot-tube measurements of thestreamwise velocity profile. By visualisations of the flow it is found that wave disturbanceson the jet has a severe effect on the flow. These waves can be predicted bylinear stability theory, which shows the presence of five convectively unstable modes.These can be divided into three types and by comparison with the experiments thetype of the visible waves is determined. These waves seem to initiate a break-up ofthe jet, which leads to strong streamwise streaks inside the jet.By flow visualisation of headbox flow of an experimental paper machine, togetherwith analysis of the resulting paper structure using the wavelet method the correspondencebetween flow disturbances and paper quality was investigated. It was shownthat the wave instability, which is present on the low Reynolds number water jet, alsocan be found in the real the headbox jet. It is shown that these waves play an importantrole in the dynamics of the headbox jet and also have an influence on the final papersheet. / QC 20100825

<p>A new model for simulating non-Brownian flexible fibres suspended in a Newtonian fluid has been developed. Special attention has been given to include realistic flow conditions found in the industrial papermaking process in the key features of the model; it is the intention of the author to employ the model in simulations of the forming section of the paper machine in future studies.</p><p>The model considers inert fibres of various shapes and finite stiffness, interacting with each other through normal, frictional and lubrication forces, and with the surrounding fluid medium through hydrodynamic forces. Fibre-fluid interactions in the non-creeping flow regime are taken into account, and the two-way coupling between the solids and the fluid phase is included by enforcing momentum conservation between phases. The incompressible three-dimensional Navier-Stokes equations are employed to model the motion of the fluid medium.</p><p>The validity of the model has been tested by comparing simulation results with experimental data from the literature. It was demonstrated that the model predicts the motion of isolated fibres in shear flow over a wide range of fibre flexibilities. It was also shown that the model predicts details of the orientation distribution of multiple straight, rigid fibres in a sheared suspension. Model predictions of the viscosity and first normal stress difference were in good agreement with experimental data found in the literature. Since the model is based solely on first-principles physics, quantitative predictions could be made without any parameter fitting.</p> / <p>En ny modell för simulering av rörelserna hos icke-Brownska böjliga fibrer dispergerade i en Newtonsk vätska har utvecklats. Eftersom det är författarens avsikt att modellen skall kunna tillämpas vid simulering av arkformning under de förhållanden som råder i en modern pappersmaskin, har särskilt omsorg givits till att inkludera motsvarande flödesvillkor i modellens giltighetsområde.</p><p>Modellen hanterar fibrer av varierande form, massa och styvhet, som växelverkar sinsemellan via normal-, friktions- och smörjkrafter. Deras växelverkan med den omgivande vätskan sker via hydrodynamiska krafter vid finita Reynolds-tal. Den så kallade tvåvägskopplingen mellan fibrerna och vätskefasen har tagits i beaktande genom att kräva att rörelsemängden bevaras vid interaktionen mellan faserna. Vidare har Navier-Stokes ekvationer för inkompressibla vätskor använts för att beskriva mediets rörelser.</p><p>Modellens giltighet har undersökts genom att jämföra resultat från simuleringar med experimentella data från litteraturen. Det har påvisats att modellen förutsäger rörelsen hos ensamma fibrer i ett skjuvflöde, för vitt skilda fiberflexibiliteter. Det visades också att modellen förutsäger detaljer hos fiberorienteringsdistributionen i suspensioner utsatta för skjuvflöde. Det kunde också konstateras att modellens förutsägelser av fibersuspensioners viskositet och första normalspänningsdifferens under skjuvning väl överensstämde med experimentella data i litteraturen. Kvantitativa förutsägelser har kunnat göras utan någon parameteranpassning, då modellen bygger uteslutande på väletablerade fysikaliska samband inom klassisk mekanik och strömningslära.</p>

A new model for simulating non-Brownian flexible fibres suspended in a Newtonian fluid has been developed. Special attention has been given to include realistic flow conditions found in the industrial papermaking process in the key features of the model; it is the intention of the author to employ the model in simulations of the forming section of the paper machine in future studies. The model considers inert fibres of various shapes and finite stiffness, interacting with each other through normal, frictional and lubrication forces, and with the surrounding fluid medium through hydrodynamic forces. Fibre-fluid interactions in the non-creeping flow regime are taken into account, and the two-way coupling between the solids and the fluid phase is included by enforcing momentum conservation between phases. The incompressible three-dimensional Navier-Stokes equations are employed to model the motion of the fluid medium. The validity of the model has been tested by comparing simulation results with experimental data from the literature. It was demonstrated that the model predicts the motion of isolated fibres in shear flow over a wide range of fibre flexibilities. It was also shown that the model predicts details of the orientation distribution of multiple straight, rigid fibres in a sheared suspension. Model predictions of the viscosity and first normal stress difference were in good agreement with experimental data found in the literature. Since the model is based solely on first-principles physics, quantitative predictions could be made without any parameter fitting. / En ny modell för simulering av rörelserna hos icke-Brownska böjliga fibrer dispergerade i en Newtonsk vätska har utvecklats. Eftersom det är författarens avsikt att modellen skall kunna tillämpas vid simulering av arkformning under de förhållanden som råder i en modern pappersmaskin, har särskilt omsorg givits till att inkludera motsvarande flödesvillkor i modellens giltighetsområde. Modellen hanterar fibrer av varierande form, massa och styvhet, som växelverkar sinsemellan via normal-, friktions- och smörjkrafter. Deras växelverkan med den omgivande vätskan sker via hydrodynamiska krafter vid finita Reynolds-tal. Den så kallade tvåvägskopplingen mellan fibrerna och vätskefasen har tagits i beaktande genom att kräva att rörelsemängden bevaras vid interaktionen mellan faserna. Vidare har Navier-Stokes ekvationer för inkompressibla vätskor använts för att beskriva mediets rörelser. Modellens giltighet har undersökts genom att jämföra resultat från simuleringar med experimentella data från litteraturen. Det har påvisats att modellen förutsäger rörelsen hos ensamma fibrer i ett skjuvflöde, för vitt skilda fiberflexibiliteter. Det visades också att modellen förutsäger detaljer hos fiberorienteringsdistributionen i suspensioner utsatta för skjuvflöde. Det kunde också konstateras att modellens förutsägelser av fibersuspensioners viskositet och första normalspänningsdifferens under skjuvning väl överensstämde med experimentella data i litteraturen. Kvantitativa förutsägelser har kunnat göras utan någon parameteranpassning, då modellen bygger uteslutande på väletablerade fysikaliska samband inom klassisk mekanik och strömningslära.