GeomBD3 is a robust Brownian dynamics simulation package designed to easily handle natural or engineered systems in diverse environments and arrangements. The software package described herein allows users to design, execute, and analyze BD simulations. The simulations use all-atom, rigid molecular models that diffuse according to overdamped Langevin dynamics and interact through electrostatic, Lennard-Jones, and ligand desolvation potentials. The program automatically calculates molecular association rates, surface residence times, and association statistics for any number of user-defined association criteria. Users can also extract molecular association pathways, diffusion coefficients, intermolecular interaction energies, intermolecular contact probability maps, and more using the provided supplementary analysis scripts. We detail the use of the package from start to finish and apply it to a protein-ligand system and a large nucleic acid biosensor. GeomBD3 provides a versatile tool for researchers from various disciplines that can aid in rational design of engineered systems or play an explanatory role as a complement to experiments. GeomBD version 3 is available on our website at http://chemcha-gpu0.ucr.edu/geombd3/ and KBbox at https://kbbox.h-its.org/toolbox/methods/molecular-simulation/geombd/.