In accordance with our Quality Policy, COMSOL maintains a library of hundreds of documented model examples that are regularly tested against the latest version of the COMSOL Multiphysics® software, including benchmark problems from ASME and NAFEMS, as well as TEAM problems.

Our Verification and Validation (V&V) test suite provides consistently accurate solutions that are compared against analytical results and established benchmark data. The documented models below are part of the COMSOL Multiphysics® software’s built-in Application Libraries. They include reference values and sources for a wide range of benchmarks, as well as step-by-step instructions to reproduce the expected results on your own computer. You can use these models not only to document your software quality assurance (SQA) and numerical code verification (NCV) efforts, but also as part of an in-house training program.

Micromechanical Model of a Piezoelectric Fiber Composite

In this example, the micromechanical properties of a piezoelectric fiber composite are studied. The homogenized electromechanical properties of the composite are derived from the individual microscopic properties of matrix and fiber. Read More

Expansion Fan

An important and interesting phenomenon with supersonic flows are expansion fans, which take place when the flow encounters a convex or expansion corner. The direction of the flow changes smoothly across the fan, while the Mach number increases. This 2D example models an expansion fan ... Read More

Lemaitre–Chaboche Viscoplastic Model

Most metals and alloys undergo viscoplastic deformation at high temperatures. In case of cyclic loading, a constitutive law with both isotropic and kinematic hardening is necessary to describe effects such as ratcheting, cyclic softening/hardening, and stress relaxation. The Lemaitre ... Read More

A Cross-Bridge Kelvin Resistor Model for the Extraction of Specific Contact Resistivity

This benchmark example builds two models of a cross-bridge Kelvin resistor used for extracting the specific contact resistivity. The first model simulates the system in 3D, using the contact resistance feature built in the Semiconductor interface. The other model is a 2D approximation of ... Read More

Failure Prediction in a Layered Shell

Laminated shells made of carbon fiber reinforced polymer (CRFP) are common in a large variety of applications due to their high strength-to-weight ratio. The widespread attention to laminated shells is due to their diverse applications in aerospace, marine, automotive, and various ... Read More

Total Internal Reflection Thin-Film Achromatic Phase Shifter (TIRTF APS)

The capability to alter the polarization of light is crucial to a wide variety of optical devices. For example, the polarization of light has a significant effect on the performance of optical isolators, attenuators, and beam splitters. By assigning a specific polarization to light, most ... Read More

Eppler Airfoil Transition

The flow around the Eppler 387 airfoil is computed with the SST turbulence model both with and without the transition model. The reuslts are compared with experimental values. Read More

Two-Stage Powder Compaction Process

Powder compaction is a key process in powder metallurgy, where it gives the flexibility to produce quality products of complex shapes for sintering. The density of the compact is a key factor to determine the overall quality of the sintered product, as regions with lower density could ... Read More

Hartmann Boundary Layer

Classical MHD benchmark problem was solved analytically by J. Hartmann. He considered laminar incompressible flow between two planes (in planar duct) in transversal imposed magnetic field under next assumptions: * fully developed flow (far from inlet) * fluid properties are constant * ... Read More

Inertial Focusing Between Two Parallel Walls

For more than 50 years, it has been known that neutrally buoyant particles in a flow channel tend to converge to specific locations in the channel cross section. For a cylindrical pipe, or two parallel planes carrying a Poiseuille flow, the equilibrium position is about 0.6 times the ... Read More

s